应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 求助：偏微分方程求解
192/2返回列表上一页12
查看: 1704  |  回复: 18
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

zhangzhiguang

至尊木虫 (文坛精英)
别太着急,就这几天,突然有点事儿要办
赞一下 回复此楼
好好学习,天天向上!
11楼2013-01-26 11:46:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangzhiguang

至尊木虫 (文坛精英)

【答案】应助回帖

★
csgt0: 金币+1, 谢谢 2013-01-28 09:35:01
clear;clc;
h=zeros(100,100);
h(:,1)=10;
h(1,=0.1;
h(:,100)=0.1;
omega=1.1;

for d=1:100
    for n=1:98
        for i=2:99
           
            h(100,n+1)=1e-5*(h(100,n))^-5*(h(99,n)-2*h(100,n)+h(99,n))+h(100,n);
            h(i,n+1)=(1-omega)* h(i,n+1)+omega*(1e-5*(h(i,n))^-5*(h(i+1,n)-2*h(i,n)+h(i-1,n))+h(i,n));
        end
    end
end


mesh(h)
赞一下 回复此楼
好好学习,天天向上!
12楼2013-01-27 13:22:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangzhiguang

至尊木虫 (文坛精英)

【答案】应助回帖

★
csgt0: 金币+1, 谢谢 2013-01-28 09:34:43
那个笑脸是:加),说明,空间步长是0.01H,时间是0.01T,取得时间长度为100个步长,这里你回去自己修改吧,如果取得较大把mesh(h)去掉，不知你计算的是什么，感觉算出来不太对劲，你那公式里h的5次方h是系数吗？我是按函数算得
赞一下(1人) 回复此楼
好好学习,天天向上!
13楼2013-01-27 13:26:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

southpark007

禁虫 (著名写手)
★
csgt0: 金币+1, 应助指数+1, 谢谢 2013-01-28 09:35:13
本帖内容被屏蔽
14楼2013-01-27 21:32:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

土老九

银虫 (小有名气)
★
fegg7502: 金币+1, 鼓励交流 2013-01-29 09:03:57
引用回帖:
14楼: Originally posted by southpark007 at 2013-01-27 21:32:44
用PDEPE来做就可以了，把t设定成无穷大，可以绕过t无穷大时候的h情况。

function pdex1
m=0;
x=linspace(0,0.14,14);
t=linspace(0,100000,100);
sol=pdepe(m,@pdex1pde,@pdex1ic,@pdex1bc,x,t);
% Extract the first solution component as u.
u=sol(:,:,1);
% A surface plot is often a good way to study a solution.
surf(x,t,u)
title('Numerical solution computed with 20 mesh points.')
xlabel('Distance x')
ylabel('Time t')
% A solution profile can also be illuminating.
figure
plot(x,u(end,)
title('Solution at t=2')
xlabel('Distance x')
ylabel('u(x,2)')
% --------------------------------------------------------------
function [c,f,s]=pdex1pde(x,t,u,DuDx)
c=u^5;
f=0.0000001*DuDx;
s=0;
% -------------为保证收敛，初始u设置成反正切函数（三角形）增加-------------------------
function u0=pdex1ic(x)
u0=10*atan(100000*x)/(pi/2);
% --------------------------------------------------------------
function [pl,ql,pr,qr]=pdex1bc(xl,ul,xr,ur,t)
pl=u;
q1=0;
pr=0;
qr=1;
赞一下 回复此楼
耕者有其田
15楼2013-01-28 08:42:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

土老九

银虫 (小有名气)
引用回帖:
14楼: Originally posted by southpark007 at 2013-01-27 21:32:44
用PDEPE来做就可以了，把t设定成无穷大，可以绕过t无穷大时候的h情况。

结果不对，请帮助看看问题出在哪里？
赞一下 回复此楼
耕者有其田
16楼2013-01-28 08:43:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

southpark007

禁虫 (著名写手)
★
fegg7502: 金币+1, 鼓励交流 2013-01-29 09:04:09
本帖内容被屏蔽
17楼2013-01-28 09:47:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

southpark007

禁虫 (著名写手)
fegg7502: 应助指数+1, 3ks 2013-01-29 09:04:22
本帖内容被屏蔽
18楼2013-01-28 12:07:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

土老九

银虫 (小有名气)
送鲜花一朵
thanks a lot
回复此楼
耕者有其田
19楼2013-01-28 13:50:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 土老九 的主题更新
192/2返回列表上一页12
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 基金啊基金 +3 longfie172 2026-08-20 3/150 2026-08-21 00:04 by wuchongjun
[基金申请] 科研孤儿太难了 +14 我4大白菜 2026-08-20 15/750 2026-08-20 22:33 by 袁向阳007
[基金申请] 言之凿凿，明天官网会通知放榜？ +15 医学老男孩 2026-08-20 17/850 2026-08-20 20:35 by CaoLG
[基金申请] 时间戳又变了 +12 wuchongjun 2026-08-20 17/850 2026-08-20 20:05 by spzhao
[基金申请] 时间戳变了，能看出什么问题？ +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 今天基金会出结果吗？20260819 +14 kkkl_v 2026-08-19 15/750 2026-08-20 08:21 by sunzitan
[基金申请] filecode，4个jtjc了 +11 ziyangfang 2026-08-19 13/650 2026-08-19 21:28 by aasahr
[基金申请] 今天放榜没戏了吧 +4 yuleib84 2026-08-19 5/250 2026-08-19 19:52 by hhs666
[基金申请] 今天放榜吗？ +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] filecode=后面第一个是大写字母 +10 wangze12014 2026-08-14 12/600 2026-08-19 16:56 by 苦难博士
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 什么时候开奖？ +6 CrisMessi 2026-08-18 6/300 2026-08-19 08:06 by Equinoxhua
[基金申请] 明天放榜？ +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] 各位道友，我要去昆明玩几天，回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
信息提示
关闭
请填处理意见
关闭