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

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

24小時(shí)熱門版塊排行榜

返回列表

亮亮1283

新蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 378
散金: 18
帖子: 156
在線: 70小時(shí)
蟲號(hào): 890575
注冊(cè): 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ù)此樓

» 猜你喜歡

0703化學(xué)調(diào)劑，求導(dǎo)師收已經(jīng)有3人回復(fù)
生物學(xué)學(xué)碩求調(diào)劑已經(jīng)有6人回復(fù)
305分求調(diào)劑（食品工程）已經(jīng)有7人回復(fù)
材料與化工考研調(diào)劑已經(jīng)有5人回復(fù)
299求調(diào)劑已經(jīng)有6人回復(fù)
材料292調(diào)劑已經(jīng)有5人回復(fù)
085404電子信息284分求調(diào)劑已經(jīng)有4人回復(fù)
收08調(diào)劑生已經(jīng)有3人回復(fù)
081700 調(diào)劑 267分已經(jīng)有9人回復(fù)
293求調(diào)劑已經(jīng)有5人回復(fù)

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

求助matlab問題已經(jīng)有5人回復(fù)
請(qǐng)教關(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ù)
【求助】麻煩大家看看我的邀請(qǐng)信有沒有問題已經(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的回帖

czxf

銀蟲 (正式寫手)

應(yīng)助: 118 (高中生)
金幣: 242.9
散金: 3033
紅花: 8
帖子: 843
在線: 364.3小時(shí)
蟲號(hào): 1254739
注冊(cè): 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í)
蟲號(hào): 890575
注冊(cè): 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的回帖

csgt0

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

彩色掛圖

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

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

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

不知道對(duì)不對(duì)，畫出來的結(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的回帖

亮亮1283

新蟲 (小有名氣)

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

引用回帖:

4樓: Originally posted by csgt0 at 2012-05-28 14:35:31
不知道對(duì)不對(duì)，畫出來的結(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的回帖

亮亮1283

新蟲 (小有名氣)

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

送鮮花一朵

引用回帖:

不好意思，我還想問一下，你知道我的方程組是一個(gè)拋物線型的非線性偏微分方程組，方程組中的這兩個(gè)方程都是拋物線型的，如果我要是再加上一個(gè)方程，但是方程的形式是雙曲型的，應(yīng)該怎么來解決呢？

贊一下

回復(fù)此樓

6樓2012-05-28 15:35:59

已閱回復(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í)
蟲號(hào): 449046
注冊(cè): 2007-11-02
專業(yè): 金屬非晶態(tài)、準(zhǔn)晶和納米晶
管轄: 計(jì)算模擬

這個(gè)我不太懂，你用hyperbolic試試

贊一下

回復(fù)此樓

showmethemoney

7樓2012-05-28 16:15:44

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

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

返回列表

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

[基金申請(qǐng)] 請(qǐng)教下大家 2026年國家基金申請(qǐng)是雙盲審嗎？ +3	lishucheng1 2026-03-22	5/250	2026-03-24 08:22 by gltch
[考研] 335求調(diào)劑 +4	yuyu宇 2026-03-23	5/250	2026-03-23 23:49 by Txy@872106
[考研] 求調(diào)劑 +7	十三加油 2026-03-21	7/350	2026-03-23 23:48 by 熱情沙漠
[考研] 一志愿國科過程所081700，274求調(diào)劑 +3	三水研0水立方 2026-03-23	3/150	2026-03-23 23:11 by MajorWen
[考研] 284求調(diào)劑 +3	yanzhixue111 2026-03-23	6/300	2026-03-23 22:58 by pswait
[考研] 303求調(diào)劑 +4	元夕元 2026-03-20	4/200	2026-03-23 19:00 by macy2011
[考研] 北科281學(xué)碩材料求調(diào)劑 +8	tcxiaoxx 2026-03-20	9/450	2026-03-23 12:16 by tcxiaoxx
[考研] 求調(diào)劑材料學(xué)碩080500，總分289分 5+3	@taotao 2026-03-19	21/1050	2026-03-23 10:17 by 冠c哥
[考研] 求調(diào)劑一志愿海大，0703化學(xué)學(xué)碩304分，有大創(chuàng)項(xiàng)目，四級(jí)已過 +6	幸運(yùn)哩哩 2026-03-22	10/500	2026-03-22 20:10 by edmund7
[考研] 303求調(diào)劑 +5	安憶靈 2026-03-22	6/300	2026-03-22 12:46 by 素顏傾城1988
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +3	邱gl 2026-03-20	7/350	2026-03-21 19:05 by 15709483992
[考研] 297求調(diào)劑 +3	喜歡還是不甘心 2026-03-20	3/150	2026-03-21 18:33 by 學(xué)員8dgXkO
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境總分308求調(diào)劑 +7	墨墨漠 2026-03-20	7/350	2026-03-21 16:36 by barlinike
[考研] 279求調(diào)劑 +5	紅衣隱官 2026-03-21	5/250	2026-03-21 14:59 by lature00
[考研] 22 350 本科985求調(diào)劑，求老登收留 +3	李軼男003 2026-03-20	3/150	2026-03-21 13:28 by 搏擊518
[考研] 303求調(diào)劑 +5	睿08 2026-03-17	7/350	2026-03-21 03:11 by JourneyLucky
[考研] 304求調(diào)劑 +7	司空. 2026-03-18	7/350	2026-03-20 23:08 by JourneyLucky
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +11	Ymlll 2026-03-18	15/750	2026-03-20 19:40 by 丁丁*
[考研] 求調(diào)劑 +3	@taotao 2026-03-20	3/150	2026-03-20 19:35 by JourneyLucky
[考研] 一志愿南理工085701環(huán)境302求調(diào)劑院校 +3	葵梓衛(wèi)隊(duì) 2026-03-20	3/150	2026-03-20 19:28 by zhukairuo

24小時(shí)熱門版塊排行榜

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

» 猜你喜歡

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

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

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

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

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

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