samp5_18.m

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

%Samp5_18
wp=1500*2*pi;ws=1000*2*pi;Rp=1;Rs=30;   %滤波器设计参数
[N,Wn]=ellipord(wp,ws,Rp,Rs,'s');   %求滤波器的最小阶数和截止频率
[b,a]=ellip(N,Rp,Rs,Wn,'high','s');  %设计高通椭圆滤波器
w=linspace(1,3000,1000)*2*pi;   %给出计算复数频率响应的频率点
H=freqs(b,a,w);   %计算给定频率的复数频率响应
magH=abs(H);phaH=unwrap(angle(H));  %求幅频响应和相频响应
figure(1)
plot(w/(2*pi),20*log10(magH));  %以频率为横轴绘幅频响应
xlabel('频率/Hz');ylabel('振幅/dB');
title('椭圆模拟高通滤波器');
hold on;plot([1000 1000],ylim, 'r'); %阻带边界
grid on
figure(2)   %给出另一个图形窗口
dt=1/10000;   %采样间隔
f1=400;f2=1600;%信号中所含频率成分
t=0:dt:0.04;%时间序列
x=sin(2*pi*f1*t)+0.5*cos(2*pi*f2*t);  %输入信号
H=[tf(b,a)];   %滤波器在MATLAB系统中的表示
[y,t1]=lsim(H,x,t);  %模拟输出
subplot(2,1,1),plot(t,x),title('输入信号') %绘输入信号
subplot(2,1,2),plot(t1,y)  %绘制输出信号
title('输出信号'),xlabel('时间/s')