sh_ge_se_ex_so.m

《MATLAB电机仿真精华50例》源代码

% 编写求解自激过程微分方程的脚本函数
% 将该脚本函数定义为sh_ge_se_ex_so(shunt_generator_self_excited_solver)
[t,iff]=ode23(@sh_ge_se_ex_ode,[0 10],[0])
% 指定全局变量（注意全局变量必须在同时使用的多个函数中同时指定）
global   a1 a2 a3 a4 Rf
% 绘图（这里采用多子图法同时绘制ea=f(if)、uf=f(if)、ea=f(t)、if=f(t)等多条曲线）
% 首先绘制ea=f(if)和uf=f(if)曲线
subplot(3,1,1),
ea=a1*iff.^3+a2*iff.^2+a3*iff.^1+a4;
plot(iff,ea)
hold on
uf=Rf.*iff;
plot(iff,uf,'r'); % 'r'表示用红颜色绘制曲线
% 下面指定纵横轴标签
xlabel('If[A]')
ylabel('Ea\Uf[V]')
% 下面绘制if=f(t)曲线
subplot(3,1,2),
plot(t,iff);
xlabel('Time[T]')
ylabel('If[A]')
% 下面绘制ea=f(t)曲线
subplot(3,1,3),
plot(t,(a1*iff.^3+a2*iff.^2+a3*iff.^1+a4));
xlabel('Time[T]')
ylabel('ea[V]')
% 下面计算励磁电流和感应电势的数值，然后通过指令窗口读取二者的稳态值
iff,a1*iff.^3+a2*iff.^2+a3*iff.^1+a4