數(shù)字信號處理 系統(tǒng)響應

1、數(shù)字信號處理實驗報告信號、系統(tǒng)及系統(tǒng)響應默認分類 2007-12-05 23:36:03 閱讀165 評論0  字號：大中小 訂閱1.實驗目的2.實驗環(huán)境（軟件、硬件及條件）   Matlab6.5及其以上版本 3.實驗方法 4.實驗分析 （1）編寫信號產(chǎn)生子程序 5.實驗步驟a.取采樣頻率fs=1kHz,觀察所得采樣xa(n)的幅頻特性|X(ejw)|.b.改變采樣頻率fs=300Hz,觀察|X(ejw)|的變化，并作記錄;進一步降低采樣頻率fs=200Hz,觀察

2、頻譜是否混疊.6.實驗結(jié)論  思考題解答：     1.數(shù)字頻率度量不同，模擬頻率相同。因為模擬頻率是連續(xù)域的傅立葉變換，與采樣頻率無關(guān)，而數(shù)字頻率與模擬頻率是T倍的關(guān)系，T 是采樣周期，所以數(shù)字頻率會隨著采樣頻率的不同而不同，而模擬頻率不會。     2.沒有差別，因為我用的是自己編寫的傅立葉變換，而沒有用FFT，假如用FFT則有差異。結(jié)論：1.當采樣頻率減小到fs<2fh時，所采樣的信號頻譜有混疊，所以采樣的頻率應該大于或等于原信號最高頻率的兩倍.2.信號

3、經(jīng)過系統(tǒng)的響應經(jīng)過系統(tǒng)函數(shù)與信號函數(shù)卷積運算可得到，在經(jīng)過傅立葉變換就可得到響應得頻域特性.3.將信號的頻域函數(shù)與系統(tǒng)的頻域函數(shù)相乘的到的函數(shù)與第二步得到的函數(shù)相同，從而驗證了卷積定理.1.子程序代碼functionx,n=signal1(A,a,q,f)n=0:49;x=A*exp(-a*n*(1/f).*sin(q*n*(1/f); function x=signal2(n,n0,n1)x=n1*(n-n0)=0; function x=signal3(n,N)x=(n-0)>=0-n-0-N>=0; function 

4、;X=dtft(x,w)X=x*exp(-j*1:length(x)'*w); 實驗第一步的代碼及結(jié)果function main()         %ʵÑéÒ»µÚÒ»²½A=444.128;a=50*sqrt(2)*pi;Q=a;f0=1000;f1=300;f2=200;x0,n=signal1(A,a,Q,f0);x1,n=signal1(A,a,Q

5、,f1);x2,n=signal1(A,a,Q,f2);w=0:0.01:2*pi;X0=dtft(x0,w);X1=dtft(x1,w);X2=dtft(x2,w); figure;subplot(3,2,1);stem(n,x0,'.');title('x(n)²¨ÐÎ');xlabel('n');ylabel('x(n)'); subplot(3,2,2);plot(w,abs(X0);title('X(w)²¨ÐÎ(

6、f=1kHz)');xlabel('w');ylabel('X(w)'); subplot(3,2,3);stem(n,x1,'.');title('x(n)²¨ÐÎ');xlabel('n');ylabel('x(n)'); subplot(3,2,4);plot(w,abs(X1);title('X(w)²¨ÐÎ(f=300Hz)');xlabel('w')

7、;ylabel('X(w)'); subplot(3,2,5);stem(n,x2,'.');title('x(n)²¨ÐÎ');xlabel('n');ylabel('x(n)'); subplot(3,2,6);plot(w,abs(X2);title('X(w)²¨ÐÎ(f=200Hz)');xlabel('w');ylabel('X(w)'); 

8、60;     2.實驗第二步代碼及結(jié)果function main0()  %ʵÑéÒ»µÚ¶þ²½n=0:10;xb=signal2(n,0,1);hb=signal2(n,0,1)+signal2(n,1,2.5)+signal2(n,2,2.5)+signal2(n,3,1);yb=conv(xb,hb);w=0:0.01:2*pi;Xb=dtft(xb,w);Hb=dtft(hb,w);

9、Yb=dtft(yb,w); figure;subplot(3,2,1);stem(n,xb,'.');title('xb(n)²¨ÐÎ');xlabel('n');ylabel('xb(n)'); subplot(3,2,2);plot(w,abs(Xb);title('Xb(w)²¨ÐÎ');xlabel('w');ylabel('Xb(w)'); subplot(3,2

10、,3);stem(n,hb,'.');title('hb(n)²¨ÐÎ');xlabel('n');ylabel('hb(n)'); subplot(3,2,4);plot(w,abs(Hb);title('Hb(w)²¨ÐÎ');xlabel('w');ylabel('Hb(w)'); subplot(3,2,5);stem(0:20,yb,'.');title(&#

11、39;yb(n)²¨ÐÎ');xlabel('n');ylabel('yb(n)'); subplot(3,2,6);plot(w,abs(Yb);title('Yb(w)²¨ÐÎ');xlabel('w');ylabel('Yb(w)'); n=0:20;xc1=signal3(n,10);ha=signal3(n,10);xc2=signal3(n,5);y1=conv(xc1,ha);y2=conv(x

12、c2,ha);Xc1=dtft(xc1,w);Xc2=dtft(xc2,w);Ha=dtft(ha,w);Y1=dtft(y1,w);Y2=dtft(y2,w); figure;subplot(3,2,1);stem(n,xc1,'.');title('xc(n)²¨ÐÎ(N=10)');xlabel('n');ylabel('xc(n)'); subplot(3,2,2);plot(w,abs(Xc1);title('Xc(w)²¨Ð

13、;Î');xlabel('w');ylabel('Xc(w)'); subplot(3,2,3);stem(n,ha,'.');title('ha(n)²¨ÐÎ');xlabel('n');ylabel('ha(n)'); subplot(3,2,4);plot(w,abs(Ha);title('Ha(w)²¨ÐÎ');xlabel('w');ylab

14、el('Ha(w)'); subplot(3,2,5);stem(0:40,y1,'.');title('y(n)²¨ÐÎ');xlabel('n');ylabel('y(n)'); subplot(3,2,6);plot(w,abs(Y1);title('Y(w)²¨ÐÎ');xlabel('w');ylabel('Y(w)'); figure;subplo

15、t(2,2,1);stem(n,xc2,'.');title('xc(n)²¨ÐÎ(N=5)');xlabel('n');ylabel('xc(n)'); subplot(2,2,2);plot(w,abs(Xc2);title('Xc(w)²¨ÐÎ');xlabel('w');ylabel('Xc(w)'); subplot(2,2,3);stem(0:40,y2,'.

16、9;);title('y(n)²¨ÐÎ');xlabel('n');ylabel('y(n)'); subplot(2,2,4);plot(w,abs(Y2);title('Y(w)²¨ÐÎ');xlabel('w');ylabel('Y(w)'); 實驗第三步代碼及結(jié)果function main00()   %ʵÑéÒ»µÚÈý²½A=1;a=0.4;Q=2.0734;f=1;n=0:49;xa,n=signal1(A,a,Q,f);hb=signal2(n,0,1)+signal2(n,1,2.5)+signal2(n,2,2.5)+signal2(n,3,1);y1=conv(xa,hb);w=0:0.01:2*pi;Xa=dtft(xa,w);Hb=dtft(hb,w);Y1=dtft(y1,w);Y2=Xa.*Hb;