| 查看: 665 | 回复: 1 | |||
novembering捐助贵宾 (小有名气)
|
[求助]
求助混凝土湿度扩散方程的matlab数值求解程序及方法
|
| 本人模拟混凝土湿度扩散问题，需要求解湿度扩散的偏微分方程，打算用matlab求解数值解，跪求偏微分方程的求解程序，主要是边界条件的给定问题，本人把自己编写的传上来，但是与实际不相符，不知道哪里有问题，希望大神们不吝赐教，金币全给！！！ |
回复此楼» 猜你喜欢
科研孤儿太难了
已经有9人回复
-
时间戳又变了
已经有15人回复
-
时间戳变了，能看出什么问题？
已经有23人回复
-
只有每年这种时候来逛逛小木虫
已经有22人回复
-
今天系统多次维护，明天很可能放榜！
已经有10人回复
-
放榜前的不淡定
已经有12人回复
-
言之凿凿，明天官网会通知放榜？
已经有16人回复
-
重要消息，中午系统在维护
已经有12人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有20人回复
-
今天基金会出结果吗？20260819
已经有15人回复
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