微分方程，求matlab 代码，不胜感激！求t=[0,20],对应的速度V
答案:2 悬赏:80 手机版
解决时间 2021-03-20 09:38
- 提问者网友：树红树绿
- 2021-03-19 10:07
微分方程，求matlab 代码，不胜感激！求t=[0,20],对应的速度V
最佳答案
- 五星知识达人网友：渡鹤影
- 2021-03-19 10:35
用ode23（）函数可以较好求得微分方程的数值解。
v0=[2.3,0];
h=0.5;
a=0;
b=20;
[t,v]=ode23(@fun,a:h:b,v0);
求解结果：
t V(t) dV(t)
0 2.3000 0
0.5000 5.2562 3.9856
1.0000 4.4588 -7.0559
1.5000 2.8150 6.5266
2.0000 5.5829 0.7497
2.5000 3.3187 -7.8804
3.0000 3.9351 7.8954
3.5000 5.4400 -2.5105
4.0000 2.4192 -3.2536
4.5000 4.9425 5.3320
5.0000 4.7890 -5.9263
5.5000 2.5304 4.4832
6.0000 5.4779 1.8393
6.5000 3.6658 -7.9648
7.0000 3.5824 7.9280
7.5000 5.4794 -1.6476
8.0000 2.5728 -4.7430
8.5000 4.7258 5.9901
9.0000 4.9100 -5.2930
9.5000 2.4462 3.2745
10.0000 5.3931 2.3321
10.5000 3.7920 -7.8113
11.0000 3.4505 7.7593
11.5000 5.4600 -1.3779
12.0000 2.6260 -5.0057
12.5000 4.6576 6.0813
13.0000 4.8927 -5.2044
13.5000 2.4588 3.1977
14.0000 5.3626 2.2714
14.5000 3.7309 -7.7293
15.0000 3.5055 7.6847
15.5000 5.4053 -1.6398
16.0000 2.5687 -4.3342
16.5000 4.7281 5.6945
17.0000 4.7518 -5.6085
17.5000 2.5561 4.1328
18.0000 5.3754 1.6801
18.5000 3.4960 -7.5912
19.0000 3.7339 7.5863
19.5000 5.2974 -2.4283
20.0000 2.4512 -2.6040
追问您好，初始条件只有速度V=2.3，没有加速度，刚刚搞错了。
能否麻烦你把定义函数的代码也一并说下，所有代码，谢谢追答用ode23（）函数可以较好求得微分方程的数值解。执行代码如下：
v0=0
h=0.5;a=0;b=20;
[t,v]=ode23(@fun,a:h:b,v0);
求解结果：
t V(t)
0 0
0.5 3.63716291516561
1 3.64158930086147
1.5 3.6428487523781
2 3.63885309188913
2.5 3.6339010775317
3 3.6342955704192
3.5 3.63569824585008
4 3.63869545971492
4.5 3.64207939462941
5 3.64279555495064
5.5 3.64076897493905
6 3.6375611660217
6.5 3.63456757976042
7 3.63417276743521
7.5 3.63640542725326
8 3.63970151119111
8.5 3.64250253471956
9 3.64247879596381
9.5 3.64003118285705
10 3.63670450855706
10.5 3.63413885787979
11 3.63455245488448
11.5 3.63719178528953
。。。。。。。
18.5 3.63397281939756
19 3.63569830040656
19.5 3.63887628526499
20 3.64241702048224
追问请问这个里面定义函数那个代码是怎么样的，谢谢
v0=[2.3,0];
h=0.5;
a=0;
b=20;
[t,v]=ode23(@fun,a:h:b,v0);
求解结果：
t V(t) dV(t)
0 2.3000 0
0.5000 5.2562 3.9856
1.0000 4.4588 -7.0559
1.5000 2.8150 6.5266
2.0000 5.5829 0.7497
2.5000 3.3187 -7.8804
3.0000 3.9351 7.8954
3.5000 5.4400 -2.5105
4.0000 2.4192 -3.2536
4.5000 4.9425 5.3320
5.0000 4.7890 -5.9263
5.5000 2.5304 4.4832
6.0000 5.4779 1.8393
6.5000 3.6658 -7.9648
7.0000 3.5824 7.9280
7.5000 5.4794 -1.6476
8.0000 2.5728 -4.7430
8.5000 4.7258 5.9901
9.0000 4.9100 -5.2930
9.5000 2.4462 3.2745
10.0000 5.3931 2.3321
10.5000 3.7920 -7.8113
11.0000 3.4505 7.7593
11.5000 5.4600 -1.3779
12.0000 2.6260 -5.0057
12.5000 4.6576 6.0813
13.0000 4.8927 -5.2044
13.5000 2.4588 3.1977
14.0000 5.3626 2.2714
14.5000 3.7309 -7.7293
15.0000 3.5055 7.6847
15.5000 5.4053 -1.6398
16.0000 2.5687 -4.3342
16.5000 4.7281 5.6945
17.0000 4.7518 -5.6085
17.5000 2.5561 4.1328
18.0000 5.3754 1.6801
18.5000 3.4960 -7.5912
19.0000 3.7339 7.5863
19.5000 5.2974 -2.4283
20.0000 2.4512 -2.6040
追问您好，初始条件只有速度V=2.3，没有加速度，刚刚搞错了。
能否麻烦你把定义函数的代码也一并说下，所有代码，谢谢追答用ode23（）函数可以较好求得微分方程的数值解。执行代码如下：
v0=0
h=0.5;a=0;b=20;
[t,v]=ode23(@fun,a:h:b,v0);
求解结果：
t V(t)
0 0
0.5 3.63716291516561
1 3.64158930086147
1.5 3.6428487523781
2 3.63885309188913
2.5 3.6339010775317
3 3.6342955704192
3.5 3.63569824585008
4 3.63869545971492
4.5 3.64207939462941
5 3.64279555495064
5.5 3.64076897493905
6 3.6375611660217
6.5 3.63456757976042
7 3.63417276743521
7.5 3.63640542725326
8 3.63970151119111
8.5 3.64250253471956
9 3.64247879596381
9.5 3.64003118285705
10 3.63670450855706
10.5 3.63413885787979
11 3.63455245488448
11.5 3.63719178528953
。。。。。。。
18.5 3.63397281939756
19 3.63569830040656
19.5 3.63887628526499
20 3.64241702048224
追问请问这个里面定义函数那个代码是怎么样的，谢谢
全部回答
- 1楼网友：话散在刀尖上
- 2021-03-19 10:55
ode45可以追问能否给个代码，谢谢
我要举报
如以上回答内容为低俗、色情、不良、暴力、侵权、涉及违法等信息，可以点下面链接进行举报！
点此我要举报以上问答信息
大家都在看
推荐资讯