funtion sp=pam(t,fc,fp,fs,tao,pha)
%該程序是用來產生脈沖調制信號（PAM）的
%參數t是所要產生脈沖幅度調制信號的時間，單位為s
%參數fc為脈內信號的頻率，單位為Hz
%參數fp為脈沖信號的頻率，單位為H
%參數fs為采樣時鐘頻率，單位為Hz
%參數tao為脈沖信號的占空比，單位為%
%參數pta為信號的初始相位，單位為rad
if nargin <=6  %    如果不輸入tao 和pha,
    tao=50;
        pha=0;
    end;
 n=0:1/fs:1/fp;
 tn=0:1/fs:t;
 m=t/(1fp);
 spt=(square(2*pi*fp*n)+1)/2;
 st=cos(2*pi*fc*n);
 sq1=st.*spt;
 sp=repmat(sq1,1,m);