关于微分方程部分←_←出问题了所以求助
mathematica吧
全部回复
仅看楼主
level 9
289728525 楼主
在 t == 0.8012366321093402` 处,步长实际上为零;可能存在奇点或者刚性系统
←_←二楼上代码
2016年01月12日 14点01分 1
level 9
289728525 楼主
NDSolveValue[{x''[t] + (1 - (x[t])^2) (x'[t]) - x[t] == 0, x[0] == 2,x'[0] == 0}, x[t], {t, 0, 2}]
Plot[%, {t, 0, 2}]
得到的图非常奇葩←_←我怀疑是溢出了←_←
然而我貌似没办法把图贴上来←_←
2016年01月12日 14点01分 2
level 9
289728525 楼主
换了个游览器←_←所以才←_←
2016年01月12日 14点01分 3
level 9
289728525 楼主
室友matlab的答案←_←
2016年01月12日 14点01分 4
吧务
level 12
你的方程和你舍友解出的结果不一致,找他确认一下你们的方程一样吧。
由x[0]==2,x'[0]==0,可知蓝色是x[t],绿色是x'[t],那从图上可以明显的看出x''[0]<0
而根据你的方程可以得到x''[t] == x[t] - x'[t] + x[t]^2 x'[t],从而x''[0]=2>0
你这个方程我试着用MATLAB解了下,结果和Mathematica一样,都是到0.8左右就开始出问题,所以我觉得多半是你方程输错了
2016年01月12日 15点01分 5
有道理;明天问下,顺便看看是不是另一个方程的解←_←
2016年01月12日 15点01分
1