ts=1.e-7;
t=-0.01:ts:0.01;
fm=1000;
fc=100000;??%frequency?of?carrier
m_sig=square（2*pi*fm*t50）;?
figure（1）;
h1=plot（tm_sig）;
xlabel（‘t（sec）‘）;ylabel（‘m（t）‘）;
axis（[-0.005?0.005?-1.5?1.5]）;
title（‘調制信號時域波形‘）;
set（h1‘linewidth‘2）;
grid?on;
%?FFT?of?m_sig
Lfft=length（t）;%the?length?of?t
Lfft=2^ceil（log2（Lfft））;%?shift?the?variable?of?‘Lfft’
M_fre=fftshift（fft（m_sigLfft）/（length（t）））;?
freqm=（-Lfft/2:Lfft/2-1）/（Lfft*ts）;%frequency?range
figure（2）;
h2=plot（freqmabs（M_fre））;
xlabel（‘f（Hz）‘）;ylabel（‘M（f）‘）;
title（‘調制信號頻譜波形‘）;
axis（[1000?1400?0?0.8]）;
set（h2‘Linewidth‘2）;
grid?on;
s_carr1=cos（2*pi*fc*t）;%carrier?signal?1
s_carr2=sin（2*pi*fc*t）;%carrier?signal?2
figure（3）;
plot（ts_carr1）;
xlabel（‘t（sec）‘）;ylabel（‘s_carr1（t）‘）;
axis（[-0.00005?0.00005?-1.5?1.5]）;
title（‘余弦載波時域波形‘）;
grid;
Figure（4）;
plot（ts_carr2）;
axis（[-0.00005?0.00005?-1.5?1.5]）;
xlabel（‘t（sec）‘）;ylabel（‘s_carr2（t）‘）;
title（‘正弦載波時域波形‘）;
grid;
S_carr1=fftshift（fft（s_carr1Lfft）/（length（ts）））;
S_carr2=fftshift（fft（s_carr2Lfft）/（length（ts）））;?
Figure（5）;
plot（freqmabs（S_carr1））;
axis（[-200000?200000?0?15000]）;
xlabel（‘f（hz）‘）;ylabel（‘S_carr1（f）‘）;
title（‘余弦載波頻譜波形‘）;grid;
Figure（6）;
plot（freqmabs（S_carr2））;
axis（[-200000?200000?0?15000]）;
xlabel（‘f（hz）‘）;ylabel（‘S_carr2（f）‘）;
title（‘正弦載波頻譜波形‘）;grid;
s_carr=cos（2*pi*fc*t）;
S_carr=fftshift（fft（s_carrLfft）/（length（ts）））;?
figure（7）;
plot（ts_carr2）;
xlabel（‘t（