[数字信号处理设计实验报告西电]数字信号处理史林.pdf

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

林

...

...

数字信号处理

设计实验报告

一、实验目的通过实验学会设计IIR和FIR数字滤波器分离多个信

号，并用matlab实现。

二、实验容

用数字信号处理技术实现两个时域重叠信号的分离，及相位检波，设

计分离和检波的方法，编写计算机程序，模拟信号处理过程，绘出时域和

频域的处理结果。

滤波器2滤波器1

滤波器2

滤波器1

采样

采样

滤波器4滤波器3

滤波器4

滤波器3

三、程序设计

模拟信号的时域波形,频谱

Fs=__;

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(s1)))/__;

S2=abs(fftshift(fft(s2)))/__;

ST=abs(fftshift(fft(st)))/__;

F=(-__:__)*0.25

figure(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])

subplot(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(F,ST);title(st频谱);

xlabel(频率F);ylabel(幅值);gridon;

axis([-1000100001])

采样信号的时域波形,频谱

Fs1=4000;

t1=0:1/Fs1:4;N=0:length(t1)-1

s1n=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(fftshift(fft(sn)))/8000;

F1=(-8000:8000)*0.25

figure(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时域波形);

xlabel(时间t);ylabel(幅度);gridon;

axis([00.025-11])

subplot(324);

plot(F1,S2N);title(S2N频谱);

xlabel(频率F);ylabel(幅值);gridon;

axis([