信号调制频率为5MHz调制方式为2PSK相移为0信道噪声为高斯白噪声请给出用PLL稳定信号频率和相位的matlab程序。你还需要什么信息

根据提供的信息，可以编写以下程序：

% 参数设置
fs = 100e6;     % 采样率
fc = 5e6;       % 载波频率
T = 1/fs;       % 采样间隔
N = 10000;      % 采样点数
SNR = 10;       % 信噪比

% 生成2PSK调制信号
bits = randi([0,1],1,N);
symbols = 2*bits - 1;   % 映射到-1和1
t = (0:N-1)*T;
s = sqrt(2/T)*real(symbols.*exp(1i*2*pi*fc*t));   % 2PSK调制

% 加入高斯白噪声
noise = awgn(s,SNR);

% PLL稳定频率和相位
f0 = fc;        % 初始频率
theta0 = 0;     % 初始相位
Kp = 1;         % 比例增益
Ki = 0.1;       % 积分增益
Kd = 0;         % 微分增益
alpha = 0.01;   % 低通滤波器带宽
beta = alpha/10;% 带通滤波器带宽

% 初始化PLL变量
theta_hat = theta0;
f_hat = f0;
phi_hat = 0;
e = zeros(1,N);
theta = zeros(1,N);
f = zeros(1,N);
phi = zeros(1,N);

% PLL算法
for n = 2:N
    % 相频检测器
    delta_phi = angle(conj(exp(1i*phi_hat))*noise(n));
    
    % 低通滤波器
    e(n) = Kp*delta_phi + Ki*T*sum(e(1:n-1));
    e_lp = filter(alpha,[1 alpha-1],e(n));
    
    % 频率估计器
    f_hat = f_hat + Kd/T*delta_phi + Ki*beta/T*e_lp;
    
    % 相位估计器
    phi_hat = phi_hat + 2*pi*(f_hat/fs) + Kp*delta_phi + Ki*beta*e_lp;
    
    % 记录PLL变量
    theta_hat = theta_hat + 2*pi*(f_hat/fs);
    theta(n) = theta_hat;
    f(n) = f_hat;
    phi(n) = phi_hat;
end

% 绘图
figure;
subplot(3,1,1);
plot(t,symbols);
title('调制信号');
subplot(3,1,2);
plot(t,noise);
title('加噪声后的信号');
subplot(3,1,3);
plot(t,real(exp(1i*phi)));
hold on;
plot(t,real(exp(1i*phi_hat)));
title('PLL稳定相位');
legend('原始相位','PLL估计相位');

需要注意的是，PLL算法中的比例增益、积分增益、微分增益和低通滤波器带宽等参数需要根据具体情况进行调整。此外，程序中的SNR也可以根据需要进行调整