function sp=pam(t,fc,fp,fs,tao,pha)
%產生脈沖幅度調制（PAM）信號
%t是PAM的時間(s)，fc是脈內信號頻率(Hz)，fp是脈沖信號頻率(Hz)
%fs是采樣時鐘頻率（Hz），tao是脈沖信號占空比(%)，pha是信號初始相位(rad)
if nargin<=6          %若不輸入tao和pha，則tao=50，pha=0
    tao=50;
    pha=0;
end;
n=0:1/fs:1/fp;
tn=0:1/fs:t;
m=t/(1/fp);
spt=(square(2*pi*fp*n)+1)/2;
st=cos(2*pi*fc*n);
sp1=st.*spt;
sp=repmat(sp1,1,m);