- 1、本文档共46页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
题目:用MATLAB求解偏微分方程主讲人:班级 :时间 :基础知识预习微分方程的MATLAB求解包含1:常微分方程的MATLAB求解(上 节课已经讲过)这里不再赘述。2:偏微分方程的MATLAB求解(本 次教学内容)偏微分方程概念偏微分方程(Partial Differential Equation,简称PDE)指含有未知函数及其偏导数的方程。描述自变量、未知函数及其偏导数之间的关系。偏微分方程分为①线性偏微分方程式与②非线性偏微分方程式,常常有几个解而且涉及额外的边界条件。常微分方程:在微分方程中,若自变量的个数只有一个的微分方程。偏微分方程:自变量的个数有两个或两个以上的微分方程。求解偏微分方程的方法求解偏微分方程的数值方法:1. 有限元法(Finite Element Method, FEM)--- hp-FEM2. 有限体积法(Finite Volume Method, FVM)3. 有限差分法(Finite Difference Method, FDM)。其它:广义有限元法(Generalized Finite Element Method, FFEM)、扩展有限元法(eXtended Finite Element Method, XFEM)、无网格有限元法(Meshfree Finite Element Method)、离散迦辽金有限元法(Discontinuous Galerkin Finite Element Method, DGFEM)等。MATLAB解偏微分方程MATLAB提供了两种方法解决PDE 问题:①pdepe()函数,它可以求解一般的PDEs,具有较大的通用性,但只支持命令行形式调用。②PDE 工具箱,可以求解特殊PDE 问题,PDEtool 有较大的局限性,比如只能求解二阶PDE 问题,并且不能解决偏微分方程组,但是它提供了GUI界面,从繁杂的编程中解脱出来了,同时还可以通过File-Save As直接生成M代码使用pdeval()直接计算某个点的函数值???一般偏微分方程组(PDEs)的MATLAB求解直接求解一般偏微分方程(组),它的调用格式为sol=pdepe(m,@pdefun,@pdeic,@pdebc,x,t)输出参数问题描述函数初值条件边界条件参数自变量【输入参数】(1)@pdefun:是PDE 的问题描述函数,它必须换成下面的标准形式PDE 就可以编写下面的入口函数[c,f,s]=pdefun(x,t,u,du)m,x,t就是对应于(式1)中相关参数和自变量,du是u的一阶导数,由给定的输入变量即可表示出出c,f,s这三个函数【输入参数】(2)@pdeic:是PDE 的初值条件,必须化为下面的形式我们使用下面的简单的函数来描述为u0=pdeic(x)【输入参数】(3)@pdebc:是PDE的边界条件描述函数,必须先化为下面的形式于是边值条件可以编写下面函数描述为[pa,qa,pb,qb]=pdebc(x,t,u,du)其中a 表示下边界,b 表示下边界【输入参数】(4)m:就是对应于(式1)中相关参数x,t:就是对应于(式1)中自变量【输出参数】sol:是一个三维数组,sol(:,:,i)表示ui的解,换句话说uk对应x(i)和t(j)时的解为sol(i,j,k)实例讲解(题目)例:初值条件边界条件实例讲解(解法)【解】第一步根据(1)对照给出的偏微分方程,则原方程可以改写为输入参数(1’)目标PDE函数%% 目标PDE函数function [c,f,s]=pdefun (x,t,u,du)c=[1;1];f=[0.024*du(1);0.17*du(2)];temp=u(1)-u(2);s=[-1;1].*(exp(5.73*temp)-exp(-11.46*temp));输入参数(2’)初值条件初值条件改写为%% 初值条件函数function u0=pdeic(x)u0=[1;0];输入参数(3’)边界条件边界条件改写为%% 边界条件函数function [pa,qa,pb,qb]=pdebc(xa,ua,xb,ub,t)%a表示左边界,b表示右边界pa=[0;ua(2)];qa=[1;0];pb=[ub(1)-1;0];qb=[0;1];(4)主调函数clcx=0:0.05:1;t=0:0.05:2;m=0;sol=pdepe(m,@pdefun,@pdeic,@pdebc,x,t);figure(numbertitle,off,name,PDE Demo——by Matlabsky)%创建个窗口,窗口名字是name后边的名字NumberTitle,off是关掉默认显示名字。subplot(211)surf(x,t,sol(:,:,1))%sol(:,
您可能关注的文档
- it项目管理课程介绍17.pptx
- 合同法律风险防范知识讲座.pptx
- 12浮力及应用解析.pptx
- 丁涛络能源管理的三大纪律八注意.pptx
- 431角的定义与表示方法.pptx
- 一偏导数的定义及其计算法.pptx
- 二年级上册音乐大鹿.pptx
- 中医眼科学——暴盲.pptx
- matlab数据和函数的可视化.pptx
- snp分型技术简介.pptx
- 自身免疫性肝炎研究进展2025.pdf
- 比亚迪E5 课件全套 项目1--7 纯电动汽车的认识与安全操作 ---- 空调、转向故障.pdf
- 原发性胆汁性胆管炎 2025.pdf
- 2024 在中国的美国企业特别报告.pdf
- 智能财税轮岗实务 课件项目四.pptx
- 新世纪放眼中国德语高级阅读教程教学课件Lektion 11_Der Chinesische Traum - Wunsch nach dem Wiederaufleben.pptx
- 新世纪放眼中国德语高级阅读教程教学课件Lektion 7_Chinas Raumfahrtwesen - Langer Marsch und Himmelpalast.pptx
- 新世纪放眼中国德语高级阅读教程教学课件Lektion 6_Wirtschaft - Made in China mit Innovation.pptx
- 新世纪放眼中国德语高级阅读教程教学课件Lektion 2_Chinesisch - Aussprache und Schriftzeichen.pptx
- 新世纪放眼中国德语高级阅读教程教学课件Lektion 12_Internationaler Austausch - Zusammenarbeit und gemeinsame Entwicklung.pptx
文档评论(0)