基于Matlab的载波同步建模与仿真科斯塔斯环.zip

基于Matlab的载波同步建模与仿真（科斯塔斯环）——超全代码和图
含对于2PSK数字调制和DSB双边带模拟调制系统的利用科斯塔斯环法进行的载波同步，含载波提取对比图和输出信号图等
代码直接可运行，已包含注释
（设计原理文档报告可参考链接：https://blog.csdn.net/weixin_41683971/article/details/103844756）

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