第九章 常微分方程数值解1_百度文库
两大类热门资源免费畅读
续费一年阅读会员，立省24元！
文档贡献者贡献于
评价文档：
151页免费72页免费50页免费33页免费29页免费 45页免费47页免费19页免费13页免费16页免费
喜欢此文档的还喜欢1页免费4页免费4页免费4页免费65页1下载券
第九章 常微分方程数值解1|常​微​分​方​程​数​值​解​的​典​型​用​法
把文档贴到Blog、BBS或个人站等：
普通尺寸(450*500pix)
较大尺寸(630*500pix)
大小：528.50KB
登录百度文库，专享文档复制特权，财富值每天免费拿！
你可能喜欢matlab解微分方程 matlab 微分方程 matlab解微分方程组 用matlab解微分方程 matla..
扫扫二维码，随身浏览文档
手机或平板扫扫即可继续访问
常微分方程初值问题的MATLAB解法
举报该文档为侵权文档。
举报该文档含有违规或不良信息。
反馈该文档无法正常浏览。
举报该文档为重复文档。
推荐理由：
将文档分享至：
分享完整地址
文档地址：
粘贴到BBS或博客
flash地址：
支持嵌入FLASH地址的网站使用
html代码：
&embed src='/DocinViewer-4.swf' width='100%' height='600' type=application/x-shockwave-flash ALLOWFULLSCREEN='true' ALLOWSCRIPTACCESS='always'&&/embed&
450px*300px480px*400px650px*490px
支持嵌入HTML代码的网站使用
您的内容已经提交成功
您所提交的内容需要审核后才能发布，请您等待！
3秒自动关闭窗口常微分方程组初值问题数值解的实现和算法分析_百度文库
两大类热门资源免费畅读
续费一年阅读会员，立省24元！
评价文档：
83页2下载券36页2下载券22页免费4页¥2.0028页免费 26页1下载券22页7下载券18页免费2页免费83页4下载券
喜欢此文档的还喜欢4页免费6页免费26页1下载券36页2下载券33页免费
常微分方程组初值问题数值解的实现和算法分析|
把文档贴到Blog、BBS或个人站等：
普通尺寸(450*500pix)
较大尺寸(630*500pix)
你可能喜欢MATLAB中，编写ode45求解单摆运动的一阶常微分方程组_百度知道
MATLAB中，编写ode45求解单摆运动的一阶常微分方程组
%函数文件（保存yfun2.m）function
exer2=yfun2(t,y);g=9.8;l=25;exer2=[x(2);(-g/l)*sin(x(1))];%程序文件tspan=[0,10];y0=[0.1475,0];[t,y]=ode45(‘yfun2’,tspan,y0);plot(t,y(:,1));请问这个程序错在哪?我们只上了四节课老师便要我们自己编写程序，所以很多不懂
我有更好的答案
按默认排序
exer2=[x(2);(-g/l)*sin(x(1))]; 敢问x(2)在这个子程序中是什么
这个是根据网上的例题跟着写的格式，就是把单摆运动的方程化为一阶微分方程组后的数& & & & & & & & & & & & & & &&
其他类似问题
ode45的相关知识
等待您来回答
下载知道APP
随时随地咨询
出门在外也不愁一阶微分方程求解，并优化参数值 - 计算模拟 - 小木虫 - 学术 科研 第一站
&& 查看话题
一阶微分方程求解，并优化参数值
紧急求助，非常感谢！！！
我有一系列y-t值，如何计算dy/dt=a*y*(b-y)/(c*d+d*y+e*(b-y)+y*(b-y)).
或者说我已知a，b，c，d，e的值如何根据时间t的变化求出Y的值，然后依据求得的y-t值来拟合y-t值求得最佳的参数值a，c，d，e。
十分感谢，现在对这个相当困惑。希望大家指导一下。
微分方程拟合可以试下1stOpt，简单易上手。有数据的话贴上来可以帮你试试。 : Originally posted by dingd at
微分方程拟合可以试下1stOpt，简单易上手。有数据的话贴上来可以帮你试试。 嗯，我现在只是这个思路，还没有数据，而且我以后还会加参数，进行拟合，会很麻烦，唉，我学生物的，做反应动力学，但是不懂计算，我在网上下到了1.5破解版，是不是没有这个功能，我以后一两个月会经常用，也不好一直让你弄吧。。。如果买1stopt4.0或5.0版需要多少钱？ 初步看了一下，你这个其实不难，用matlab，mathtype，或者是自己写一个小程序（具体你可以参考任何一本数值分析的书，如果不需要太高的精度的话，你就用欧拉方法）就可以解决你的第一个问题。 想将y-t值做数值微分得到dy/dt-t的曲线，然后用最小二乘法拟合该曲线，得到abcde。 : Originally posted by wjx0162 at
想将y-t值做数值微分得到dy/dt-t的曲线，然后用最小二乘法拟合该曲线，得到abcde。 嗯，谢谢，版主已经帮我了，谢谢。 : Originally posted by 六本木 at
初步看了一下，你这个其实不难，用matlab，mathtype，或者是自己写一个小程序（具体你可以参考任何一本数值分析的书，如果不需要太高的精度的话，你就用欧拉方法）就可以解决你的第一个问题。 嗯，谢谢，版主已经帮我解决了，谢谢 : Originally posted by dingd at
微分方程拟合可以试下1stOpt，简单易上手。有数据的话贴上来可以帮你试试。 你好，你能帮我跑一下程序吗？非常感谢！！！
方程dy/dt = V*y*(2-y)/(Kia*Kmb+Kmb*y+Kma*(2-y)+y*(2-y))
0& && &&&0.
0.167& && &&&0.
0.333& && &&&0.
0.5& && &&&0.
0.667& && &&&0.
0.833& && &&&0.
1& && &&&0.
1.167& && &&&0.
1.333& && &&&0.
1.5& && &&&0.
1.667& && &&&0.
1.833& && &&&0.
2& && &&&0.
2.167& && &&&0.
2.333& && &&&0.
2.5& && &&&0.
2.667& && &&&0.
2.833& && &&&0.1085209
3& && &&&0.
3.167& && &&&0.
3.333& && &&&0.
3.5& && &&&0.
3.667& && &&&0.
3.833& && &&&0.
4& && &&&0.
4.167& && &&&0.
4.333& && &&&0.1
4.5& && &&&0.
4.667& && &&&0.
4.833& && &&&0.
5& && &&&0.
5.167& && &&&0.
5.333& && &&&0.
5.5& && &&&0.
5.667& && &&&0.
5.833& && &&&0.
6& && &&&0.
6.167& && &&&0.
6.333& && &&&0.
6.5& && &&&0.
6.667& && &&&0.
6.833& && &&&0.
7& && &&&0.
7.167& && &&&0.
7.333& && &&&0.
7.5& && &&&0.
7.667& && &&&0.
7.833& && &&&0.
8& && &&&0.
8.167& && &&&0.
8.333& && &&&0.
8.5& && &&&0.
8.667& && &&&0.
8.833& && &&&0.
9& && &&&0.
9.167& && &&&0.
9.333& && &&&0.
9.5& && &&&0.
9.667& && &&&0.
9.833& && &&&0.
10& && &&&0.
10.167& && &&&0.
10.333& && &&&0.
10.5& && &&&0.
10.667& && &&&0.
10.833& && &&&0.
11& && &&&0.
11.167& && &&&0.
11.333& && &&&0.
11.5& && &&&0.
11.667& && &&&0.
11.833& && &&&0.
12& && &&&0.
12.167& && &&&0.
12.333& && &&&0.
12.5& && &&&0.
12.667& && &&&0.
12.833& && &&&0.
13& && &&&0.
13.167& && &&&0.
13.333& && &&&0.
13.5& && &&&0.
13.667& && &&&0.
13.833& && &&&0.
14& && &&&0.
14.167& && &&&0.
14.333& && &&&0.
14.5& && &&&0.
14.667& && &&&0.0585209
14.833& && &&&0.
15& && &&&0.
15.167& && &&&0.
15.333& && &&&0.
15.5& && &&&0.
15.667& && &&&0.
15.833& && &&&0.
16& && &&&0.
16.167& && &&&0.
16.333& && &&&0.
16.5& && &&&0.
16.667& && &&&0.
16.833& && &&&0.
17& && &&&0.
17.167& && &&&0.
17.333& && &&&0.
17.5& && &&&0.
17.667& && &&&0.
17.833& && &&&0.05
18& && &&&0.
18.167& && &&&0.
18.333& && &&&0.
18.5& && &&&0.
18.667& && &&&0.
18.833& && &&&0.
19& && &&&0.
19.167& && &&&0.
19.333& && &&&0.
19.5& && &&&0.
19.667& && &&&0.
19.833& && &&&0.
20& && &&&0.
20.167& && &&&0.
20.333& && &&&0.
20.5& && &&&0.
20.667& && &&&0.
20.833& && &&&0.
21& && &&&0.
21.167& && &&&0.
21.333& && &&&0.
21.5& && &&&0.
21.667& && &&&0.
21.833& && &&&0.
22& && &&&0.
22.167& && &&&0.
22.333& && &&&0.
22.5& && &&&0.
22.667& && &&&0.
22.833& && &&&0.
23& && &&&0.
23.167& && &&&0.
23.333& && &&&0.
23.5& && &&&0.
23.667& && &&&0.
23.833& && &&&0.
24& && &&&0.
24.167& && &&&0.
24.333& && &&&0.
24.5& && &&&0.
24.667& && &&&0.
24.833& && &&&0.
25& && &&&0.
25.167& && &&&0.
25.333& && &&&0.
25.5& && &&&0.
25.667& && &&&0.
25.833& && &&&0.
26& && &&&0.
26.167& && &&&0.
26.333& && &&&0.
26.5& && &&&0.
26.667& && &&&0.
26.833& && &&&0.
27& && &&&0.
27.167& && &&&0.
27.333& && &&&0.
27.5& && &&&0.
27.667& && &&&0.
27.833& && &&&0.
28& && &&&0.
28.167& && &&&0.
28.333& && &&&0.
28.5& && &&&0.
28.667& && &&&0.
28.833& && &&&0.
29& && &&&0.
29.167& && &&&0.
29.333& && &&&0.
29.5& && &&&0.
29.667& && &&&0.
29.833& && &&&0.
30& && &&&0.
30.167& && &&&0.
30.333& && &&&0.
30.5& && &&&0.
30.667& && &&&0.
30.833& && &&&0.
31& && &&&0.
31.167& && &&&0.
31.333& && &&&0.
31.5& && &&&0.
31.667& && &&&0.
31.833& && &&&0.
32& && &&&0.
32.167& && &&&0.
32.333& && &&&0.
32.5& && &&&0.
32.667& && &&&0.
32.833& && &&&0.
33& && &&&0.
33.167& && &&&0.
33.333& && &&&0.
33.5& && &&&0.
33.667& && &&&0.
33.833& && &&&0.
34& && &&&0.
34.167& && &&&0.
34.333& && &&&0.
34.5& && &&&0.
34.667& && &&&0.
34.833& && &&&0.
35& && &&&0.
35.167& && &&&0.
35.333& && &&&0.
35.5& && &&&0.
35.667& && &&&0.
35.833& && &&&0.
36& && &&&0.
36.167& && &&&0.
36.333& && &&&0.
36.5& && &&&0.
36.667& && &&&0.
36.833& && &&&0.
37& && &&&0.
37.167& && &&&0.
37.333& && &&&0.
37.5& && &&&0.
37.667& && &&&0.
37.833& && &&&0.
38& && &&&0.
38.167& && &&&0.
38.333& && &&&0.
38.5& && &&&0.
38.667& && &&&0.
38.833& && &&&0.
39& && &&&0.
39.167& && &&&0.
39.333& && &&&0.
39.5& && &&&0.
39.667& && &&&0.
39.833& && &&&0.
40& && &&&0.
0& && &&&0.
0.167& && &&&0.
0.333& && &&&0.
0.5& && &&&0.
0.667& && &&&0.
0.833& && &&&0.
1& && &&&0.
1.167& && &&&0.
1.333& && &&&0.
1.5& && &&&0.
1.667& && &&&0.
1.833& && &&&0.
2& && &&&0.
2.167& && &&&0.
2.333& && &&&0.
2.5& && &&&0.
2.667& && &&&0.
2.833& && &&&0.
3& && &&&0.
3.167& && &&&0.
3.333& && &&&0.
3.5& && &&&0.
3.667& && &&&0.
3.833& && &&&0.
4& && &&&0.
4.167& && &&&0.
4.333& && &&&0.
4.5& && &&&0.
4.667& && &&&0.
4.833& && &&&0.
5& && &&&0.
5.167& && &&&0.
5.333& && &&&0.
5.5& && &&&0.
5.667& && &&&0.
5.833& && &&&0.
6& && &&&0.
6.167& && &&&0.
6.333& && &&&0.
6.5& && &&&0.
6.667& && &&&0.
6.833& && &&&0.
7& && &&&0.
7.167& && &&&0.
7.333& && &&&0.
7.5& && &&&0.
7.667& && &&&0.
7.833& && &&&0.
8& && &&&0.
8.167& && &&&0.
8.333& && &&&0.
8.5& && &&&0.
8.667& && &&&0.
8.833& && &&&0.
9& && &&&0.
9.167& && &&&0.
9.333& && &&&0.
9.5& && &&&0.
9.667& && &&&0.
9.833& && &&&0.
10& && &&&0.
10.167& && &&&0.
10.333& && &&&0.
10.5& && &&&0.
10.667& && &&&0.
10.833& && &&&0.
11& && &&&0.
11.167& && &&&0.
11.333& && &&&0.
11.5& && &&&0.
11.667& && &&&0.
11.833& && &&&0.
12& && &&&0.
12.167& && &&&0.
12.333& && &&&0.
12.5& && &&&0.
12.667& && &&&0.
12.833& && &&&0.
13& && &&&0.
13.167& && &&&0.
13.333& && &&&0.
13.5& && &&&0.
13.667& && &&&0.
13.833& && &&&0.
14& && &&&0.
14.167& && &&&0.
14.333& && &&&0.
14.5& && &&&0.
14.667& && &&&0.
14.833& && &&&0.
15& && &&&0.
15.167& && &&&0.
15.333& && &&&0.
15.5& && &&&0.
15.667& && &&&0.
15.833& && &&&0.
16& && &&&0.
16.167& && &&&0.
16.333& && &&&0.
16.5& && &&&0.
16.667& && &&&0.
16.833& && &&&0.
17& && &&&0.
17.167& && &&&0.
17.333& && &&&0.
17.5& && &&&0.
17.667& && &&&0.
17.833& && &&&0.
18& && &&&0.
18.167& && &&&0.
18.333& && &&&0.
18.5& && &&&0.
18.667& && &&&0.
18.833& && &&&0.
19& && &&&0.
19.167& && &&&0.
19.333& && &&&0.
19.5& && &&&0.
19.667& && &&&0.
19.833& && &&&0.
20& && &&&0.
20.167& && &&&0.0085209
20.333& && &&&0.
20.5& && &&&0.
20.667& && &&&0.0085209
20.833& && &&&0.0085209
21& && &&&0.0085209
21.167& && &&&0.
21.333& && &&&0.
21.5& && &&&0.
21.667& && &&&0.
21.833& && &&&0.
22& && &&&0.
22.167& && &&&0.
22.333& && &&&0.
22.5& && &&&0.
22.667& && &&&0.
22.833& && &&&0.
23& && &&&0.
23.167& && &&&0.
23.333& && &&&0.
23.5& && &&&0.
23.667& && &&&0.
23.833& && &&&0.
24& && &&&0.
24.167& && &&&0.
24.333& && &&&0.
24.5& && &&&0.
24.667& && &&&0.
24.833& && &&&0.
25& && &&&0.
25.167& && &&&0.
25.333& && &&&0.
25.5& && &&&0.
25.667& && &&&0.
25.833& && &&&0.
26& && &&&0.
26.167& && &&&0.
26.333& && &&&0.
26.5& && &&&0.
26.667& && &&&0.
26.833& && &&&0.
27& && &&&0.
27.167& && &&&0.
27.333& && &&&0.
27.5& && &&&0.
27.667& && &&&0.
27.833& && &&&0.
28& && &&&0.
28.167& && &&&0.
28.333& && &&&0.
28.5& && &&&0.
28.667& && &&&0.
28.833& && &&&0.
29& && &&&0.
29.167& && &&&0.
29.333& && &&&0.
29.5& && &&&0.
29.667& && &&&0.
29.833& && &&&0.
30& && &&&0. : Originally posted by dingd at
微分方程拟合可以试下1stOpt，简单易上手。有数据的话贴上来可以帮你试试。 恳请帮助，谢谢 : Originally posted by dingd at
微分方程拟合可以试下1stOpt，简单易上手。有数据的话贴上来可以帮你试试。 V=,Kia=,kmb=,kma=
这是约束条件。
恳请帮助，谢谢 供参考：
均方差(RMSE): 0.51791
残差平方和(SSE): 0.157
相关系数(R): 0.062
相关系数之平方(R^2): 0.637
参数& && && & & & & & 最佳估算
--------------------& & & & -------------
v& & & & -0.017
kia& & & & 0.211
kmb& & & & 19.
kma& & & & 4.5E-5 : Originally posted by dingd at
均方差(RMSE): 0.51791
残差平方和(SSE): 0.157
相关系数(R): 0.062
相关系数之平方(R^2): 0.637
参数& && && & & & & & 最佳估算
----------------- ... 谢谢，谢谢 这个用不着数值解，因为可得到解析解。经整理可得到：
{1+/*dy=a*dt
令/ ≡ α/y+β/(y-b)
整理后经比较系数后可得：α=e+c*d/b; β=-d-c*d/b
故有：｛1+/y-/(y-b)}dy=a*dt
故：&&y+*lny-*ln(y-b)=a*t+ A ,&&A为积分常数。
令δ（t）=y+*lny-*ln(y-b)-a*t-A ,&&
l令Δ=Σ^2,&&(1=1~n),n为实验此次数。
令 Δ分别对a、b、c、d、e、A求偏导数并令其为零，得到五个方程，求出五个未知数后代入原方程中即得结果。
解题完毕。
如果楼主觉得没用就算了。若觉得还行，还有用的话，就请送朵红花吧。 可以很方便地由分离变量法解出：
&&y^2/2+(c*d+b*e)/b*ln(y)-(1+c/b)*d*ln(y-b)=a*t+C
其中C为积分常数，由方程的附加条件确定。
令f(t)=y^2/2+(c*d+b*e)/b*ln(y)-(1+c/b)*d*ln(y-b)-a*t-C。设共测有n组（ti,yi)数据，i=1~n。再令Δ(t)=Σ^2 ,( i=1~n)
由Δ(t)分别对a、b、c、d求偏导数并令为零，得到四个方程。求解这个非线性方程组得到的a、b、c、d的值，这就是所要求的拟合值。 : Originally posted by peterflyer at
可以很方便地由分离变量法解出：
&&y^2/2+(c*d+b*e)/b*ln(y)-(1+c/b)*d*ln(y-b)=a*t+C
其中C为积分常数，由方程的附加条件确定。
令f(t)=y^2/2+(c*d+b*e)/b*ln(y)-(1+c/b)*d*ln(y-b)-a*t-C。设共测有n组（ti,yi) ... 楼上方程可化简为&&y^2/2+c*ln(y)-d*ln(y-b)=a*t+f;
对第一组数据，有
a = -6.77E-12
c = 0.6797
f = 0.0755
对第二组数据，有
a = 1.67E-11
c = 0.358757
var cpro_id = 'u1216994';
欢迎监督和反馈：本帖内容由
提供，小木虫为个人免费站点，仅提供交流平台，不对该内容负责。欢迎协助我们监督管理，共同维护互联网健康，如果您对该内容有异议，请立即发邮件到
联系通知管理员，也可以通过QQ周知，我们的QQ号为：8835100
我们保证在1个工作日内给予处理和答复，谢谢您的监督。
小木虫，学术科研第一站，为中国学术科研研究提供免费动力
欢迎监督，发现不妥请立即
E-mail： & QQ：8835100