逆变器的电感电路双环PI控制,最后matlab求出相位裕度103度,

[复制链接]
3352|2
 楼主| snakeemail 发表于 2008-1-18 23:31 | 显示全部楼层 |阅读模式
 楼主| snakeemail 发表于 2008-1-21 16:04 | 显示全部楼层

求电流环的Kip

%%%%%%%%%%%&nbsp;求电流环的Kip&nbsp;%%%%%%%%%%%%%%%%%%%%%%%%<br />PI&nbsp;=&nbsp;3.1415926;<br />Vout&nbsp;=&nbsp;45;<br />Iout&nbsp;=&nbsp;0.7;<br />Alpha&nbsp;=&nbsp;(3.225/2)&nbsp;/&nbsp;Iout;&nbsp;%&nbsp;输出&nbsp;=&nbsp;Alpha&nbsp;*&nbsp;输入<br />Beta&nbsp;=&nbsp;(3.225/2)&nbsp;/&nbsp;Vout;<br />Kpwm&nbsp;=&nbsp;48;<br />R&nbsp;=&nbsp;70;<br />C&nbsp;=&nbsp;10.8&nbsp;*&nbsp;10^(-6);<br />L&nbsp;=&nbsp;1.2&nbsp;*&nbsp;10^(-3);<br />fn&nbsp;=&nbsp;1&nbsp;/&nbsp;(2&nbsp;*&nbsp;PI&nbsp;*&nbsp;((L&nbsp;*&nbsp;C)&nbsp;^&nbsp;0.5));<br />%fc&nbsp;=&nbsp;fn&nbsp;/&nbsp;16;&nbsp;%穿越频率<br />fc&nbsp;=&nbsp;(1/10)&nbsp;*20000;<br />s&nbsp;=&nbsp;2&nbsp;*&nbsp;PI&nbsp;*&nbsp;fc&nbsp;*&nbsp;i;<br />%求得电流环的Kip<br />Gio&nbsp;=&nbsp;Alpha&nbsp;*&nbsp;Kpwm&nbsp;*&nbsp;(R*C*s&nbsp;+&nbsp;1)&nbsp;/(R*C*L*s^2&nbsp;+&nbsp;L*s&nbsp;+&nbsp;R);<br />Gio&nbsp;=&nbsp;abs(Gio);<br />Kip&nbsp;=&nbsp;1/Gio<br />%电流环开环单位阶越响应<br />SysGio&nbsp;=&nbsp;tf([Gio*Alpha*Kip*Kpwm*&nbsp;R*C&nbsp;Alpha*Kip*Kpwm],&nbsp;[R*L*C&nbsp;L&nbsp;R]);<br />step(SysGio)<br />%margin(SysGio);<br />%&nbsp;求电流环闭环传递函数<br />SysGI&nbsp;=&nbsp;tf([Kip*Kpwm*R*C&nbsp;Kip*Kpwm],&nbsp;[R*C*L&nbsp;L&nbsp;R]);<br />SysHI&nbsp;=&nbsp;tf([Alpha],&nbsp;[1]);<br />SysGic&nbsp;=&nbsp;feedback(SysGI,&nbsp;SysHI,&nbsp;-1)%电流闭环函数<br />%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<br />Kip&nbsp;=<br />0.0708<br />电流闭环函数为:<br />Transfer&nbsp;function:<br />0.00257&nbsp;s&nbsp;+&nbsp;3.399<br />-----------------------------------<br />9.072e-007&nbsp;s^2&nbsp;+&nbsp;0.007119&nbsp;s&nbsp;+&nbsp;77.83
 楼主| snakeemail 发表于 2008-1-21 21:12 | 显示全部楼层

电压环Kvi, Kvp参数求解

%电压环Kvi,&nbsp;Kvp参数求解<br /><br />Fvcut&nbsp;=&nbsp;250;&nbsp;%电压环的穿越频率<br /><br />zeta&nbsp;=&nbsp;2*PI*&nbsp;Fvcut*i;<br /><br />Fvzero&nbsp;=&nbsp;300;&nbsp;%电压环PI的零点<br /><br />&nbsp;<br /><br />GvNum&nbsp;=&nbsp;(zeta&nbsp;+&nbsp;2*PI*Fvzero)&nbsp;*R&nbsp;*&nbsp;Beta&nbsp;*&nbsp;(0.00257&nbsp;*zeta&nbsp;+&nbsp;3.399);<br /><br />GvDen&nbsp;=&nbsp;zeta&nbsp;*&nbsp;(R*C*zeta&nbsp;+1)&nbsp;*&nbsp;(9.072e-007&nbsp;*zeta^2&nbsp;+&nbsp;&nbsp;0.007119&nbsp;*zeta&nbsp;+&nbsp;77.83);<br /><br />Gv&nbsp;=&nbsp;abs(GvNum/GvDen);<br /><br />Kvp&nbsp;=&nbsp;1/&nbsp;Gv<br /><br />Kvi&nbsp;=&nbsp;2&nbsp;*PI*&nbsp;Fvzero&nbsp;*Kvp<br /><br />&nbsp;<br /><br />%系统的开环伯德图<br /><br />Gvpi&nbsp;=&nbsp;tf([Kvp&nbsp;Kvi],&nbsp;[1&nbsp;0]);<br /><br />Gr&nbsp;=&nbsp;tf([Beta&nbsp;*R],&nbsp;[R*C&nbsp;1]);<br /><br />%Gbeta&nbsp;=&nbsp;tf([Beta],&nbsp;[1]);<br /><br />Gvo&nbsp;=&nbsp;series(Gvpi,&nbsp;series(SysGic,&nbsp;Gr))<br /><br />margin(Gvo)
您需要登录后才可以回帖 登录 | 注册

本版积分规则

547

主题

1909

帖子

8

粉丝
快速回复 在线客服 返回列表 返回顶部