| 查看: 666 | 回复: 1 | |||
novembering捐助贵宾 (小有名气)
|
[求助]
求助混凝土湿度扩散方程的matlab数值求解程序及方法
|
| 本人模拟混凝土湿度扩散问题，需要求解湿度扩散的偏微分方程，打算用matlab求解数值解，跪求偏微分方程的求解程序，主要是边界条件的给定问题，本人把自己编写的传上来，但是与实际不相符，不知道哪里有问题，希望大神们不吝赐教，金币全给！！！ |
回复此楼» 猜你喜欢
让我中一个面上吧！
已经有13人回复
-
看来今天不会放榜了？
已经有6人回复
-
人气不行了
已经有7人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有18人回复
-
感觉是下周放榜了
已经有12人回复
-
我面上完蛋了
已经有8人回复
-
只有每年这种时候来逛逛小木虫
已经有25人回复
-
什么时候开奖？
已经有7人回复
-
时间戳又变了
已经有18人回复
-
时间戳又变了8-15
已经有27人回复
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-06 11:40:48
40