资源描述：
《数字信号处理设计实验报告》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、数字信号处理设计实验报告一、实验目的通过实验学会设计IIR和FIR数字滤波器分离多个信号，并用matlab实现。二、实验内容用数字信号处理技术实现两个时域重叠信号的分离，及相位检波，设计分离和检波的方法，编写计算机程序，模拟信号处理过程，绘出时域和频域的处理结果。滤波器2滤波器1采样滤波器4滤波器3三、程序设计模拟信号的时域波形,频谱Fs=40000;t=0:1/Fs:4;s1=cos(2*pi*30*t).*cos(2*pi*100*t);s2=cos(2*pi*70*t).*cos(2*pi*700*t);st=s1+s2;S1=abs(fftshift(fft(
2、s1)))/80000;S2=abs(fftshift(fft(s2)))/80000;ST=abs(fftshift(fft(st)))/80000;F=(-80000:80000)*0.25figure(1)subplot(321);plot(t,s1);title('s1时域波形');xlabel('时间t');ylabel('幅度');gridon;axis([00.1-11])subplot(322);plot(F,S1);title('s1频谱');xlabel('频率F');ylabel('幅值');gridon;axis([-1000100001])s
3、ubplot(323);plot(t,s2);title('s2时域波形');xlabel('时间t');ylabel('幅度');gridon;axis([00.05-11])subplot(324);plot(F,S2);title('s2频谱');xlabel('频率F');ylabel('幅值');gridon;axis([-1000100001])subplot(325);plot(t,st);title('st时域波形');xlabel('时间t');ylabel('幅度');gridon;axis([00.05-11])subplot(326);plot
4、(F,ST);title('st频谱');xlabel('频率F');ylabel('幅值');gridon;axis([-1000100001])采样信号的时域波形,频谱Fs1=4000;t1=0:1/Fs1:4;N=0:length(t1)-1s1n=cos(2*pi*30*N/Fs1).*cos(2*pi*100*N/Fs1);s2n=cos(2*pi*70*N/Fs1).*cos(2*pi*700*N/Fs1);sn=s1n+s2n;S1N=abs(fftshift(fft(s1n)))/8000;S2N=abs(fftshift(fft(s2n)))/800
5、0;SN=abs(fftshift(fft(sn)))/8000;F1=(-8000:8000)*0.25figure(2)subplot(321);stem(t1,s1n);title('s1n时域波形');xlabel('时间t');ylabel('幅度');gridon;axis([00.05-11])subplot(322);plot(F1,S1N);title('S1N频谱');xlabel('频率F');ylabel('幅值');gridon;axis([-1000100001])subplot(323);stem(t1,s2n);title('s2n时域
6、波形');xlabel('时间t');ylabel('幅度');gridon;axis([00.025-11])subplot(324);plot(F1,S2N);title('S2N频谱');xlabel('频率F');ylabel('幅值');gridon;axis([-1000100001])subplot(325);stem(t1,sn);title('sn时域波形');xlabel('时间t');ylabel('幅度');gridon;axis([00.025-11])subplot(326);plot(F1,SN);title('SN频谱');xlabel
7、('频率F');ylabel('幅值');gridon;axis([-1000100001])通过前级滤波器的波形fp1=300;fs1=400;Rp=1;Rs=40Wp1=2*fp1/Fs1;Ws1=2*fs1/Fs1;%%滤波器1[M1,Wc1]=buttord(Wp1,Ws1,Rp,Rs);[Bz1,Az1]=butter(M1,Wc1,'low');fp2=500;fs2=400;Rp=1;Rs=40Wp2=2*fp2/Fs1;Ws2=2*fs2/Fs1;%%滤波器2[M2,Wc2]=buttord(Wp2,Ws2,Rp,Rs);[Bz2,A
