samp6_6.m

MATLAB7.x数字信号处理 光盘内容

%Samp6_6
wp=0.2*pi;ws=0.3*pi;Rp=1;Rs=15;   %数字滤波器截止频率、通带波纹和阻带衰减
T=0.01;Nn=128;   %采样间隔
Wp=wp/T;Ws=ws/T;   %得到模拟滤波器的频率—采用脉冲响应不变法的频率转换形式
[N,Wn]=buttord(Wp,Ws,Rp,Rs,'s');   %计算模拟滤波器的最小阶数
[z,p,k]=buttap(N);   %设计低通原型数字滤波器
[Bap,Aap]=zp2tf(z,p,k);   %零点极点增益形式转换为传递函数形式
[b,a]=lp2lp(Bap,Aap,Wn);  %低通滤波器频率转换
[bz,az]=impinvar(b,a,1/T);  %脉冲响应不变法设计数字滤波器传递函数
figure(1)
[H,f]=freqz(bz,az,Nn,1/T);  %输出幅频响应和相频响应
subplot(2,1,1),plot(f,20*log10(abs(H)));
xlabel('频率/Hz');ylabel('振幅/dB');grid on;
subplot(2,1,2),plot(f,180/pi*unwrap(angle(H)))
xlabel('频率/Hz');ylabel('相位/^o');grid on;
figure(2)
f1=5;f2=30;   %输入信号含有的频率
N=100;   %数据点数
n=0:N-1;t=n*T;   %时间序列
x=sin(2*pi*f1*t)+0.5*cos(2*pi*f2*t);  %输入信号
subplot(2,1,1),plot(t,x),title('输入信号')
y=filtfilt(bz,az,x);   %对信号进行滤波
subplot(2,1,2),plot(t,y),title('输出信号'),xlabel('时间/s')