ms11r25.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;
    if i<=2 c=0,else c=0.2,end
    fs=(4-i+c)*fm;Ts=1/fs;  %确定采样频率和周期      
    n=-4:Ts:4;      
    f=sinc(n);              %生成抽样信号
    st=length(n);           %求时间轴上采样点数
    %设计契比雪夫Ⅰ型低通滤波器
    Wp=2*pi*fm;             
    dw=(2*Wp)/(st-1);
    w=0:dw:2*Wp;
    [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