| 查看: 1091 | 回复: 11 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
哪位大仙帮我看下程序，我用的matlab7.1，运行提示du未定义
|
||
|
function [c,f,s]=pdefun(x,t,u,du) c=[1;1]; f=[0.024*du(1);0.17*du(2)]; temp=u(1)-u(2); s=[-1;1].*(exp(5.73*temp)-exp(-11.46*temp)); function [pa,qa,pb,qb]=pdebc(xa,ua,xb,ub,t) pa=[0;ua(2)]; qa=[1;0]; pb=[ub(1)-1;0]; qb=[0;1]; function u0=pdeic(x) u0=[1;0]; clc x=0:0.05:1; t=0:0.05:2; m=0; sol=pdepe(m,@pdefun,@pdeic,@pdebc,x,t); figure('numbertitle','off','name','PDE Demo——by Matlabsky') subplot(211) surf(x,t,sol(:,:,1)) title('The Solution of u_1') xlabel('X') ylabel('T') zlabel('U') subplot(212) surf(x,t,sol(:,:,2)) title('The Solution of u_2') xlabel('X') ylabel('T') zlabel('U') |
回复此楼» 猜你喜欢
寻求一种能扛住强氧化性腐蚀性的容器密封件
已经有7人回复
-
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有8人回复
-
申请2026年博士
已经有6人回复
-
请问哪里可以有青B申请的本子可以借鉴一下。
已经有5人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有5人回复
-
2025冷门绝学什么时候出结果
已经有7人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有6人回复
-
Bioresource Technology期刊，第一次返修的时候被退回好几次了
已经有7人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
Kuiyanlong
铁虫 (初入文坛)
- 应助: 4 (幼儿园)
- 金币: 1135.5
- 帖子: 25
- 在线: 27.1小时
- 虫号: 2166508
- 注册: 2012-12-04
- 性别: GG
- 专业: 机械动力学
3楼2012-12-05 00:52:47
Kuiyanlong
铁虫 (初入文坛)
- 应助: 4 (幼儿园)
- 金币: 1135.5
- 帖子: 25
- 在线: 27.1小时
- 虫号: 2166508
- 注册: 2012-12-04
- 性别: GG
- 专业: 机械动力学
【答案】应助回帖
★ ★
感谢参与，应助指数 +1
xiegangmai: 金币+2, 谢谢参与 2012-12-05 19:55:09
感谢参与，应助指数 +1
xiegangmai: 金币+2, 谢谢参与 2012-12-05 19:55:09
|
换一下表达就可以了，将 sol函数提前 function ABC x=0:0.05:1; t=0:0.05:2; m=0; sol=pdepe(m,@pdefun,@pdeic,@pdebc,x,t); figure('numbertitle','off','name','PDE Demo——by Matlabsky') subplot(211) surf(x,t,sol(:,:,1)) title('The Solution of u_1') xlabel('X') ylabel('T') zlabel('U') subplot(212) surf(x,t,sol(:,:,2)) title('The Solution of u_2') xlabel('X') ylabel('T') zlabel('U') function [c,f,s]=pdefun(x,t,u,du) c=[1;1]; f=[0.024*du(1);0.17*du(2)]; temp=u(1)-u(2); s=[-1;1].*(exp(5.73*temp)-exp(-11.46*temp)); function [pa,qa,pb,qb]=pdebc(xa,ua,xb,ub,t) pa=[0;ua(2)]; qa=[1;0]; pb=[ub(1)-1;0]; qb=[0;1]; function u0=pdeic(x) u0=[1;0]; clc |
2楼2012-12-05 00:47:40
4楼2012-12-05 16:47:55
5楼2012-12-05 16:58:32