數字信號處理實驗.doc_第1頁
數字信號處理實驗.doc_第2頁
數字信號處理實驗.doc_第3頁
數字信號處理實驗.doc_第4頁
全文預覽已結束

下載本文檔

版權說明：本文檔由用戶提供并上傳，收益歸屬內容提供方，若內容存在侵權，請進行舉報或認領

文檔簡介

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)系上傳者。文件的所有權益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙，網頁內容里面會有圖紙預覽，若沒有圖紙預覽就沒有圖紙。
  • 4. 未經權益所有人同意不得將文件中的內容挪作商業(yè)或盈利用途。
  • 5. 人人文庫網僅提供信息存儲空間，僅對用戶上傳內容的表現方式做保護處理，對用戶上傳分享的文檔內容本身不做任何修改或編輯，并不能對任何下載內容負責。
  • 6. 下載文件中如有侵權或不適當內容，請與我們聯(lián)系，我們立即糾正。
  • 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

最新文檔

評論

0/150

提交評論