samp6_11.m

数字信号处理的Matlab实现的随书光盘

%Samp6_11
Fs=1000;   %采样频率
ws=[100 200]*2/Fs;wp=[50 250]*2/Fs;   %通带和阻带边界频率（归一化频率）
Rp=1;Rs=50;Nn=512; %通带波纹和阻带衰减以及绘制频率特性的数据点数
[N,Wn]=ellipord(wp,ws,Rp,Rs);  %求取数字滤波器的最小阶数和归一化截止频率
[b,a]=ellip(N,Rp,Rs,Wn,'stop'); %按最小阶数、截止频率、通带波纹和阻带衰减设计滤波器
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=50;f2=200;f3=270;   %输入信号的频率成分
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)+0.5*sin(2*pi*f3*t); %输入信号
subplot(2,1,1),plot(t,x),title('输入信号') %绘制输入信号
y=filtfilt(b,a,x);  %对输入信号滤波
subplot(2,1,2),plot(t,y) %绘制输出信号
title('输出信号'),xlabel('时间/s')