版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

zhangzhiguang

至尊木蟲 (文壇精英)

應(yīng)助: 121 (高中生)
金幣: 53088.8
紅花: 24
沙發(fā): 222
帖子: 26519
在線: 881.4小時
蟲號: 492270
注冊: 2008-01-08
性別: GG
專業(yè): 基礎(chǔ)物理學(xué)

別太著急,就這幾天,突然有點事兒要辦

贊一下

回復(fù)此樓

好好學(xué)習(xí),天天向上!

11樓2013-01-26 11:46:51

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zhangzhiguang

至尊木蟲 (文壇精英)

應(yīng)助: 121 (高中生)
金幣: 53088.8
紅花: 24
沙發(fā): 222
帖子: 26519
在線: 881.4小時
蟲號: 492270
注冊: 2008-01-08
性別: GG
專業(yè): 基礎(chǔ)物理學(xué)

【答案】應(yīng)助回帖

★
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)

贊一下

回復(fù)此樓

好好學(xué)習(xí),天天向上!

12樓2013-01-27 13:22:00

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zhangzhiguang

至尊木蟲 (文壇精英)

應(yīng)助: 121 (高中生)
金幣: 53088.8
紅花: 24
沙發(fā): 222
帖子: 26519
在線: 881.4小時
蟲號: 492270
注冊: 2008-01-08
性別: GG
專業(yè): 基礎(chǔ)物理學(xué)

【答案】應(yīng)助回帖

★
csgt0: 金幣+1, 謝謝 2013-01-28 09:34:43

那個笑臉是:加),說明,空間步長是0.01H,時間是0.01T,取得時間長度為100個步長,這里你回去自己修改吧,如果取得較大把mesh(h)去掉，不知你計算的是什么，感覺算出來不太對勁，你那公式里h的5次方h是系數(shù)嗎？我是按函數(shù)算得

贊一下(1人)

回復(fù)此樓

好好學(xué)習(xí),天天向上!

13樓2013-01-27 13:26:08

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

southpark007

禁蟲 (著名寫手)

★
csgt0: 金幣+1, 應(yīng)助指數(shù)+1, 謝謝 2013-01-28 09:35:13

本帖內(nèi)容被屏蔽

14樓2013-01-27 21:32:44

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

土老九

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1357.3
散金: 236
紅花: 1
帖子: 126
在線: 39小時
蟲號: 841055
注冊: 2009-09-05
專業(yè): 作物栽培與耕作學(xué)

★
fegg7502: 金幣+1, 鼓勵交流 2013-01-29 09:03:57

引用回帖:

14樓: Originally posted by southpark007 at 2013-01-27 21:32:44
用PDEPE來做就可以了，把t設(shè)定成無窮大，可以繞過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設(shè)置成反正切函數(shù)（三角形）增加-------------------------
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;

贊一下

回復(fù)此樓

耕者有其田

15樓2013-01-28 08:42:44

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

土老九

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1357.3
散金: 236
紅花: 1
帖子: 126
在線: 39小時
蟲號: 841055
注冊: 2009-09-05
專業(yè): 作物栽培與耕作學(xué)

引用回帖:

14樓: Originally posted by southpark007 at 2013-01-27 21:32:44
用PDEPE來做就可以了，把t設(shè)定成無窮大，可以繞過t無窮大時候的h情況。

結(jié)果不對，請幫助看看問題出在哪里？

贊一下

回復(fù)此樓

耕者有其田

16樓2013-01-28 08:43:49

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

southpark007

禁蟲 (著名寫手)

★
fegg7502: 金幣+1, 鼓勵交流 2013-01-29 09:04:09

本帖內(nèi)容被屏蔽

17樓2013-01-28 09:47:58

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

southpark007

禁蟲 (著名寫手)

fegg7502: 應(yīng)助指數(shù)+1, 3ks 2013-01-29 09:04:22

本帖內(nèi)容被屏蔽

18樓2013-01-28 12:07:55

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

土老九

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1357.3
散金: 236
紅花: 1
帖子: 126
在線: 39小時
蟲號: 841055
注冊: 2009-09-05
專業(yè): 作物栽培與耕作學(xué)

送鮮花一朵

thanks a lot

回復(fù)此樓

耕者有其田

19樓2013-01-28 13:50:03

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主土老九的主題更新

返回列表

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 275求調(diào)劑 +5	jjjjjjjjjjl 2026-03-27	5/250	2026-03-27 12:52 by lature00
[材料工程] 一志愿C9材料與化工專業(yè)總分300求調(diào)劑 +6	曼111 2026-03-24	7/350	2026-03-27 11:58 by 不吃魚的貓
[考研] 材料292調(diào)劑 +11	橘頌思美人 2026-03-23	11/550	2026-03-27 11:25 by wangjy2002
[考研] 調(diào)劑推薦 +5	清酒714 2026-03-26	6/300	2026-03-27 11:12 by 不吃魚的貓
[考研] 324求調(diào)劑 +5	hanamiko 2026-03-26	5/250	2026-03-27 10:33 by wangjy2002
[考研] 一志愿哈爾濱工業(yè)大學(xué)材料與化工方向336分 +3	辰沐5211314 2026-03-26	3/150	2026-03-27 10:24 by 盡舜堯1
[考研] 求調(diào)劑323材料與化工 +7	1124361 2026-03-24	7/350	2026-03-27 10:22 by wangjy2002
[考研] 安徽大學(xué)專碩生物與醫(yī)藥專業(yè)(086000)324分，英語已過四六級，六級521，求調(diào)劑 +3	美味可樂雞翅 2026-03-26	3/150	2026-03-27 08:46 by hxsm
[考研] 316求調(diào)劑 +5	江辭666 2026-03-26	5/250	2026-03-27 08:08 by hypershenger
[考研] 304材料求調(diào)劑 +4	鐘llll 2026-03-26	4/200	2026-03-27 03:42 by wxiongid
[考研] 303求調(diào)劑 +7	元夕元 2026-03-20	8/400	2026-03-26 20:38 by 不吃魚的貓
[考研] 08開頭275求調(diào)劑 +3	拉誰不重要 2026-03-26	3/150	2026-03-26 20:22 by barlinike
[考研] 271求調(diào)劑 +6	生如夏花… 2026-03-22	6/300	2026-03-26 16:48 by 張凱十八號
[考研] 尋找調(diào)劑 +5	倔強芒? 2026-03-21	8/400	2026-03-26 13:25 by 0906ljy
[考研] 考研調(diào)劑 +6	來好運來來來 2026-03-21	7/350	2026-03-25 22:43 by 418490947
[考研] 290分調(diào)劑求助 +3	吉祥止止陳 2026-03-25	3/150	2026-03-25 19:58 by barlinike
[考研] 材料專碩331求調(diào)劑 +4	鮮當(dāng)牛 2026-03-24	4/200	2026-03-24 15:58 by JourneyLucky
[考博] 26申博自薦 +3	whh869393 2026-03-24	3/150	2026-03-24 09:55 by 21018060
[考研] 求調(diào)劑 +7	十三加油 2026-03-21	7/350	2026-03-23 23:48 by 熱情沙漠
[論文投稿] 急發(fā)核心期刊論文 +3	賢達問津 2026-03-23	5/250	2026-03-23 17:13 by 妹子不好惹