appl6_2.m

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

%Appl6_2
load kerl030916ud.txt
Fs=50;    %采样频率
wp=1.5*2/Fs;ws=2.5*2/Fs; %根据采样频率将滤波器边界频率进行转换
Rp=1;Rs=30;    %通带波纹和阻带衰减
Nn=128;
[N,Wn]=buttord(wp,ws,Rp,Rs);  %求得滤波器的最小阶数和截止频率（归一化频率）
[b,a]=butter(N,Wn);  %设计Butterworth低通数字滤波器
figure(1)
[H,f]=freqz(b,a,Nn,Fs);    %用Nn点求得频率特性
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;
dt=1/Fs;   %采样间隔
x=kerl030916ud'; %输入信号
t=[0:(length(x)-1)]*dt;  %时间序列
figure(2)
subplot(2,1,1),plot(t,x); title('输入信号')  %绘制输入信号
y=filter(b,a,x);   %对输入信号进行滤波
subplot(2,1,2),plot(t,y),title('输出信号') %绘制输出信号
xlabel('时间/s')