chugezhuang

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

[求助] 哪位大仙幫我看下程序，我用的matlab7.1，運(yùn)行提示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ù)此樓

» 猜你喜歡

288資源與環(huán)境專(zhuān)碩求調(diào)劑，不限專(zhuān)業(yè)，有學(xué)上就行已經(jīng)有20人回復(fù)
311（085601）求調(diào)劑已經(jīng)有11人回復(fù)
材料調(diào)劑已經(jīng)有12人回復(fù)
一志愿南京航空航天大學(xué) ，080500材料科學(xué)與工程學(xué)碩已經(jīng)有3人回復(fù)
0856材料與化工調(diào)劑，339 已經(jīng)有6人回復(fù)
08工科275分求調(diào)劑已經(jīng)有6人回復(fù)
334分一志愿武理-080500 材料求調(diào)劑已經(jīng)有10人回復(fù)
324分 085600材料與化工已經(jīng)有7人回復(fù)
各位老師好，我的一志愿為北京科技大學(xué)085601材料專(zhuān)碩已經(jīng)有11人回復(fù)
309求調(diào)劑已經(jīng)有6人回復(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

鐵蟲(chóng) (初入文壇)

應(yīng)助: 4 (幼兒園)
金幣: 1135.5
帖子: 25
在線(xiàn): 27.1小時(shí)
蟲(chóng)號(hào): 2166508
注冊(cè): 2012-12-04
性別: GG
專(zhuān)業(yè): 機(jī)械動(dòng)力學(xué)

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

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

換一下表達(dá)就可以了，將 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

鐵蟲(chóng) (初入文壇)

應(yīng)助: 4 (幼兒園)
金幣: 1135.5
帖子: 25
在線(xiàn): 27.1小時(shí)
蟲(chóng)號(hào): 2166508
注冊(cè): 2012-12-04
性別: GG
專(zhuān)業(yè): 機(jī)械動(dòng)力學(xué)

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

太可惜了運(yùn)行出來(lái)的圖片放不進(jìn)去

贊一下

回復(fù)此樓

Bebornthisway

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

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

chugezhuang

鐵蟲(chóng) (初入文壇)

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

引用回帖:

3樓: Originally posted by Kuiyanlong at 2012-12-05 00:52:47
太可惜了運(yùn)行出來(lái)的圖片放不進(jìn)去

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

贊一下

回復(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

鐵蟲(chóng) (初入文壇)

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

引用回帖:

2樓: Originally posted by Kuiyanlong at 2012-12-05 00:47:40
換一下表達(dá)就可以了，將 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 ...

你好，再麻煩你下，能不能幫我把這個(gè)解的程序編出來(lái)，我想要解的圖，可我不太會(huì)弄。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

鐵蟲(chóng) (初入文壇)

應(yīng)助: 4 (幼兒園)
金幣: 1135.5
帖子: 25
在線(xiàn): 27.1小時(shí)
蟲(chóng)號(hào): 2166508
注冊(cè): 2012-12-04
性別: GG
專(zhuān)業(yè): 機(jī)械動(dòng)力學(xué)

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

t在哪里？？？？

贊一下

回復(fù)此樓

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

chugezhuang

Bebornthisway

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

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

chugezhuang

鐵蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 46
帖子: 14
在線(xiàn): 9小時(shí)
蟲(chóng)號(hào): 1407185
注冊(cè): 2011-09-19
性別: GG
專(zhuān)業(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

鐵蟲(chóng) (初入文壇)

應(yīng)助: 4 (幼兒園)
金幣: 1135.5
帖子: 25
在線(xiàn): 27.1小時(shí)
蟲(chóng)號(hào): 2166508
注冊(cè): 2012-12-04
性別: GG
專(zhuān)業(yè): 機(jī)械動(dòng)力學(xué)

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

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

我試了一下，while循環(huán)結(jié)果不是很理想，用for循環(huán)驗(yà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

鐵蟲(chóng) (初入文壇)

應(yīng)助: 4 (幼兒園)
金幣: 1135.5
帖子: 25
在線(xiàn): 27.1小時(shí)
蟲(chóng)號(hào): 2166508
注冊(cè): 2012-12-04
性別: GG
專(zhuān)業(yè): 機(jī)械動(dòng)力學(xué)

【答案】應(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

還是上面說(shuō)的問(wèn)題，那個(gè)求和的部分太小

贊一下

回復(fù)此樓

Bebornthisway

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

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

Kuiyanlong

鐵蟲(chóng) (初入文壇)

應(yīng)助: 4 (幼兒園)
金幣: 1135.5
帖子: 25
在線(xiàn): 27.1小時(shí)
蟲(chóng)號(hào): 2166508
注冊(cè): 2012-12-04
性別: GG
專(zhuān)業(yè): 機(jī)械動(dòng)力學(xué)

【答案】應(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的回帖

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

[考研] 化學(xué)308分調(diào)劑 +11	你好明天你好 2026-03-30	12/600	2026-03-31 19:14 by 楊金金金
[考研] 375求調(diào)劑 +7	雨夏整夜 2026-03-29	7/350	2026-03-31 18:52 by xhai2011
[考研] 311求調(diào)劑一志愿合肥工業(yè)大學(xué) +11	秋二十二 2026-03-30	11/550	2026-03-31 18:09 by 253863592
[考研] 本2一志愿C9-333分，材料科學(xué)與工程，求調(diào)劑 +9	升升不降 2026-03-31	9/450	2026-03-31 18:01 by 無(wú)際的草原
[考研] 322求調(diào)劑 +3	熹僖XX 2026-03-31	3/150	2026-03-31 16:32 by 記事本2026
[考研] 理學(xué)07化學(xué) 303求調(diào)劑 +8	睿08 2026-03-27	8/400	2026-03-31 16:31 by 690616278
[考研] 調(diào)劑310 +13	溫柔的晚安 2026-03-25	14/700	2026-03-31 13:03 by 記事本2026
[考研] 調(diào)劑求院校招收 +7	鶴鯨鴿 2026-03-28	7/350	2026-03-31 11:21 by oooqiao
[考研] 求調(diào)劑 +4	研s. 2026-03-25	4/200	2026-03-31 10:05 by fmesaito
[考研] 291求調(diào)劑 +12	Y-cap 2026-03-29	15/750	2026-03-31 09:25 by Huaxue_Wang
[考研] 一志愿中海洋320化學(xué)工程與技術(shù)學(xué)碩求調(diào)劑 +8	披星河 2026-03-30	8/400	2026-03-31 08:53 by lbsjt
[考研] 085701環(huán)境工程求調(diào)劑 +11	多久上課 2026-03-27	12/600	2026-03-30 21:21 by 研究僧導(dǎo)導(dǎo)
[考研] 085602 化學(xué)工程專(zhuān)碩 340分求調(diào)劑 +4	qianbai11 2026-03-29	4/200	2026-03-30 11:34 by 唐沐兒
[考研] 311求調(diào)劑 +10	lin0039 2026-03-26	10/500	2026-03-30 10:26 by herarysara
[考研] 290求調(diào)劑 +3	dfffsar 2026-03-29	3/150	2026-03-29 22:38 by 毛毛毛阿莫2
[考研] 340求調(diào)劑 +6	Amber00 2026-03-26	6/300	2026-03-29 12:06 by 無(wú)際的草原
[考研] 356求調(diào)劑 +3	gysy?s?a 2026-03-28	3/150	2026-03-29 00:33 by 544594351
[考研] 求調(diào)劑 +7	爭(zhēng)取九點(diǎn)睡 2026-03-28	8/400	2026-03-28 21:07 by 爭(zhēng)取九點(diǎn)睡
[考研] 藥學(xué)105500求調(diào)劑 +3	Ssun。。 2026-03-28	3/150	2026-03-28 11:24 by lxf170613
[考研] 調(diào)劑 +4	柚柚yoyo 2026-03-26	4/200	2026-03-26 20:43 by fmesaito