| 查看: 3728 | 回复: 5 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
goleo木虫 (正式写手)
|
[求助]
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 ] |
回复此楼» 猜你喜欢
如何从铁电相到顺电相。
已经有1人回复
-
各位大佬，求一份最新Wien2k版本作为学习用途
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有67人回复
-
南科大活性流体和软物质课题组诚招2027级博士、硕士研究生和博后
已经有42人回复
-
瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向
已经有24人回复
-
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
PVA溶解
已经有2人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 变系数偏微分方程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人回复
goleo
木虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 3708.7
- 散金: 63
- 红花: 1
- 帖子: 401
- 在线: 129.8小时
- 虫号: 741341
- 注册: 2009-04-06
- 性别: GG
- 专业: 分离过程
6楼2014-02-24 17:36:16
arising2010
铜虫 (正式写手)
- 应助: 32 (小学生)
- 金币: 1433.4
- 红花: 11
- 帖子: 996
- 在线: 412.6小时
- 虫号: 2464689
- 注册: 2013-05-15
- 专业: 可再生与替代能源利用中的
4楼2014-02-23 01:06:07
goleo
木虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 3708.7
- 散金: 63
- 红花: 1
- 帖子: 401
- 在线: 129.8小时
- 虫号: 741341
- 注册: 2009-04-06
- 性别: GG
- 专业: 分离过程
|
能否帮我看一下如下解偏微分方程的代码，这是采用的什么数值解法，有限差分？但看着公式不像 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; |
5楼2014-02-24 17:33:56
100