ZKX's LAB

全非线性抛物型偏微分方程 如何用matlab解二维的非线性偏微分方程组, 其中每个方程是抛物线型的

2020-07-21知识9

数学题的分类 数理逻辑与数学基础演绎逻辑学亦称符号逻辑学证明论亦称元数学递归论模型论公理集合论数学基础数理逻辑与数学基础其他学科数论初等。抛物型偏微分方程的抛物方程 。二阶线性偏微分方程(6)在区域Q内称为是抛物型的,如果存在常数α>;0,使得对于任意ξ∈Rn,(x1,x2,…,xn,t)∈Q 有。的形式。(7)称为具有散度形式的抛物型方程,(6)称为非散度形式的抛物型方程。时,(6)与(7)是有区别的,不能互推。如果方程(6)、(7)中的系数和右端还依赖于u,墷u,则(6)和(7)称为拟线性抛物型方程。抛物型方程和椭圆型方程的研究有相似的地方,它们互相影响、互为借鉴。椭圆型方程理论很多结果在抛物型方程中都有相应的定理,例如先验估计、极值原理等。怎样判断微分方程的线性与非线性 对于线性2113微分方程,其中只能出现函数本身,5261以及函数4102的任何阶次的导函数;函数本身跟所有的导1653函数之间除了加减之外,不可以有任何运算;函数本身跟本身、各阶导函数本身跟本身,都不可以有任何加减之外的运算;不允许对函数本身、各阶导函数做任何形式的复合运算,例如:siny、cosy、tany、lny、lgx、y2、y3。若一个微分方程不符合上面的条件,就是非线性微分方程。扩展资料线性方程:在代数方程中,仅含未知数的一次幂的方程称为线性方程。这种方程的函数图象为一条直线,所以称为线性方程。可以理解为:即方程的最高次项是一次的,允许有0次项,但不能超过一次。比如ax+by+c=0,此处c为关于x或y的0次项。微分方程:含有自变量、未知函数和未知函数的导数的方程称为微分方程。如果一个微分方程中仅含有未知函数及其各阶导数作为整体的一次幂,则称它为线性微分方程。可以理解为此微分方程中的未知函数y是不超过一次的,且此方程中y的各阶导数也应该是不超过一次的。参考资料-线性微分方程抛物型偏微分方程的拟线性蜕化 考虑在绝热过程中气体通抄过多孔介质的流动,这个过程可由下述方程来刻画:,式中m>;1,u是气体密度,通常研究它的非负解。由于当u=0时方程蜕化,因此它是一个拟线性蜕化抛物型方程。对于袭这个问题的系统理论研究是从 1957年开始的。解u的支集的边界是一条自由边界,通过自由边u一般不连续,因此这个方程知一般只存在在索伯列夫意义下的广义解,而且由于当u=0时方程蜕化为一阶方程,因此与热传导方程不同,扰动的传播速道度是有限的。请问具体如何区分,抛物型偏微分方程,双曲型偏微分方程,椭圆型偏微分方程? 依次是椭圆型,双曲型,双曲型AUxx+BUxy+CUyy+.=0Δ=B^2-4ACΔ=0:抛物型Δ>;0:双曲型Δ“数学”是一门什么样的学科 数学是研究数量、结构、变化、空间以及信息等概念的一门学科,从某种角度看属于形式科学的一种。数学家和哲学家对数学的确切范围和定义有一系列的看法。。怎样判断微分方程的线性与非线性 区别线性微分方程和非线性微分方程如下:1.微分方程中的线性,指的是y及其导数y'都是一次方。如y'=2xy。2.非线性,就是除了线性的。如y'=2xy^2。扩展资料微分方程指含有未知函数及其导数的关系式。解微分方程就是找出未知函数。微分方程是伴随着微积分学一起发展起来的。微积分学的奠基人Newton和Leibniz的著作中都处理过与微分方程有关的问题。微分方程的应用十分广泛,可以解决许多与导数有关的问题。参考资料:-微分方程如何用matlab解二维的非线性偏微分方程组, 其中每个方程是抛物线型的 MATLAB提供了两种方法解决PDE问题:一是pdepe()函数,它可以求解一般的PDEs,据用较大的通用性,但只支持命令行形式调用。二是PDE工具箱,可以求解特殊PDE问题,PDEtool有较大的局限性,比如只能求解二阶PDE问题,并且不能解决偏微分方程组,但是它提供了GUI界面,从繁杂的编程中解脱出来了,同时还可以通过File->;Save As直接生成M代码MATLAB语言提供了pdepe()函数,可以直接求解一般偏微分方程(组),它的调用格式为sol=pdepe(m,@pdefun,@pdeic,@pdebc,x,t)【输入参数】pdefun:是PDE的问题描述函数,它必须换成下面的标准形式这样,PDE就可以编写下面的入口函数[c,f,s]=pdefun(x,t,u,du)m,x,t就是对应于(式1)中相关参数,du是u的一阶导数,由给定的输入变量即可表示出出c,f,s这三个函数pdebc:是PDE的边界条件描述函数,必须先化为下面的形式于是边值条件可以编写下面函数描述为[pa,qa,pb,qb]=pdebc(x,t,u,du)其中a表示下边界,b表示下边界pdeic:是PDE的初值条件,必须化为下面的形式股我们使用下面的简单的函数来描述为u0=pdeic(x)m,x,t:就是对应于(式1)中相关参数【输出参数】sol:是一个三维数组,sol(:,:,i)表示ui的解,换句话说uk对应x(i)。求解二维抛物线型偏微分方程matlab程序 function[u,x,y,t]=TDE(A,D,T,ixy0,bxyt,Mx,My,N)解方程 u_t=c(u_xx+u_yy)for D(1)(2),D(3)(4),0初值:u(x,y,0)=ixy0(x,y)边界条件:u(x,y,t)=bxyt(x,y,t)for(x,y)cBMx/My:x轴和y轴的等分段数N:t 轴的等分段数dx=(D(2)-D(1))/Mx;x=D(1)+[0:Mx]*dx;dy=(D(4)-D(3))/My;y=D(3)+[0:My]'*dy;dt=T/N;t=[0:N]*dt;初始化ufor i=1:Mx+1for j=1:My+1u(i,j)=ixy0(x(i),y(j));endendrx=A*dt/(dx*dx);rx1=1+2*rx;rx2=1-2*rx;ry=A*dt/(dy*dy);ry1=1+2*ry;ry2=1-2*ry;for i=1:Mx-1%(11.2.21a)P(i,i)=ry1;if i>;1P(i-1,i)=-ry;P(i,i-1)=-ry;endendfor j=1:My-1%(11.2.21b)Q(j,j)=rx1;if j>;1Q(j-1,j)=-rx;Q(j,j-1)=-rx;endendfor k=1:Nu_1=u;t=k*dt;for i=1:Mx+1%边界条件u(i,1)=feval(bxyt,x(i),y(1),t);u(i,My+1)=feval(bxyt,x(i),y(My+1),t);endfor j=1:My+1u(1,j)=feval(bxyt,x(1),y(j),t);u(Mx+1,j)=feval(bxyt,x(Mx+1),y(j),t);endif mod(k,2)=0for i=2:Mxj=2:My;bx=[ry*u(i,1)zeros(1,Mx-3)ry*u(i,My+1)]+rx*(u_1(i-1,j)+u_1(i+1,j))+rx2*u_1(i,j);u(i,j)=linsolve(P,bx');(11.2.21a)endelsefor j=2:Myi=2:Mx;by=[rx*u(1,j);zeros(My-3,1);rx*u。如何用matlab解二维的非线性偏微分方程组, 其中每个方程是抛物线型的 如何用matlab解二维的非线性偏微分方程组,其中每个方程是抛物线型的 MATLAB提供两种解决PDE问题:pdepe()函数求解般PDEs据用较通用性支持命令行形式调用 二PDE工具箱求解。

#线性#数学#二维#matlab函数#matlab

随机阅读

qrcode
访问手机版