| 查看: 641 | 回复: 2 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
novembering捐助贵宾 (小有名气)
|
[求助]
求助混凝土湿度扩散方程的matlab数值求解程序及方法
|
|
| 本人模拟混凝土湿度扩散问题，需要求解湿度扩散的偏微分方程，打算用matlab求解数值解，跪求偏微分方程的求解程序，主要是边界条件的给定问题，本人把自己编写的传上来，但是与实际不相符，不知道哪里有问题，希望大神们不吝赐教，金币全给！！！ |
回复此楼» 猜你喜欢
求调剂
已经有4人回复
-
070300化学求调剂
已经有11人回复
-
抱歉
已经有5人回复
-
环境285分，过六级，求调剂
已经有9人回复
-
求材料调剂，一志愿郑州大学289分
已经有11人回复
-
化学调剂
已经有9人回复
-
303求调剂
已经有4人回复
-
化工调剂303分，过四级
已经有25人回复
-
081200-11408-276学硕求调剂
已经有6人回复
-
298分 070300求调剂
已经有11人回复
novembering
捐助贵宾 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 117.7
- 散金: 120
- 红花: 1
- 帖子: 79
- 在线: 32.8小时
- 虫号: 2711702
- 注册: 2013-10-10
- 专业: 结构工程
3楼2016-02-07 10:46:01
novembering
捐助贵宾 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 117.7
- 散金: 120
- 红花: 1
- 帖子: 79
- 在线: 32.8小时
- 虫号: 2711702
- 注册: 2013-10-10
- 专业: 结构工程
|
function pde x=0:1:30; t=0:50:300; m=0; sol=pdepe(m,@mpde,@mpic,@mpbc,x,t); u=sol(:,:,1); surf(x,t,u) figure; plot(x,u(end,) function [c,f,s]=mpde(x,t,u,du) c=1; f=8.2e-8+(1.56e-6/(1+(5-5*u)^15)); s=0; function [pa,qa,pb,qb]=mpbc(xa,ua,xb,ub,t) pa=0; qa=7.354*(0.6-ua)/8.2e-8+(1.56e-6/(1+(5-5*u)^15)); pb=30; qb=7.354*(0.6-ub)/8.2e-8+(1.56e-6/(1+(5-5*u)^15)); function u0=mpic(x) u0=1; |
2楼2016-02-07 10:44:59
