资源描述：
《数字信号处理实验源程序》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、数字信号处理实验源程序实验一源程序%实验一%根据差分方程用filter函数求响应closeall;clearallA=[1,-0.9];B=[0.05,0.05];x1n=[11111111zeros(1,50)];x2n=ones(1,128);hn=impz(B,A,58);subplot(2,2,1);y='h(n)';stem(hn)title('(a)系统单位脉冲响应h(n)')y1n=filter(B,A,x1n);subplot(2,2,2);y='y1(n)';stem(y1n);title('(b)系统对R8(n)的响应y1(n)')y2n=filter(B,A,x
2、2n);subplot(2,2,3);y='y2(n)';stem(y2n);title('(c)系统对u(n)的响应y2(n)')%实验一%利用单位脉冲响应用conv函数求响应closeall;clearallx1n=[11111111];h1n=[ones(1,10)zeros(1,10)];h2n=[12.52.51zeros(1,10)];y21n=conv(h1n,x1n);y22n=conv(h2n,x1n);subplot(2,2,1);y='h21(n)';stem(h1n)title('(a)系统单位脉冲响应h1(n)')subplot(2,2,2);y='y21(
3、n)';stem(y21n);title('(b)h1(n)与R8(n)的卷积y21(n)')subplot(2,2,3);y='h2(n)';stem(h2n);title('(c)系统单位脉冲响应h2(n)')subplot(2,2,4);y='y22(n)';stem(y22n);title('(d)h2(n)与R8(n)的卷积y22(n)')实验二源程序%%实验二%采样频率为1kHzTp=64/1000;Fs=1000;T=1/Fs;M=Tp*Fs;n=0:M-1;A=444.128;alph=pi*50*2^0.5;omega=pi*50*2^0.5;xnt=A*exp(-
4、alph*n*T).*sin(omega*n*T);Xk=T*fft(xnt,M);yn='xa(nT)';subplot(3,2,1);stem(xnt)boxon;title('(a)Fs=1000HZ')k=0:M-1;fk=k/Tp;subplot(3,2,2);plot(fk,abs(Xk));title('(a)T*FT[xa](nT),Fs=1000Hz')xlabel('f(Hz)');ylabel('幅度');axis([0,Fs,0,1.2*max(abs(Xk))])%采样频率为300HzTp=64/1000;Fs=300;T=1/Fs;M=ceil(Tp*Fs
5、);n=0:M-1;A=444.128;alph=pi*50*2^0.5;omega=pi*50*2^0.5;xnt=A*exp(-alph*n*T).*sin(omega*n*T);Xk=T*fft(xnt,M);yn='xa(nT)';subplot(3,2,3);stem(xnt)boxon;title('(a)Fs=300HZ')k=0:M-1;fk=k/Tp;subplot(3,2,4);plot(fk,abs(Xk));title('(a)T*FT[xa](nT),Fs=300Hz')xlabel('f(Hz)');ylabel('幅度');axis([0,Fs,0,1.
6、2*max(abs(Xk))])%采样频率为200HzTp=64/1000;Fs=200;T=1/Fs;M=ceil(Tp*Fs);n=0:M-1;A=444.128;alph=pi*50*2^0.5;omega=pi*50*2^0.5;xnt=A*exp(-alph*n*T).*sin(omega*n*T);Xk=T*fft(xnt,M);yn='xa(nT)';subplot(3,2,5);stem(xnt)boxon;title('(a)Fs=200HZ')k=0:M-1;fk=k/Tp;subplot(3,2,6);plot(fk,abs(Xk));title('(a)T*F
7、T[xa](nT),Fs=200Hz')xlabel('f(Hz)');ylabel('幅度');axis([0,Fs,0,1.2*max(abs(Xk))])实验三源程序%实验三%利用FFT对信号进行频谱分析closeall;clearallx1n=[1111];M=8;xa=1:(M/2);xb=(M/2):-1:1;x2n=[xa,xb];%产生长度为8的三角波序列x2(n)x3n=[xb,xa];X1k8=fft(x1n,8);%计算x1n的8点D
