| 查看: 603 | 回复: 1 | |||
novembering捐助贵宾 (小有名气)
|
[求助]
求助混凝土湿度扩散方程的matlab数值求解程序及方法
|
| 本人模拟混凝土湿度扩散问题，需要求解湿度扩散的偏微分方程，打算用matlab求解数值解，跪求偏微分方程的求解程序，主要是边界条件的给定问题，本人把自己编写的传上来，但是与实际不相符，不知道哪里有问题，希望大神们不吝赐教，金币全给！！！ |
回复此楼» 猜你喜欢
过年走亲戚时感受到了所开私家车的鄙视链
已经有10人回复
-
今年春晚有几个节目很不错，点赞！
已经有12人回复
-
情人节自我反思：在爱情中有过遗憾吗？
已经有13人回复
-
体制内长辈说体制内绝大部分一辈子在底层，如同你们一样大部分普通教师忙且收入低
已经有12人回复
-
基金正文30页指的是报告正文还是整个申请书
已经有5人回复
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