t0=0.15;    %信號保持時間
ts=0.001;   %采樣時間間隔
fc=250;     %載波頻率
fs=1/ts;    %采樣頻率
df=0.3;     %頻率分辨力
a=0.8;      %調制系數
t=[0:ts:t0];%時間矢量
m=[ones(1,t0/(3*ts)),-2*ones(1,t0/(3*ts)),zeros(1,t0/(3*ts)+1)];  
    %定義信號序列
c=cos(2*pi*fc.*t);                       %載波信號
m_n=m/max(abs(m));                                  %調制信號
[M,m,dfl]=fft_seq(m,ts,df);              %傅立葉變換
M=M/fs;
u=(1+a*m_n).*c;                         %調制信號載波
[U,u,dfl]=fft_seq(u,ts,df);          
U=U/fs;
f=[0:dfl:dfl*(length(m)-1)]-fs/2;        %頻率矢量

clf
subplot(2,2,1);plot(t,m(1:length(t)));   %未調制信號
title('未調制信號');
subplot(2,2,2);plot(t,u(1:length(t)));   %已調制信號
title('已調制信號');
subplot(2,2,3);plot(f,abs(fftshift(M))); %未調制信號頻譜
title('未調制信號頻譜');
subplot(2,2,4);plot(f,abs(fftshift(U))); %已調制信號頻譜
title('已調制信號頻譜');