 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 2940 | 回復: 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]; |
回復此樓» 猜你喜歡
（081700）化學工程與技術-298分求調劑
已經(jīng)有4人回復
-
求調劑
已經(jīng)有5人回復
-
一志愿中科大材料與化工。353分能調劑到哪？
已經(jīng)有7人回復
-
337一志愿華南理工0805材料求調劑
已經(jīng)有5人回復
-
伙伴們，祝我生日快樂吧
已經(jīng)有4人回復
-
301求調劑
已經(jīng)有3人回復
-
材料與化工304求B區(qū)調劑
已經(jīng)有4人回復
-
讀博申請
已經(jīng)有7人回復
-
一志愿江南大學085701環(huán)境工程專碩總分287求調劑
已經(jīng)有5人回復
-
材料工程330分求調劑，一志愿985
已經(jīng)有5人回復
|
你先把該合并的項合并了，該化簡的常數(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ā)表 | |
|---|---|---|---|---|
|
[考研] 274環(huán)境工程求調劑 +8 | 扶柳盈江 2026-03-05 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 考研材料與化工，求調劑 +7 | 戲精丹丹丹 2026-03-09 | 7/350 |
|
|
|
[考研] 材料工程307，求調劑 +7 | 我要燃燒你的夢 2026-03-08 | 7/350 |
|
|
|
[考研] 086000生物與醫(yī)藥319分求調劑 +4 | Tolkien 2026-03-07 | 8/400 |
|
|
|
[考研] 332材料求調劑 +6 | zjy101327 2026-03-05 | 8/400 |
|
|
|
[考研] 求調劑！ +3 | 朔朔話 2026-03-09 | 3/150 |
|
|
|
[考研] 0703化學調劑 +3 | 三dd. 2026-03-10 | 3/150 |
|
|
|
[考研] A區(qū)一本交叉課題組，低分調劑，招收機械電子信息通信等交叉方向 +56 | lisimayy 2026-03-04 | 73/3650 |
|
|
|
[考研] 研究生招生 +3 | 徐海濤11 2026-03-10 | 6/300 |
|
|
|
[考研] 考研調劑，一志愿山東大學材料與化工，328分，政治51 +5 | 關你西紅柿929 2026-03-08 | 6/300 |
|
|
|
[考研] 求0856調劑 +5 | 小力氣珂珂 2026-03-08 | 5/250 |
|
|
|
[考研] 一志愿211 材料與化工 280求調劑 +3 | Sanity蔣 2026-03-08 | 3/150 |
|
|
|
[考博] 26年博士申請 +4 | 科研狗111 2026-03-07 | 4/200 |
|
|
|
[考研] 求調劑 +3 | 拾柒12。 2026-03-08 | 3/150 |
|
|
|
[考研] 346分材料求調劑 +5 | snow_反季節(jié)版 2026-03-07 | 5/250 |
|
|
|
[考研] 第一志愿上海大學，專業(yè)化學工程與技術，總分288，求調劑 +3 | 1829197082 2026-03-07 | 3/150 |
|
|
|
[考研] 一志愿211 085600 280數(shù)二英二求調劑 +3 | 月山斜 2026-03-06 | 3/150 |
|
|
|
[考研] 287求調劑 +3 | 看看我. 2026-03-05 | 6/300 |
|
|
|
[考博] 2026申博自薦 六級440電催化方向 +4 | 櫻落成影花成雙 2026-03-05 | 4/200 |
|
|
|
[考研] 316求調劑 +3 | 林小星發(fā)大財 2026-03-04 | 3/150 |
|