具体函数如下所示，function [wpos,ypos]=bd_asymp(G,w)G1=zpk(G); Gtf=tf(G);if nargin==1,w=freqint2(Gtf.num{1},Gtf.den{1},100);end zer=G1.z{1}; pol=G1.p{1}; gain=G1.k;wpos=[]; pos1=[];for i=1:length(zer);if isreal(zer(i))wpos=[wpos, abs(zer(...