Matlab编程解一阶非线性常微分方程组,谁知道程序怎么修改?
来源:学生作业帮 编辑:搜搜做题作业网作业帮 分类:综合作业 时间:2024/07/04 15:51:41
Matlab编程解一阶非线性常微分方程组,谁知道程序怎么修改?
![](http://img.wesiedu.com/upload/c/55/c5550ede7c99ce183a19498d635e1023.jpg)
function Yd = power_tra(t,YY)
global y0 Tj Pt E U X1 YY
Yd=[(YY(2)-1)*y0;(Pt-(E*U/X1)*sin(YY(1)))/Tj;];
global w0 Tj Pt E U X1
w0=2*pi*50;Tj=11.28;Pt=1;E=1.47;U=1;
X1=2.82;
tspan=[0.0 0.3];
y1=[31.54*pi/180;1]
[t,YY]=ode45('power_tra',tspan,y1);
x=YY(:,1);
y=YY(:,2);
plot(t,x*180/pi);
xlabel('t/s');
ylabel('delta/deg');
grid on
![](http://img.wesiedu.com/upload/c/55/c5550ede7c99ce183a19498d635e1023.jpg)
function Yd = power_tra(t,YY)
global y0 Tj Pt E U X1 YY
Yd=[(YY(2)-1)*y0;(Pt-(E*U/X1)*sin(YY(1)))/Tj;];
global w0 Tj Pt E U X1
w0=2*pi*50;Tj=11.28;Pt=1;E=1.47;U=1;
X1=2.82;
tspan=[0.0 0.3];
y1=[31.54*pi/180;1]
[t,YY]=ode45('power_tra',tspan,y1);
x=YY(:,1);
y=YY(:,2);
plot(t,x*180/pi);
xlabel('t/s');
ylabel('delta/deg');
grid on
![Matlab编程解一阶非线性常微分方程组,谁知道程序怎么修改?](/uploads/image/z/4752810-18-0.jpg?t=Matlab%E7%BC%96%E7%A8%8B%E8%A7%A3%E4%B8%80%E9%98%B6%E9%9D%9E%E7%BA%BF%E6%80%A7%E5%B8%B8%E5%BE%AE%E5%88%86%E6%96%B9%E7%A8%8B%E7%BB%84%2C%E8%B0%81%E7%9F%A5%E9%81%93%E7%A8%8B%E5%BA%8F%E6%80%8E%E4%B9%88%E4%BF%AE%E6%94%B9%3F)
clc
Yd=inline('[(YY(2)-1)*2*pi*50;(1-(1.47*1/2.82)*sin(YY(1)))/11.28]','t','YY');
[t,YY]=ode45(Yd,0.3,[31.54*pi/180;1])
x=YY(:,1);
y=YY(:,2);
plot(t,x*180/pi);
xlabel('t/s'),ylabel('delta/deg')
![](http://img.wesiedu.com/upload/0/a8/0a87890fabed6026591cc96f9004bde8.jpg)
Yd=inline('[(YY(2)-1)*2*pi*50;(1-(1.47*1/2.82)*sin(YY(1)))/11.28]','t','YY');
[t,YY]=ode45(Yd,0.3,[31.54*pi/180;1])
x=YY(:,1);
y=YY(:,2);
plot(t,x*180/pi);
xlabel('t/s'),ylabel('delta/deg')
![](http://img.wesiedu.com/upload/0/a8/0a87890fabed6026591cc96f9004bde8.jpg)