samp6_12.m

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

%Samp6_12
Order=10;   %滤波器的阶数
f=0:0.1:1;   %归一化频率点
m=[0 0 1 1 0 0  1 1 1 0 0];  %幅度点
[b,a]=yulewalk(Order,f,m);  %设计滤波器
[h,w]=freqz(b,a,128);   %计算128个点的频率特性
figure(1)
plot(f,m,'b-',w/pi,abs(h),'k:');  %绘制理想滤波器和设计滤波器的幅频特性
xlabel('归一化频率');ylabel('振幅');
title('运用yulewalk方法设计IIR滤波器');
legend('理想特性','实际设计',1);  %给定图例
figure(2)
Fs=50;   %信号采样频率
f1=6;f2=17;  %信号的频率成分
N=100; %数据点数
dt=1/Fs;n=0:N-1;t=n*dt; %时间序列
x=sin(2*pi*f1*t)+0.5*cos(2*pi*f2*t);  %输入信号
subplot(2,1,1),plot(t,x),title('输入信号') %绘制输入信号
y=filtfilt(b,a,x);  %对信号进行滤波
subplot(2,1,2),plot(t,y)   %绘制输出信号
title('输出信号'),xlabel('时间/s')