%?AM_amplitude_modulation_test.m
clc;?
close?all;
clear?all;
%--參數
%--采樣參數
fs?=10e6;?%--數字采樣速率?fs?>=?2（fc+fm+0.5*Bm）?這里取?fs?=?10?MHz
Ts?=?1/fs;
ps?=?1e6;
N_number?=?200;
N=N_number*fs/ps;??%--采樣點個數?N?>?fix（2*fs/fm）;?%--至少一個周期內采兩個點
n?=0:N-1;?%--采樣序列
t?=n/fs;??%--采樣時間序列

%--調制信號
Am?=1;?%--歸一化幅值
fm?=0.1e6;?%--調制信號的頻率?這里取?fm?=?0.1MHz
Bm?=?0;%--帶寬這里取為單頻信號，所以?Bm=0
%-----------------------
%--調制信號表達式
%----------------------
sm?=?Am*cos（2*pi*fm*t）;
%?sm?=?sawtooth（2*pi*10000*t0.5）;??%%三角波
%?sm?=?square（2*pi*10000*t75）;
%--載波信號
Ac?=1;?%--歸一化幅值
fc?=1e6;?%--載波頻率?一般?fc?>?fm?這里取?fc?=?1?MHz
%-----------------------
%--載波信號表達式
%----------------------
sc?=?Ac*cos（2*pi*fc*t）;
sc1?=??Ac*sin（2*pi*fc*t）;

%--調制度
mf?=?0.5;%--mf?取值在?0?和?1?之間.?mf?=?0?表示沒有調制；mf?=1?是過調制的邊界

%--普通幅度調制：載波+雙邊帶?%--?點乘：.*??兩個相等長度的矢量對應點相乘
%?s_am?=?（1+mf*Am*cos（2*pi*fm*t））.*（Ac*cos（2*pi*fc*t））;
%----------------------------
%--普通幅度調制?調幅波表達式
%----------------------------
s_am_general?=?（1+mf*sm）.*sc;?%--

%--雙邊帶調制：抑制載波?%--?點乘：.*??兩個相等長度的矢量對應點相乘
%?s_am_DSB?=?mf*Am*cos（2*pi*fm*t）.*（Ac*cos（2*pi*fc*t））;
%----------------------------
%--雙邊帶調制?調幅波表達式
%----------------------------
s_am_DSB?=?mf*sm.*sc;?%--
s_am_DSB_1?=?mf*sm.*sc1;?%--

Signal_Channel?=?s_am_DSB?+?1i?*?s_am_DSB_1;
T?=?1/ps;
Sample_Length?=?length（t）;
Simulation_Length?=?200;
decimator?=?1;
symbol_rate?=?ps;
I_data1?=?sm;
Q_data1?=?zeros（1length（t））;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%以下為COSTAS環載波恢復處理過程
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%載波同步參數清零和初始化
Signal_PLL=zeros（1Sample_Length）;
Sample_Length1=Simulation_Length*round（T/（Ts*decimator））;
Signal_PLL_D=zeros（1Sample_Length1）;
I_PLL_D=zeros（1Sample_Length1）;
Q_PLL_D=zeros（1Sample_Length1）;
NCO_Phase?=?zeros（1Sample_Length1）;
Discriminator_Out=zeros（1Sample_Length1）;
Freq_Control=zeros（1Sample_Length1）;
PLL_Phase_Part=zeros（1Sample_Length1）;
PLL_Freq_Part=zeros（1Sample_Length1）;

Ko=1;???????????????????????????%壓控振蕩器增益???????????????
Kd=1;???????????????????????????%鑒相器增益??
K=Ko*Kd;
sigma=0.707;????????????????????%環路阻尼系數
BL=0.98*symbol_rate;????????????%環路等效噪聲帶寬
Wn=8*sigma*BL/（1+4*sigma^2）;????%環路自由震蕩角頻率
T_nco=Ts*decimator;?????????????%壓控振蕩器NCO頻率字更新周期
K1=（2*sigma*Wn*T_nco）/（K）;??????%環路濾波器系數K1
K2=（（T_nco*Wn）^2）/（K）;??????????%環路濾波器系數K2??%環路系數

kuai=2*sigma*Wn;
%K1=2*0.