 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 691 | 回復(fù): 2 | ||
小木蟲zb木蟲 (正式寫手)
|
[求助]
求助matlab解一個(gè)復(fù)雜非線性方程組
|
|
解一個(gè)偏微分方程組，用有限差分法離散出了下列復(fù)雜的非線性方程組，不知道用matlab怎么求解 kf=200;k1=0.50427;n1=1/0.29823;k2=0.17356;n2=1/0.44799; n11=21.9;n12=3118.9;k12=0.3;n22=-0.5;n21=-1;k21=0.1; c0=7.03;a=0.363;u=0.361032/a;c01=7.03;c02=7.03;L=8.31e-2;p=436.8/(1-0.363); m=20;dz=L/m;t=[20:20:1000]/60;t_end=t(end);n=20;dt=t_end/n; A=u*dt/(4*dz);B=(1-a)*p/a;D=dt*kf/p; for i=1:n%時(shí)間節(jié)點(diǎn) %-----------定義平衡方程------------------------------ q1(i,1:m)=k1*ce1(i,1:m).^(1/n1+1/n11)/(ce1(i,1:m).^(1/n11)+k12*ce2(i,1:m).^(1/n12)); q2(i,1:m)=k2*ce2(i,1:m).^(1/n2+1/n22)/(ce2(i,1:m).^(1/n22)+k21*ce1(i,1:m).^(1/n21)); %---------定義吸附速率方程----------------------------- q1(i+1,1:m)=D*(c1(i,1:m)-ce1(i,1:m))+q1(i,1:m); q2(i+1,1:m)=D*(c2(i,1:m)-ce2(i,1:m))+q2(i,1:m); %---------定義床層物料衡算方程------------------------- c1(i,1)=A*(c1(i,2)-2*c01+c1(i+1,2))+c1(i+1,1)+B*(q1(i+1,1)-q1(i,1));%邊界條件j=1 c2(i,1)=A*(c2(i,2)-2*c02+c2(i+1,2))+c2(i+1,1)+B*(q2(i+1,1)-q2(i,1));%邊界條件j=1 c1(i,m)=c1(i+1,m)+B*(q1(i+1,m)-q1(i,m));%j=m c2(i,m)=c2(i+1,m)+B*(q2(i+1,m)-q2(i,m));%j=m c1(i,2:m-1)=A*(c1(i,3:m)-c1(i,1:m-2)+c1(i+1,3:m)-c1(i+1,1:m-2))+c1(i+1,2:m-1)+B*(q1(i+1,2:m-1)-q1(i,2:m-1)); c2(i,2:m-1)=A*(c2(i,3:m)-c2(i,1:m-2)+c2(i+1,3:m)-c2(i+1,1:m-2))+c2(i+1,2:m-1)+B*(q2(i+1,2:m-1)-q2(i,2:m-1)); end 想求c1、c2,用disp(c1)嘗試了下，結(jié)果不理想，求大俠賜教 [ Last edited by 小木蟲zb on 2013-5-21 at 16:46 ] |
回復(fù)此樓» 猜你喜歡
2026年第五屆電氣、電子與信息工程國際會(huì)議（ISEEIE 2026）
已經(jīng)有0人回復(fù)
-
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑
已經(jīng)有6人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費(fèi)?
已經(jīng)有210人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑
已經(jīng)有5人回復(fù)
-
求調(diào)劑
已經(jīng)有0人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有0人回復(fù)
-
請問還有沒有用Latex寫文章的小伙伴們？
已經(jīng)有0人回復(fù)
-
光學(xué)工程學(xué)碩調(diào)劑信息
已經(jīng)有30人回復(fù)
-
歡迎加入課題組
已經(jīng)有0人回復(fù)
-
散金幣，求好運(yùn)，祝面上順利！
已經(jīng)有63人回復(fù)
-
華東師范大學(xué)芯片設(shè)計(jì)徐瓏真實(shí)水平咋樣
已經(jīng)有1人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:
- matlab非線性擬合求助
已經(jīng)有7人回復(fù)
- MATLAB求解非線性方程組
已經(jīng)有5人回復(fù)
- 求Matlab解非線性方程代碼
已經(jīng)有12人回復(fù)
- 求助matlab解一個(gè)非線性方程組
已經(jīng)有8人回復(fù)
- matlab解非線性方程組解不開
已經(jīng)有3人回復(fù)
- 擬牛頓法求解復(fù)雜（帶根號(hào)的）非線性方程組
已經(jīng)有17人回復(fù)
- 求助matlab 非線性方程組全部解！
已經(jīng)有6人回復(fù)
- matlab求解非線性方程組，錯(cuò)誤提示怎么解決
已經(jīng)有5人回復(fù)
- 求助 MATLAB解方程組-fslove
已經(jīng)有7人回復(fù)
- 求助matlab---fsolve解非線性方程組
已經(jīng)有6人回復(fù)
- 用不動(dòng)點(diǎn)迭代求非線性方程組，求得的不是想要的解，怎么辦
已經(jīng)有12人回復(fù)
- matlab求解非線性方程組，求助！
已經(jīng)有6人回復(fù)
- 求助！matlab用fsolve函數(shù)求解非線性方程組的問題！
已經(jīng)有19人回復(fù)
- matlab求解非線性方程組
已經(jīng)有16人回復(fù)
- 求助有關(guān)Matlab有約束非線性最優(yōu)化問題
已經(jīng)有9人回復(fù)
- 求matlab求解一個(gè)方程組的問題
已經(jīng)有6人回復(fù)
- 求高人指點(diǎn)用matlab求解非線性方程組，解決了追加100金幣；
已經(jīng)有11人回復(fù)
- matlab的fsove 命令求解非線性方程組
已經(jīng)有6人回復(fù)
- 【求助】向各位大俠求助matlab求解微分方程組遇到的一個(gè)問題
已經(jīng)有21人回復(fù)
- 【求助】用mathematica 5.0求解一個(gè)非線性方程組失敗，特發(fā)帖求助！
已經(jīng)有5人回復(fù)
- 【求助】求教matlab解非線性方程組
已經(jīng)有9人回復(fù)
thermalE
木蟲 (正式寫手)
- 應(yīng)助: 24 (小學(xué)生)
- 金幣: 1999.2
- 紅花: 8
- 帖子: 864
- 在線: 357.2小時(shí)
- 蟲號(hào): 2057299
- 注冊: 2012-10-13
- 性別: GG
- 專業(yè): 傳遞過程
2樓2013-05-22 14:56:33
小木蟲zb
木蟲 (正式寫手)
- 應(yīng)助: 11 (小學(xué)生)
- 金幣: 1560.9
- 紅花: 5
- 帖子: 345
- 在線: 142.1小時(shí)
- 蟲號(hào): 1227944
- 注冊: 2011-03-10
- 性別: GG
- 專業(yè): 分離過程
3樓2013-05-23 19:40:56
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 小木蟲zb 的主題更新
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085700資源與環(huán)境308求調(diào)劑 +7 | 墨墨漠 2026-03-18 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 081700化工學(xué)碩調(diào)劑 +3 | 【1】 2026-03-16 | 3/150 |
|
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-18 | 3/150 |
|
|
|
[考研] 求調(diào)劑，一志愿:南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分 +3 | @taotao 2026-03-19 | 3/150 |
|
|
|
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +7 | 小材化本科 2026-03-18 | 7/350 |
|
|
|
[教師之家] 焦慮 +9 | 水冰月月野兔 2026-03-13 | 13/650 |
|
|
|
[考研] 本科鄭州大學(xué)物理學(xué)院，一志愿華科070200學(xué)碩，346求調(diào)劑 +4 | 我不是一根蔥 2026-03-18 | 4/200 |
|
|
|
[考研] 材料專碩英一數(shù)二306 +5 | z1z2z3879 2026-03-18 | 5/250 |
|
|
|
[考研] 材料工程專碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 6/300 |
|
|
|
[考研] 一志愿華中科技大學(xué)，080502，354分求調(diào)劑 +4 | 守候夕陽CF 2026-03-18 | 4/200 |
|
|
|
[考研] 303求調(diào)劑 +4 | 睿08 2026-03-17 | 6/300 |
|
|
|
[考研] 材料，紡織，生物（0856、0710），化學(xué)招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|
|
|
[考研] 268求調(diào)劑 +6 | 簡單點(diǎn)0 2026-03-17 | 6/300 |
|
|
|
[考研] 293求調(diào)劑 +11 | zjl的號(hào) 2026-03-16 | 16/800 |
|
|
|
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級(jí)食品/生物方向碩士 +3 | 1喜春8 2026-03-17 | 5/250 |
|
|
|
[考研] 290求調(diào)劑 +3 | p asserby. 2026-03-15 | 4/200 |
|
|
|
[考研] 283求調(diào)劑 +3 | 聽風(fēng)就是雨； 2026-03-16 | 3/150 |
|
|
|
[考研] 中科大材料與化工319求調(diào)劑 +3 | 孟鑫材料 2026-03-14 | 3/150 |
|
|
|
[考研] 材料與化工 323 英一+數(shù)二+物化，一志愿：哈工大 本人本科雙一流 +4 | 自由的_飛翔 2026-03-13 | 5/250 |
|
|
|
[考研] 297一志愿上交085600求調(diào)劑 +5 | 指尖八千里 2026-03-14 | 5/250 |
|