作业帮matlab微分方程求解 用ode23函数
来源：学生作业帮 编辑：作业帮 分类：综合作业 时间：2024/10/06 01:50:01
matlab微分方程求解 用ode23函数
就是这样...
最后要做出图像
就是这样...
最后要做出图像
比如我给初始值：初始值：θ(0)=0.5; θ'(0)=0;
下面用y代表θ,然后,函数程序：
function yprime=verderpol(t,y)
global g;
global l;
g=9.8;
l=10;
yprime=[y(2);-g/l*sin(y(1))];
end
主程序：
clc; clear all
global g;
global l;
g=9.8;
l=10;
y0=[0.5;0];%初始值：y(0)=0.5;y'(0)=0;
hold on
[t,y]=ode23('verderpol',[0,3*pi],y0);
plot(t,y(:,1),'r','linewidth',2);
plot(t,y(:,2),'g','linewidth',2);
legend('θ-t图','θ''-t图')
grid on
hold on
画出的图像：
下面用y代表θ,然后,函数程序：
function yprime=verderpol(t,y)
global g;
global l;
g=9.8;
l=10;
yprime=[y(2);-g/l*sin(y(1))];
end
主程序：
clc; clear all
global g;
global l;
g=9.8;
l=10;
y0=[0.5;0];%初始值：y(0)=0.5;y'(0)=0;
hold on
[t,y]=ode23('verderpol',[0,3*pi],y0);
plot(t,y(:,1),'r','linewidth',2);
plot(t,y(:,2),'g','linewidth',2);
legend('θ-t图','θ''-t图')
grid on
hold on
画出的图像：