我有一个颂歌解决者
[T, Y] = ode23(@(t,y)dydt(t,y),t,y0);
返回T列表和Y列表。
我还尝试绘制变量v
,它是Y
对T
的函数。在特定情况下,v
接近1,例如,0.999999950378838,但随T变化。
但是,T
会返回非常接近的时间,例如
0.336349531356886
0.336405066970296
0.336460602583705
0.336516138197115
0.336571661388434
0.336627184579754
0.336682707771074
0.336738230962393
0.336793745468554
0.336849259974715
0.336904774480876
0.336960288987037
0.337015784154427
0.337071279321817
0.337126774489207
0.337182269656597
0.337237743635592
0.337293217614587
0.337348691593582
0.337404165572577
0.337459629339774
0.337515093106971
0.337570556874168
0.337626020641365
0.337681473992671
0.337736927343978
0.337792380695284
0.337847834046591
0.337903265963968
以10e-5
为单位。当我绘制v
vs T
时,而不是细线,它是一个厚带:
我怀疑是因为单位10e-5
太小或者颂歌解决方案来回解决。但输出T
应该是所有成功的步骤。我怎么处理这个?例如,更改ode选项并获取细线。 v
不完全是1.它可以像0.999999950378838,0.999999905129327和0.999999885006476并随T
一样变化。如何用这个绘制细线并显示精度呢?