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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 求助：偏微分方程求解
北京石油化工学院2026年研究生招生接收调剂公告
192/2返回列表上一页12
查看: 1503  |  回复: 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
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 327求调剂 +3 Xxjc1107. 2026-04-06 3/150 2026-04-06 20:15 by lijunpoly
[考研] 一志愿南科大生物学297分，求调剂推荐 +5 Y-yyusx 2026-04-06 6/300 2026-04-06 19:53 by lijunpoly
[考研] 调剂 +9 不逢春 2026-04-05 10/500 2026-04-06 11:58 by 考研学校招点人
[考研] 材料专硕283求调剂 +17 试试看呗 2026-04-04 18/900 2026-04-06 09:24 by 286640313
[考研] 材料专硕322分 +10 哈哈哈吼吼吼哈 2026-04-04 10/500 2026-04-05 21:22 by 学员8dgXkO
[考研] 一志愿苏州大学材料工程（085601）专硕有科研经历三项国奖两个实用型专利一项省级立项 +10 大火山小火山 2026-04-05 10/500 2026-04-05 19:58 by 小蜗牛小蜗牛
[考研] 301求调剂 +3 XYPLR 2026-04-05 4/200 2026-04-05 19:07 by XYPLR
[考研] 275求调剂 +16 waltzh 2026-04-01 16/800 2026-04-05 17:14 by Hdyxbekcb
[考研] 求调剂 +3 小沢 2026-04-03 3/150 2026-04-05 09:10 by sihailian3
[考研] 306求调剂 +3 hyb上名工 2026-04-02 3/150 2026-04-04 18:12 by 热情沙漠
[考研] 考研求调剂 +3 木心想继续深造 2026-04-03 3/150 2026-04-03 21:56 by 啵啵啵0119
[考研] 283求调剂 +3 jiouuu 2026-04-03 4/200 2026-04-03 13:28 by jiouuu
[考研] 081200-11408-276学硕求调剂 +6 崔wj 2026-04-02 6/300 2026-04-03 10:19 by 蓝云思雨
[考研] 362求调剂 +14 西南交材料专硕3 2026-03-31 14/700 2026-04-02 17:50 by yunlongyang
[考研] 一志愿北京科技大学085601材料工程英一数二初试总分335求调剂 +8 双马尾痞老板2 2026-04-02 9/450 2026-04-02 14:45 by 5896
[考研] 350求调剂 +7 阿佳～ 2026-03-31 7/350 2026-04-01 16:12 by yanflower7133
[考研] 085600，321分求调剂 +13 大馋小子 2026-03-31 13/650 2026-04-01 12:35 by chemdavid
[考研] 一志愿北交材料工程总分358 +5 cs0106 2026-04-01 7/350 2026-04-01 11:45 by wangjy2002
[考研] 考研材料工程351分调剂 +5 整个好的 2026-03-31 5/250 2026-04-01 09:36 by topgun2009
[考研] 一志愿大连理工大学，机械工程学硕，341 +3 西瓜田的守望者 2026-03-30 3/150 2026-03-31 11:08 by asdfzly
信息提示
关闭
请填处理意见
关闭