% 《數字信號處理教程——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