1：问题
常微分方程的初值问题的标准数学表述为：y'=f(t,y),a<=t<=b,y(a)=y(0)；
我们要求解的任何高阶常微分方程都可以用替换法化为上式所示的一阶形式，其中 y 为向量，
yo 为初始值。
2：Matlab 中解决以上问题的步骤
（1）：化方程组为标准形式。
例如：y'''-3y''-y’y=0,y(0)=0,y'(0)=1,y''(0)=-1.
把微分方程的高阶导数写为低阶导数的算式，即：
y'''=3y''+y'y,设：y1=y,y2=y',y3=y'',则原方程化为下列等价的
方程组：
满足初值条件： 已把该方程化成了标准形式。
其中：y'->（y1’,y2’,y3’）,a->(0,0,0),y0->(0,1,-1),f(t,y)->(y2,y3,3y3+y2y1).
(2):把微分方程组编成 m 函数文件。
如：function dy=F(t,y)
dy=[y(2);y(3);3*y(3)+y(2)*y(1)];
注意：A:在函数文件里，虽然写微分方程时并不同时包含参数 t 和 y，但第一行必须包含这
两个输入变量。B:向量 dy 必须为列向量。
（3）：调用一个微分方程的求解函数求解。
[T,Y]=solver(‘F’,tspan,y0);
其中：solver:求解函数名；
F:包含微分方程的 m 文件;
tspan 为积分的数据范围，其格式为：[t0,tfinal];
y0 为 t0 时刻的初值列向量。
输出参数 T 和 Y 为列向量
T 为时刻向量。
Y 表是不同时刻的函数值。
3：（例）一个求解常微分方程初值问题的完整过程。
问题：求解方程 y’’-3(1-y^2)y’+y=0 在初值 y’(0)=3,y(0)=2 的解。
1 化成标准形式：
设 y1=y,y2=y’,则： 初值为：
2 编写函数文件 ode.m,内容为：
function dy=ode(t,y)
资源评论
pppptttpp6774
- 粉丝: 0
最新资源
- 数据库基础理论—体系结构ppt课件.ppt
- 王志良jian养老服务物联网.pptx
- 中国软件投融资高峰论坛演讲稿.pptx
- 互联网产品经理的九道思考题.pdf
- 第七章项目管理班子的人员配备、素质及管理经验.doc
- 用户综合布线手册.doc
- 软件详细设计专题讲座.pptx
- 因特网信息交流与网络安全摘要.ppt
- 生化药物和基因工程药物分析概念专题知识专家讲座.pptx
- 生物学网络考试系统市公开课金奖市赛课一等奖课件.pptx
- 网站运营数据分析之内容指标.docx
- 网络营销基础知识培训.pptx
- 汽车制造企业的OTD管理与信息化ERP建设.doc
- 数字图像处理-试卷.doc
- 计算机体系结构期末考试试题及答案.doc
- actionscript程序结构与语法.pptx
资源上传下载、课程学习等过程中有任何疑问或建议，欢迎提出宝贵意见哦~我们会及时处理！
点击此处反馈
