| 5 | 1/1 | 返回列表 |
| 查看: 698 | 回復: 2 | ||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||
小木蟲zb木蟲 (正式寫手)
|
[求助]
求助matlab解一個復雜非線性方程組
|
|
|
解一個偏微分方程組,用有限差分法離散出了下列復雜的非線性方程組,不知道用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%時間節(jié)點 %-----------定義平衡方程------------------------------ 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ā)表 | |
|---|---|---|---|---|
|
[考研] 材料科學與工程339求調(diào)劑 +3 | hyz0119 2026-03-31 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 復試調(diào)劑 +5 | 雙馬尾痞老板2 2026-03-31 | 5/250 |
|
|
[考研] 08開頭看過來。! +3 | wwwwffffff 2026-03-31 | 5/250 |
|
|
[考研] 生物學296求調(diào)劑 +8 | 湯圓包 2026-03-29 | 12/600 |
|
|
[考研] 322求調(diào)劑:一志愿湖南大學 材料與化工(085600),已過六級。 +10 | XX小鄧 2026-03-29 | 10/500 |
|
|
[考研] 考研調(diào)劑求助 +7 | 13287130938 2026-03-31 | 7/350 |
|
|
[考研] 0856 335分 | 封裝or激光加工 老師看看我!! +9 | cccchenso 2026-03-29 | 9/450 |
|
|
[考研]
|
nnnnnnn5 2026-03-25 | 12/600 |
|
|
[考研] 343求調(diào)劑 +8 | 愛羈絆 2026-03-28 | 8/400 |
|
|
[考研] 求收留 +8 | 1943443204 2026-03-28 | 8/400 |
|
|
[考研] 320分,材料與化工專業(yè),求調(diào)劑 +10 | 一定上岸aaa 2026-03-27 | 14/700 |
|
|
[考研] 08工科求調(diào)劑286 +5 | tgs_001 2026-03-28 | 5/250 |
|
|
[考研] 303求調(diào)劑 +7 | DLkz1314. 2026-03-30 | 7/350 |
|
|
[考研] 359求調(diào)劑 +5 | 王了個楠 2026-03-25 | 5/250 |
|
|
[考研] 284求調(diào)劑 +14 | junqihahaha 2026-03-26 | 15/750 |
|
|
[考研] 311求調(diào)劑 +6 | 冬十三 2026-03-24 | 6/300 |
|
|
[考研] 071000生物學求調(diào)劑,初試成績343 +7 | 小小甜面團 2026-03-25 | 7/350 |
|
|
[考研] 復試調(diào)劑 +3 | raojunqi0129 2026-03-28 | 3/150 |
|
|
[考研] 一志愿南京航空航天大學材料學碩求調(diào)劑 +3 | @taotao 2026-03-28 | 3/150 |
|
|
[考研] 285求調(diào)劑 +4 | AZMK 2026-03-27 | 7/350 |
|