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

%--調(diào)制信號
Am?=1;?%--歸一化幅值
fm?=0.1e6;?%--調(diào)制信號的頻率?這里取?fm?=?0.1MHz
Bm?=?0;%--帶寬這里取為單頻信號，所以?Bm=0
%-----------------------
%--調(diào)制信號表達(dá)式
%----------------------
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
%-----------------------
%--載波信號表達(dá)式
%----------------------
sc?=?Ac*cos（2*pi*fc*t）;
sc1?=??Ac*sin（2*pi*fc*t）;

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

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

%--雙邊帶調(diào)制：抑制載波?%--?點(diǎn)乘：.*??兩個相等長度的矢量對應(yīng)點(diǎn)相乘
%?s_am_DSB?=?mf*Am*cos（2*pi*fm*t）.*（Ac*cos（2*pi*fc*t））;
%----------------------------
%--雙邊帶調(diào)制?調(diào)幅波表達(dá)式
%----------------------------
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環(huán)載波恢復(fù)處理過程
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%載波同步參數(shù)清零和初始化
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;????????????????????%環(huán)路阻尼系數(shù)
BL=0.98*symbol_rate;????????????%環(huán)路等效噪聲帶寬
Wn=8*sigma*BL/（1+4*sigma^2）;????%環(huán)路自由震蕩角頻率
T_nco=Ts*decimator;?????????????%壓控振蕩器NCO頻率字更新周期
K1=（2*sigma*Wn*T_nco）/（K）;??????%環(huán)路濾波器系數(shù)K1
K2=（（T_nco*Wn）^2）/（K）;??????????%環(huán)路濾波器系數(shù)K2??%環(huán)路系數(shù)

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