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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 求助：偏微分方程求解
北京石油化工学院2026年研究生招生接收调剂公告
192/2返回列表上一页12
查看: 1502  |  回复: 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
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂 +4 熊二想上岸 2026-04-06 4/200 2026-04-06 17:17 by dongzh2009
[考研] 277工科求调剂 +8 1915668 2026-04-04 9/450 2026-04-06 15:52 by lin-da
[考研] 一志愿郑州大学材料与化工085600，求调剂 +27 吃的不少 2026-04-02 27/1350 2026-04-06 15:19 by 495374996
[考研] 301求调剂 +13 121. 2026-04-04 13/650 2026-04-06 14:20 by lqwchd
[考研] 求调剂 一志愿西南交通大学085701环境工程 282分 +9 多多爱吃汉堡 2026-04-04 9/450 2026-04-06 14:14 by lqwchd
[考研] 0703化学调剂325分 +12 15771691647 2026-04-04 13/650 2026-04-06 12:00 by lijunpoly
[考研] 081200-11408-276学硕求调剂 +4 崔wj 2026-04-04 5/250 2026-04-05 14:06 by imissbao
[考研] 313求调剂 +3 海日海日 2026-04-04 3/150 2026-04-05 07:48 by 544594351
[考研] 296材料专硕求调剂 +21 202451007219 2026-04-02 22/1100 2026-04-04 21:48 by hemengdong
[考研] 342求调剂 +3 Liang7111 2026-04-04 5/250 2026-04-04 19:47 by dongzh2009
[考研] 085602 找调剂 +4 逆时针快乐 2026-04-02 4/200 2026-04-04 19:32 by 蓝云思雨
[考研] 求调剂 +3 usbdndj 2026-04-03 3/150 2026-04-03 14:10 by dxiaoxin
[考博] 申博求助 +3 Reee1Llll 2026-04-01 3/150 2026-04-02 22:29 by 这是一个无聊的�
[考研] 求调剂推荐 +3 南山南@ 2026-04-01 3/150 2026-04-02 12:09 by xiaoranmu
[考研] 085900土木水利336分求调剂 +4 Zhangjiangj 2026-03-31 6/300 2026-04-02 11:40 by 1753564080
[考研] 372求调剂 +3 jj涌77 2026-04-02 3/150 2026-04-02 09:57 by olim
[考研] 0710生物学，325求调剂 +3 mkkkkkl 2026-04-01 3/150 2026-04-02 09:48 by Jaylen.
[考研] 310分求调剂 +4 成功上岸wang 2026-04-01 4/200 2026-04-01 20:35 by liu823948201
[考研] 08工科275求调剂，可跨考。 +5 AaAa7420 2026-03-31 5/250 2026-04-01 15:21 by 159357hjz
[考研] 085404 22408 315分 +5 zhuangyan123 2026-03-31 6/300 2026-03-31 13:48 by limeifeng
信息提示
关闭
请填处理意见
关闭