數(shù)字信號處理實(shí)驗(yàn).doc_第1頁
數(shù)字信號處理實(shí)驗(yàn).doc_第2頁
數(shù)字信號處理實(shí)驗(yàn).doc_第3頁
數(shù)字信號處理實(shí)驗(yàn).doc_第4頁
全文預(yù)覽已結(jié)束

下載本文檔

版權(quán)說明：本文檔由用戶提供并上傳，收益歸屬內(nèi)容提供方，若內(nèi)容存在侵權(quán)，請進(jìn)行舉報(bào)或認(rèn)領(lǐng)

文檔簡介

function y=circonvt(x1,x2,N)if length(x1)N error(N =x1 )endif length(x2)N error(N =x2 )endx1=x1 zeros(1,N-length(x1);x2=x2 zeros(1,N-length(x2);m=0:1:N-1;x2=x2(mod(-m,N)+1);H=zeros(N,N); for n=1:1:N H(n,:)=cirshift(x2,n-1,N);end y=x1*H;function y=cirshift(x,m,N)if length(x)N error(N =x )endx=x zeros(1,N-length(x);n=0:1:N-1;n=mod(n-m,N);y=x(n+1);clc;clear alln=0:1:11;m=0:1:5;g1=1 2 3 4 5;g2=2 2 0 1 1;y1n=conv(g1,g2);ycn=circonvt(g1,g2,5);ny1=0:1:length(y1n)-1;ny2=0:1:length(ycn)-1;subplot(2,1,1);stem(ny1,y1n);subplot(2,1,2);stem(ny2,ycn);function y=circonvt(x1,x2,N)if length(x1)N error(N =x1 )endif length(x2)N error(N =x2 )endx1=x1 zeros(1,N-length(x1);x2=x2 zeros(1,N-length(x2);m=0:1:N-1;x2=x2(mod(-m,N)+1);H=zeros(N,N); for n=1:1:N H(n,:)=cirshift(x2,n-1,N);end y=x1*H;function y=cirshift(x,m,N)if length(x)N error(N =x )endx=x zeros(1,N-length(x);n=0:1:N-1;n=mod(n-m,N);y=x(n+1);clc;clear alln=1:2:17;m=0:1:5;x1=1 3 5 7 9 11 13 15 17;x2=1 -2 3 -2 1;y1n=conv(x1,x2);ycn=circonvt(x1,x2,9);ny1=0:1:length(y1n)-1;ny2=0:1:length(ycn)-1;subplot(2,1,1);stem(ny1,y1n);subplot(2,1,2);stem(ny2,ycn);function Y=overlapsav(x,h,N)Lx=length(x);M=length(h);M1=M-1;L=N-M1; h=h,zeros(1,N-M); x=zeros(1,M1),x,zeros(1,N-1);a=floor(Lx+M1-1)/(L)+1;Y=zeros(1,N);for k=0:a-1 xk=x(k*L+1:k*L+N); b=fft(xk,N); C=fft(h,N); Z=b.*C; Y(k+1,:)=ifft(Z,N);endY=Y(:,M:N);Y=(Y(:)clc;clear allx=1 2 3 4 5 6 7 8 9 10;h=1 0 -1;N=6;Y=overlapsav(x,h,N);Y = Columns 1 through 8 1.0000 2.0000 2.0000 2.0000 2.0000 2.0000 2.0000 2.0000 Columns 9 through 122.0000 2.0000 -9.0000 -10.0000 Y=fftfilt(x,h,N)Y = 1.0000 2.0000 2.0000clc;clear allx1=1 3 5 7 9 11 13 15 17;x2=1 -2 3 -2 1;y=conv(x1,x2);stem (y);clc;clear allx1=1 3 5 7 9 11 13 15 17;x2=1 -2 3 -2 1;N=length(x1) + length(x2)-1;n=0:N-1;x1=x1,zeros(1,N-length(x1);x2=x2,zeros(1,N-length(x2);X1=fft(x1);X2=fft(x2);Y=X1.*X2;y=ifft(Y);subplot(221),stem(n,x1,.);xlabel(n);ylabel(x1(n);title( x1(n);grid;subplot(222),stem(n,x2,.);xlabel(n);ylabel(x2(n);title( x2(n);grid;subplot(223),stem(n,y,.);xlabel(n);ylabel(y(n);title( y(n);grid;clc;clear allx1=1 3 5 7 9 11 13 15 17;x2=1 -2 3 -2 1;N=length(x1) + length(x2)-1;n=0:N-1;x1=x1,zeros(1,N-length(x1);x2=x2,zeros(1,N-length(x2);y0=conv(x1,x2);Y=x1.*x2;y1=ifft(Y);subplot(121),stem(y0,.);subplot(122),stem(n,y1,.);clc;clear allx1=1 3 5 7 9 11 13 15 17;x2=1 -2 3 -2 1;N=length(x1) + length(x2)-1;n=0:N-1;x1=x1,zeros(1

溫馨提示

  • 1. 本站所有資源如無特殊說明，都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等，如果需要附件，請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙，網(wǎng)頁內(nèi)容里面會有圖紙預(yù)覽，若沒有圖紙預(yù)覽就沒有圖紙。
  • 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 人人文庫網(wǎng)僅提供信息存儲空間，僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理，對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯，并不能對任何下載內(nèi)容負(fù)責(zé)。
  • 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容，請與我們聯(lián)系，我們立即糾正。
  • 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

評論

0/150

提交評論