求助，Matlab编程解一阶非线性常微分方程组，给高分。谁知道程序怎么修改？
functionYd=power_tra(t,YY)globaly0TjPtEUX1YYYd=[(YY(2)-1)*y0;(Pt-(E*U/X1)*sin(YY(1)))...
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 展开
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 展开
展开全部
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')
本回答由提问者推荐
已赞过
已踩过<
评论
收起
你对这个回答的评价是？
Sievers分析仪
2024-10-13 广告
是的。传统上，对于符合要求的内毒素检测，最终用户必须从标准内毒素库存瓶中构建至少一式两份三点标准曲线；必须有重复的阴性控制；每个样品和PPC必须一式两份。有了Sievers Eclipse内毒素检测仪，这些步骤可以通过使用预嵌入的内毒素标准...
点击进入详情页
本回答由Sievers分析仪提供
推荐律师服务：
若未解决您的问题，请您详细描述您的问题，通过百度律临进行免费专业咨询
