数字信号处理设计实验报告 西电

《数字信号处理设计实验报告 西电》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、数字信号处理设计实验报告西电数字信号处理设计实验报告一、实验目的通过实验学会设计IIR和FIR数字滤波器分离多个信号，并用matlab实现。二、实验内容14用数字信号处理技术实现两个时域重叠信号的分离，及相位检波，设计分离和检波的方法，编写计算机程序，模拟信号处理过程，绘出时域和频域的处理结果。滤波器1s(t)?s1(t)?s2(t)s1(n)a(n)滤波器2s(n)采样s2(n)滤波器3b(n)滤波器144三、程序设计模拟信号的时域波形,频谱Fs=40000;t=0:1/Fs:4;s1=cos(2*pi*30*t).*cos(2*pi*

2、100*t);s2=cos(2*pi*70*t).*cos(2*pi*700*t);st=s1+s2;S1=abs(fftshift(fft(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('幅度');grid

3、on;axis([00.1-11])subplot(322);plot(F,S1);title('s1频谱');xlabel('频率F');ylabel('幅值');gridon;axis([-1000100001])14subplot(323);plot(t,s2);title('s2时域波形');xlabel('时间t');ylabel('幅度');gridon;axis([00.05-11])subplot(324);plot(F,S2)

4、;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(F,ST);title('st频谱');xlabel('频率F'

5、);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)))/8000;SN=abs(f

6、ftshift(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('频率14F');ylabel('幅值');gridon;ax

7、is([-1000100001])subplot(323);stem(t1,s2n);title('s2n时域波形');xlabel('时间t');ylabel('幅度');grid14on;axis([00.025-11])subplot(324);plot(F1,S2N);title('S2N频谱');xlabel('频率F');ylabel('幅值');gridon;axis([-1000100001])subplot(325);stem

8、(t1,sn);title('sn时域波形');xlabel('时间t');ylabel('幅度');gridon;axis([00.025-1