阶变系统的开环传递函数.doc

《阶变系统的开环传递函数.doc》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、阶变系统的开环传递函数阶变系统的开环传递函数clear此资料由网络收集而来，如有侵权请告知上传者立即删除。资料共分享，我们负责传递知识。all;Ap=1.68e-2;In=0.03;ps=4e6;pL=2*ps/3;Ki=188.6;Vt=2.873e-3;Kf=1;bate=6900e5;m=35000;Wh=sqrt(4*bate*Ap/(m*Vt))zuni1=0.3;sys1=tf(1/Ap,[1/Wh2*zuni1/Wh10])Wsv=157;zuni2=0.7;Ksv=1.96e-3;sys2=tf(Ksv,[1/Wsv2*zun

2、i1/Wsv此资料由网络收集而来，如有侵权请告知上传者立即删除。资料共分享，我们负责传递知识。1])%系统的开环传递函数sys_open=Ki*sys1*sys2sysclose=feedback(sys_open,1);figure;%绘制nyquist曲线subplot(121);pzmap(sys_open);gridon;xlabel(‘实轴’);ylabel(‘虚轴’);title(‘零极点图’);此资料由网络收集而来，如有侵权请告知上传者立即删除。资料共分享，我们负责传递知识。subplot(122);nyquist(sys_open)

3、;gridon;xlabel(‘实轴’);ylabel(‘虚轴’);title(‘Nyquist图’);figure;%时域分析subplot(121);step(sysclose);gridon;xlabel(‘时间’);ylabel(‘振幅’);title(‘阶跃响应’);subplot(122);impulse(sysclose);gridon;xlabel(‘时间’);ylabel(‘振幅’);title(‘脉冲图响应’);figure;此资料由网络收集而来，如有侵权请告知上传者立即删除。资料共分享，我们负责传递知识。%绘制Bode图及其参

4、数求解w=logspace(-1,2);gridon;margin(sys_open);xlabel(‘频率’);title(‘Bode图’);[mag,phase,W]=bode(sysclose,w);[l,c]=size(mag);mag1=zeros(c,l);for此资料由网络收集而来，如有侵权请告知上传者立即删除。资料共分享，我们负责传递知识。i=1:cmag1(i)=10*log10(mag(1,1,i));enddisp(‘-3dB频率：’);W_3dB=abs(interp1(mag1,W,-3,’spline’));gridon

5、;Kv=Ki*Ksv/Ap;ef=0.02*0.003/Ki*Kf;vm=2.2e-2;er=vm/Kv;e=ef+er;[W_3dBe]此资料由网络收集而来，如有侵权请告知上传者立即删除。资料共分享，我们负责传递知识。