samp6_10.m

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

%Samp6_10
Fs=1000;   %采样频率
wp=[100 200]*2/Fs;ws=[50 250]*2/Fs; %通带和阻带边界频率（归一化频率）（6-20式）
Rp=1;Rs=30;Nn=128;   %通带波纹和阻带衰减以及绘制频率特性的数据点数
[N,Wn]=cheb2ord(wp,ws,Rp,Rs);  %求取数字滤波器的最小阶数和归一化截止频率
[b,a]=cheby2(N,Rs,Wn);  %按最小阶数截止频率和阻带衰减设计数字滤波器
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;
figure(2)
f1=10;f2=100;   %输入信号的频率成分
N=500;   %输入信号的数据点数
dt=1/Fs;n=0:N-1;t=n*dt;  %时间序列
x=sin(2*pi*f1*t)+0.8*cos(2*pi*f2*t); %输入信号
subplot(2,1,1),plot(t,x),title('输入信号') %绘制输入信号
y=filtfilt(b,a,x);   %对输入信号进行滤波，输出为y.
subplot(2,1,2),plot(t,y)   %绘制输出信号波形
title('输出信号'),xlabel('时间/s')