版塊導航: 正在加載中...

調(diào)劑小程序

登錄注冊

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

24小時熱門版塊排行榜

>論壇更新日志 (6119)
>考研 (1487)
>導師招生 (1255)
>文獻求助 (464)
>蟲友互識 (274)
>考博 (244)
>碩博家園 (133)
>博后之家 (111)
>基金申請 (108)
>招聘信息布告欄 (89)
>論文投稿 (85)
>休閑灌水 (68)
>第一性原理 (65)
>教師之家 (64)
>公派出國 (41)
>外文書籍求助 (32)

返回列表

亮亮1283

新蟲 (小有名氣)

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

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

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

%目標函數(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];

回復此樓

» 猜你喜歡

材料學碩333求調(diào)劑已經(jīng)有3人回復
289求調(diào)劑已經(jīng)有7人回復
收08調(diào)劑生已經(jīng)有6人回復
材料專碩找調(diào)劑已經(jīng)有5人回復
085602 289分求調(diào)劑已經(jīng)有5人回復
化工專碩求調(diào)劑已經(jīng)有3人回復
274求調(diào)劑已經(jīng)有5人回復
調(diào)劑已經(jīng)有3人回復
08工學調(diào)劑已經(jīng)有15人回復
招08考數(shù)學已經(jīng)有16人回復

» 本主題相關價值貼推薦，對您同樣有幫助:

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

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

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

czxf

銀蟲 (正式寫手)

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

【答案】應助回帖

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

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

贊一下

回復此樓

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

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

亮亮1283

新蟲 (小有名氣)

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

引用回帖:

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

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

贊一下

回復此樓

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

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

csgt0

榮譽版主 (著名寫手)

彩色掛圖

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

【答案】應助回帖

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

不知道對不對，畫出來的結(jié)果好像都是1*10^16
原程序主要有兩個問題，你的那些參數(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')

%目標函數(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];

贊一下

回復此樓

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

亮亮1283

showmethemoney

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

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

亮亮1283

新蟲 (小有名氣)

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

引用回帖:

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

function pdex1()

clc
c ...

問題解決了

回復此樓

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

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

亮亮1283

新蟲 (小有名氣)

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

送鮮花一朵

引用回帖:

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

贊一下

回復此樓

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

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

csgt0

榮譽版主 (著名寫手)

彩色掛圖

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

這個我不太懂，你用hyperbolic試試

贊一下

回復此樓

showmethemoney

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

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

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

返回列表

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

[考研] 0703化學調(diào)劑，求導師收 +6	天天好運來上岸?/a> 2026-03-24	6/300	2026-03-24 17:19 by dick_runner
[考研] 材料學碩，求調(diào)劑 6+3	糖葫蘆888ll 2026-03-22	7/350	2026-03-24 17:11 by hello七七
[考研] 一志愿北京化工大學材料與化工 264分各科過A區(qū)國家線 +3	哈哈157349 2026-03-21	3/150	2026-03-24 14:11 by zhyzzh
[考研] 344求調(diào)劑 +3	desto 2026-03-24	3/150	2026-03-24 10:09 by 搏擊518
[考研] 一志愿哈工大，085400，320，求調(diào)劑 +3	gdlf9999 2026-03-24	3/150	2026-03-24 10:08 by 搏擊518
[考研] 07化學280分求調(diào)劑 +4	722865 2026-03-23	4/200	2026-03-24 00:01 by chixmc
[考研] 材料/農(nóng)業(yè)專業(yè)，07/08開頭均可，過線就行 +3	呵唔哦豁 2026-03-23	4/200	2026-03-23 22:30 by 汪！？！
[考研] 328求調(diào)劑 +4	LHHL66 2026-03-23	4/200	2026-03-23 14:55 by lbsjt
[考研] 287求調(diào)劑 +8	晨昏線與星海 2026-03-19	9/450	2026-03-22 17:01 by i_cooler
[考研] 319求調(diào)劑 +4	小力氣珂珂 2026-03-20	4/200	2026-03-22 15:53 by ColorlessPI
[考研] 求調(diào)劑 +5	Zhangbod 2026-03-21	7/350	2026-03-22 13:13 by Zhangbod
[考研] 材料學碩301分求調(diào)劑 +7	Liyouyumairs 2026-03-21	7/350	2026-03-21 22:31 by peike
[考研] 求調(diào)劑 +3	13341 2026-03-20	3/150	2026-03-21 18:28 by 學員8dgXkO
[考研] 311求調(diào)劑 +3	勇敢的小吳 2026-03-20	3/150	2026-03-21 17:40 by ColorlessPI
[考研] 294求調(diào)劑材料與化工專碩 +15	陌の森林 2026-03-18	15/750	2026-03-20 23:28 by JourneyLucky
[考研] 288求調(diào)劑 +16	于海海海海 2026-03-19	16/800	2026-03-20 22:28 by JourneyLucky
[考研] 295復試調(diào)劑 +8	簡木ChuFront 2026-03-19	8/400	2026-03-20 20:44 by zhukairuo
[考研] 求調(diào)劑 +3	@taotao 2026-03-20	3/150	2026-03-20 19:35 by JourneyLucky
[考研] 085410人工智能專碩317求調(diào)劑（0854都可以） +4	xbxudjdn 2026-03-18	4/200	2026-03-20 09:07 by 不168
[考研] 288求調(diào)劑，一志愿華南理工大學071005 +5	ioodiiij 2026-03-17	5/250	2026-03-19 18:22 by zcl123

24小時熱門版塊排行榜

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

» 猜你喜歡

» 本主題相關價值貼推薦，對您同樣有幫助:

【答案】應助回帖

【答案】應助回帖

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

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

» 本主題相關價值貼推薦，對您同樣有幫助: