ZKX's LAB

二维抛物型方程matlab程序 抛物型方程的数值求解,用C语言或MATLAB(不使用工具箱)

2020-07-19知识7

如何用matlab解二维的非线性偏微分方程组, 其中每个方程是抛物线型的 如何用matlab解二维的非线性偏微分方程组,其中每个方程是抛物线型的 MATLAB提供两种解决PDE问题:pdepe()函数求解般PDEs据用较通用性支持命令行形式调用 二PDE工具箱求解。用matlab求解抛物型方程,急啊!!用最简隐格式(向后差分格式)求解抛物型方程 你的精确定绝对有问题。你自己将精确解代入那个泛定方程,或者初值都不符的。一维热传导方程的差分格式k=1/16;xleft=0;xright=1;tend=0.2;时间终值dx=0.1;dt=0.05;n=(xright-xleft)/dx;x=xleft:dx:xright;beta=k*dt/dx/dx;A=diag((1+2*beta*ones(n+1,1)))+diag(-beta*ones(n,1),1)+diag(-beta*ones(n,1),-1);Q=dt/gou/c*ones(n+1,1);边界条件A(1,1)=1;A(1,2)=0;A(end,end)=1;A(end,end-1)=0;T0=25*log(2*pi*x(:));Tseriers=T0;leg_info{1}='t=0';T=T0;i=1;for t=0:dt:tendi=i+1;right=T+Q;边界条件right(1)=0;right(end)=0;T=A\\right;Tseriers=[Tseriers,T];leg_info{i}=['t=',num2str(t)];endplot(x,Tseriers)legend(leg_info)plot(x,T,x,2*exp(-pi*tend/4)*sin(2*pi*x),'r*')legend({['T=',num2str(tend)],'精确解'})求解二维抛物线型偏微分方程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)如何用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解偏微分方程组该方程组由两个抛物型偏微分方程组成 这个没有自带的函数,需要把插分格式写出来以后自己编程。帮帮忙求解这个抛物方程的matlab程序 什么意思?急求matlab二阶抛物型方程的有限元程序 1,不一定百有效果,因为pdetool具体编程是不知道的,如果解决小问题两者的结果一样说明不了度什麽问题,尤其对于偏微分方程。2有限元的边知界必须固定,从数理方程上讲静态有道限元问题就是边值问题,回如果边界变化的话,初始一下别的专业有限元软件,答比如anasys,adima等。跪求MATLAB解抛物型偏微分方程的程序 1,不一定有效果,因为pdetool具体编程是不知道的,如果解决小问题两者的结果一样说明不了什麽问题,尤其对于偏微分方程。2有限元的边界必须固定,从数理方程上讲静态有限元问题就是边值问题,如果边界变化的话,初始一下别的专业有限元软件,比如anasys,adima等。

#matlab函数#matlab#有限元#偏微分方程

随机阅读

qrcode
访问手机版