三阶锁相环 matlab,matlab仿真实现三阶锁相环路跟踪高动态多普勒变化出现问题...

[/color]

[color=#323456]clc;

clear all;

Chip_rate = 6.3*10^6;          % 扩频码速率(chips/s)

Spreadcode_len = 63;           % 扩频码长度(chips)

Data_len = 240;                % 数据长度(symbols)

Freq_carrier = 4*10^10;        % 射频载波频率(Hz)

Freq_sample = 6.3*10^6;        % 接收机采样频率(Hz)

Symbol_rate = 1*10^5;          % 符号速率(sym/s)

SNR = -5;                     % 信噪比(dB)

load GoldSequence.mat;

PNcode = 2*GoldSequence(1:Data_len,1:63)-1;       % 伪码长度：数据长度*63

n_chip = 63;

g = 9.8;                       % 重力加速度(m/s^2)

c = 3*10^8;                    % 光速(m/s)

ts = 1/1890;                    % JPL的采样时间(s),所设置的值与扩码序列的长度有关

V0 = -500;                     % 初始速度(m/s)

Theta0 = 2*pi*randn(1);        % 初始相位(rad)

Time = [3 0.5 2 0.5 2];

Total_le