早教吧作业答案频道 -->其他-->
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 微分方程数值解...的网友还看了以下:
感应电动势的应用E=n△φ/△t是感生电动势的公式,因为B的变化才出现了感应电动势,为什么不写成E= 2020-03-30 …
衰变公式N(t)=N0e^(-λt)中N表示的是原子核的个数还是核子(质子和中子)的总个数? 2020-04-05 …
PI的 初始值为什么是pi=1#includemain(){int s;float n,t,pi; 2020-05-16 …
云一剐,玉一梭,澹澹衫儿薄薄罗.这里的澹、薄应该怎么读呢?薄拼音:[báo][bó][bò]澹拼音 2020-05-22 …
讨论:关于如何求卷积x(t)*h(-t)的积分表达式?以前信号与系统里学过了x(t)*h(t)的表 2020-06-06 …
设{a下n}是公比大于1的等比数列,s下n为其前n项和,已知s下3=7,且a下1+3,3a下2,a 2020-07-30 …
对于数列{an},若∀m,n∈N*(m≠n),都有am-anm-n≥t(t为常数)成立,则称数列{ 2020-08-01 …
拼读下列音标并写出相应单词./reə//kaɪts//bedz//ka:dz//hændz//'tʃ 2020-12-04 …
有一幅磁通量随时间变化的正弦图像,磁通量最大值是a(这是第一个波峰),这是时间为t1;一开始时间t0 2020-12-15 …
高中数学1、已知等比数列{An}前n项和Sn=2*3^n=t,求t的值及数列{An}通向公式2、一个 2020-12-31 …