appl6_4.m

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

%Appl6_4
Fcp=4;Fcs=6;   %设置边界频率为4Hz，6Hz
load iir_signaldata.txt; %调入数据文件
Xt=iir_signaldata; %原始波形数据
Fs=50; %采样频率50Hz
wp=Fcp*2/Fs;ws=Fcs*2/Fs; %将通带临界频率Fcp和阻带临界频率Fcs转换为归一化频率
Rp=1;Rs=25; %通带衰减和阻带衰减
Nn=128; %绘频谱图所用点数
[N,Wn]=buttord(wp,ws,Rp,Rs);  %求滤波器的最小阶数和临界频率（归一化频率）
[b,a]=butter(N,Wn);  %求Butterworth数字滤波器传递函数的分子b和分母a
figure(1); %图形(一)
[H,f]=freqz(b,a,Nn,Fs); %滤波器特性图
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;
n=0:length(Xt)-1; t=n/Fs; %转换为时间序列
figure(2); %图形(二)
subplot(2,1,1),plot(t,Xt);title('输入信号'); %绘制输入信号
xlabel('时间/s');ylabel('振幅'); %坐标轴标识
Yt=filter(b,a,Xt); %对输入信号进行滤波
subplot(2,1,2),plot(t,Yt),title('输出信号');%绘制输出信号
xlabel('时间/s');ylabel('振幅'); %坐标轴标识