ms15r240.m

《数字信号处理实验(MATLAB版)》这

%用模拟滤波器重建信号
dt=0.1;
t=-4:dt:4;
f=sinc(t); 
subplot(4,1,1),plot(t,f,'k');grid
title('原信号与重建的信号');
axis([min(t) max(t) 1.1*min(f) 1.1*max(f)]);
%对fs<2fm.fs=2fm.fs>2fm三种抽样信号进行滤波
fm=1;N=6;rp=1;as=20;
for i=1:3;
    fs=i*fm;Ts=1/fs;  %确定采样频率和周期      
    n=-4:Ts:4;      
    f=sinc(n);              %生成抽样信号
    %设计契比雪夫Ⅰ型低通滤波器
    Wp=2*pi*fm;             
    [b,a]=cheby1(N,rp,Wp,'s');
    %[b,a]=butter(N,Wp,'s');
    %[b,a]=cheby2(N,as,Wp,'s');
    y=lsim(b,a,f,n);
    subplot(4,1,i+1),plot(n,y,'k');grid
    axis([min(n) max(n) 1.1*min(y) 1.1*max(y)]);
end