数字信号处理实验源程序.docVIP

数字信号处理实验源程序 实验一源程序 %实验一 %根据差分方程用filter函数求响应 close all;clear all A=[1,-0.9];B=[0.05,0.05]; x1n=[1 1 1 1 1 1 1 1 zeros(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,x2n); subplot(2,2,3);y=y2(n);stem(y2n); title((c)系统对u(n)的响应y2(n)) %实验一 %利用单位脉冲响应用conv函数求响应 close all;clear all x1n=[1 1 1 1 1 1 1 1 ]; h1n=[ones(1,10) zeros(1,10)]; h2n=[1 2.5 2.5 1 zeros(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(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)) 实验二源程序 %%实验二 %采样频率为1kHz Tp=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(-alph*n*T).*sin(omega*n*T); Xk=T*fft(xnt,M); yn=xa(nT);subplot(3,2,1); stem(xnt) box on;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))]) %采样频率为300Hz Tp=64/1000; Fs=300;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,3); stem(xnt) box on;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.2*max(abs(Xk))]) %采样频率为200Hz Tp=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) box on;title((a)Fs=200HZ) k=0:M-1;fk=k/Tp; subplot(3,2,6);plot(fk,abs(Xk));title((a)T*FT[xa](nT),Fs=200Hz) xlabel(f(Hz));ylabel(幅度);axis([0,Fs,0,1.2*max(abs(Xk))]) 实验三源程序 %实验三 %利用FFT对信号进行频谱分析 close all;clear all x1n=[1 1 1 1]; M=8;xa=1:(M/2);xb=(M/2):-1:1;x2n=[xa,xb]; %