如何用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等。
随机阅读
- 水稻土游离氧化铁 水稻土的划分
- 合肥租房网在哪?
- 陕西省兴平市赵村镇邮编 陕西省咸阳市兴平市赵村镇邮编是什么?
- 腾创控科技有限公司 广州市腾创节能科技有限公司怎么样?
- 2016牙克石谁是公安局长 免渡河镇镇长
- 柴厂屯是属于永乐店镇的吗 永乐店镇南堤寺西村
- 贵州凉米线图片 陈村过桥米线
- 满天星草莓花束韩式 满天星加百合的花语·就是百合满天星的花束·是什么意思
- 武陵源到张家界大峡谷怎么走 从武陵源有车到凤凰吗 想明天去完大峡谷直接去凤凰古城 武林源到张家界多久
- 二氧化碳的性质和用途 自然氧化塘工艺
- 地方政府与自己控股的国有企业合作 政府和社会资本合作模式项目需要具备什么样 的条件? 企业社会资本与企业绩效
- 从杭州如家精选河坊街店到武林码头 和家人要去杭州旅游 除了河坊街我想知道杭州还有哪里晚上比较热闹的 不要商场还有西湖
- 上海食品工业年鉴 中国现在有多少人吃素? 世界上大概有多少人吃素? 吃素会成为一种潮流吗?
- 节妇吟寄东平李司空师道古诗 表达爱意的诗词
- 南如皋港红星花苑 深圳市区的电话有2开头和8开头的,代表什么意思的啊?
- 立案后多久通知被告 民和县松树乡派出所长电话
- 请问尿素的主要成分是什么?其主要作用又是什么? 尿素主要成分是什么
- 渡江战役主要战场 下列地区既是淮海战役战场又是渡江战役战场的是(
- 反渗透水处理电气设备控制系统说明书 二级反渗透装置使用说明书
- 试验筛r40 3是什么意思 HIV抗体初筛试验是什么意思啊