matlab 微分方程数值解小弟用ode45函数求解α=1的情况失败,以下是函数function dt=eq1(r,t)dt(2)=-2/r*t(2)+8.5e-28*t(1)-1e14*cos(pi*r/2); %因为不会写导数,我自己手动展开了一下dt(1)=t(2);dt=dt';r0=0;rf=1;[r,t]=od

来源:学生作业帮助网 编辑:六六作业网 时间:2024/12/25 02:43:43
matlab微分方程数值解小弟用ode45函数求解α=1的情况失败,以下是函数functiondt=eq1(r,t)dt(2)=-2/r*t(2)+8.5e-28*t(1)-1e14*cos(pi*r

matlab 微分方程数值解小弟用ode45函数求解α=1的情况失败,以下是函数function dt=eq1(r,t)dt(2)=-2/r*t(2)+8.5e-28*t(1)-1e14*cos(pi*r/2); %因为不会写导数,我自己手动展开了一下dt(1)=t(2);dt=dt';r0=0;rf=1;[r,t]=od
matlab 微分方程数值解

小弟用ode45函数求解α=1的情况失败,以下是函数

function dt=eq1(r,t)
dt(2)=-2/r*t(2)+8.5e-28*t(1)-1e14*cos(pi*r/2);  %因为不会写导数,我自己手动展开了一下
dt(1)=t(2);
dt=dt';


r0=0;
rf=1;
[r,t]=ode145(@eq1,[r0,rf],[0;0]);


错误原因:t算出来都是NaN.

matlab 微分方程数值解小弟用ode45函数求解α=1的情况失败,以下是函数function dt=eq1(r,t)dt(2)=-2/r*t(2)+8.5e-28*t(1)-1e14*cos(pi*r/2); %因为不会写导数,我自己手动展开了一下dt(1)=t(2);dt=dt';r0=0;rf=1;[r,t]=od
改一下初值