二階環路濾波器的matlab 設計

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*90
Wn=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)/K
c2=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);