【正文】 （3） 采用MATLAB仿真工具建立了電機(jī)的線性數(shù)學(xué)模型，對(duì)SRG線性數(shù)學(xué)模型進(jìn)行了一系列定性分析，獲得了相應(yīng)結(jié)論。（2） 闡釋了開關(guān)磁阻發(fā)電機(jī)的運(yùn)行和控制機(jī)理，指出了開關(guān)磁阻電機(jī)作為發(fā)電機(jī)所具有的獨(dú)特性。仿真結(jié)果驗(yàn)證了開關(guān)磁阻發(fā)電機(jī)的輸出電壓的控制系統(tǒng)的正確性，同時(shí)也證明了不對(duì)稱功率變換器的正確性、有效性和可行性。詳細(xì)分析了SRG線性模型。 end end end仿真波形如圖下所示。 end end else if(i_m=0) u=s。 else if(deta_im=*i_ref) u=s。 else u=v_q_m。(thetatheta_off)) if(deta_im=*i_ref) u=s。 if ((thetatheta_on)amp。 else V=0。 end end end end end % function V=voltage(v_q_m,v_uc) global s if(v_q_m==s) V=v_uc。 else if (theta) dotL=(L_maxL_min)/theta_rise*180/pi。 else if (theta) dotL=(L_maxL_min)/theta_rise*180/pi。 else L=L_min。 else if (theta) L=L_max。 end end % function L=inductance(theta) global L_max L_min theta_rise if (theta) L=L_min。 end theta_C=theta_A+theta_pole/3。 end figure(1) subplot(411) plot(data_time,data_inductance_a,data_time,data_inductance_b,data_time,data_inductance_c) subplot(412) plot(data_time,data_current_a,data_time,data_current_b,data_time,data_current_c) subplot(413)plot(data_time,data_current_ic,data_time,data_current_dc,data_time,data_current_l) subplot(414) plot(data_time,data_voltage_a,data_time,data_voltage_b,data_time,data_voltage_c) figure(2) subplot(411) plot(data_time,data_torque_a,data_time,data_torque_b,data_time,data_torque_c)subplot(412) plot(data_time,data_voltage_q_a,data_time,data_voltage_q_b,data_time,data_voltage_q_c)subplot(413)plot(data_time,data_voltage_uc) end% function [theta_B,theta_C]=other_theta(theta_A) global theta_pole theta_B=theta_A+theta_pole*2/3。 data_torque_b(num+1,1)=T_b。 % data_current_l(num+1,1)=i_l。 % data_current_dc(num+1,1)=i_dc。 data_voltage_q_c(num+1,1)=v_q_c。 data_voltage_q_a(num+1,1)=v_q_a。 data_voltage_b(num+1,1)=v_b。 data_current_c(num+1,1)=i_c。 data_current_a(num+1,1)=i_a。 data_inductance_b(num+1,1)=L_b。 T_c=*dotinduct(theta_c)*i_c*i_c。% T_a=*dotinduct(theta_a)*i_a*i_a。 v_q_b=zhihuan(i_b,v_q_b,theta_b)。 if(v_uc0) v_uc=0。 dotv_uc=i_ic/C*T_timer。 % i_l=v_uc/R_l。 end [theta_b,theta_c]=other_theta(theta_a)。 end theta_a=theta_a+w_act*T_timer1*180/pi。 i_c=i_c+doti_c*T_timer1。 dotL_c=dotinduct(theta_c)。 if(i_b0) i_b=0。 doti_b=(v_bR_m*i_bdotL_b*w_act*i_b)/L_b。 end L_b=inductance(theta_b)。 i_a=i_a+doti_a*T_timer1。 dotL_a=dotinduct(theta_a)。 end N=0。 i_ref=dotu*k_p+i_2。 i_1=dotu*k_i*T_timer*25。 N=N+1。 v_b=voltage(v_q_b,v_uc)。 data_torque_c(1)=T_c。% data_torque_a(1)=T_a。 data_voltage_q_c(1)=v_q_c。 data_voltage_q_a(1)=v_q_a。 data_voltage_b(1)=v_b。 data_current_c(1)=i_c。% data_current_a(1)=i_a。% data_current_ic(1)=i_ic。 data_inductance_c(1)=L_c。 data_inductance_a(1)=L_a。 data_torque_b=zeros(num_1,1)。% data_voltage_uc=zeros(num_1,1)。 data_voltage_q_b=zeros(num_1,1)。 data_voltage_c=zeros(num_1,1)。 data_voltage_a=zeros(num_1,1)。 data_current_b=zeros(num_1,1)。% data_current_l=zeros(num_1,1)。 data_current_dc=zeros(num_1,1)。 data_inductance_b=zeros(num_1,1)。 num_2=T_timer/T_timer1。 data_time=[0:T_timer:T_sum]39。 T_timer=。 T_c=0。 T_a=*dotinduct(theta_a)*i_a*i_a。 L_b=inductance(theta_b)。 v_q_c=0。% v_q_a=0。 v_c=0。 v_a=0。 i_ic=0。 i_c=0。 i_a=0。 theta_a=35。 n_act=800。 i_2=0。 k_p=。 i_ref=0。 theta_on=。 s=1。 theta_rise=14。 L_min=。 R_l=300。 end N=0。 i_ref=dotu*k_p+i_2。 i_1=dotu*k_i*T_timer*25。其中使用的PI調(diào)節(jié)器的仿真程序?yàn)椋? N=N+1。 v_q_b=zhihuan(i_b,v_q_b,theta_b)。 if(v_uc0) v_uc=0。 dotv_uc=i_ic/C*T_timer。 % i_l=v_uc/R_l。 end endEnd 不對(duì)稱功率變換器仿真模型在不對(duì)稱橋式功率變換器的自勵(lì)模式中，由自勵(lì)模式的拓?fù)淇梢缘玫街绷髂妇€電流idc和直流母線電壓Uc的表達(dá)式： （51） （52）式中： ik1—?jiǎng)?lì)磁階段第k相繞組的電流； ik2—發(fā)電階段第k相繞組的電流； RL—負(fù)載電阻； iL—負(fù)載電流。 else if(v_q_m==s) V=v_uc。 else u=0。 else u=v_q_m。 end end elseif(theta) if(deta_im=*i_ref) u=0。 else if(deta_im=*i_ref) u=s。amp。滯環(huán)控制的仿真程序如下所示： function u=zhihuan(i_m,v_q_m,theta) global theta_on theta_off i_ref s deta_im=i_refi_m。并且當(dāng)導(dǎo)通角較小，轉(zhuǎn)速較低的時(shí)候，在它的發(fā)電階段的運(yùn)動(dòng)電動(dòng)勢(shì)小于它的發(fā)電電壓，這使得它的相電流不能夠增加，會(huì)影響它的發(fā)電能力。 else L=L_min。 else if (theta) L=L_max。它的仿真程序如下所示：function L=inductance(theta) global L_max L_min theta_rise if (theta) L=L_min。θa=45176。θ3=176。在一個(gè)周期內(nèi)，相電感的模型可以劃分為五段來(lái)表示，其中，θ1=176。則對(duì)于三相12/8極的開關(guān)磁阻發(fā)電機(jī)來(lái)說(shuō)，轉(zhuǎn)子極距是45176。 開關(guān)磁阻發(fā)電機(jī)線性模型仿真以下用A相為例子說(shuō)明，輸入的量有電機(jī)的轉(zhuǎn)速n、勵(lì)磁電壓Us；輸出的量由電機(jī)的相電流ia、相轉(zhuǎn)矩Ta還有相電感La。表51 仿真參數(shù)電機(jī)相數(shù)3定子極數(shù)12轉(zhuǎn)子極數(shù)8電感最小值30mH電感最大值300mH相繞組電阻定子極弧14186。 仿真參數(shù)的設(shè)置本次的畢業(yè)設(shè)計(jì)以一臺(tái)三相12/8極的開關(guān)磁阻發(fā)電機(jī)作為研究對(duì)象，它的額定電壓為220V，額定轉(zhuǎn)速為800r/min，勵(lì)磁電源為50V的直流電源。經(jīng)過(guò)三十多年的發(fā)展，目前已經(jīng)成為科學(xué)計(jì)算和系統(tǒng)仿真的首先軟件工具。大致的描述了開關(guān)磁阻發(fā)電機(jī)的雙閉環(huán)的調(diào)節(jié)的過(guò)程，并且詳細(xì)的描述了其中的穩(wěn)壓裝置、滯環(huán)控制器和電壓發(fā)生指示器的部分。還有關(guān)于開關(guān)磁阻發(fā)電機(jī)的各種控制方式，分別為APC方式、CCC方式和PWM方式，它們?cè)诓煌膱?chǎng)合有各自的優(yōu)點(diǎn)，而對(duì)于本文的研究，采用電流斬波控制方式有很