早教吧作业答案频道 -->其他-->
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 微分方程数值解...的网友还看了以下：
对某式两边取积分是不是这么理解的(不懂装懂的同学请绕道啊别害人啊)比如对i(t)=Cdv(t)/d 2020-05-13 …
函数分解一个函数x=f（t）,在原点附近是连续可导的在原点附近可分解为x=t*df/dt+R(t) 2020-05-14 …
matlab 微分方程数值解小弟用ode45函数求解α=1的情况失败,以下是函数function 2020-05-16 …
拉氏变换的问题这个题的拉氏变换有几个步骤不清楚,求详解：已知系统微分方程为d²x0（t）/dt²+ 2020-05-23 …
设f（x），g（x）在[a，b]上连续，且满足∫xaf(t)dt≥∫xag(t)dt，x∈[a，b 2020-06-12 …
关於e(x)常微分方程的问题dx(t)/dt=exp(x(t))如何求x(t)之前有试过,dx/d 2020-07-23 …
当x∈[0,4]求定积分∫1/(1+√x)dx设√x=t则t∈[0,2]dx=2tdt原式=2∫t 2020-07-23 …
f(x)=e^x+∫tf(t)dt-x∫f(t)dt解f'(x)=e^x+xf(x)-∫f(t)d 2020-07-31 …
高数变上限积分求导问题!题目是这样的：曲线y=∫sin(x-t)dt(下限为0,上限为x)在点x= 2020-07-31 …
一个简单多变函数的积分证明..f(x)=∫h(x,t)dt=h(x,g(x))g(x)+∫hx(x 2020-08-02 …