call_caiyang2.m

数字信号处理的matlab实现代码 对电子信息专业同学的学习很有帮助

% 《数字信号处理教程——MATLAB释义与实现》第四章采样演示程序
% 电子工业出版社出版  陈怀琛编著 2004年9月
% caiyang   演示采样定理，时域采样，频谱周期延拓
% 同时演示采样频率小于2*fc时，产生的混叠效应
%
% call_caiyang2  在频域绘制周期延拓波形
function []=call_caiyang2(h)
fs=get(gcbo,'string');
fs=str2num(fs);
axes(h);hold off;
fs=fs/1000;
n=floor(1/fs);
x=[0:0.1:1,0.9:-0.1:0];
k=-10:10;
df=0.01;
fc=0.1;
w=df*k;
for i=0:n
    f1=w+i*fs;
    f2=w-i*fs;
    plot(f1,x);hold on
    plot(f2,x);hold on
    if fs<2*fc
        n1=fs/df;
        x1=(n1-10)*df+i*fs; x3=(n1-10)*df-i*fs;
        x2=10*df+i*fs;  x4=10*df-i*fs; 
        y1=(20-n1)*0.1; y3=(20-n1)*0.1;
        plot([x1,x2],[y1,y1],'color','r','linestyle','-.');
        plot([x3,x4],[y3,y3],'color','r','linestyle','-.');
        hold on
    end
    axis([-1,1,0,1.5]);pause(0.1);hold on
end