 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 692 | 回復(fù): 2 | ||
小木蟲(chóng)zb木蟲(chóng) (正式寫(xiě)手)
|
[求助]
求助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 小木蟲(chóng)zb on 2013-5-21 at 16:46 ] |
回復(fù)此樓» 猜你喜歡
2026年第五屆電氣、電子與信息工程國(guó)際會(huì)議（ISEEIE 2026）
已經(jīng)有0人回復(fù)
-
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑
已經(jīng)有6人回復(fù)
-
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有154人回復(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ù)
-
請(qǐng)問(wèn)還有沒(méi)有用Latex寫(xiě)文章的小伙伴們？
已經(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à)值貼推薦，對(duì)您同樣有幫助:
- matlab非線性擬合求助
已經(jīng)有7人回復(fù)
- MATLAB求解非線性方程組
已經(jīng)有5人回復(fù)
- 求Matlab解非線性方程代碼
已經(jīng)有12人回復(fù)
- 求助matlab解一個(gè)非線性方程組
已經(jīng)有8人回復(fù)
- matlab解非線性方程組解不開(kāi)
已經(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ù)求解非線性方程組的問(wèn)題！
已經(jīng)有19人回復(fù)
- matlab求解非線性方程組
已經(jīng)有16人回復(fù)
- 求助有關(guān)Matlab有約束非線性最優(yōu)化問(wèn)題
已經(jīng)有9人回復(fù)
- 求matlab求解一個(gè)方程組的問(wèn)題
已經(jīng)有6人回復(fù)
- 求高人指點(diǎn)用matlab求解非線性方程組，解決了追加100金幣；
已經(jīng)有11人回復(fù)
- matlab的fsove 命令求解非線性方程組
已經(jīng)有6人回復(fù)
- 【求助】向各位大俠求助matlab求解微分方程組遇到的一個(gè)問(wèn)題
已經(jīng)有21人回復(fù)
- 【求助】用mathematica 5.0求解一個(gè)非線性方程組失敗，特發(fā)帖求助！
已經(jīng)有5人回復(fù)
- 【求助】求教matlab解非線性方程組
已經(jīng)有9人回復(fù)
thermalE
木蟲(chóng) (正式寫(xiě)手)
- 應(yīng)助: 24 (小學(xué)生)
- 金幣: 1999.2
- 紅花: 8
- 帖子: 864
- 在線: 357.2小時(shí)
- 蟲(chóng)號(hào): 2057299
- 注冊(cè): 2012-10-13
- 性別: GG
- 專(zhuān)業(yè): 傳遞過(guò)程
2樓2013-05-22 14:56:33
小木蟲(chóng)zb
木蟲(chóng) (正式寫(xiě)手)
- 應(yīng)助: 11 (小學(xué)生)
- 金幣: 1560.9
- 紅花: 5
- 帖子: 345
- 在線: 142.1小時(shí)
- 蟲(chóng)號(hào): 1227944
- 注冊(cè): 2011-03-10
- 性別: GG
- 專(zhuān)業(yè): 分離過(guò)程
3樓2013-05-23 19:40:56
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 小木蟲(chóng)zb 的主題更新
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 304求調(diào)劑 +4 | 曼殊2266 2026-03-18 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 一志愿蘇州大學(xué)材料求調(diào)劑，總分315（英一） +3 | sbdksD 2026-03-19 | 3/150 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +10 | 妮妮ninicgb 2026-03-15 | 14/700 |
|
|
|
[考研] 307求調(diào)劑 +9 | 冷笙123 2026-03-17 | 9/450 |
|
|
|
[考研] 281求調(diào)劑（0805） +9 | 煙汐憶海 2026-03-16 | 19/950 |
|
|
|
[考研] 【同濟(jì)軟件】軟件（085405）考研求調(diào)劑 +3 | 2026eternal 2026-03-18 | 3/150 |
|
|
|
[考研] 一志愿武理材料305分求調(diào)劑 +5 | 想上岸的鯉魚(yú) 2026-03-18 | 6/300 |
|
|
|
[考研] 311求調(diào)劑 +6 | 26研0 2026-03-15 | 6/300 |
|
|
|
[考研]
 0854，計(jì)算機(jī)類(lèi)招收調(diào)劑
+3
|
胡辣湯放糖 2026-03-15 | 6/300 |
|
|
|
[考研] 280求調(diào)劑 +6 | 咕嚕曉曉 2026-03-18 | 7/350 |
|
|
|
[考研] 334求調(diào)劑 +3 | 志存高遠(yuǎn)意在機(jī)?/a> 2026-03-16 | 3/150 |
|
|
|
[考研] 301求調(diào)劑 +4 | A_JiXing 2026-03-16 | 4/200 |
|
|
|
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級(jí)食品/生物方向碩士 +3 | 1喜春8 2026-03-17 | 5/250 |
|
|
|
[考研] 332求調(diào)劑 +6 | Zz版 2026-03-13 | 6/300 |
|
|
|
[考研] 一志愿蘇州大學(xué)材料工程（085601）專(zhuān)碩有科研經(jīng)歷三項(xiàng)國(guó)獎(jiǎng)兩個(gè)實(shí)用型專(zhuān)利一項(xiàng)省級(jí)立項(xiàng) +6 | 大火山小火山 2026-03-16 | 8/400 |
|
|
|
[考博] 26申博 +4 | 八6八68 2026-03-16 | 4/200 |
|
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|
|
|
[基金申請(qǐng)]
今年的國(guó)基金是打分制嗎？
 50+3
|
zhanghaozhu 2026-03-14 | 3/150 |
|
|
|
[考研] 266求調(diào)劑 +4 | 學(xué)員97LZgn 2026-03-13 | 4/200 |
|