當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

亮亮1283

新蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 378
散金: 18
帖子: 156
在線: 70小時(shí)
蟲號: 890575
注冊: 2009-11-01
專業(yè): 等離子體物理

[求助] matlab程序問題，總是出錯(cuò)，麻煩大家看看問題在那里？

function pdex1()

clc
clear
Up =400;
Un =8000;
Kb =1.38*10^(-23);
h =6.63*10^(-34);
e =1.6*10^(-19);
r =0.35;
y =0.6;
n0 =1*10^16;
p0 =1*10^16;
tn =1*10^(-9);
tp =1*10^(-9);
p1 =p0;
n1 =n0;
T =300;
a =1.65*10^4;
v =(3/532)*10^17;
X =10;
P=10;
%主函數(shù)----------------------------------------------------------------
m=0;
x=0:0.1:10;
t=0:20:1000;
sol = pdepe(m,@pdex1pde,@pdex1ic,@pdex1bc,x,t);
figure('numbertitle','off','name','PDE Demo-by matlab')
subplot(211)
surf(x,t,sol(:,:,1))
title('The Solution of p')
xlable('x')
ylable('t')
zlable('p')
subplot(212)
surf(x,t,sol(:,:,2))
title('The Solution of n')
xlable('x')
ylable('t')
zlable('n')

%目標(biāo)函數(shù)----------------------------------------------------------------
function [c,f,s] = pdex1pde(x,t,u,DuDx)
c= [1;1];
f= [(Up*Kb*T/e)*DuDx(1); (Un*Kb*T/e)*DuDx(2)];
G= a*y*(1-r)*P*exp(-a*x)/(h*v);
R1= u(1)*u(2)-p0*n0;
R2=tn*(u(1)+p1)+tp*(u(2)+n1);
R=R1/R2;
s= [1;1].*(G-R);

%初值條件函數(shù)------------------------------------------------------------
function u0 = pdex1ic(x)
u0= [p0;n0];

%邊界條件函數(shù)------------------------------------------------------------
function [pa,qa,pb,qb] = pdex1bc(xa,ua,xb,ub,t)    %a為下邊界， b為上邊界
pa =[0; 0];
qa =[1; 1];
pb =[0; 0];
qb =[1; 1];

回復(fù)此樓

» 猜你喜歡

288求調(diào)劑一志愿哈工大材料與化工已經(jīng)有9人回復(fù)
化學(xué)工程085602 305分求調(diào)劑已經(jīng)有31人回復(fù)
070300化學(xué)專業(yè)279調(diào)劑已經(jīng)有3人回復(fù)
304求調(diào)劑已經(jīng)有6人回復(fù)
085404 22408 315分已經(jīng)有5人回復(fù)
一志愿a區(qū)211，085601-307分求調(diào)劑已經(jīng)有10人回復(fù)
求調(diào)劑已經(jīng)有7人回復(fù)
290分調(diào)劑求助已經(jīng)有7人回復(fù)
一志愿華東理工大學(xué)，080500學(xué)碩，317分，求調(diào)劑已經(jīng)有3人回復(fù)
考研調(diào)劑求助已經(jīng)有4人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:

求助matlab問題已經(jīng)有5人回復(fù)
請教關(guān)于matlab矩陣轉(zhuǎn)換問題已經(jīng)有6人回復(fù)
關(guān)于Matlab作圖的問題已經(jīng)有14人回復(fù)
matlab調(diào)試問題已經(jīng)有3人回復(fù)
求matlab求解一個(gè)方程組的問題已經(jīng)有6人回復(fù)
matlab非線性參數(shù)擬合問題已經(jīng)有7人回復(fù)
【求助】麻煩大家看看我的邀請信有沒有問題已經(jīng)有9人回復(fù)
【求助】大家看看，matlab安裝過程的截圖，找出錯(cuò) 已經(jīng)有14人回復(fù)
【求助】MATLAB多個(gè)mat中數(shù)據(jù)畫圖問題已經(jīng)有7人回復(fù)
【求助】大家?guī)臀铱纯催@個(gè)MATLAB程序？已經(jīng)有6人回復(fù)

1樓 2012-05-25 17:16:04

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

csgt0

榮譽(yù)版主 (著名寫手)

彩色掛圖

專家經(jīng)驗(yàn): +2
應(yīng)助: 367 (碩士)
貴賓: 0.457
金幣: 7807
散金: 8897
紅花: 20
沙發(fā): 7
帖子: 2809
在線: 515.8小時(shí)
蟲號: 449046
注冊: 2007-11-02
專業(yè): 金屬非晶態(tài)、準(zhǔn)晶和納米晶
管轄: 計(jì)算模擬

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

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
亮亮1283: 金幣+20, ★★★★★最佳答案, 哈哈，就是這個(gè)原因。太謝謝你了！ 2012-05-28 15:17:02

不知道對不對，畫出來的結(jié)果好像都是1*10^16
原程序主要有兩個(gè)問題，你的那些參數(shù)需要用全局變量，子函數(shù)不能直接用主函數(shù)的值的
還有xlabel，ylabel，zlabel，不是xlable等

function pdex1()

clc
clear
global Up Kb T P a y r h v Un e tp p0 n0 n1 p1 G tn
Up =400;
Un =8000;
Kb =1.38*10^(-23);
h =6.63*10^(-34);
e =1.6*10^(-19);
r =0.35;
y =0.6;
n0 =1*10^16;
p0 =1*10^16;
tn =1*10^(-9);
tp =1*10^(-9);
p1 =p0;
n1 =n0;
T =300;
a =1.65*10^4;
v =(3/532)*10^17;
X =10;
P=10;

%主函數(shù)----------------------------------------------------------------
m=0;
x=0:0.1:10;
t=0:20:1000;
sol = pdepe(m,@pdex1pde,@pdex1ic,@pdex1bc,x,t);
% u =[ sol(:,:,1), sol(:,:,2)];
figure('numbertitle','off','name','PDE Demo-by matlab')
subplot(211)
surf(x,t,sol(:,:,1))
title('The Solution of p')
xlabel('x')
ylabel('t')
zlabel('p')
subplot(212)
surf(x,t,sol(:,:,2))
title('The Solution of n')
xlabel('x')
ylabel('t')
zlabel('n')

%目標(biāo)函數(shù)----------------------------------------------------------------
function [c,f,s] = pdex1pde(x,t,u,DuDx)
   global Up Kb T P a y r h v Un e tp p0 n0 n1 p1 G tn
c= [1;1];
f= [(Up*Kb*T/e)*DuDx(1); (Un*Kb*T/e)*DuDx(2)];
G= a*y*(1-r)*P*exp(-a*x)/(h*v);
R1= u(1)*u(2)-p0*n0;
R2=tn*(u(1)+p1)+tp*(u(2)+n1);
R=R1/R2;
s= [1;1].*(G-R);

%初值條件函數(shù)------------------------------------------------------------
function u0 = pdex1ic(x)
global p0 n0
u0= [p0;n0];

%邊界條件函數(shù)------------------------------------------------------------
function [pa,qa,pb,qb] = pdex1bc(xa,ua,xb,ub,t)    %a為下邊界， b為上邊界
pa =[0; 0];
qa =[1; 1];
pb =[0; 0];
qb =[1; 1];

贊一下

回復(fù)此樓

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

亮亮1283

showmethemoney

4樓2012-05-28 14:35:31

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

查看全部 7 個(gè)回答

czxf

銀蟲 (正式寫手)

應(yīng)助: 118 (高中生)
金幣: 242.9
散金: 3033
紅花: 8
帖子: 843
在線: 364.3小時(shí)
蟲號: 1254739
注冊: 2011-04-04
專業(yè): 電磁場與波

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

感謝參與，應(yīng)助指數(shù) +1

樓主，你的函數(shù)文件能正常調(diào)用嗎？有些格式都是錯(cuò)的。注意FUNCTION申明的是函數(shù)文件，而不是c中的函數(shù)

贊一下

回復(fù)此樓

2樓2012-05-26 21:43:32

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

亮亮1283

新蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 378
散金: 18
帖子: 156
在線: 70小時(shí)
蟲號: 890575
注冊: 2009-11-01
專業(yè): 等離子體物理

引用回帖:

2樓: Originally posted by czxf at 2012-05-26 21:43:32
樓主，你的函數(shù)文件能正常調(diào)用嗎？有些格式都是錯(cuò)的。注意FUNCTION申明的是函數(shù)文件，而不是c中的函數(shù)

實(shí)在是不好意思啊，我不是很懂你說的什么意思，這個(gè)程序是我參考matlab幫助里面一個(gè)專門求解這樣的方程組的例子寫的。能不能麻煩你把錯(cuò)誤的地方給我指一下，我好再看看。

贊一下

回復(fù)此樓

3樓2012-05-27 12:52:04

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

亮亮1283

新蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 378
散金: 18
帖子: 156
在線: 70小時(shí)
蟲號: 890575
注冊: 2009-11-01
專業(yè): 等離子體物理

引用回帖:

4樓: Originally posted by csgt0 at 2012-05-28 14:35:31
不知道對不對，畫出來的結(jié)果好像都是1*10^16
原程序主要有兩個(gè)問題，你的那些參數(shù)需要用全局變量，子函數(shù)不能直接用主函數(shù)的值的
還有xlabel，ylabel，zlabel，不是xlable等

function pdex1()

clc
c ...

問題解決了

回復(fù)此樓

5樓2012-05-28 15:20:24

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

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

[考研] 一志愿a區(qū)211，085601-307分求調(diào)劑 +5	黨嘉豪 2026-03-31	10/500	2026-03-31 13:19 by Dyhoer
[考研] 本科211生物醫(yī)學(xué)工程085409求調(diào)劑339分 +5	里子木yy 2026-03-29	5/250	2026-03-31 10:44 by WWM123
[考研] 282求調(diào)劑 +3	呼吸都是減肥 2026-03-31	3/150	2026-03-31 10:00 by wangjy2002
[考研] 哈爾濱工業(yè)大學(xué)材料與化工專碩378求調(diào)劑 +3	塔比烏斯 2026-03-30	3/150	2026-03-30 22:55 by 無際的草原
[考研] 286分調(diào)劑 +10	Faune 2026-03-30	11/550	2026-03-30 20:35 by 啊李999
[考研] 一志愿華中師范化學(xué)332分求調(diào)劑 +3	Lyy930824@ 2026-03-29	3/150	2026-03-30 20:15 by DHUSHUAI
[考研] 26考研-291分-廈門大學(xué)（085601）-柔性電子學(xué)院材料工程專業(yè)求調(diào)劑 +5	min3 2026-03-24	6/300	2026-03-30 18:42 by 544594351
[考研] 求調(diào)劑 +10	家佳佳佳佳佳 2026-03-29	10/500	2026-03-30 18:34 by 544594351
[考研] 292求調(diào)劑 +13	是妍子也是研子 2026-03-30	13/650	2026-03-30 18:01 by 小徐0109
[考研] 283求調(diào)劑（080500） +14	A child 2026-03-27	14/700	2026-03-30 12:06 by 探123
[考研] 332求調(diào)劑 +6	Lyy930824@ 2026-03-29	6/300	2026-03-30 10:53 by 1172367218
[考研] 356求調(diào)劑 +4	gysy?s?a 2026-03-28	4/200	2026-03-29 10:32 by 唐沐兒
[考研] 調(diào)劑 +3	好好讀書。 2026-03-28	3/150	2026-03-28 12:04 by 王保杰33
[考研] 081200-314 +3	LILIQQ 2026-03-27	4/200	2026-03-28 09:41 by 保護(hù)地球你我做?/a>
[考研] 331環(huán)境科學(xué)與工程求調(diào)劑 +3	熠然好運(yùn)氣 2026-03-27	3/150	2026-03-28 04:11 by fmesaito
[考研] 315分求調(diào)劑 +7	26考研上岸版26 2026-03-26	7/350	2026-03-28 04:05 by fmesaito
[考研] 070300化學(xué)求調(diào)劑 +4	起個(gè)名咋這么難 2026-03-27	4/200	2026-03-27 21:39 by 83503孫老師
[考研] 一志愿上海理工能源動(dòng)力（085800）310分求調(diào)劑 +3	zhangmingc 2026-03-27	4/200	2026-03-27 19:01 by 給你你注意休息
[考研] 0856調(diào)劑 +5	求求讓我有書讀?/a> 2026-03-26	6/300	2026-03-27 15:12 by caszguilin
[考研] 材料求調(diào)劑 +5	.m.. 2026-03-25	5/250	2026-03-27 11:08 by 不吃魚的貓