ch3example23prg1.m

通信程序源代码和模型通信程序源代码和模型通信程序源代码和模型

% ch3example23prg1.m
kc=150e3;              % Hz/V  VCO控制灵敏度
omega_n=2*pi*15e3/2.5; % PLL自然角频率
K=2*pi*(0.5*kc);       % 估算环路增益
zeta=1;                % 临界阻尼
tau_1=K/((omega_n).^2);
tau_2=2*zeta/omega_n-1/K;
freq=0:10:100e3;  % 计算频率范围0到100KHz
s=j*2*pi*freq;
G_s=(1+tau_2*s)./(1+tau_1*s);       % 环路滤波器传递函数
figure(1);semilogx(freq,(abs(G_s)));% 作出环路滤波器的频率响应
title('环路滤波器幅频响应');xlabel('Hz');ylabel('|G(s)|');
grid on;
b=[tau_2,1]; % 环路滤波器分子系数向量
a=[tau_1,1]; % 环路滤波器分母系数向量
H_s=(G_s*K./s)./(1+G_s*K./s);
figure(2);semilogx(freq,20*log10(abs(H_s)));% 作出闭环频率响应
title('PLL线性相位模型闭环频率响应');xlabel('Hz');ylabel('20log|H(s)|(dB)');
grid on;