版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、实验六 用matlab求解常微分方程1微分方程的概念未知的函数以及它的某些阶的导数连同自变量都由一方程联系在一起的方程称为微分方程。如果未知函数是一元函数,称为常微分方程。常微分方程的一般形式为F(t,y,y',W, ,y(n) 0如果未知函数是多元函数,成为偏微分方程。联系一些未知函数的一组微分方程组称为微分方程组。微分方程中出现的未知函数的导数的最高阶解数称为微分方程的阶。假设方程中未知函数及其各阶导数都是一次的,称为线性常微分方程,一般表示为y(n) a!(t)y(n 1) an i(t)y' an(t)y b(t)假设上式中的系数 ai(t),i1,2,n均与t无关,称
2、之为常系数。2 常微分方程的解析解dy 彳y 1有些微分方程可直接通过积分求解例如,一解常系数常微分方程dt可化为虫dtty 1,两边积分可得通解为y ce 1.其中c为任意常数.有些常微分方程可用一些技巧,如别离变量法,积分因子法,常数变异法,降阶法等可化为可积分的方程而求得解析解线性常微分方程的解满足叠加原理,从而他们的求解可归结为求一个特解和相应齐次微 分方程的通解一阶变系数线性微分方程总可用这一思路求得显式解。高阶线性常系数微分 方程可用特征根法求得相应齐次微分方程的根本解,再用常数变异法求特解。一阶常微分方程与高阶微分方程可以互化,已给一个n阶方程y(n)f(t,y',y&q
3、uot;,y(n1)(n 1)设y1y,y2y, ,yn y,可将上式化为一阶方程组%' y2y2 y y'n 1 ynyn' f(t,y1,y2,yn)反过来,在许多情况下,一阶微分方程组也可化为高阶方程。所以一阶微分方程组与高 阶常微分方程的理论与方法在许多方面是相通的,一阶常系数线性微分方程组也可用特征根法求解。3 微分方程的数值解法除常系数线性微分方程可用特征根法求解,少数特殊方程可用初等积分法求解外,大局部微分方程无限世界,应用中主要依靠数值解法。考虑一阶常微分方程初值问题y'(t)f(t,y(t),t° t tfy(to)yo其中 y(yi
4、,y2,ym)',f(fl,f2,fm)',y0(y10 ,y20 ,ym0)'.所谓数值解法,就是寻求y(t)在一系列离散节点to tltn tf上的近似值yk,k 0,1, ,n称hktk 1 tk为步长,通常取为常量 h。最简单的数值解法是Euler法。Euler法的思路极其简单:在节点出用差商近似代替导数y(tk 1) y(tk)y (tk)h这样导出计算公式称为 Euler格式yk i yk hf(tk,yk),k o,1,2,他能求解各种形式的微分方程。Euler法也称折线法。Euler方法只有一阶精度,改进方法有二阶Runge-Kutta法、四阶Runge
5、-Kutta法、五阶Runge-Kutta-Felhberg法和先行多步法等,这些方法可用于解高阶常微分方程组初值问题。边值问题采用不同方法,如差分法、有限元法等。数值算法的主要缺点是它缺乏物理理解。4.解微分方程的 MATLAB 命令MATLAB中主要用dsolve求符号解析解,ode45,ode23,ode15s求数值解。s=dsolve(方程1'方程2初始条件1','初始条件2',自变量) 用字符串方程表示,自变量缺省值为t。导数用D表示,2阶导数用D2表示,以此类推。S返回解析解。在方程组情形,s为一个符号结构。tout,yout=ode45( ypri
6、me ' ,t采用变步长四阶Runge-Kutta 法和五阶Runge-Kutta-Felhberg法求数值解,yprime是用以表示f(t,y)的M文 件名,t0表示自变量的初始值,tf表示自变量的终值,y0表示初始向量值。 输出向量tout表示节点(t0,t1,tn)T输出矩阵yout表示数值解,每一列对 应y的一个分量。假设无输出参数,那么自动作出图形。ode45是最常用的求解微分方程数值解的命令,对于刚性方程组不宜采用。ode23与ode45类似,只是精度低一些。ode12s用来求解刚性方程组,是用格式同ode45。可以用helpdsolve, help ode45查阅有关这些
7、命令的详细信息例1求以下微分方程的解析解们 y' ay b2 y'' sin(2x) y,y(0) 0,y'(0) 13 f' f g,g' g f, f'(0) 1,g'(0) 1 方程 1求解的 MATLAB 代码为:>>clear;>>s=dsolve('Dy=a*y+b') 结果为s =-b/a+exp(a*t)*C1 方程 2求解的 MATLAB 代码为:>>clear; >>s=dsolve('D2y=sin(2*x)-y','y(
8、0)=0','Dy(0)=1','x') >>simplify(s) % 以最简形式显示 s 结果为s =(-1/6*cos(3*x)-1/2*cos(x)*sin(x)+(-1/2*sin(x)+1/6*sin(3*x)*cos(x)+5/3*sin(x) ans =-2/3*sin(x)*cos(x)+5/3*sin(x) 方程 3求解的 MATLAB 代码为:>>clear;>>s=dsolve('Df=f+g','Dg=g-f','f(0)=1','g(0
9、)=1')>>simplify(s.f) %s 是一个结构 >>simplify(s.g)结果为ans =exp(t)*cos(t)+exp(t)*sin(t)ans =-exp(t)*sin(t)+exp(t)*cos(t)例2求解微分方程y' y t 1, y(0) 1,先求解析解,再求数值解,并进行比拟。由>>clear;>>s=dsolve('Dy=-y+t+1','y(0)=1','t') >>simplify(s)t 可得解析解为 y t e 。下面再求其数值
10、解,先编写 M 文件%M函数function f=fun8(t,y) f=-y+t+1;再用命令>>clear; close; t=0:0.1:1;>>y=t+exp(-t); plot(t,y); %化解析解的图形>>hold on; %保存已经画好的图形,如果下面再画图,两个图形和并在一起>>t,y=ode45('fun8',0,1,1);>>plot(t,y,'ro'); %画数值解图形,用红色小圈画>>xlabel('t'),ylabel('y')结果
11、见图1.41.351.31.251.151.11.051t图解析解与数值解y 1.2由图可见,解析解和数值解吻合得很好。 例3求方程ml "mgsin , (0)0, '(0) 0的数值解不妨取11,g9.8, (0)15 .那么上面方程可化为MATLAB代码"9.8si n , (0)15, '(0)0运行先看看有没有解析解>>clear;>>s=dsolve('D2y=9.8*si n(y)','y(0)=15','Dy(0)=0','t')>>simpl
12、ify(s)知原方程没有解析解下面求数值解令可将原方程化为如下方程组y2y 9.8si n(y1)y1(0)15, y2(0)0建立M文件如下%M文件fun cti on f=fun 9(t,y)f=y(2), 9.8*s in (y(1)' %f向量必须为一列向量运行MATLAB代码>>clear; close;>>t,y=ode45('fu n9',0,10,15,0);>>plot(t,y(:,1); %画随时间变化图,y(:2) 那么表示'的值>>xlabel('t'),ylabel('y1')结果见图由图可见,15067图7.2数值解随时间t周期变化。习题16-61 求以下微分方程的解析解(1) T +召亠一3 = 2产鈕 j*+= sinx (o >0)(4) 2 -1 = 0 h血+ 2(b-W妙=0,=13(5) y + / + y=cosx
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 八下历史第三课的教育课件
- 5年中考3年模拟试卷初中道德与法治八年级下册02期末素养综合测试
- 不定时工作制实施方案范本(四篇)
- 四年级下册音乐教案
- 古韵悠长 第一课《雕梁画栋》课件 2024-2025学年人美版(2024)初中美术七年级上册
- (统考版)2023版高考化学一轮复习课时作业10碳、硅及其化合物
- 工艺品交易桥梁
- 互联网数据中心运输协议
- 塑料原料购销运输合同范本
- 办公室改造追加协议模板
- 我国农产品品牌价值及品牌战略管理研究
- 技能认证变配电运行值班员中级考试(习题卷7)
- 特殊学生心理健康基本档案表
- 泌尿系结石的诊治课件
- 自然资源学概论期末考试试卷
- 小学语文-部编本三年级上册《道德与法治》第10课《父母多爱我》教学设计学情分析教材分析课后反思
- 河南省高校毕业生基层就业国家助学贷款代偿申请表
- 全国住户收支调查业务知识考试复习题库(含答案)
- 市中医医院推行6S精益管理实施方案
- rmk1608电阻型号命名规则
- 一种40CrNiMoA钢棒的热处理方法及应用与流程
评论
0/150
提交评论