matlab微分方程求解 用ode23函数就是这样...最后要做出图像

来源:学生作业帮助网 编辑:作业帮 时间:2024/11/28 23:07:54
matlab微分方程求解 用ode23函数就是这样...最后要做出图像
xVnE~UP*zgfflЋJ&Bk/RTMR@m RhJMk6AxxƻkI#3眙o/6iQmOޚ~>@9j\in\=\y&lms4Xs'wrYݘɰVpˣ?7'OWɳ{E^. 4e ʍG=-JXȺ'I-M!Z.rZ5b%eAn3V-U& OOϒ P#6$Qk|gMRs$<Jd$8.'`@b,fCGg-\>qPd¥“6^wBr|E$R>t]fRTbֳI/NC)xW7W6j=t^t築N˩4n7,;I< ndN*c,('W? :럑R|hl4'BVq,fȆX01#(7-

matlab微分方程求解 用ode23函数就是这样...最后要做出图像
matlab微分方程求解 用ode23函数

就是这样...

最后要做出图像

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
画出的图像:

给的条件不够,只好自己假设了

 l = 1m, g=9.8  θ0=π/12


%存为fun.m

function xdot=fun(t,theta)

g=9.8; l=1;

xdot=[theta(2);-g*sin(theta(1))/l];


%存为main.m

theta0=[pi/12,0];

[t, theta]=ode23('fun',[0,10],theta0);

plot(t,theta(:,1));


g=10;

l=1;

fun=@(t,x) [x(2),;-g/l*sin(x(1))];

ts=[0,20];

x0=[pi/20,0];

[T,X]=ode23(fun,ts,x0);

t=T';

theta=X(:,1)';

plot(t,theta);

xlabel('t');

ylabel('\theta');