fangzhen.m

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

function fangzhen 
%RIGIDODE  Euler equations of a rigid body without external forces
tspan = [0:.1:1];
y0 = [0;0;0;0;0;0];
% Solve the problem using ode45
[t,y]=ode23tb('tongbufadianji',tspan,y0);
%ic=cos(t+2*pi/3).*y(:,1)-sin(t+2*pi/3).*y(:,2);
plot(t,y(:,4))

% 编写求解自激过程微分方程的脚本函数
% 将该脚本函数定义为sh_ge_se_ex_so(shunt_generator_self_excited_solver)
%[t,iff]=ode23(@sh_ge_se_ex_ode,[0 10],[0])
% 指定全局变量（注意全局变量必须在同时使用的多个函数中同时指定）
% 绘图（这里采用多子图法同时绘制ea=f(if)、uf=f(if)、ea=f(t)、if=f(t)等多条曲线）
% 首先绘制ea=f(if)和uf=f(if)曲线
%subplot(3,1,1),
% 下面指定纵横轴标签
xlabel('If[A]')
ylabel('Ea\Uf[V]')
% 下面绘制if=f(t)曲线