clc;Rs=10*10^6;      %码速率es=0.707;       %阻尼系数fs=8*10^6;     %采样频率、系统时钟频率Bit_Loop_out=27;%环路滤波器输入输出位宽N=31;           %NCO频率字位宽Tdds=8;        %NCO频率字更新所需的系统时钟周期个数K=2*pi*Tdds*2^(Bit_Loop_out-2)/2^N    %环路总增益K=K*90Wn=73.867*10^3;     %环路固有谐振频率(rad/s)c1=((8*es*Wn/fs)+4*Wn*Wn/fs/fs)/(4+4*es*Wn/fs+Wn*Wn/fs/fs)/Kc2=4*Wn*Wn/fs/fs/(4+4*es*Wn/fs+Wn*Wn/fs/fs)/K% c1=2^(-6)% c2=2^(-13)b=[0 K*c1 K*c2-K*c1];a=[1 K*c1-2 1+K*c2-K*c1];roots(a)%画出零极点图zplane(b,a);