 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 2941 | 回復: 3 | |||
[求助]
Matlab求解非線性偏微分方程組 已有1人參與
|
|
在課題研究中建立了一個非線性偏微分方程組，用Matlab中pdepe函數(shù)求解時，無法求出結果，Matlab一直正忙，我不是數(shù)學專業(yè)的，一直找不出問題，懇請各位指導一下。謝謝！ 程序如下 function pdexx11 m = 0; x = linspace(0,0.02,3); t = linspace(0,1,3); sol = pdepe(m,@pdefun,@icfun,@bcfun,x,t); u1 = sol(:,:,1); u2 = sol(:,:,2); u3 = sol(:,:,3); figure; surf(x,t,u1) title('u1(x,t)') xlabel('Distance x') ylabel('Time t') figure; surf(x,t,u2) title('u2(x,t)') xlabel('Distance x') ylabel('Time t') figure; surf(x,t,u3) title('u3(x,t)') xlabel('Distance x') ylabel('Time t') function [ c,f,s ] = pdefun( x,t,u,DuDx ) c = [1;1;1]; f = [1.25e-9; 1.828e-9; 1.856e-9].*DuDx+[0.59*u(1);0.59*u(2);0.59*u(3)]; s = [-1.03e-4/60*1.9*((u(1)*u(2)*u(3)/(5.51e-14))^(1/3)-1); -1.03e-4/60*1.9*((u(1)*u(2)*u(3)/(5.51e-14))^(1/3)-1); -1.03e-4/60*1.9*((u(1)*u(2)*u(3)/(5.51e-14))^(1/3)-1)]; function u0 = icfun(x) u0 = [0;0;0]; function [ pl,ql,pr,qr ] = bcfun( xl,ul,xr,ur,t ) pl = [ul(1)-0.0493; ul(2)-0.07142; ul(3)-0.0387]; ql = [0;0;0]; pr = [ur(1)-0.0106;ur(2)-0.03272;ur(3)]; qr = [0;0;0]; |
回復此樓» 猜你喜歡
一志愿天津大學材料與化工275求調(diào)劑
已經(jīng)有22人回復
-
0703求調(diào)劑
已經(jīng)有5人回復
-
0857 資源與環(huán)境 285分
已經(jīng)有6人回復
-
材料工程,326分，求調(diào)劑
已經(jīng)有5人回復
-
301求調(diào)劑
已經(jīng)有5人回復
-
334求調(diào)劑
已經(jīng)有3人回復
-
材料與化工085600調(diào)劑求老師收留
已經(jīng)有5人回復
-
一志愿華科071000生物學 338求調(diào)劑
已經(jīng)有4人回復
-
290求調(diào)劑
已經(jīng)有4人回復
-
308求調(diào)劑
已經(jīng)有4人回復
|
你先把該合并的項合并了，該化簡的常數(shù)化簡了，再試試吧。比如K_sp值很小，做了除數(shù)，又開方，讓程序算很容易掛掉，你該把它提出來，和k Ac合并起來算。再比如D的量級是10的-9次方，后面系數(shù)量級為1，你可以先把空間偏導項忽略試試 發(fā)自小木蟲Android客戶端 |
2樓2019-03-01 20:18:56
3樓2019-03-01 20:25:27
ofwhy
至尊木蟲 (著名寫手)
- 應助: 3 (幼兒園)
- 金幣: 13429.6
- 散金: 93
- 紅花: 7
- 帖子: 2309
- 在線: 177.9小時
- 蟲號: 2777064
- 注冊: 2013-11-04
- 性別: GG
- 專業(yè): 工程熱力學
【答案】應助回帖
|
可能公式推導有錯誤，如果按照你導出的公式計算，會計算出虛值，然后不停的死循環(huán)。如將pu3/pt的最后k*Ac前面的負號改為正號則有解。從你的邊界條件來看，似乎這個分量也應該跟其他兩個應該是不一樣的。 function [ output_args ] = PDEs1( input_args ) m=0; x=linspace(0,0.02,3); t=linspace(0,1,3); sol=pdepe(m,@pdes1pde,@pdes1ic,@pdes1bc,x,t); u1=sol(:,:,1); u2=sol(:,:,2); u3=sol(:,:,3); figure surf(x,t,u1) title('u1(x,t)') xlabel('Distance x') ylabel('Time t') figure surf(x,t,u2) title('u2(x,t)') xlabel('Distance x') ylabel('Time t') figure surf(x,t,u3) title('u3(x,t)') xlabel('Distance x') ylabel('Time t') function[c,f,s]=pdes1pde(x,t,u,DuDx) c=[1;1;1]; D1=1.25*10^(-9); D2=1.838*10^(-9); D3=1.856*10^(-9); V=0.59; k=1.03*10^(-4)/60; Ac=1.9; Ksp=5.51*10^(-14); f=[D1/(1+V);D2/(1+V);D3/(1+V)].*DuDx; s=[-k*Ac/(1+V)*((u(1)*u(2)*u(3)/Ksp)^(1/3)-1);-k*Ac/(1+V)*((u(1)*u(2)*u(3)/Ksp)^(1/3)-1);k*Ac/(1+V)*((u(1)*u(2)*u(3)/Ksp)^(1/3)-1);]; end function u0 = pdes1ic(x) u0 = [0;0;0]; end function [pl,ql,pr,qr] = pdes1bc(xl,ul,xr,ur,t) pl = [ul(1)-0.0493; ul(2)-0.07142;ul(3)-0.0387]; ql = [0;0;0]; pr = [ur(1)-0.0106;ur(2)-0.03272;ur(3)]; qr = [0;0;0]; end end |
4樓2019-03-31 12:08:01
相關版塊跳轉
我要訂閱樓主 iMitsuha 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 291分獸醫(yī)專碩 一志愿華中農(nóng) +3 | 種大大太陽 2026-03-05 | 3/150 |
|
|---|---|---|---|---|
|
|
[碩博家園] 深圳大學碩士招生（2026秋，傳感器方向，僅錄取第一志愿） +4 | xujiaoszu 2026-03-11 | 6/300 |
|
|
|
[考研] 調(diào)劑 +5 | 調(diào)劑的考研學生 2026-03-09 | 5/250 |
|
|
|
[考研] 調(diào)劑 +5 | 呵唔哦豁 2026-03-10 | 5/250 |
|
|
[考研] 085602化工求調(diào)劑 +7 | 董boxing 2026-03-10 | 7/350 |
|
|
|
[考研] 0817學碩華東區(qū)求調(diào)劑 +3 | 30660438 2026-03-08 | 3/150 |
|
|
|
[考研] 材料與化工求調(diào)劑 +11 | 與冬清寧 2026-03-07 | 14/700 |
|
|
|
[考博] 找博導 +4 | 小呆呆熊 2026-03-07 | 4/200 |
|
|
|
[考研] 家人們 調(diào)劑不迷路 看這里 +8 | likeihood 2026-03-09 | 13/650 |
|
|
|
[考研] 070300化學求調(diào)劑 +5 | 撲風鈴的貓 2026-03-08 | 10/500 |
|
|
|
[考研] 中科大材料299求調(diào)劑 +10 | DAIjiayo 2026-03-05 | 16/800 |
|
|
|
[考研] 290 材料與化工求調(diào)劑 +7 | Nebulala 2026-03-08 | 7/350 |
|
|
|
[考研] 297求調(diào)劑 +3 | 胡達靈 2026-03-05 | 5/250 |
|
|
|
[考研] 理學287分求調(diào)劑 +4 | 看看我. 2026-03-05 | 4/200 |
|
|
|
[考研] 301求調(diào)劑 +11 | 朝天椒的雙馬尾 2026-03-05 | 11/550 |
|
|
|
[考研] 085701環(huán)境工程專業(yè)，初試305，均過國家A區(qū)線 +7 | 卡卡來了@ 2026-03-07 | 8/400 |
|
|
|
[考研] 08工科求調(diào)劑 +3 | 隆LLL 2026-03-06 | 4/200 |
|
|
|
[考研] 材料與化工354調(diào)劑 +4 | Lucy-xiao 2026-03-06 | 7/350 |
|
|
|
[考研] 334求調(diào)劑 +6 | Trying] 2026-03-05 | 8/400 |
|
|
|
[考研] 紡織、生物、化學、材料等專業(yè) +3 | Eember. 2026-03-05 | 7/350 |
|