24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

goleo

木虫 (正式写手)

应助: 1 (幼儿园)
金币: 3372.7
散金: 63
红花: 1
帖子: 401
在线: 129.6小时
虫号: 741341
注册: 2009-04-06
性别: GG
专业: 分离过程

[求助] MATLAB求解偏微分方程组

本人菜鸟，有两个问题想请教
1.用pde求解偏微分方程组，所有的pdefun都可以化为标准格式吗？
2.有如下方程组（附件），请高人帮忙，可以用pdepe命令求数值解吗？想模拟出一条c关于t的曲线，与实验测得的c-t曲线对比
其中，q,c,cp,r和t为变量

另外，mathtype编辑的公式，感觉第二行以下有点像悬挂缩进，怎么解决

孔扩散.gif

[ Last edited by goleo on 2014-2-19 at 11:02 ]

回复此楼

» 猜你喜欢

反铁磁体中的磁性切换：两种不同的机制已成功可视化已经有0人回复
求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有265人回复
新西兰Robinson研究所招收全奖PhD 已经有0人回复
石墨烯转移--二氧化硅衬底石墨烯已经有0人回复
笼目材料中量子自旋液体基态的证据已经有0人回复
数学教学论硕士可以读数学物理博士吗？已经有0人回复
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生已经有4人回复
澳门大学应用物理及材料工程研究院潘晖教授课题组诚招博士后已经有11人回复

» 本主题相关价值贴推荐，对您同样有帮助:

变系数偏微分方程matlab求解已经有6人回复
一个偏微分方程组的求解问题已经有5人回复
matlab 非线性微分方程求解已经有3人回复
用matlab求解方程已经有6人回复
偏微分方程求解！已经有4人回复
matlab如何求解一个非线性微分方程组已经有8人回复
如何用Matlab求时变微分方程的解？已经有5人回复
matlab解偏微分方程求助已经有12人回复
matlab求解高阶方程组已经有3人回复
求解一个二阶偏微分方程~~有结果，但是不知道中间过程（完美解决者答谢100￥）~ 已经有4人回复
有偿求助：matlab怎么解联立方程组？已经有11人回复
求高手解答Matlab求解微分方程的初始值和边界值问题已经有5人回复
Matlab求解二阶偏微分方程组，希望能给出相应的求解代码已经有10人回复
有限差分离散二阶偏微分方程后如何编译matlab程序解方程组？已经有3人回复
三元偏微分方程求解已经有4人回复
matlab有限差分法求解复杂的偏微分方程组已经有19人回复
求助MATLAB解方程组已经有6人回复
Matlab求解偏微分方程组已经有22人回复
matlab数值求解边界条件微分方程组已经有7人回复
高金求助matlab解微分方程组已经有12人回复
【求助】matlab 求解微分方程中的未知参数已经有20人回复
【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题已经有21人回复
【求助】MATLAB 有限差分法（FDM）求解偏微分方程已经有22人回复
【求助】积分微分方程matlab求解已经有6人回复
【求助】matlab怎么求解偏微分方程组啊，先谢谢了已经有13人回复
【求助】用MATLAB求解偏微分方程的数值解悬赏20 已经有8人回复

路漫漫其修远兮，吾将上下而求索

1楼 2014-02-19 11:01:20

已阅关注TA 给TA发消息送TA红花 TA的回帖

goleo

木虫 (正式写手)

应助: 1 (幼儿园)
金币: 3372.7
散金: 63
红花: 1
帖子: 401
在线: 129.6小时
虫号: 741341
注册: 2009-04-06
性别: GG
专业: 分离过程

引用回帖:

4楼: Originally posted by arising2010 at 2014-02-23 01:06:07
你把q带入第一个偏微分方程，整理成完全关于cm的方程
matlab里好像有直接求解这类方程的命令。
这个解出来了，后面两个就好求了

能否帮我看一下如下解偏微分方程的代码，这是采用的什么数值解法，有限差分？但看着公式不像
NN=20;
hk=1/NN;
RT=[(1-1/sqrt(3))/2 (1+1/sqrt(3))/2]';
H=[(1-RT).^2.*(1+2*RT) RT.*(1-RT).^2 RT.^2.*(3-2*RT) RT.^2.*(RT-1)];
A=[6*RT.*(RT-1) (1-RT).*(1-3*RT) -6*RT.*(RT-1) RT.*(3*RT-2)];
B=[12*RT-6 6*RT-4 6-12*RT 6*RT-2];
E=diag([1,hk,1,hk]);
H=H*E; A=A*E; B=B*E;
MatrixM=zeros(2*NN,2*NN);
MatrixH=zeros(2*NN+2,2*NN+2);
MatrixA=zeros(2*NN+2,2*NN+2);
MatrixB=zeros(2*NN+2,2*NN+2);
for i=1:NN
MatrixM(2*i-1:2*i,2*i-1:2*i+2)=H(:,

;
MatrixH(2*i:2*i+1,2*i-1:2*i+2)=H(:,

;
MatrixA(2*i:2*i+1,2*i-1:2*i+2)=A(:,

;
MatrixB(2*i:2*i+1,2*i-1:2*i+2)=B(:,

;
end
%MatrixH(2*NN+1:2*NN+2,2*NN:2*NN+2)=H(:,[1 2 3]);
%MatrixA(2*NN+1:2*NN+2,2*NN:2*NN+2)=A(:,[1 2 3]);
%MatrixB(2*NN+1:2*NN+2,2*NN:2*NN+2)=B(:,[1 2 3]);
MatrixH(1,1)=1;
MatrixH(1,2)=-1/Peclet; %Peclet就是一个具体数
MatrixH(2*NN+2,2*NN+2)=1;
MatrixH;
MatrixA;
MatrixB;
%--------------------------------------------------------------------------
A=MatrixA/hk;