早教吧 育儿知识 作业答案 考试题库 百科 知识分享

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.
▼优质解答
答案和解析
改一下初值