【正文】 ,[1 2*pi*kv*kd])] %error transmittance given in the handout[y12 t] = lsim(g12,u,t)figureplot(t,y12)xlabel(39。TIME IN SECONDS39。)ylabel(39。AMPLITUDE39。)title(39。一階鎖相環(huán)誤差傳輸函數(shù)kd*kv=439。)%****************************************************************a = zeta = sqrt((pi*kv*kd)/(2*a))omegan = sqrt(2*pi*kv*kd*a)g21 = [tf([2*zeta*omegan omegan^2],[1 2*zeta*omegan omegan^2])] [y21 t] = lsim(g21,u,t)figureplot(t,y21)xlabel(39。TIME IN SECONDS39。)ylabel(39。AMPLITUDE39。)title(39。二階鎖相環(huán)傳輸函數(shù)kd*kv=4439。)%****************************************************************g22 = [tf([1 0 0],[1 2*zeta*omegan omegan^2])] [y22 t] = lsim(g22,u,t)figureplot(t,y22)xlabel(39。TIME IN SECONDS39。)ylabel(39。AMPLITUDE39。)title(39。二階鎖相環(huán)誤差傳輸函數(shù)kd*kv=439。)%****************************************************************附錄B一種新的快速鎖定電荷泵鎖相環(huán)：分析與設計Vahideh Sadat Sadeghi Hossein MiarNaimi收稿日期：2012年11月5日/修訂：2012年12月16日/接受日期：2012年12月19日/在線發(fā)布時間：2013年1月6日