| 查看: 929 | 回复: 6 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】 求解偏微分方程组的程序 最好时用matlab编写的
|
|||
|
求助 如题 谢谢 是非线性偏微分方程组 [ Last edited by 121049115 on 2009-4-10 at 14:00 ] |
回复此楼» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有46人回复
-
各位道友，我要去昆明玩几天，回来见。
已经有4人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有16人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有11人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有4人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有8人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有37人回复
-
好奇怪的filecode
已经有7人回复
-
filecode
已经有9人回复
-
请教兼职经验
已经有3人回复
hitzhang
木虫 (正式写手)
- 仿真EPI: 1
- 应助: 0 (幼儿园)
- 贵宾: 2.15
- 金币: 1356.7
- 散金: 969
- 红花: 8
- 帖子: 863
- 在线: 226.5小时
- 虫号: 390575
- 注册: 2007-06-02
- 性别: GG
- 专业: 无机非金属类电介质与电解
★ ★ ★ ★
sunxiao(金币+4,VIP+0):谢谢参与，欢迎常来仿真版交流 3-10 23:50
sunxiao(金币+4,VIP+0):谢谢参与，欢迎常来仿真版交流 3-10 23:50
|
function pdex4 m = 0; x = [0 0.005 0.01 0.05 0.1 0.2 0.5 0.7 0.9 0.95 0.99 0.995 1]; t = [0 0.005 0.01 0.05 0.1 0.5 1 1.5 2]; sol = pdepe(m,@pdex4pde,@pdex4ic,@pdex4bc,x,t); u1 = sol(:,:,1); u2 = sol(:,:,2); figure surf(x,t,u1) title('u1(x,t)') xlabel('Distance x') ylabel('Time t') figure surf(x,t,u2) title('u2(x,t)') xlabel('Distance x') ylabel('Time t') % -------------------------------------------------------------- function [c,f,s] = pdex4pde(x,t,u,DuDx) c = [1; 1]; f = [0.024; 0.17] .* DuDx; y = u(1) - u(2); F = exp(5.73*y)-exp(-11.47*y); s = [-F; F]; % -------------------------------------------------------------- function u0 = pdex4ic(x); u0 = [1; 0]; % -------------------------------------------------------------- function [pl,ql,pr,qr] = pdex4bc(xl,ul,xr,ur,t) pl = [0; ul(2)]; ql = [1; 0]; pr = [ur(1)-1; 0]; qr = [0; 1]; [ Last edited by hitzhang on 2009-3-10 at 22:48 ] |
4楼2009-03-10 22:46:55
woshilsh
荣誉版主 (职业作家)
- 应助: 43 (小学生)
- 贵宾: 10.33
- 金币: 21115.1
- 散金: 737
- 红花: 22
- 帖子: 3777
- 在线: 639.7小时
- 虫号: 293038
- 注册: 2006-11-04
- 专业: 计算机应用技术
- 管辖: 仿真模拟
2楼2008-12-11 22:17:38
3楼2008-12-12 14:38:43
kleinkinder
铁杆木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 8530.9
- 散金: 1
- 红花: 2
- 帖子: 569
- 在线: 117.6小时
- 虫号: 280499
- 注册: 2006-09-17
- 专业: 分离过程
★
sunxiao(金币+1,VIP+0):thanks a lot 3-24 09:07
sunxiao(金币+1,VIP+0):thanks a lot 3-24 09:07
|
推荐你一本书，上面有相关问题。 http://rs495.rapidshare.com/file ... .www.netbks.com.rar |
5楼2009-03-19 13:10:47
