ms15_5.m

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

f0=1;T0=1/f0;           %输入基波的频率.周期
fm=5*f0;Tm=1/fm;        %最高频率为基波的5倍频
t=-2:0.1:2;
f=sin(2*pi*f0*t)+1/3*sin(6*pi*f0*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三种抽样信号进行滤波
N=6;rp=1;as=20;
for i=1:3;
    fs=i*fm;Ts=1/fs;     %确定采样频率和周期      
    n=-2:Ts:2;      
    f=sin(2*pi*f0*n)+1/3*sin(6*pi*f0*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