24小时热门版块排行榜    

查看: 676  |  回复: 2
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

novembering

捐助贵宾 (小有名气)

[求助] 求助混凝土湿度扩散方程的matlab数值求解程序及方法

本人模拟混凝土湿度扩散问题,需要求解湿度扩散的偏微分方程,打算用matlab求解数值解,跪求偏微分方程的求解程序,主要是边界条件的给定问题,本人把自己编写的传上来,但是与实际不相符,不知道哪里有问题,希望大神们不吝赐教,金币全给!!!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

novembering

捐助贵宾 (小有名气)

引用回帖:
2楼: Originally posted by novembering at 2016-02-07 10:44:59
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 =mpde(x,t,u,du)
c=1;
f=8.2e-8+(1.56e-6/(1+(5-5*u)^15)) ...

这个程序是从网上找的,我感觉我的边界条件有问题,导致出的结果与实际不符
3楼2016-02-07 10:46:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 3 个回答

novembering

捐助贵宾 (小有名气)

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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +16 医学老男孩 2026-08-13 39/1950 2026-08-14 12:01 by zxy7909
[论文投稿] 职称评审,求友友推荐见刊最快的期刊 +5 工厂打螺丝 2026-08-08 5/250 2026-08-14 11:03 by 玖戈弋
[基金申请] FileCode能看出啥? +10 要乐观耀哥 2026-08-10 32/1600 2026-08-14 09:37 by 要乐观耀哥
[基金申请] 关于Filecode分析方法 +9 majunge000 2026-08-10 12/600 2026-08-13 23:20 by iwuli
[文学芳草园] 阿姨 +4 汪汪锅 2026-08-09 4/200 2026-08-13 19:43 by arzu_hma
[硕博家园] 一作与独作在应聘高校教师时区别大吗 +3 mbygzh 2026-08-08 4/200 2026-08-13 19:31 by 龙-樱
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +4 majunge000 2026-08-11 6/300 2026-08-13 16:50 by 四季常青藤
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +22 Lanmanbaby 2026-08-09 36/1800 2026-08-12 22:50 by sdfapple719
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[基金申请] 是这周出结果还是下周出结果? +3 yuleib84 2026-08-11 3/150 2026-08-11 21:53 by jnhyjjm
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 关于代码变化问题,想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
信息提示
请填处理意见