clc;?clear;?close?all;?
%仿真數據長度?
L=2000;?
%構造數字基帶信號
I_Data=（randi（2L1）-2）*2+1;?
Q_Data=zeros（L11）;
Signal_Source=I_Data?+?j*Q_Data;?
%載波信號?
Freq_Sample=2400;???????????????????????%采樣率，Hz?
Delta_Freq=60;??????????????????????????%載波頻率?
Time_Sample=1/Freq_Sample;??????????????%采樣間隔
Delta_Phase=rand（1）*2*pi;???????????????%隨機初相，rad?
Carrier=exp（j*（Delta_Freq*Time_Sample*（1:L）+Delta_Phase））;??????%構造載波信號
%調制處理?
Signal_Channel=Signal_Source.*Carrier‘;?

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%?
%以下為鎖相環處理過程?
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%?
%參數清零及初始化
Signal_PLL=zeros（L1）;??????????????????%鎖相環鎖定及穩定后的數據
NCO_Phase?=?zeros（L1）;?????????????????%鎖定的相位
Discriminator_Out=zeros（L1）;???????????%鑒相器輸出
Freq_Control=zeros（L1）;????????????????%頻率控制
PLL_Phase_Part=zeros（L1）;??????????????%鎖相環相位響應函數
PLL_Freq_Part=zeros（L1）;???????????????%鎖相環頻率響應函數
I_PLL?=?zeros（L1）;?
Q_PLL?=?zeros（L1）;?
%環路處理?
C1=0.022013;????????????????????%環路濾波器系數C1
C2=0.00024722;??????????????????%環路濾波器系數C2??

for?i=2:L?
????Signal_PLL（i）=Signal_Channel（i）*exp（-j*mod（NCO_Phase（i-1）2*pi））;???%得到環路濾波器前的相乘器的輸入
????I_PLL（i）=real（Signal_PLL（i））;???????????????????????????????????????%環路濾波器前的相乘器的I路輸入信息數據
????Q_PLL（i）=imag（Signal_PLL（i））;???????????????????????????????????????%環路濾波器前的相乘器的Q路輸入信息數據
????Discriminator_Out（i）=sign（I_PLL（i））*Q_PLL（i）/abs（Signal_PLL（i））;????%鑒相器的輸出誤差電壓信號
????PLL_Phase_Part（i）=Discriminator_Out（i）*C1;??????????????????????????%環路濾波器對鑒相器輸出的誤差電壓信號處理后得到鎖相環相位響應函數
????Freq_Control（i）=PLL_Phase_Part（i）+PLL_Freq_Part（i-1）;???????????????%控制壓控振蕩器的輸出信號頻率
????PLL_Freq_Part（i）=Discriminator_Out（i）*C2+PLL_Freq_Part（i-1）;????????%環路濾波器對鑒相器輸出的誤差電壓信號處理后得到鎖相環頻率響應函數
????NCO_Phase（i）=NCO_Phase（i-1）+Freq_Control（i）;????????????????????????%壓控振蕩器進行相位調整
end?

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%畫圖顯示結果?
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%畫出鎖相環提取的載波和發射載波
figure（1）
plot（cos（NCO_Phase）‘r‘）;grid?on????????%鎖相環提取的載波
hold?on?
plot（real（Carrier））?????????????????????%發射載波
legend（‘鎖相環提取的載波‘‘發射載波‘）

%畫出鎖相環頻率響應曲線和相位響應曲線
figure（2）
subplot（211）?
plot（-PLL_Freq_Part（2:L）*Freq_Sample）;?grid?on;?
title（‘鎖相環頻率響應曲線‘）;?
axis（[1?L?-100?100]）;?
subplot（212）?
plot（PLL_Phase_Part（2:L）*180/pi）;?
title（‘鎖相環相位響應曲線‘）;?
axis（[1?L?-2?2]）;?grid?on;?

%畫出調制信號（調制信號）和解調后輸出的信號
figure（3）
%設定顯示范圍?
Show_D=300;?%起始位置?
Show_U=350;?%終止位置?
Show_Length=Show_U-Show_D;?
subplot（221）?
plot（I_Data（Show_D:Show_U））;?grid?on;?
title（‘I路信息數據（調制信號）‘）;?
axis（[1?Show_Length?-2?2]）;?
subplot（222）?
plot（Q_Data（Show_D:Show_U））;?grid?on;?
title（‘Q路信息數據‘）;?
axis（[1?Show_Length?-2?2]）;?
subplot（223）?
plot（I_PLL（Show_D:Show_U））;?grid?on;?
title（‘鎖相環輸出I路信息數據（解調信號）‘）;?
axis（[1?Show_Length?-2?2]）;?
subplot（224）?
plot（Q_PLL（Show_D:Show_U））;?grid?on;?
title（‘鎖相環輸出Q路信息數據‘）;?
axis（[1?Show_Length?-2?2]）;?

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2020-01-05?16:40??基于Matlab的載波同步建模與仿真（科斯塔斯環）\
?????文件????????3349??2019-07-15?00:33??基于Matlab的載波同步建模與仿真（科斯塔斯環）\costas_2PSK.m
?????文件????????4618??2019-07-15?00:00??基于Matlab的載波同步建模與仿真（科斯塔斯環）\costas_DSB.m
?????目錄???????????0??2020-01-05?16:39??基于Matlab的載波同步建模與仿真（科斯塔斯環）\仿真圖\
?????目錄???????????0??2020-01-05?16:39??基于Matlab的載波同步建模與仿真（科斯塔斯環）\仿真圖\2PSK信號載波同步\
?????文件???????55048??2019-07-14?22:16??基于Matlab的載波同步建模與仿真（科斯塔斯環）\仿真圖\2PSK信號載波同步\IQ信號輸入輸出.fig
?????文件??????118341??2019-07-11?15:25??基于Matlab的載波同步建模與仿真（科斯塔斯環）\仿真圖\2PSK信號載波同步\載波提取.fig
?????文件??????114958??2019-07-14?22:16??基于Matlab的載波同步建模與仿真（科斯塔斯環）\仿真圖\2PSK信號載波同步\鎖相環頻率_相位響應.fig
?????目錄???????????0??2020-01-05?16:39??基于Matlab的載波同步建模與仿真（科斯塔斯環）\仿真圖\DSB信號載波同步\
?????文件??????561381??2019-07-14?23:57??基于Matlab的載波同步建模與仿真（科斯塔斯環）\仿真圖\DSB信號載波同步\I路Q路輸入輸出.fig
?????文件??????289331??2019-07-14?23:58??基于Matlab的載波同步建模與仿真（科斯塔斯環）\仿真圖\DSB信號載波同步\提取載波_發射載波.fig
?????文件??????340801??2019-07-14?23:58??基于Matlab的載波同步建模與仿真（科斯塔斯環）\仿真圖\DSB信號載波同步\鎖相環頻率_相位響應.fig