【正文】 )。 dotL_a=dotinduct(theta_a)。 doti_a=(v_aR_m*i_adotL_a*w_act*i_a)/L_a。 i_a=i_a+doti_a*T_timer1。 if(i_a0) i_a=0。 end L_b=inductance(theta_b)。 dotL_b=dotinduct(theta_b)。 doti_b=(v_bR_m*i_bdotL_b*w_act*i_b)/L_b。 i_b=i_b+doti_b*T_timer1。 if(i_b0) i_b=0。 end L_c=inductance(theta_c)。 dotL_c=dotinduct(theta_c)。 doti_c=(v_cR_m*i_cdotL_c*w_act*i_c)/L_c。 i_c=i_c+doti_c*T_timer1。 if(i_c0) i_c=0。 end theta_a=theta_a+w_act*T_timer1*180/pi。 if(theta_atheta_pole) theta_a=theta_atheta_pole。 end [theta_b,theta_c]=other_theta(theta_a)。 i_dc=(v_q_a*i_a+v_q_b*i_b+v_q_c*i_c)。 % i_l=v_uc/R_l。 i_ic=i_dci_l。 dotv_uc=i_ic/C*T_timer。 v_uc=v_uc+dotv_uc。 if(v_uc0) v_uc=0。 end end v_q_a=zhihuan(i_a,v_q_a,theta_a)。 v_q_b=zhihuan(i_b,v_q_b,theta_b)。 v_q_c=zhihuan(i_c,v_q_c,theta_c)。% T_a=*dotinduct(theta_a)*i_a*i_a。 T_b=*dotinduct(theta_b)*i_b*i_b。 T_c=*dotinduct(theta_c)*i_c*i_c。 data_inductance_a(num+1,1)=L_a。 data_inductance_b(num+1,1)=L_b。 data_inductance_c(num+1,1)=L_c。 data_current_a(num+1,1)=i_a。 data_current_b(num+1,1)=i_b。 data_current_c(num+1,1)=i_c。 data_voltage_a(num+1,1)=v_a。 data_voltage_b(num+1,1)=v_b。 data_voltage_c(num+1,1)=v_c。 data_voltage_q_a(num+1,1)=v_q_a。 data_voltage_q_b(num+1,1)=v_q_b。 data_voltage_q_c(num+1,1)=v_q_c。 data_voltage_uc(num+1,1)=v_uc。 % data_current_dc(num+1,1)=i_dc。 % data_current_ic(num+1,1)=i_ic。 % data_current_l(num+1,1)=i_l。 data_torque_a(num+1,1)=T_a。 data_torque_b(num+1,1)=T_b。 data_torque_c(num+1,1)=T_c。 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。 if(theta_Btheta_pole) theta_B=theta_Btheta_pole。 end theta_C=theta_A+theta_pole/3。 if(theta_Ctheta_pole) theta_C=theta_Ctheta_pole。 end end % function L=inductance(theta) global L_max L_min theta_rise if (theta) L=L_min。 else if (theta) L=L_min+(L_maxL_min)/theta_rise*()。 else if (theta) L=L_max。 else if (theta) L=L_max(L_maxL_min)/theta_rise*()。 else L=L_min。 end end end end end % function dotL=dotinduct(theta) global L_max L_min theta_rise if (theta) dotL=0。 else if (theta) dotL=(L_maxL_min)/theta_rise*180/pi。 else if (theta) dotL=0。 else if (theta) dotL=(L_maxL_min)/theta_rise*180/pi。 else dotL=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(v_q_m==s) V=v_uc。 else V=0。 end end end function u=zhihuan(i_m,v_q_m,theta) global theta_on theta_off i_ref s deta_im=i_refi_m。 if ((thetatheta_on)amp。amp。(thetatheta_off)) if(deta_im=*i_ref) u=s。 else if(deta_im=*i_ref) u=s。 else u=v_q_m。 end end elseif(theta) if(deta_im=*i_ref) u=0。 else if(deta_im=*i_ref) u=s。 else u=v_q_m。 end end else if(i_m=0) u=s。 else u=0。 end end end仿真波形如圖下所示。以下為相電感、相電流、相電壓、相轉(zhuǎn)矩的仿真波形圖：圖51 CCC方式下的相電感、相電流、相電壓、相轉(zhuǎn)矩加入閉環(huán)后的系統(tǒng)仿真圖：圖52 加閉環(huán)后的系統(tǒng)仿真圖初始電壓為300V時的仿真波形圖：圖53 初始電壓為300V的系統(tǒng)仿真圖R=500時的仿真波形圖：圖54 R=500時的系統(tǒng)仿真圖R突然增大時的仿真波形：圖55 R突然增大時的仿真波形圖 本章小結(jié)本章使用MATLAB仿真軟件，搭建了開關(guān)磁阻發(fā)電機的仿真模型，給出了開關(guān)磁阻發(fā)電機模型的仿真結(jié)果。詳細分析了SRG線性模型。建立了典型的開關(guān)磁阻風力發(fā)電系統(tǒng)以及基于不對稱功率變換器的開關(guān)磁阻發(fā)電機穩(wěn)壓系統(tǒng)的仿真模型。仿真結(jié)果驗證了開關(guān)磁阻發(fā)電機的輸出電壓的控制系統(tǒng)的正確性，同時也證明了不對稱功率變換器的正確性、有效性和可行性。 結(jié)論結(jié)論本文對國內(nèi)外開關(guān)磁阻電機近期的研究成果，尤其是開關(guān)磁阻電機作為發(fā)電運行的研究成果進行分析總結(jié)的基礎上，以三相12/8極的開關(guān)磁阻發(fā)電機為對象，從發(fā)電原理、控制方案研究等角度深入的進行了研究，取得了一些具有理論和實用價值的研究成果，這方面的工作主要有:（1） 全面總結(jié)和分析了當前開關(guān)磁阻發(fā)電控制系統(tǒng)的發(fā)展現(xiàn)狀和前景，提出了開關(guān)磁阻發(fā)電機的幾大優(yōu)勢，進一步明確了研究開關(guān)磁阻電機發(fā)電系統(tǒng)的意義。（2） 闡釋了開關(guān)磁阻發(fā)電機的運行和控制機理，指出了開關(guān)磁阻電機作為發(fā)電機所具有的獨特性。詳細推導了SRG電機的線性模型。（3） 采用MATLAB仿真工具建立了電機的線性數(shù)學模型，對SRG線性數(shù)學模型進行了一系列定性分析，獲得了相應結(jié)論。（4） 開關(guān)磁阻發(fā)電機在自勵模式下的穩(wěn)態(tài)輸出電壓與起勵電壓無關(guān)，但起勵電壓越高，輸出電