《matlab阶跃响应频域求导三种方法拟合检验》由会员分享，可在线阅读，更多相关《matlab阶跃响应频域求导三种方法拟合检验（2页珍藏版）》请在装配图网上搜索。
1、%针对非最小相位系统，读者可自己修改GO,下面其他程序不用修改G0=tf(-4*1 -1.25,conv(5 1,0.4 1)G_step=zhuanhua_step(GO)G_WC=zhuanhua_WC(GO)G_qiudao=zhuanhua_qiudao(GO)figure (3)step(GO,y)hold onstep(G_step,b-)hold onstep(G_WC,r)hold onstep(G_qiudao,g)title(CxiBjiaii6i6)gridse t ( gcf, Pos i t ion ,1 O O 1OO 56O 42O);set (gca, Posi
2、tion ,O O 1 1);figure_FontSize=8; set(get(gca,XLabel),FontSize,figure_FontSize,Vertical,top); set(get(gca,YLabel),FontSize,figure_FontSize,Vertical,middle); set(findobj(FontSize,1O),FontSize,figure_FontSize);set(findobj(get(gca,Children),LineWidth,O.5),LineWidth,2);以下为三个子函数function G_step=zhuanhua_s
3、tep(G)子函数一，用阶跃响应的办法，对G函数进行拟合,x,t);% G=tf(3.O67*O 1,conv(25 1,25 1) y t=step(G);fun=inline(x(1)*(1-exp(-(t-x(2)/x(3) x=lsqcurvefit(fun,1 1 1,t,y);K=dcgain(G);L=x(2);T=x(3); G_step=tf(K,T 1,InputDelay,L);function G_WC=zhuanhua_WC(G)子函数二，用频域的方法Kc,Pm,wc,wcp=margin(G);%uOUp6dIi6|pAuEikey=0;L=1.6*pi/(3*wc
4、);K=dcgain(G);T=0.5*Kc*K*L; if isfinite(Kc),x0=L;T;while ikey=0,u=wc*x0(1);v=wc*x0(2);FF=K*Kc*(cos(u)-v*sin(u)+1+vA2;sin(u)+v*cos(u);J=-K*Kc*wc*sin(u)-K*Kc*wc*cos(u),-K*Kc*wc*sin(u)+2*wc*v;wc*cos(u)-wc*v *sin(u),wc*cos(u);x1=x0-inv(J)*FF;if norm(x1-x0)1e-8,ikey=1;else,x0=x1;end,endL=x0(1); T=x0(2);e
5、ndG_WC=tf(K,T 1,iodelay,L);function G_qiudao=zhuanhua_qiudao(G)% global fun fx1 子函数三，用求导的方法，求出带纯延迟的一阶惯性传递函数模型 global K T Lsyms xxishu.fenzi=G.num1;xishu.fenmu=G.den1; xishu.delay=G.InputDelay;fx=poly2sym(xishu.fenzi)/poly2sym(xishu.fenmu)*exp(-xishu.delay*x); K=subs(fx,0);K1=diff(fx);K2=diff(fx,2);K1_0=subs(K1,0);K2_0=subs(K2,0);TAR=-K1_0/K;T=sqrt(K2_0/K-TAR2);L=TAR-T;G_qiudao=tf(K,T 1,InputDelay,L);
matlab阶跃响应频域求导三种方法拟合检验
