早教吧作业答案频道 -->其他-->
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.
小弟用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 微分方程数值解...的网友还看了以下:
f(x)=∫(0,2x)f(t/2)dt+ln2,显然f(0)=ln2两边求导f'(x)=f(2x/ 2020-03-31 …
matlab 微分方程数值解小弟用ode45函数求解α=1的情况失败,以下是函数function 2020-05-16 …
matlab求解微分方程组dy(1)/dt=(e12-y(1)/c13)/r6-A0*y(3)/I 2020-05-16 …
积分:(sin t)^2 dt=∫(sin t)^2 dt= 2020-05-16 …
求不定积分∫dx/[(x^2+1)^3]^(1/2)∫dx/[(x^2+1)^3]^(1/2)今天 2020-07-03 …
定积分,求解.求极限limx→0∫上限x下限0tˆ½eˆ(tˆ3)dt/∫上限xˆ2/3下限0eˆ 2020-07-09 …
∫[(sinx)^3-(sinx)^5]dx∫x^3(1-x^2)^1/2dx设x=sint,(1 2020-07-20 …
理论力学中动能定理的证明问题,和求导有关(大学)我在书上看到一个动能定理的证明,其中有一个步骤有疑 2020-07-29 …
请教大家一个变上限积分变量替换的问题复习全书29页,例1.46里F(x)=∫(0,x^2)f[(x 2020-08-02 …
ds^2=-dt^2+(ds-v^2f)(r^2)dt)^2+dy^2+dz^2这个关于相对论的方程 2020-12-23 …