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

登錄注冊

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

24小時熱門版塊排行榜

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

返回列表

chugezhuang

鐵蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 46
帖子: 14
在線: 9小時
蟲號: 1407185
注冊: 2011-09-19
性別: GG
專業(yè): 巖土與基礎(chǔ)工程

[求助] 哪位大仙幫我看下程序，我用的matlab7.1，運行提示du未定義

function [c,f,s]=pdefun(x,t,u,du)
c=[1;1];
f=[0.024*du(1);0.17*du(2)];
temp=u(1)-u(2);
s=[-1;1].*(exp(5.73*temp)-exp(-11.46*temp));
function [pa,qa,pb,qb]=pdebc(xa,ua,xb,ub,t)
pa=[0;ua(2)];
qa=[1;0];
pb=[ub(1)-1;0];
qb=[0;1];
function u0=pdeic(x)
u0=[1;0];
clc
x=0:0.05:1;
t=0:0.05:2;
m=0;
sol=pdepe(m,@pdefun,@pdeic,@pdebc,x,t);
figure('numbertitle','off','name','PDE Demo——by Matlabsky')
subplot(211)
surf(x,t,sol(:,:,1))
title('The Solution of u_1')
xlabel('X')
ylabel('T')
zlabel('U')
subplot(212)
surf(x,t,sol(:,:,2))
title('The Solution of u_2')
xlabel('X')
ylabel('T')
zlabel('U')

回復(fù)此樓

» 猜你喜歡

308求調(diào)劑已經(jīng)有10人回復(fù)
285求調(diào)劑已經(jīng)有3人回復(fù)
070300一志愿211，312分求調(diào)劑院校已經(jīng)有11人回復(fù)
考研生物與醫(yī)藥調(diào)劑已經(jīng)有6人回復(fù)
環(huán)境工程 085701，267求調(diào)劑已經(jīng)有13人回復(fù)
08開頭看過來！�。� 已經(jīng)有5人回復(fù)
086000調(diào)劑已經(jīng)有5人回復(fù)
070300化學專業(yè)279調(diào)劑已經(jīng)有7人回復(fù)
本2一志愿C9-333分，材料科學與工程，求調(diào)劑已經(jīng)有8人回復(fù)
328求調(diào)劑已經(jīng)有5人回復(fù)

Life isn't about waiting for the storm to pass, it's about learning to dance in the rain.

1樓 2012-12-03 10:11:23

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

Kuiyanlong

鐵蟲 (初入文壇)

應(yīng)助: 4 (幼兒園)
金幣: 1135.5
帖子: 25
在線: 27.1小時
蟲號: 2166508
注冊: 2012-12-04
性別: GG
專業(yè): 機械動力學

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

★ ★
感謝參與，應(yīng)助指數(shù) +1
xiegangmai: 金幣+2, 謝謝參與 2012-12-05 19:55:09

換一下表達就可以了，將 sol函數(shù)提前
function ABC
x=0:0.05:1;
t=0:0.05:2;
m=0;
sol=pdepe(m,@pdefun,@pdeic,@pdebc,x,t);
figure('numbertitle','off','name','PDE Demo——by Matlabsky')
subplot(211)
surf(x,t,sol(:,:,1))
title('The Solution of u_1')
xlabel('X')
ylabel('T')
zlabel('U')
subplot(212)
surf(x,t,sol(:,:,2))
title('The Solution of u_2')
xlabel('X')
ylabel('T')
zlabel('U')

function [c,f,s]=pdefun(x,t,u,du)
c=[1;1];
f=[0.024*du(1);0.17*du(2)];
temp=u(1)-u(2);
s=[-1;1].*(exp(5.73*temp)-exp(-11.46*temp));
function [pa,qa,pb,qb]=pdebc(xa,ua,xb,ub,t)
pa=[0;ua(2)];
qa=[1;0];
pb=[ub(1)-1;0];
qb=[0;1];
function u0=pdeic(x)
u0=[1;0];
clc

贊一下

回復(fù)此樓

Bebornthisway

2樓2012-12-05 00:47:40

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

Kuiyanlong

鐵蟲 (初入文壇)

應(yīng)助: 4 (幼兒園)
金幣: 1135.5
帖子: 25
在線: 27.1小時
蟲號: 2166508
注冊: 2012-12-04
性別: GG
專業(yè): 機械動力學

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

太可惜了運行出來的圖片放不進去

贊一下

回復(fù)此樓

Bebornthisway

3樓2012-12-05 00:52:47

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

chugezhuang

鐵蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 46
帖子: 14
在線: 9小時
蟲號: 1407185
注冊: 2011-09-19
性別: GG
專業(yè): 巖土與基礎(chǔ)工程

引用回帖:

3樓: Originally posted by Kuiyanlong at 2012-12-05 00:52:47
太可惜了運行出來的圖片放不進去

已經(jīng)解決了，不過還是非常感謝

贊一下

回復(fù)此樓

Life isn't about waiting for the storm to pass, it's about learning to dance in the rain.

4樓2012-12-05 16:47:55

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

chugezhuang

鐵蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 46
帖子: 14
在線: 9小時
蟲號: 1407185
注冊: 2011-09-19
性別: GG
專業(yè): 巖土與基礎(chǔ)工程

引用回帖:

2樓: Originally posted by Kuiyanlong at 2012-12-05 00:47:40
換一下表達就可以了，將 sol函數(shù)提前
function ABC
x=0:0.05:1;
t=0:0.05:2;
m=0;
sol=pdepe(m,@pdefun,@pdeic,@pdebc,x,t);
figure('numbertitle','off','name','PDE Demo——by Matlabsky')
subplot ...

你好，再麻煩你下，能不能幫我把這個解的程序編出來，我想要解的圖，可我不太會弄。L,Tv,p0都是常數(shù)

TC[OK9ZXP3BTFH{S825$[49.jpg

贊一下

回復(fù)此樓

Life isn't about waiting for the storm to pass, it's about learning to dance in the rain.

5樓2012-12-05 16:58:32

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

Kuiyanlong

鐵蟲 (初入文壇)

應(yīng)助: 4 (幼兒園)
金幣: 1135.5
帖子: 25
在線: 27.1小時
蟲號: 2166508
注冊: 2012-12-04
性別: GG
專業(yè): 機械動力學

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

t在哪里？？？？

贊一下

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

chugezhuang

Bebornthisway

6樓2012-12-05 21:06:53

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

chugezhuang

鐵蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 46
帖子: 14
在線: 9小時
蟲號: 1407185
注冊: 2011-09-19
性別: GG
專業(yè): 巖土與基礎(chǔ)工程

送鮮花一朵

引用回帖:

6樓: Originally posted by Kuiyanlong at 2012-12-05 21:06:53
t在哪里？？？？

Tv=Cv*t，Cv是常數(shù)。非常感謝。我看的頭大了這幾天，呵呵

贊一下

回復(fù)此樓

Life isn't about waiting for the storm to pass, it's about learning to dance in the rain.

7樓2012-12-06 10:38:46

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

Kuiyanlong

鐵蟲 (初入文壇)

應(yīng)助: 4 (幼兒園)
金幣: 1135.5
帖子: 25
在線: 27.1小時
蟲號: 2166508
注冊: 2012-12-04
性別: GG
專業(yè): 機械動力學

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

★
xiegangmai: 金幣+1, 謝謝參與 2012-12-06 19:45:17

我試了一下，while循環(huán)結(jié)果不是很理想，用for循環(huán)驗算，發(fā)現(xiàn)u的結(jié)果太小
ep=1e-6;
p0=0.25;
t=0.02;
Cv=0.04;
Tv=Cv*t;
L=0.5;
z=3;
n=0;
c0=2.*L.*p0./pi;
u=0;
for n=0:100
m=2.*n+1
u=u+1/m*sin(m/(2*L))*pi*z*eps(-Tv*m^2*pi*pi/4)
end

贊一下

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

chugezhuang

Bebornthisway

8樓2012-12-06 18:50:32

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

Kuiyanlong

鐵蟲 (初入文壇)

應(yīng)助: 4 (幼兒園)
金幣: 1135.5
帖子: 25
在線: 27.1小時
蟲號: 2166508
注冊: 2012-12-04
性別: GG
專業(yè): 機械動力學

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

★
xiegangmai: 金幣+1, 謝謝參與 2012-12-06 19:45:07

eps=1e-6;
p0=0.25;
t=0.02;
Cv=0.04;
Tv=Cv*t;
L=0.5;
z=3;
n=1;
c0=2.*L.*p0./pi;
u=0;
u=inline('1/m*sin(m/(2*L))*pi*z*eps(-Tv*m^2*pi*pi/4)','m','L','Tv','z');
u=u(n,L,Tv,z);
eu=0;
while eu>=eps
m=2*n+1
eu=u(m,L,Tv,z)
u=u+eu;
eu=co*eu;
end
u=c0*u-p0

還是上面說的問題，那個求和的部分太小

贊一下

回復(fù)此樓

Bebornthisway

9樓2012-12-06 19:12:05

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

Kuiyanlong

鐵蟲 (初入文壇)

應(yīng)助: 4 (幼兒園)
金幣: 1135.5
帖子: 25
在線: 27.1小時
蟲號: 2166508
注冊: 2012-12-04
性別: GG
專業(yè): 機械動力學

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

★
xiegangmai: 金幣+1, 謝謝參與 2012-12-06 19:45:00

eps=1e-6;
p0=0.25;
t=0.02;
Cv=0.04;
Tv=Cv*t;
L=0.5;
z=3;
n=1;
c0=2.*L.*p0./pi;
c=inline('1/m*sin(m/(2*L))*pi*z*eps(-Tv*m^2*pi*pi/4)','m','L','Tv','z');
u=c(n,L,Tv,z);
eu=1;
while abs(eu)>=eps
m=2*n+1;
eu=c(m,L,Tv,z);
u=u+eu;
n=n+1;
eu=c0*eu;
end
u=c0*u-p0;

贊一下

回復(fù)此樓

Bebornthisway

10樓2012-12-06 19:27:34

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

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

返回列表

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

[考研] 一志愿武理材料工程302調(diào)劑環(huán)化或化工 +6	Doleres 2026-03-31	6/300	2026-03-31 17:38 by 2026材料調(diào)劑
[考研] 080500-315分復(fù)試調(diào)劑 +9	上岸3821 2026-03-31	9/450	2026-03-31 17:29 by 唐沐兒
[考研] 286求調(diào)劑 +5	丟掉懶惰 2026-03-27	8/400	2026-03-31 11:27 by Delta2012
[考研] 一志愿食品科學與工程083200求調(diào)劑 +4	XQTJZ 2026-03-30	4/200	2026-03-31 04:10 by fmesaito
[考研] 合肥區(qū)域性重點一本招收調(diào)劑 +4	6266jl 2026-03-30	4/200	2026-03-30 21:48 by zhuangyan123
[考研] 279求調(diào)劑 +12	j的立方 2026-03-29	12/600	2026-03-30 20:30 by dick_runner
[考研] 求調(diào)劑 +10	家佳佳佳佳佳 2026-03-29	10/500	2026-03-30 18:34 by 544594351
[考研] 328求調(diào)劑 +8	嗯滴的基本都 2026-03-27	8/400	2026-03-30 17:20 by Wang200018
[考研] 342求調(diào)劑 +4	加油a李zs 2026-03-26	4/200	2026-03-30 16:39 by 晶體之美
[考研] 329求調(diào)劑 +8	星野? 2026-03-26	8/400	2026-03-30 13:41 by chemdavid
[考研] 318一志愿吉林大學生物與醫(yī)藥求調(diào)劑 +5	篤行致遠. 2026-03-28	5/250	2026-03-30 06:56 by ilovexiaobin
[考研] 305求調(diào)劑 +8	RuiFairyrui 2026-03-28	8/400	2026-03-29 08:22 by fmesaito
[考研] 081200-11408-276學碩求調(diào)劑 +6	崔wj 2026-03-26	6/300	2026-03-29 01:11 by hanserlol
[考研] 0856，材料與化工321分求調(diào)劑 +12	大饞小子 2026-03-27	13/650	2026-03-28 10:56 by self2008
[考研] 275求調(diào)劑 +10	jjjjjjjjjjl 2026-03-27	10/500	2026-03-27 23:47 by barnett0632
[考研] 266分求材料化工冶金礦業(yè)等專業(yè)的調(diào)劑 +4	哇呼哼呼哼 2026-03-26	4/200	2026-03-27 17:02 by zhyzzh
[碩博家園] 北京林業(yè)大學碩導(dǎo)招生廣告 +6	kongweilin 2026-03-26	8/400	2026-03-27 10:18 by FF_16
[考研] 機械學碩310分，數(shù)一英一，一志愿211本科雙非找調(diào)劑信息 +3	@357 2026-03-25	3/150	2026-03-26 16:34 by by.MENG
[考研] 打過很多競賽，085406控制工程300分，求調(diào)劑 +3	askeladz 2026-03-26	3/150	2026-03-26 09:08 by 給你你注意休息
[考研] 【2026考研調(diào)劑】制藥工程 284分求相關(guān)專業(yè)調(diào)劑名額 +4	袁奐奐 2026-03-25	8/400	2026-03-25 14:32 by lbsjt

24小時熱門版塊排行榜

[求助] 哪位大仙幫我看下程序，我用的matlab7.1，運行提示du未定義

» 猜你喜歡

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

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

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

» 本帖已獲得的紅花（最新10朵）

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

» 本帖已獲得的紅花（最新10朵）

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

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

[求助] 哪位大仙幫我看下程序，我用的matlab7.1，運行提示du未定義