 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1541 | 回復(fù): 6 | ||||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||||
亮亮1283新蟲 (小有名氣)
|
[求助]
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ù)此樓» 猜你喜歡
085600 295分求調(diào)劑
已經(jīng)有14人回復(fù)
-
環(huán)境工程297分求調(diào)劑一志愿杭高院
已經(jīng)有3人回復(fù)
-
理學(xué)07化學(xué) 303求調(diào)劑
已經(jīng)有7人回復(fù)
-
290分調(diào)劑求助
已經(jīng)有6人回復(fù)
-
362求調(diào)劑
已經(jīng)有6人回復(fù)
-
求調(diào)劑推薦 材料 304
已經(jīng)有17人回復(fù)
-
304求調(diào)劑
已經(jīng)有4人回復(fù)
-
張芳銘-中國農(nóng)業(yè)大學(xué)-環(huán)境工程專碩-298
已經(jīng)有8人回復(fù)
-
070300化學(xué)354求調(diào)劑
已經(jīng)有12人回復(fù)
-
08工科，295，接受跨專業(yè)調(diào)劑
已經(jīng)有3人回復(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ù)
亮亮1283
新蟲 (小有名氣)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 378
- 散金: 18
- 帖子: 156
- 在線: 70小時(shí)
- 蟲號(hào): 890575
- 注冊: 2009-11-01
- 專業(yè): 等離子體物理
3樓2012-05-27 12:52:04
2樓2012-05-26 21:43:32
csgt0
榮譽(yù)版主 (著名寫手)
彩色掛圖
-
專家經(jīng)驗(yàn): +2 - 應(yīng)助: 367 (碩士)
- 貴賓: 0.457
- 金幣: 7807
- 散金: 8897
- 紅花: 20
- 沙發(fā): 7
- 帖子: 2809
- 在線: 515.8小時(shí)
- 蟲號(hào): 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
感謝參與，應(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]; |
» 本帖已獲得的紅花（最新10朵）
亮亮1283
4樓2012-05-28 14:35:31
亮亮1283
新蟲 (小有名氣)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 378
- 散金: 18
- 帖子: 156
- 在線: 70小時(shí)
- 蟲號(hào): 890575
- 注冊: 2009-11-01
- 專業(yè): 等離子體物理
5樓2012-05-28 15:20:24
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 290分調(diào)劑求助 +6 | 吉祥止止陳 2026-03-25 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 362求調(diào)劑 +6 | 西南交材料專碩3 2026-03-31 | 6/300 |
|
|
|
[考研] 學(xué)碩274求調(diào)劑 +16 | Li李魚 2026-03-26 | 16/800 |
|
|
|
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +6 | hanamiko 2026-03-30 | 6/300 |
|
|
|
[考研] 哈爾濱工業(yè)大學(xué)材料與化工專碩378求調(diào)劑 +3 | 塔比烏斯 2026-03-30 | 3/150 |
|
|
|
[考研] 332求調(diào)劑 +14 | 032500 2026-03-25 | 14/700 |
|
|
|
[考研] 抱歉 +3 | 田洪有 2026-03-30 | 3/150 |
|
|
|
[考研] 334分 一志愿武理 材料求調(diào)劑 +16 | 李李不服輸 2026-03-26 | 16/800 |
|
|
|
[考研] 317分 一志愿南理工材料工程 本科湖工大 求調(diào)劑 +12 | 芋泥小鈴鐺 2026-03-28 | 12/600 |
|
|
|
[考研] 375求調(diào)劑 +6 | 雨夏整夜 2026-03-29 | 6/300 |
|
|
|
[考研] 考研調(diào)劑 +7 | 小蠟新筆 2026-03-29 | 7/350 |
|
|
|
[考研] 291求調(diào)劑 +5 | Y-cap 2026-03-29 | 6/300 |
|
|
|
[考研] 086000生物與醫(yī)藥調(diào)劑 +5 | Feisty。 2026-03-28 | 9/450 |
|
|
|
[考研] 085405 考的11408求各位老師帶走 +3 | Qiu學(xué)ing 2026-03-28 | 3/150 |
|
|
|
[考研] 材料求調(diào)劑一志愿哈工大324 +7 | 閆旭東 2026-03-28 | 9/450 |
|
|
|
[考研] 265求調(diào)劑11408 +3 | 劉小鹿lu 2026-03-27 | 3/150 |
|
|
|
[考研] 08開頭275求調(diào)劑 +4 | 拉誰不重要 2026-03-26 | 4/200 |
|
|
|
[考研] 324求調(diào)劑 +5 | hanamiko 2026-03-26 | 5/250 |
|
|
|
[考研] 303求調(diào)劑 +6 | 藍(lán)山月 2026-03-25 | 6/300 |
|
|
|
[考研] 347求調(diào)劑 +4 | L when 2026-03-25 | 4/200 |
|