| 查看: 903 | 回复: 6 | |||
| 当前主题已经存档。 | |||
[交流]
【求助】 求解偏微分方程组的程序 最好时用matlab编写的
|
|||
|
求助 如题 谢谢 是非线性偏微分方程组 [ Last edited by 121049115 on 2009-4-10 at 14:00 ] |
» 猜你喜欢
评审有感
已经有6人回复
面上本子正文33页,违规吗?会被低分嘛?
已经有8人回复
今年审到国自然15份,谈谈感受
已经有16人回复
国自然上会要求
已经有5人回复
上海大学实验技术岗位非升即走
已经有8人回复
考博自荐
已经有6人回复
青C资助名额大幅增加!
已经有16人回复
重磅!青年科学基金项目(C类)资助增幅预计超过50%
已经有10人回复
我在等一个没有答案的答案
已经有3人回复
半夜喝咖啡
已经有5人回复
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
hitzhang
木虫 (正式写手)
- 仿真EPI: 1
- 应助: 0 (幼儿园)
- 贵宾: 2.15
- 金币: 1376.7
- 散金: 969
- 红花: 8
- 帖子: 862
- 在线: 226.4小时
- 虫号: 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
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
![]() ![]() |
6楼2009-03-24 04:24:14
7楼2009-12-10 13:47:26












回复此楼

20