clc;
clear?all;

SimTime?=?0.1;????%?總的仿真時間（單位：秒）
DispTime?=?0.001;?%?顯示的時間范圍（單位：秒）
fs?=?450000;??????%?系統總的采樣頻率（單位：Hz）
Ts?=?1/fs;????????%?采樣時間間隔
t?=?0:Ts:SimTime;?%?時間矢量
N?=?length（t）;????%?總采樣點數
A=1;??????????????%調制信號的幅度因子
m=0.5;????????????%調制指數

df?=?fs/N;????????%?計算頻率分辨率（事實上df=1/SimTime）
f?=?-fs/2:df:fs/2-df;?%?定義頻率矢量（頻譜圖的橫坐標）

%被調制信號m（t）=sin6000πt
mt?=?sin（6000*pi*t）;??
figure（1）;
plot（tmt）;grid?on;
axis（[0?DispTime?-1.1?1.1]）;
xlabel（‘時間（s）‘）;ylabel（‘電壓值（V）‘）;
title（‘被調制信號m（t）=sin6000πt‘）;

figure（2）;
fmt?=?fft（mt）;
fmt?=?fftshift（fmt）;
fmt?=?abs（fmt）;
plot（ffmt）;grid?on;
minF?=?min（abs（fmt））;
maxF?=?max（abs（fmt））;?
axis（[-4000?4000?minF?maxF]）;
xlabel（‘頻率（Hz）‘）;ylabel（‘頻譜幅度值‘）;
title（‘和信號幅頻特性‘）;

%載波信號c（t）=cos90000πt
ct?=?cos（90000*pi*t）;??
figure（3）;
plot（tct）;grid?on;
axis（[0?DispTime?-1.1?1.1]）;
xlabel（‘時間（s）‘）;ylabel（‘電壓值（V）‘）;
title（‘載波信號m（t）=cos90000πt‘）;

%AM已調制信號Sam（t）=?A*[1+m*m（t）]*c（t）
samt=?A*（1+m*mt）.*ct;??
figure（4）