数字信号处理第一次上机实验报告.docVIP

数字信号处理实验报告 实验一：信号、系统及系统响应 1, 采样后的频谱变化 %用1000Hz采样频率： T1=1/1000;?????????????????????????? n=0:64/(1000*T1);??????????????????? A=444.128;?a=50*sqrt(2.0);w0=50*sqrt(2.0)*pi;? xn1=A*exp(-a*n*T1).*sin(w0*n*T1);????? Xk1=fft(xn1,1024);???????? subplot(322);stem(n,xn1,.);?grid?on? xlabel(n);ylabel(x(n));title(1000Hz采样x(n));? k=0:1023;wk=2*k/1024;? subplot(321);plot(wk,abs(Xk1));grid?on;? xlabel(\omega/\pi);ylabel(|X(e^j^\omega)|);title(FT[x(n)]);? ? %用300Hz采样频率 ? T2=1/300;?n2=0:64/(1000*T2);? A=444.128;?a=50*sqrt(2.0);w0=50*sqrt(2.0)*pi;? xn2=A*exp(-a*n2*T2).*sin(w0*n2*T2);???? Xk2=fft(xn2,1024);?? subplot(324);stem(n2,xn2,.);grid?on title(?300Hz采样x(n2));xlabel(n);ylabel(x(n2));? k=0:1023;wk=2*k/1024;subplot(323);? plot(wk,abs(Xk2));grid?on title((a)FT[x(n2)]);?xlabel(\omega/\pi);ylabel(|X(e^j^\omega)|); %用200采样频率 ? T3=1/200;?n3=0:64/(1000*T3);?A=444.128;?a=50*sqrt(2.0);? w0=50*sqrt(2.0)*pi;xn3=A*exp(-a*n3*T3).*sin(w0*n3*T3);?????? Xk3=fft(xn3,1024);?subplot(326);stem(n3,xn3,.);grid?on title(200Hz采样?x(n3));? xlabel(n);ylabel(x(n3));k=0:1023;wk=2*k/1024;? subplot(325);plot(wk,abs(Xk3));grid?on title((a)FT[x(n3)]);xlabel(\omega/\pi);ylabel(|X(e^j^\omega)|); 可知200Hz的图形失真最严重。原因是因为采样频率必须大于两倍的fc，如果小于会真。且当NM时发生混叠失真。1000Hz采样效果最好。 2,离散信号、系统和系统响应分析 n=1:50;? hb=zeros(1,50);? hb(1)=1;hb(2)=2.5;? hb(3)=2.5;hb(4)=1;? close?all;?subplot(3,1,1);??stem(hb);? title(系统?h_{b}(n));? m=1:50;??T=0.001; A=444.128;? a=50*sqrt(2.0)*pi;? w0=50*sqrt(2.0)*pi;? x=A*exp(-a*m*T).*sin(w0*m*T);? subplot(3,1,2);??stem(x);? title(输入信号?x(n));? y=conv(x,hb);? subplot(3,1,3);? stem(y);? title(输出?y(n)); 图： 3，证明卷积合 n=1:50;? xc=zeros(1,50);? xc(1)=1;xc(2)=1;xc(3)=1;xc(4)=1;xc(5)=1;xc(6)=1;xc(7)=1;xc(8)=1;xc(9)=1;xc(10)= 1;? ha=zeros(1,50);? ha(1)=1;ha(2)=1;ha(3)=1;ha(4)=1;ha(5)=1;ha(6)=1;ha(7)=1;ha(8)=1;ha(9)=1;ha(10)= 1;? yn=conv(xc,ha);? subplot(321);stem(xc,.,r);?title(响应xc(n)的时域波形);? subplot(323);stem(ha,.,r);?title(响应ha(n)的时域波形);? subplot(325);stem(yn,.,g);?title(xc(n)与ha(n)卷积的响应y(n)的时域波形); n2=1