matlab电机仿真精华50例

目录前言 (1)1 异步电动机动态数学模型 (2)1.1电压方程 (2)1.2磁链方程 (3)1.3转矩方程 (5)1.4运动方程 (6)2 坐标变化和变换矩阵 (8)2.1三相--两相变换（3/2变换） (8)3 异步电动机仿真 (9)3.1异步电机仿真框图及参数 (9)3.2异步电动机的仿真模型 (11)4 仿真结果 (15)5 结论 (16)参考文献 (17)前言随着电力电子技术与交流电动机的调速和控制理论的迅速发展,使得异步电动机越来越广泛地应用于各个领域的工业生产。

异步电动机的仿真运行状况和用计算机来解决异步电动机控制直接转矩和电机故障分析具有重要意义。

它能显示理论上的变化,当异步电动机正在运行时,提供了直接理论基础的电机直接转矩控制(DTC),并且准确的分析了电气故障。

在过去,通过研究的异步电动机的电机模型建立了三相静止不动的框架。

研究了电压、转矩方程在该模型的功能,同相轴之间的定子、转子的线圈的角度。

θ是时间函数、电压、转矩方程是时变方程这些变量都在这个运动模型中。

这使得很难建立在αβ两相异步电动机的固定框架相关的数学模型。

但是通过坐标变换,建立在αβ两相感应电动机模型框架可以使得固定电压、转矩方程,使数学模型变得简单。

在本篇论文中,我们建立的异步电机仿真模型在固定框架αβ两相同步旋转坐标系下,并给出了仿真结果,表明该模型更加准确地反映了运行中的电动机的实际情况。

1 异步电动机动态数学模型在研究三相异步电动机数学模型时，通常做如下假设 1) 三相绕组对称，磁势沿气隙圆周正弦分布；2) 忽略磁路饱和影响，各绕组的自感和互感都是线性的； 3) 忽略铁芯损耗4) 不考虑温度和频率对电阻的影响异步电机的数学模型由下述电压方程、磁链方程、转矩方程和运动方程组成。

1.1 电压方程三相定子绕组的电压平衡方程为（1-1）与此相应，三相转子绕组折算到定子侧后的电压方程为（1-2）式中 A u , B u , C u , a u , b u ,c u —定子和转子相电压的瞬时值；A i ,B i ,C i , a i , b i ,c i —定子和转子相电流的瞬时值；A ψ,B ψ,C ψ, a ψ, b ψ,c ψ—各相绕组的全磁链； Rs, Rr —定子和转子绕组电阻上述各量都已折算到定子侧，为了简单起见，表示折算的上角标“ ’”均省略，以下同此。

学院：机电学院班级：机自09-2班姓名：谢伟学号：03090943报告关于电机互锁的仿真报告我们要先明确什么叫互锁互锁，说的是几个回路之间，利用某一回路的辅助触点，去控制对方的线圈回路，进行状态保持或功能限制。

一般对象是对其他回路的控制.例如：上图是一个很典型的互锁电路图，实现的功能是电机的正反转，当实现正传的时候，关闭断路器Qs，按下正向按钮SB3，SB3的常闭触点断开，是KM2不工作，SB3的常开触点闭合，是电磁铁KM1得电，线圈KM1吸合实现自锁，电机实现正向转动，同时线圈的常闭触点断开，保证KM2不工作，防止电路短路导致的电路故障，当实现反向转动时，按下反向按钮SB2，SB2的常闭触点断开，KM1不得电，电机无法正向转动，SB2的常开触点闭合，KM2得电，线圈KM2吸合实现自锁，电机实现反向转动，同时线圈KM2 的常闭触点断开，保证KM1不工作。

以上当电机正转时保证电机反转电路断开，反转保证正转电路断开的方法叫作互锁。

利用的是两个常闭辅助触点来实现其功能的。

在实际工作中，由于坏境是变化的，负载也随着环境的变化而变化，为了测试电动机在不同负载环境下的运转情况，我们必须测试一系列不同的或者双向的负载加在电机轴上的电机反应。

此外，一个理想的负载应该回归到发电机从电网上吸收的能量。

这样的负载应该有一个有四象限的DC2或者DC4的直流电产生，这测试的时候应该通过刚性轴这两个电动机进行机械耦合连接。

因此，这个仿真实验需要两个电动机模型，AC4电动机和DC2电动机。

其中，AC4电动机是一个直流三相感应异步电动机，DC2是一个直流单向整流电动机。

在这样的系统里，一个电动机负责速度参数的设定，一个电动机负责力矩参数的设定.同时每个电动机可以充当马达或者发电机。

DC2的额定参数一般是3马力，240伏特，1800转每分钟。

AC4电动机的额定参数是3马力，380伏特，60赫兹，4级。

上图为直流电动机的工作原理图上图为桥式整流电路图注意：对于两个电动机的信号控制可以使用用户界面的最下方的机械输入菜单，下图显示了再马达—发电机结构中如何去模拟刚性轴的互相连接。

MATLAB仿真实验报告班级：：学号：转速闭环一电机参数的设置励磁电流取1A，励磁和电枢互感（Field-armature mutual inductance）取0.673H时，电机转速U N=3000r/min。

其余参数如下：二性能指标超调量δ<=5%最大电流Imax=2I N=175A三模型的设计转速反馈系数：α=10/3000=0.0033v.min/r额定转矩计算：Ceφ=（U-I N*Ra）/n=0.0708 T L=9.55CeφIa=59N.M输出限幅：（-10,10）触发角：30～90°则触发角函数为：（90-6*u ）电源线电压：U=220*/2.34cos(30°)=110V电流限幅175A 采用设置死区来限幅，死区围（-160,160），反馈系数取0.4 平波电抗器取0.01H 调节器选用P 调节 完整模型如下：nifiav+-Voltage Measurement4v +-Voltage Measurement2v +-Voltage Measurement1v +-Voltage MeasurementgABC+-Universal BridgeNABCThree-Phase Programmable Voltage SourceTealpha_deg AB BC CA BlockpulsesSynchronized 6-Pulse GeneratorStepSeries RLC BranchScope4Saturation0.4Gain40.0033Gain3100Gain29.55Gain1(90-6*u)Fcn Dead ZoneDC Voltage SourceT LmA +F +A -F -dcDC Machine10Constant20Constant1<Speed wm (rad/s)><Armature current ia (A)><Field current if (A)><Electrical torque Te (n m)>四仿真结果与调试①放大倍数取100，额定负载启动，突加负载扰动时仿真结果：转速：电枢电流：②放大倍数取100，额定负载启动，突加电源扰动时仿真结果转速：电枢电流：③放大倍数为100，额定负载启动，控制输入电压变化时仿真结果：转速：电枢电流：五实验结论由于控制器采用P调节，转速无法实现无静差。

直流电机调速matlab仿真报告以直流电机调速Matlab仿真报告为标题引言：直流电机是一种常见的电动机，广泛应用于工业、交通、家电等领域。

在实际应用中，电机的调速控制是一项关键技术，可以使电机在不同工况下实现恒定转速或变速运行。

本文将利用Matlab软件进行直流电机调速的仿真实验，旨在通过仿真结果分析不同调速控制策略的优劣，并提供一种基于Matlab的直流电机调速方法。

一、直流电机调速原理直流电机的调速原理基于电压与转速之间的关系。

电机的转速与输入电压成正比，即在给定电压下，电机转速可以通过调整电压大小来实现调速。

常用的直流电机调速方法有电压调速、电流调速和PWM调速等。

二、Matlab仿真实验设置本次仿真实验将以直流电机调速为目标，基于Matlab软件进行实验设置。

首先，需要建立电机的数学模型，包括电机的转速、电流和电压等参数。

其次，选择合适的调速控制策略，如PID控制、模糊控制或神经网络控制等。

最后，通过调节电压输入，观察电机的转速响应和稳定性。

三、PID控制调速实验1. 实验目的本实验旨在通过PID控制器对直流电机进行调速控制，并分析不同PID参数对控制效果的影响。

2. 实验步骤(1) 建立直流电机的数学模型；(2) 设计PID控制器，包括比例系数Kp、积分系数Ki和微分系数Kd；(3) 利用Matlab软件进行仿真，设定电机的目标转速和初始转速；(4) 通过调节PID参数，观察电机的转速响应和稳定性。

3. 实验结果与分析根据实验设置，我们分别对比了不同PID参数值下的电机转速响应曲线。

结果显示，在合适的PID参数设置下，电机能够实现快速响应和稳定控制。

但是，过大或过小的PID参数值都会导致转速超调或调速不稳定的问题。

四、模糊控制调速实验1. 实验目的本实验旨在通过模糊控制器对直流电机进行调速控制，并分析不同模糊规则和输入输出的影响。

2. 实验步骤(1) 建立直流电机的数学模型；(2) 设计模糊控制器，包括模糊规则、输入变量和输出变量；(3) 利用Matlab软件进行仿真，设定电机的目标转速和初始转速；(4) 通过调节模糊规则和输入输出变量，观察电机的转速响应和稳定性。

MATLAB电机控制综合仿真实验一、他励直流电机单闭环调速仿真实验要求：利用Simpowersystem里面自带的DC电机模块，完成他励直流电机单闭环调速仿真，速度调节用PI控制方法，要求封装PI模块，给定速度100rad/s，负载由空载到1s时跳变到20N。

调节不同的PI参数，观察仿真结果总结速度波形、转矩波形的变化规律（PI参数和超调量、稳定时间、稳态误差、振荡次数）。

另外要求将scope图中的4条曲线参数导出到工作空间，并用subplot和plot 函数画在同一个窗口中，每个子图加上对应的标题。

电机相关参数的设置图：仿真原理图：在仿真试验中需要按照实验要求对PI控制器子系统进行封装，然后更改Kp、Ki参数值的大小。

封装PI模块图如下：Plot绘图程序：>>subplot(411)>> plot(t,W,'r'),title('转速')>> subplot(412)>> plot(t,Ia,'b'),title('电枢电流')>> subplot(413)>> plot(t,Te,'g'),title('转矩')>> subplot(414)>> plot(t,If,'y'),title('励磁电流')速度调节用PI控制方法，给定速度100rad/s，负载由空载到1s 时跳变到20N，调节不同的PI参数，从PI模块封装中调节，修改不同的参数Ki 、Kp观察仿真结果。

Ki=100, Kp=5;050100w (r a d /s )00.51 1.52 2.53 3.54 4.55-2000200I a (A )-202I f (A )-1000100T e (N .m )Ki=2, Kp=1;w (r a d /s)I a (A)00.51 1.52 2.53 3.54 4.55I f (A)00.51 1.52 2.53 3.54 4.55T e (N .m )二、 他励直流电机闭环调速系统仿真实验要求：利用Simulink 基本模块搭建他励直流电机闭环调速系统直流电机子模块，根据以下电机数学模型搭建：电磁转矩公式：e M a T C I =Φ 动力学平衡方程：e L m d T T B J dtωω--=电机模块要求封装，参数20.05kg m J =⋅，0.02N m s m B =⋅⋅，165m C =，0.01Wb f Φ=，恒定负载T L =20N 点击封装模块时输入。

三相异步电动机Matlab仿真..中国石油大学胜利学院综合课程设计总结报告题目：三相异步电机直接启动特性实验模型学生姓名：潘伟鹏系别：机械与电气工程系专业年级： 2012级电气工程专业专升本2班指导教师：王铭2013年 6 月 27日一、设计任务与要求普通异步电动机直接起动电流达到额定电流的6--7倍，起动转矩能达到额定转矩的1.25倍以上。

过高的温度、过快的加热速度、过大的温度梯度和电磁力，产生了极大的破坏力，缩短了定子线圈和转子铜条的使用寿命。

但在电网条件和工艺条件允许的情况下，异步电动机也可以直接启动。

本次课程设计通过MATLAB软件建模模拟三相异步电动机直接启动时的各个元器件上的电量变化。

参考：电力系统matlab仿真类书籍电机类教材二、方案设计与论证三相异步电动机直接起动就是利用开关或接触器将电动机的定子绕组直接接到具有额定电压的电网上。

由《电机学》知三相异步电动机的电磁转矩M与直流电动机的电磁转矩有相似的表达形式。

它们都与电机结构（表现为转矩常数）和每级下磁通有关，只不过在三相异步电动机中不再是通过电枢的全部电流，而是点数电流的有功分量。

三相异步电机电磁转矩的表达式为：（1-1）式中——转矩常数——每级下磁通——转子功率因数式(1-1)表明，转子通入电流后，与气隙磁场相互作用产生电磁力，因此，反映了电机中电流、磁场和作用力之间符合左手定则的物理关系，故称为机械特性的物理表达式。

该表达式在分析电磁转矩与磁通、电流之间的关系时非常方便。

从三相异步电动机的转子等值电路可知，（1-2）（1-3）将式(1-2)、(1-3)代入(1-1)得：（1-4）一：我们做如下分析：1.当s=0时，，M=0，说明电动机的理想空载转速为同步转速。

2.当s很小时，有，，说明电磁转矩T近似与s呈线性关系，即随着M的增加，略有下降。

因而，类似直流电动机的机械特性，是一条下倾的直线。

3.当s很大时，有，，说明电磁转矩M近似与s成反比，即M增加时n反而升高。

目录摘要 (1关键词 (11.引言 (12.单相半波可控整流电路 (1 2.1实验目的 (12.2实验原理 (12.3实验仿真 (23.单相桥式全控整流电路 (8 3.1实验目的 (83.2实验原理 (83.3实验仿真 (94.三相半波可控整流电路 (10 4.1实验目的 (104.2实验原理 (114.3实验仿真 (125. 三相半波有源逆变电路 (14 5.1实验目的 (145.2实验原理 (145.3实验仿真 (156.三相桥式半控整流电路 (176.1 实验目的 (176.2实验原理 (17`6.3 实验仿真 (177.小结 (19致谢 (19电力电子技术应用实例的MATLAB 仿真摘要本文是用MATLAB/SIMULINK 实现电力电子有关电路的计算机仿真的毕业设计。

论文给出了单相半波可控整流电路、单相桥式全控整流电路、三相半波可控整流电路、三相半波有源逆变电路、三相桥式全控整流电路的实验原理图、MATLAB 系统模型图、及仿真结果图。

实验过程和结果都表明:MATLAB 在电力电子有关电路计算机仿真上的应用是十分广泛的。

尤其是电力系统工具箱-Power System Blockset (PSB 使得电力系统的仿真更加方便。

关键词 MATLAB SIMULINK PSB 电力电子相关电路1.引言MATLAB 是由Math Works 公司出版发行的数学计算软件,为了准确建立系统模型和进行仿真分析,Math Works 在MATLAB 中提供了系统模型图形输入与仿真工具一SIMULINK 。

其有两个明显功能:仿真与连接,即通过鼠标在模型窗口画出所系统的模型,然后可直接对系统仿真。

这种做法使一个复杂系统模型建立和仿真变得十分容易。

[4][2]在1998年,MathWoIks 推出了电力系统仿真的电力系统工具箱-Power System Blockset (PSB 。

其中包括了电路仿真所需的各种元件模型,包括有电源模块、基础电路模块、电力电子模块、电机模块、连线器模块、检测模块以及附加功率模块等七种模块库。

【仿真实例0-1】解：用M语言编写计算励磁电流程序如下：%磁路计算求解励磁电流问题clc%清除主程序窗口clear%清除变量空间的变量A=0.8*1e-3;%已知铁心截面积m2，1e-3表示10-3kFe=0.94;%已知铁心叠片系数Ph=1*1e-3;%需产生的磁通量Wbu0=4*pi*1e-7;%已知空气磁导率H/m，1e-7表示10-7 l1=0.08;l2=0.1;l3=0.034;l4=0.04;l5=0.1;%已知各段磁路长度mN=2000;%已知励磁绕组匝数d=0.006;%已知气隙长度mAk=kFe*A;%计算净截面积m2B=Ph/Ak;%计算铁心磁通密度uFe=1900*u0;%计算铁心磁导率Hc=B/uFe;%计算铁心磁场强度Fc=Hc*(l1+l2+l3+l4+l5);%计算铁心的磁压降Ha=Ph/u0/A;%计算气隙磁场强度Fa=Ha*d;%计算气隙的磁压降F=Fc+Fa;%计算总磁压降i=F/N;%计算励磁电流s=num2str(i);%将数字转换成字符串s1='励磁电流为：';%定义字符串s=strcat(s1,s,'A');%合并字符串disp(s);%显示计算结果程序运行结果为：励磁电流为：3.0827A【仿真实例0-2】解：计算电感系数的公式为L= μ0μr AN2/l，用M语言编写计算电感系数和绘制电感系数与相对磁导率之间关系曲线程序如下：%求解电感系数和绘制L=f（μr）曲线问题clc%清除主程序窗口clear%清除变量空间的变量u0=pi*4e-7;%已知空气磁导率H/m，1e-7表示10-7 N=400;l=0.3;d=0.005;A=8e-4;A0=8e-4;%已知匝数、铁心长度、气隙、截面积R0=d/(u0*A0);%求气隙磁阻for n=1:80;%for循环语句ur(n)=100+(10000-100)*(n-1)/100;%求相对磁导率Rm(n)=l/(ur(n)*u0*A);%求铁心磁阻R=R0+Rm(n);%求计算总磁阻L(n)=N^2/R;%求电感系数end%for循环语句的结束plot(ur,L)%绘制L=f（μr）曲线title('{\itL=f}（{\it\mu}_r）关系曲线')%标题'L=f（μr）关系曲线'xlabel('铁心相对磁导率{\it\mu}_r')%设置x坐标标签“铁心相对磁导率μr”ylabel('电感系数{\itL}[H]')%设置y坐标标签“电感系数L[H]”【仿真实例0-3】解：用M语言编写拟合磁化曲线的MATLAB程序如下：%拟合与绘制磁化曲线问题clcclearHdata=[1.38,1.40,1.42,1.44,1.46,1.48,1.50,1.52,1.54,1.56,...1.58,1.60,1.62,1.64,1.66,1.69,1.71,1.74,1.76,1.78,...1.81,1.84,1.86,1.89,1.91,1.94,1.97,2.00,2.03,2.06,...2.10,2.13,2.16,2.20,2.24,2.28,2.32,2.36,2.40,2.45,...2.50,2.55,2.60,2.65,2.70,2.76,2.81,2.87,2.93,2.99,...3.06,3.13,3.19,3.26,3.33,3.41,3.49,3.57,3.65,3.74,...3.83,3.92,4.01,4.11,4.22,4.33,4.44,4.56,4.67,4.80,...4.93,5.07,5.21,5.36,5.52,5.68,5.84,6.00,6.16,6.33,...6.52,6.72,6.94,7.16,7.38,7.62,7.86,8.10,8.36,8.62,...8.90,9.20,9.50,9.80,10.1,10.5,10.9,11.3,11.7,12.1,...12.6,13.1,13.6,14.2,14.8,15.5,16.3,17.1,18.1,19.1,...20.1,21.2,22.4,23.7,25.0,26.7,28.5,30.4,32.6,35.1,...37.8,40.7,43.7,46.8,50.0,53.4,56.8,60.4,64.0,67.8];%磁场强度数据Bdata=0.40:0.01:1.69;%磁感应强度数据len=length(Hdata);%计算数组长度Hmax=Hdata(len);%提取数组最大值a=polyfit(Hdata,Bdata,13);%计算多项式拟合系数a for n=1:151%提取拟合数据Hfit(n)=Hmax*(n-1)/150;Bfit(n)=a(1)*Hfit(n)^13+a(2)*Hfit(n)^12+a(3)*Hfit(n)^11+...a(4)*Hfit(n)^10+a(5)*Hfit(n)^9+a(6)*Hfit(n)^8+a(7)*Hfit(n)^7+...a(8)*Hfit(n)^6+a(9)*Hfit(n)^5+a(10)*Hfit(n)^4+a(11)*Hfit(n)^3+...a(12)*Hfit(n)^2+a(13)*Hfit(n)+a(14);endplot(Hdata,Bdata,'*')%对原数据描点绘图hold on%保持当前坐标轴和图形plot(Hfit,Bfit)%绘制多项式拟合曲线hold on%保持当前坐标轴和图形title('磁化曲线')%标题'磁化曲线'xlabel('{\itH}[A/cm]')%x坐标标签“H[A/cm]”ylabel('{\itB}[T]')%y坐标标签“B[T]”【仿真实例0-4】解：用M语言编写绘制磁化曲线的MATLAB程序如下：%绘制磁化曲线问题clcclearHdata=[43,49,55,61,67,72.5,78,84,90,100,110,...123,137,155,173,192,210,240,300,395,...520,643,800,920,1100,1300,1800];%磁场强度H值Bdata=0.2:0.05:1.50;%磁感应强度B值ydata=0:0.001:1.6;%y坐标0~1.6xdata=interp1(Bdata,Hdata,ydata,'spline');%采用样条插值的方法分析数据plot(Hdata,Bdata,'*');%用'*'描点绘制磁化曲线hold on%保持当前坐标轴和图形plot(xdata,ydata);%绘制x，y坐标hold on%保持当前坐标轴和图形title('磁化曲线')%标题'磁化曲线'xlabel('{\itH}(A/m)')%x坐标标签'H(A/m)'ylabel('{\itB}(T)')%y坐标标签'B(T)'ylim([0,1.80])%y坐标标注0~1.8【仿真实例0-5】解：用M语言编写绘制磁滞回线的MATLAB程序如下：%绘制磁滞回线问题clcclearBdata=[0,0.2,0.4,0.6,0.7,0.8,0.9,1.0,0.95,0.9,0.8,...0.7,0.6,0.4,0.2,0,-0.2,-0.4,-0.6,-0.7,-0.8,...-0.9,-1.0,-0.95,-0.9,-0.8,-0.7,-0.6,-0.4,-0.2,0];%磁感应强度基本数据Hdata=[48,52,58,73,85,103,135,193,80,42,2,-18,...-29,-40,-45,-48,-52,-58,-73,-85,-103,-135,...-193,-80,-42,-2,18,29,40,45,48];%磁场强度基本数据plot(Hdata,Bdata)%绘制磁滞回线hold on%保持当前坐标轴和图形title('磁滞回线')%标题'磁滞回线'plot([-250,250],[0,0],'r-')%绘制红色x轴坐标线hold on%保持当前坐标轴和图形plot([0,0],[-1.5,1.5],'r-')%绘制红色y轴坐标线xlabel('{\itH}(A/m)')%x坐标标签'H(A/m)'ylabel('{\itB}(T)')%y坐标标签'B(T)'grid on%显示网格【仿真实例0-6】解：用M语言编写绘制基于【仿真实例0-5】磁滞回线的平均曲线的MATLAB程序如下：%绘制磁滞回线的平均曲线问题clcclearBdata=[-1.0,-0.95,-0.9,-0.8,-0.7,-0.6,-0.4,-0.2,...0,0.2,0.4,0.6,0.7,0.8,0.9,1.0,0.95,0.9,0.8,...0.7,0.6,0.4,0.2,0,-0.2,-0.4,-0.6,-0.7,-0.8,-0.9,-1.0];%磁感应强度基本数据Hdata=[-193,-80,-42,-2,18,29,40,45,48,52,58,73,85,...103,135,193,80,42,2,-18,-29,-40,-45,...-48,-52,-58,-73,-85,-103,-135,-193];%磁场强度基本数据plot(Hdata,Bdata,'linewidth',1.5)%绘制磁滞回线hold on%保持当前坐标轴和图形Bdata1=[-1.0,-0.95,-0.9,-0.8,-0.7,-0.6,-0.4,-0.2,0,0.2,0.4,0.6,0.7,0.8,0.9,1.0];Hdata1=[-193,-80,-42,-2,18,29,40,45,48,52,58,73,85,103,135,193];Bdata2=[-1.0,-0.9,-0.8,-0.7,-0.6,-0.4,-0.2,0,0.2,0.4,0.6,0.7,0.8,0.9,0.95,1.0];Hdata2=[-193,-135,-103,-85,-73,-58,-52,-48,-45,-40,-29,-18,2,42,80,193];%将数据分两组If=-1.0:.05:1.0;%重新设置磁感应强度参数Ean1=spline(Bdata1,Hdata1,If);Ean2=spline(Bdata2,Hdata2,If);%根据两组数据求If对应的样条差值plot((Ean1+Ean2)/2,If,':','linewidth',1.5)%绘制平均值曲线hold on%保持当前坐标轴和图形title('磁滞回线的平均曲线')%标题'磁滞回线的平均曲线'plot([-250,250],[0,0],'r-')%绘制红色x轴坐标线hold on%保持当前坐标轴和图形plot([0,0],[-1.5,1.5],'r-')%绘制红色y轴坐标线xlabel('{\itH}(A/m)')%x坐标标签'H(A/m)'ylabel('{\itB}(T)')%y坐标标签'B(T)'grid on%显示网格【仿真实例1-1】解：用M语言编写绘制基本磁化曲线和磁路未饱和的磁化电流曲线的MATLAB程序如下：clcclearBdata=[0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1,1.1,...%磁通密度基本数据1.2,1.3,1.34,1.38,1.41,1.43,1.45,1.47,...1.48,1.49,1.495,1.5,1.55,1.6];Hdata=[5,10,20,25,40,50,60,70,90,110,150,...%磁场强度基本数据200,300,400,500,600,700,800,900,...1000,1100,1200,1700,5000,10000];B=[-Bdata,0,Bdata];%将磁通密度数据扩展为正负值H=[-Hdata,0,Hdata];%将磁场强度数据扩展为正负值B=sort(B);%对磁通密度数据进行排序H=sort(H);%对磁场强度数据进行排序subplot(2,2,1);%将窗口划分2行2列使用第1子窗口Bx=0:0.01:2*pi;%正弦函数自变量从0~2π，间隔0.01 Bsin=1.1*sin(Bx);%计算正弦磁通密度值，幅值为1.1 plot(Bx,Bsin);%绘制磁通密度正弦曲线grid on;%显示网格线xlim([0,2*pi]);%限定横坐标范围为0~2πxlabel('角度{\it\omegat}/rad');%在横坐标上标注'角度ωt/rad'ylabel('磁通密度{\itB}(T)');%在纵坐标上标注'磁通密度B(T)'subplot(2,2,2);%将窗口划分2行2列使用第2子窗口hold on%保持图形plot(H,B,'ro');%用红色'o'绘制原始数据点grid on;%显示网格线xlabel('磁场强度{\itH}/(A/m)');%在横坐标上标注'磁场强度H/(A/m)' ylabel('磁通密度{\itB}(T)');%在纵坐标上标注'磁通密度B(T)'mymodel=fittype('a*sinh(b*x)');%选择sinh为拟合模型opts=fitoptions(mymodel);%初始化设置set(opts,'Robust','LAR','Normalize','Off');%设置使用线性最小二乘法拟合Fit=fit(B',H',mymodel,opts);%拟合bt=B;%拟合曲线临时磁通密度数据ht=Fit.a.*sinh(Fit.b.*bt);%拟合曲线临时磁场强度数据plot(ht,bt);%在原始数据窗口绘制拟合曲线subplot(2,2,4);%将窗口划分2行2列使用第4子窗口hold on;%保持图形MX=1.6;%磁通密度最大值BI1=sin((-MX:0.01:MX)./MX.*pi).*MX;%磁通正弦变化HI1=Fit.a.*Fit.b.*BI1;%拟合曲线映射后的磁场强度XI1=1:length(BI1);%初始化x轴刻度XI1=XI1/length(BI1)*2*pi;%折算到0~2πplot(XI1,HI1);%绘制磁化电流曲线grid on;%显示网格线xlim([0,2*pi]);%限定横坐标显示范围xlabel('角度{\it\omegat}/rad');%在横坐标上标注'角度ωt/rad'ylabel('磁化电流{\itI}/安匝');%在纵坐标上标注'磁化电流I/安匝'【仿真实例1-2】解：用M语言编写绘制基本磁化曲线和磁路未饱和的磁化电流曲线的MATLAB程序如下：clcclearBdata=[0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1,1.1,...%磁通密度基本数据1.2,1.3,1.34,1.38,1.41,1.43,1.45,1.47,...1.48,1.49,1.495,1.5,1.55,1.6];Hdata=[5,10,20,25,40,50,60,70,90,110,150,...%磁场强度基本数据200,300,400,500,600,700,800,900,...1000,1100,1200,1700,5000,10000];B=[-Bdata,0,Bdata];%将磁通密度数据扩展为正负值H=[-Hdata,0,Hdata];%将磁场强度数据扩展为正负值B=sort(B);%对磁通密度数据进行排序H=sort(H);%对磁场强度数据进行排序subplot(2,2,1);%将窗口划分2行2列使用第1子窗口Bx=0:0.01:2*pi;%正弦函数自变量从0~2π，间隔0.01 Bsin=1.5*sin(Bx);%计算正弦磁通密度值，幅值为1.5 plot(Bx,Bsin);%绘制磁通密度正弦曲线grid on;%显示网格线xlim([0,2*pi]);%限定横坐标范围为0~2πxlabel('角度{\it\omegat}/rad');%在横坐标上标注'角度ωt/rad'ylabel('磁通密度{\itB}(T)');%在纵坐标上标注'磁通密度B(T)'subplot(2,2,2);%将窗口划分2行2列使用第2子窗口hold on%保持图形plot(H,B,'ro');%用红色'o'绘制原始数据点grid on;%显示网格线xlabel('磁场强度{\itH}/(A/m)');%在横坐标上标注'磁场强度H/(A/m)' ylabel('磁通密度{\itB}(T)');%在纵坐标上标注'磁通密度B(T)' mymodel=fittype('a*sinh(b*x)');%选择sinh为拟合模型opts=fitoptions(mymodel);%初始化设置set(opts,'Robust','LAR','Normalize','Off');%设置使用线性最小二乘法拟合Fit=fit(B',H',mymodel,opts);%拟合bt=B;%拟合曲线临时磁通密度数据ht=Fit.a.*sinh(Fit.b.*bt);%拟合曲线临时磁场强度数据plot(ht,bt);%在原始数据窗口绘制拟合曲线subplot(2,2,4);%将窗口划分2行2列使用第4子窗口hold on;%保持图形MX=1.6;%磁通密度最大值BI1=sin((-MX:0.01:MX)./MX.*pi).*MX;%磁通正弦变化HI1=Fit.a.*sinh(Fit.b.*BI1);%拟合曲线映射后的磁场强度XI1=1:length(BI1);%初始化x轴刻度XI1=XI1/length(BI1)*2*pi;%折算到0~2πplot(XI1,HI1);%绘制磁化电流曲线grid on;%显示网格线xlim([0,2*pi]);%限定横坐标显示范围xlabel('角度{\it\omegat}/rad');%在横坐标上标注'角度ωt/rad'ylabel('磁化电流{\itI}/安匝');%在纵坐标上标注'磁化电流I/安匝'【仿真实例1-3】解：用M语言编写绘制磁滞回线的MATLAB程序如下：clcclearBdata=[0,0.2,0.4,0.6,0.7,0.8,0.9,1.0,0.95,0.9,0.8,...0.7,0.6,0.4,0.2,0,-0.2,-0.4,-0.6,-0.7,-0.8,...-0.9,-1.0,-0.95,-0.9,-0.8,-0.7,-0.6,-0.4,-0.2,0];%输入磁感应强度基本数据Hdata=[0.24,0.26,0.29,0.365,0.425,0.515,0.675,...0.965,0.4,0.21,0.01,-0.09,-0.145,-0.2,...-0.225,-0.24,-0.26,-0.29,-0.365,-0.425,...-0.515,-0.675,-0.965,-0.4,-0.21,-0.01,...0.09,0.145,0.2,0.225,0.24];%输入磁场强度基本数据MX=1.39;%磁通密度最大值H=Hdata;B=Bdata;%磁滞回线数据H=H.*10;B=B.*MX;%调整B和H的显示比例subplot(2,2,1);%将窗口划分2行2列使用第1子窗口Bx=-pi/2:0.01:3*pi/2;%磁通密度数据范围-π/2~3π/2Bsin=1.5*sin(Bx);%计算正弦值扩大1.5倍以适应磁滞回线plot(Bx,Bsin);%绘制磁通密度正弦曲线hold on;%保持图形plot([-pi/23*pi/2],[00],'r-')%用红色实线绘制磁密曲线的x坐标轴plot([00],[-22],'r-')%用红色实线绘制磁密曲线的y坐标轴grid on;%显示网格线xlim([-pi/2,3*pi/2]);%限定横坐标显示范围-π/2~3π/2 xlabel('角度{\it\omegat}/rad');%在横坐标上标注'ωt'ylabel('磁通密度{\itB}/(T)');%在纵坐标上标注'磁通密度B/(T)' subplot(2,2,2);%将窗口划分2行2列使用第2子窗口hold on;%保持图形plot(H,B,'ro');%用红色'o'绘制磁滞回线B1=B(1:8);H1=H(1:8);%第1象限数据B>0数据B2=B(8:15);H2=H(8:15);%第2象限数据B>0数据B3=B(16:23);H3=H(16:23);%第3象限数据B>0数据B4=B(23:30);H4=H(23:30);%第4象限数据B>0数据BB1=[B4,B1];HH1=[H4,H1];%磁滞回线的下分支XI1=-10:0.1:10;%设置横坐标值YI1=interp1(HH1,BB1,XI1,'spline');%用3次样条插值计算plot(XI1,YI1);%绘制插值后磁滞回线下分支BB2=[B2,B3];HH2=[H2,H3];%磁滞回线的上分支XI2=-10:0.1:10;%设置横坐标值YI2=interp1(HH2,BB2,XI2,'spline');%用3次样条插值计算plot(XI2,YI2);%绘制插值后磁滞回线上分支grid on;%显示网格线xlabel('磁场强度{\itH}/(A/m)');%在横坐标上标注'磁场强度H/(A/m)' ylabel('磁通密度{\itB}/(T)');%在纵坐标上标注'磁通密度B/(T)'xlim([-1515])%设置横坐标范围subplot(2,2,4);%将窗口划分2行2列使用第4子窗口XI1=sin((0:0.01:MX)./MX.*pi./2).*MX;%磁通密度正弦变化（0~π/2）YI1=interp1(B1,H1,XI1,'spline');%用3次样条插值计算lengthX=length(XI1)-1;%计算先前波形的横坐标长度X1=(0:lengthX);%设置横坐标范围（0~π/2）X1=X1/lengthX*pi./2;plot(X1,YI1);%绘制0~π/2上的曲线hold on;%保持图形XI2=sin((MX:0.01:2*MX)./MX.*pi./2).*MX;%磁通密度正弦变化（π/2~π）YI2=interp1(B(8:16),H(8:16),XI2,'spline');%用3次样条插值计算lengthX=length(XI2)-1;%计算先前波形的横坐标长度X2=(lengthX:2*lengthX);%设置横坐标范围（π/2~π）X2=X2/lengthX*pi./2;plot(X2,YI2);%绘制π/2~π上的曲线XI3=sin((MX:0.01:1.5*MX)./MX.*pi).*MX;%磁通密度正弦变化（π~3π/2）YI3=interp1(B3,H3,XI3,'spline');%用3次样条插值计算lengthX=length(XI3)-1;%计算先前波形的横坐标长度X3=(2*lengthX:3*lengthX);%设置横坐标范围（π~3π/2）X3=X3/lengthX*pi./2;plot(X3,YI3);%绘制π~3π/2上的曲线XI4=sin((1.5*MX:0.01:2*MX)./MX.*pi).*MX;%磁通密度正弦变化（3π/2~2π）YI4=interp1(B4,H4,XI4,'spline');%用3次样条插值计算lengthX=length(XI4)-1;%计算先前波形的横坐标长度X4=(3*lengthX:4*lengthX);%设置横坐标范围（3π/2~2π）X4=X4/lengthX*pi./2;plot(X4,YI4);%绘制3π/2~2π上的曲线grid on;%显示网格线plot([0,2*pi],[00],'r-')%用红色实线绘制磁化电流曲线的x坐标轴xlim([0,2*pi]);%限定横坐标显示范围0~2πxlabel('角度{\it\omegat}/rad');%在横坐标上标注'ωt'ylabel('磁化电流{\itI}/安匝');%在纵坐标上标注'磁化电流I/安匝'【仿真实例2-1】解：用M语言编写绘制【例2-1】功角特性曲线的MATLAB程序如下：clcclearU1=6000/sqrt(3);%定子相电压I1=57.8;%定子相电流xd=64.2;%直轴电抗xq=40.8;%交轴电抗cosfain=0.8;%cosϕN=0.8fain=acos(cosfain);%求ϕNsinfain=sin(fain);%求sinϕNpsi=atan((U1*sinfain+xq)/(U1*cosfain));%求内功率因数角ψId=I1*sin(psi)*exp(1i*pi/2-psi);%计算直轴电流分量Iq=I1*cos(psi)*exp(-psi);%计算交轴电流分量E0=abs(U1*exp(-1i*fain)-1i*Id*xd-1i*Iq*xq);%计算励磁电动势th=0:0.01:pi;%x坐标轴0~π变化x=th./pi.*180;%换算成角度Pem1=3*U1*E0/xd.*sin(th);%求基本电磁功率Pem2=1.5*U1^2*(1/xq-1/xd).*sin(2.*th);%求附加电磁功率Pem=Pem1+Pem2;%凸极电磁功率hold on;%保持当前坐标轴和图形plot(x,Pem1,'b-',x,Pem2,'b-',x,Pem,'b-')%绘制各功角特性plot([0180],[00],'r-')%用红色实现绘制x坐标轴plot([9090],[-2e510e5],'r--')%用红色虚线绘制90°对应值text(50,9e5,'Pem')%在相应位置标注“Pem”text(130,7e5,'Pem1')%在相应位置标注“Pem1”text(70,1.5e5,'Pem2')%在相应位置标注“Pem2”xlabel('Angle[°]');%横坐标标注“Angle[°]”ylabel('Power[kW]');%纵坐标标注“Power[kW]”title('三相凸极同步电动机功角特性');%标题“三相凸极同步电动机功角特性”【仿真实例2-2】解：用M语言编写绘制【例2-2】矩角特性曲线的MATLAB程序如下：clcclearU1=3464;%定子相电压E0=6378;%励磁电动势xc=64;%同步电抗n1=1000;%同步转速th=0:.01:pi;%x坐标轴0~π变化x=th./pi.*180;%换算成角度T=90*U1*E0/(xc*pi*n1).*sin(th);%求电磁转矩plot(x,T)%绘制矩角特性xlabel('Angle[°]');%横坐标标注“Angle[°]”ylabel('Torque[N\cdotm]');%纵坐标标注“Torque[N ·m]”title('三相隐极同步电动机矩角特性');%标题“三相隐极同步电动机矩角特性”【仿真实例3-1】解：根据例题3-1的解得该电动机固有机械特性方程为1230123061972.ss ..s s s s T T mm m +=+=编制绘制该电动机固有机械特性的M 文件程序如下：clc clear n1=1500;%输入同步转速s=0:0.005:1;%转差率变化范围0~1，间隔0.005T=197.6./(0.123./s+s./0.123);%计算电磁转矩Subplot(1,2,1)%按照1行2列绘制子图1plot(s,T,'k-');%绘制机械特性xlabel('转差率{\its}');%横坐标标注“转差率s ”ylabel('电磁转矩{\itT}/(N\cdotm)');%纵坐标标注“电磁转矩T /(N·m)”title('异步电动机固有机械特性{\itT}={\itf}({\its})')%标注标题“T =f (s )”n=n1.*(1-s);%计算转速Subplot(1,2,2)%按照1行2列绘制子图2plot(T,n,'k-');%绘制机械特性xlabel('电磁转矩{\itT}/(N\cdotm)');%横坐标标注“电磁转矩T /(N·m)”ylabel('转速{\itn}/(rpm)');%纵坐标标注“转速n /(rpm)”title('异步电动机固有机械特性{\itn}={\itf}({\itT})')%标注标题“n =f (T )”hold on;%保持当前坐标轴和图形【仿真实例3-2】解：用M 语言编写绘制【例3-2】人为机械特性曲线的MATLAB 程序如下：clc clear m1=3;%异步电动机相数U1=220;%定子相电压n1=1000;%输入同步转速p=2;%极对数f1=50;%电源频率r1=2.08;%定子绕组电阻r2=1.53;%转子绕组电阻折算值w1=2*pi*f1/p;%计算同步角速度，单位：rad/sx1=3.12;%定子漏电抗x2=4.25;%转子漏电抗折算值s=0:0.005:1;%设定转差率变化范围：0~1，间隔0.005n=n1.*(1-s);%计算转速T=(m1*p*U1^2*r2)./s./(w1.*((r1+r2./s).^2+(x1+x2)^2));%计算电磁转矩subplot(2,2,1)%按2行2列绘制子图plot(s,T,'k');%绘制固有机械特性曲线xlabel('转差率{\its}');%横坐标标注“转差率s”ylabel('电磁转矩{\itT}/(N\cdotm)');%纵坐标标注“电磁转矩T/(N·m)”str_x=0.25;%标注字符的横坐标text(str_x,max(T),strcat('U1=',num2str(int16(U1)),'V'),'FontSize',9,'Color','black'); %标注固有机械特性曲线的电压值title('降低定子电压的人为机械特性')%标题“降低定子电压的人为机械特性”hold on;%保持当前坐标轴和图形subplot(2,2,2)%按2行2列绘制子图2plot(T,n,'k');%绘制固有机械特性曲线hold on;xlabel('电磁转矩{\itT}/(N\cdotm)');%横坐标标注“电磁转矩T/(N·m)”ylabel('转速{\itn}/(rpm)');%纵坐标标注“转速n/(rpm)”text(max(T)-20,500,strcat('U1=',num2str(int16(U1)),'V'),'FontSize',9,'Color','black'); %标注固有机械特性曲线的电压值title('降低定子电压的人为机械特性')%标题“降低定子电压的人为机械特性”str_y=500;%设定字符串的纵坐标初值for coef=0.75:-0.25:0.25;%设定降低定子电压的范围U1p=U1*coef;%降低定子电压T1=(m1*p*U1p^2*r2)./s./(w1.*((r1+r2./s).^2+(x1+x2)^2));%计算电磁转矩subplot(2,2,1)%按2行2列绘制子图1plot(s,T1,'k-');%绘制降低定子电压人为机械特性str=strcat('U1=',num2str(int16(U1p)),'V');%创建标注字符串text(str_x,max(T1)+5,str,'FontSize',9,'Color','black');%标记各曲线的电压值subplot(2,2,2)%按2行2列绘制子图2plot(T1,n,'k-');%绘制降低定子电压人为机械特性str=strcat('U1=',num2str(int16(U1p)),'V');%创建标注字符串str_y=str_y+100;%修改字符串的纵坐标值text(max(T1)-5,str_y,str,'FontSize',9,'Color','black');%标记各曲线的电压值endsubplot(2,2,3)%按2行2列绘制子图3plot(s,T,'k-');%绘制固有机械特性曲线xlabel('转差率{\its}');%横坐标标注“转差率s”ylabel('电磁转矩{\itT}/(N\cdotm)');%纵坐标标注“电磁转矩T/(N·m)”str_x=0.75;%标注字符的横坐标text(str_x,max(T)-45,strcat('r2=',num2str(r2),'\Omega'),'FontSize',9,'Color','black'); %标注固有机械特性曲线的电阻值title('改变转子电阻的人为机械特性')%标题“改变转子电阻的人为机械特性”hold on;%保持当前坐标轴和图形subplot(2,2,4)%按2行2列绘制子图4plot(T,n,'k-');%绘制固有机械特性曲线hold on;xlabel('电磁转矩{\itT}/(N\cdotm)');%横坐标标注“电磁转矩T/(N·m)”ylabel('转速{\itn}/(rpm)');%纵坐标标注“转速n/(rpm)”text(50,400,strcat('r2=',num2str(r2),'\Omega'),'FontSize',9,'Color','black');%标注固有机械特性曲线的电阻值title('改变转子电阻的人为机械特性')%标题“改变转子电阻的人为机械特性”r2p=r2;%设定改变的转子电阻的初值str_y1=400;%设定字符串的纵坐标初值for coef=1:4;%设定改变转子电阻的范围r2p=r2p+0.75;%改变转子电阻T1=(m1*p*U1^2*r2p)./s./(w1.*((r1+r2p./s).^2+(x1+x2)^2));%计算电磁转矩subplot(2,2,3)%按2行2列绘制子图3plot(s,T1,'k-');%绘制改变转子电阻人为机械特性str=strcat('r2=',num2str(r2p),'\Omega');%创建标注字符串str_y=T1(length(T1))-5;%标注字符串的纵坐标值text(str_x,str_y+11,str,'FontSize',9,'Color','black');%标记各曲线的电阻值subplot(2,2,4)%按2行2列绘制子图4plot(T1,n,'k-');%绘制改变转子电阻人为机械特性str=strcat('r2=',num2str(r2p),'\Omega');%创建标注字符串str_y1=str_y1-90;%修改字符串的纵坐标值text(50+coef*8,str_y1,str,'FontSize',9,'Color','black');%标记各曲线的电阻值end【仿真实例4-1】解：用M语言编写绘制直流发电机空载特性曲线的MATLAB程序如下：%绘制直流发电机空载特性曲线问题clcclearIfdata1=[0.0,0.5,0.7,0.8,1.0,1.2,1.5];Ifdata2=[0.0,0.5,0.7,0.8,1.0,1.2,1.5];%励磁电流I f值U0data1=[2,75,93.5,99,106,111.5,117];U0data2=[18,83.6,97.5,102.5,109.5,114,117];%空载电压U0值xdata=0:.1:1.5;%y坐标0~120ydata1=interp1(Ifdata1,U0data1,xdata,'spline');ydata2=interp1(Ifdata2,U0data2,xdata,'spline');%采用样条插值的方法分析数据plot(Ifdata1,U0data1,'*')%用'*'描点绘制空载特性hold on;%保持当前坐标轴和图形plot(Ifdata2,U0data2,'*')%绘制I f,U0坐标hold on;plot(xdata,ydata1);%绘制x,y坐标hold on;%保持当前坐标轴和图形plot(xdata,ydata2);%绘制x,y坐标hold on;%保持当前坐标轴和图形title('直流发电机空载特性')%标题'直流发电机空载特性' xlabel('{\itI}_f(A)')%x坐标标签'I f(A)'ylabel('{\itU}_0(V)')%y坐标标签'U0(V)'axis([0,2,0,120])【仿真实例4-2】解：用M语言编写绘制他励直流发电机外特性曲线的MATLAB程序如下：%绘制他励直流发电机外特性曲线问题clcclearIdata=[0,4.3,7,8.8,13,15,17];%负载电流I值Udata=[118,116,114,112.5,110,108.5,107];%负载电压U值xdata=0:1:17;%x坐标0~17ydata=interp1(Idata,Udata,xdata,'spline');%采用样条插值的方法分析数据plot(Idata,Udata,'*');%用'*'描点绘制外特性hold on%保持当前坐标轴和图形plot(xdata,ydata);%绘制x,y坐标hold on%保持当前坐标轴和图形title('他励直流发电机外特性')%标题'他励直流发电机外特性' xlabel('{\itI}(A)')%x坐标标签'I(A)'ylabel('{\itU}(V)')%y坐标标签'U(V)'axis([0,20,0,150])%设置轴线数据范围【仿真实例4-3】解：用M语言编写绘制并励直流发电机外特性曲线的MATLAB程序如下：%绘制并励直流发电机外特性曲线问题clcclearnN=1500;IfN=1.4;Rf=75;Ra=2.3;%输入发电机基本数据Ifdata=[0,0.5,0.7,0.8,1.0,1.2,1.5];%输入励磁电流U0data=[5,75,95,100,108,110,117];%输入空载电压If=0:0.01:12;%重新设置励磁参数P=polyfit(Ifdata,U0data,3);U0=polyval(P,If);%空载特性曲线拟合U=If*Rf;%计算电枢端电压Ia=(U0-U)/Ra;%计算电枢电流I=Ia-If;%计算负载电流plot(I,U)%绘制并励直流发电机外特性hold on%保持当前坐标轴和图形title('并励直流发电机外特性')%标题'并励直流发电机外特性'axis([0,20,0,120])%设置轴线数据范围xlabel('{\itI}(A)')%x坐标标签'I(A)'ylabel('{\itU}(V)')%y坐标标签'U(V)'【仿真实例4-4】解：用M语言编写仿真并励直流发电机自励过程和求稳态电压的MATLAB程序如下：%仿真并励直流发电机自励过程和求稳态电压问题%sh_ge_se_ex_ode函数的M文件function dydt=sh_ge_se_ex_ode(~,iff)%定义该M—函数为sh_ge_se_ex_ode，即%shunt_generator_self_excited_odeglobal a1a2a3a4Rf%指定全局变量Lf=40;Rf=75;%输入发电机基本数据Ifdata=[0,0.5,0.7,0.8,1.0,1.2,1.5];%输入励磁电流实验数据U0data=[5,75,95,100,108,110,117];%输入感应电动势数据a=polyfit(Ifdata,U0data,3);%曲线拟合dydt=1/Lf*(a(1)*iff^3+a(2)*iff^2+a(3)*iff^1+a(4))-Rf/Lf*iff;%列写标准形式微分方程并将感应电动势用%励磁电流的拟合多项式函数表示a1=a(1);a2=a(2);a3=a(3);a4=a(4);%给全局变量赋值%脚本M文件[t,iff]=ode23(@sh_ge_se_ex_ode,[010],[0]);%选择微分方程解算指令ode23()global a1a2a3a4Rf%指定全局变量subplot(3,1,1),%定义绘制仿真曲线1ea=a1*iff.^3+a2*iff.^2+a3*iff.^1+a4;%求感应电动势plot(iff,ea)%绘制e a=f(I f)曲线hold on%保持当前坐标轴和图形uf=Rf.*iff;%计算u fplot(iff,uf,'g')%用绿颜色绘制u f=f(i f)场阻线title('并励直流发电机自励过程')%标题'并励直流发电机自励过程'xlabel('{\itI}_f(A)')%曲线1x坐标标签'I f(A)'ylabel('{\itE}_a/{\itU}_f(V)')%曲线1y坐标标签'E a/U f(V)'axis([0,1.6,0,150])%设置x、y坐标值subplot(3,1,2)%定义绘制曲线2plot(t,iff);%绘制曲线2即i f=f(t)xlabel('{\itt}(s)')%曲线2x坐标标签't(s)'ylabel('{\itI}_f(A)')%曲线2y坐标标签'I f(A)'axis([0,5,0,2])subplot(3,1,3),%定义绘制曲线3plot(t,(a1*iff.^3+a2*iff.^2+a3*iff.^1+a4));%绘制e a曲线xlabel('{\itt}(s)')%曲线3x坐标标签't(s)'ylabel('{\ite}_a(V)')%曲线3y坐标标签'e a(V)'iff,a1*iff.^3+a2*iff.^2+a3*iff.^1+a4;%计算励磁电流和感应电动势的数值axis([0,5,0,150])【仿真实例4-5】解：用M语言编写绘制他励直流电动机工作特性的MATLAB程序如下：%绘制他励直流电动机工作特性问题clcclearUN=220;PN=22;IaN=115;Nn=1500;%输入电动机参数Ra=0.21;%输入电枢电阻CePhiN=(UN-Ra*IaN)/Nn;%计算电动势常数C eΦNCTPhiN=9.55*CePhiN;%计算电磁转矩常数C TΦNIa=0:IaN;%电枢电流从0~额定电流I aNn=UN/CePhiN-Ra/(CePhiN)*Ia;%计算转速TN=CTPhiN*Ia;%计算电磁转矩TNP=TN*10;%为清楚起见，将电磁转矩扩大十倍显示plot(Ia,n,'b.-',Ia,TNP,'r.-');%绘制转速特性和转矩特性曲线xlabel('电枢电流{\itI}_a/A')%横坐标标签'电枢电流I a/A'ylabel('转速{\itn}/rpm，电磁转矩{\itT}/N.m')%纵坐标'转速n/rpm，电磁转矩T/N.m' text(30,1500,'转速{\itn}');%标记转速曲线text(50,500,'电磁转矩{\itT}（X10）');%标记转矩曲线【仿真实例4-6】解：用M语言编写绘制串励直流电动机工作特性的MATLAB程序如下：%绘制串励直流电动机工作特性问题clcclearUN=220;PN=22;IaN=115;Nn=1500;%输入电动机参数Ra=0.18;%输入电枢电阻Rf=0.31;%输入励磁电阻k=0.01;%输入比例常数CePhiN=(UN-(Ra+Rf)*IaN)/Nn;%计算电动势常数C eΦNCe=CePhiN/k/IaN;%计算电动势常数C eIa=0:IaN;%电枢电流从0~额定电流I aNn=UN./(Ce*k.*Ia)-(Ra+Rf)/(Ce*k);%计算转速start_p=30;%设置显示的起始点Ia_p=Ia(start_p:length(Ia));%截取电流显示区间n_p=n(start_p:length(n));%截取转速显示区间CTPhiN=9.55*CePhiN;%计算电磁转矩常数C TΦNCT=CTPhiN/k/IaN;%计算电磁转矩常数C TTN=CT*k*Ia.*Ia;%计算电磁转矩TNP=TN*30;%为清楚起见，将电磁转矩扩大三十倍显示plot(Ia_p,n_p,'b.-',Ia,TNP,'r.-');%绘制转速特性和转矩特性曲线xlabel('电枢电流{\itI}_a/A')%横坐标标签'电枢电流I a/A'ylabel('转速{\itn}/rpm,电磁转矩{\itT}/N.m')%纵坐标标签'转速n/rpm，电磁转矩T/N.m' text(40,5500,'转速\itn');%标记转速曲线text(15,1000,'电磁转矩{\itT}（X30）');%标记转矩曲线【仿真实例5-1】解：用M语言编写绘制他励直流电动机机械特性的MATLAB程序如下：clcclearPN=22,UN=220,IN=115,nN=1500,Ra=0.21;%输入铭牌数据IaN=IN;%计算电枢电流CePhiN=(UN-Ra*IaN)/nN;%计算电动势常数C eΦNCTPhiN=9.55*CePhiN;%计算电磁转矩常数C TΦNIa=0:IaN;%建立电枢电流数组n=UN/CePhiN-Ra/(CePhiN)*Ia;%计算转速T=CTPhiN*Ia;%计算电磁转矩figure(1);%建立1号图形窗口plot(T,n,'.-');%绘制固有机械特性曲线title('固有机械特性');%标题'固有机械特性'xlabel('电磁转矩{\itT}/N\cdotm');%横轴标注'电磁转矩T/N•m'ylabel('转速{\itn}/rpm');%纵轴标注'转速n/rpm'ylim([0,1800]);%限制纵轴显示范围figure(2);%建立2号图形窗口plot(T,n,'rd');%绘制固有机械特性曲线title('降低电枢电源电压的人为机械特性');%标题'降低电枢电源电压的人为机械特性' xlabel('电磁转矩{\itT}/N\cdot m');%横轴标注'电磁转矩T/N•m'ylabel('转速{\itn}/rpmin');%纵轴标注'转速n/rpm'hold on;%保持当前坐标轴和图形for jy=1:-0.25:0.25;U=UN*jy;%改变电枢电源电压n=U/CePhiN-Ra/(CePhiN*CTPhiN)*T;%计算对应不同电枢电源电压的转速plot(T,n,'-');%绘制改变电枢电源电压的人为机械特性str=strcat('{\it U}=',num2str(U),'V');%显示字符串处理y=1700*jy;%显示字符串纵坐标text(60,y,str);%给曲线标注电压值endfigure(3);%建立3号图形窗口Rc=0;%临时变量n=UN/CePhiN-(Ra+Rc)/(CePhiN*CTPhiN)*T;%计算转速plot(T,n,'rd');%绘制固有机械特性曲线title('电枢回路串电阻的人为机械特性');%标题'电枢回路串电阻的人为机械特性' xlabel('电磁转矩{\itT}/N\cdot m');%横轴标注'电磁转矩T/N•m'ylabel('转速{\itn}/rpmin');%纵轴标注'转速n/rpm'hold on;%保持当前坐标轴和图形Rc=0.02;%电枢串电阻值for Rc=0:0.5:1.9;n=UN/CePhiN-(Ra+Rc)/(CePhiN*CTPhiN)*T;%计算转速plot(T,n,'-');%绘制电枢回路串电阻的人为机械特性str=strcat('{\it R}=',num2str(Ra+Rc),'\Omega');%字符串处理y=400*(4-Rc*1.8);%显示字符串的纵坐标text(120,y,str);%给各曲线标记电阻值endylim([0,1700]);%限制纵轴显示范围figure(4);%建立4号图形窗口n=UN/CePhiN-Ra/(CePhiN*CTPhiN)*T;%计算转速plot(T,n,'rd');%绘制固有机械特性曲线title('减弱磁通的人为机械特性');%标题'减弱磁通的人为机械特性' xlabel('电磁转矩{\itT}/N\cdot m');%横轴标注'电磁转矩T/N•m'ylabel('转速{\itn}/rpmin');%纵轴标注'转速n/rpm'hold on;%保持当前坐标轴和图形for ct=0.5:0.25:1.3;CePhi=CePhiN*ct;CTPhi=CTPhiN*ct;%改变磁通值n=UN/CePhi-Ra/(CePhi*CTPhi)*T;%计算转速plot(T,n,'-');%绘制改变磁通时的人为机械特性str=strcat('{\it\phi}=',num2str(ct),'*\phi_N');%显示字符串处理y=3600-1850*ct;%显示字符串纵坐标text(120,y,str);%给各曲线标记磁通endylim([0,3600]);%限制纵坐标的显示范围。