数字信号处理综合实验11

%产生信号

n=0:999;

x1=sin(0.03*pi*n); %低频信号

x2=sin(0.1*pi*n); %中频信号

x3=sin(0.2*pi*n); %高频信号

x=x1+x2+x3; %混合信号

%信号频域分析

x1f=fft(x1);

x2f=fft(x2);

x3f=fft(x3);

xf=fft(x);

%滤波器设计

%低通滤波器设计

A1p=1;

A1s=20;

w1p=0.04;

w1s=0.08;

[N1 Wn1]=buttord(w1p,w1s,A1p,A1s);

[b1 a1]=butter(N1,Wn1,'low')

%带通滤波器设计

A2p=1;

A2s=20;

w2p1=0.08;

w2p2=0.12;

w2s1=0.06;

w2s2=0.14;

[N2 Wn2]=buttord([w2p1,w2p2],[w2s1,w2s2],A2p,A2s); [b2 a2]=butter(N2,Wn2,'bandpass')

%高通滤波器设计

A3p=1;

A3s=20;

w3p=0.18;

w3s=0.12;

[N3 Wn3]=buttord(w3p,w3s,A3p,A3s);

[b3 a3]=butter(N3,Wn3,'high')

%滤波

x1p=filter(b1,a1,x); %低通滤波

x2p=filter(b2,a2,x); %带通滤波

x3p=filter(b3,a3,x); %高通滤波

%画图

%信号时域图

figure(1)

subplot(4,1,1);plot(x1);axis([0 1000 -3 3]);

title('Signal x1');xlabel('sample');ylabel('magnitude');

相关文档
相关主题
热门文档