| 查看: 707 | 回复: 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%时间节点 %-----------定义平衡方程------------------------------ 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)尝试了下，结果不理想，求大侠赐教 [ Last edited by 小木虫zb on 2013-5-21 at 16:46 ] |
回复此楼» 猜你喜欢
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入机器视觉与3D光学成像课题组。
已经有0人回复
-
**
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有71人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有3人回复
-
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有2人回复
-
B区学生调剂-兰州交通大学材料科学与工程学院
已经有11人回复
-
山西大同大学物理学专业还有调剂名额，欢迎调剂！
已经有12人回复
-
桂林理工大学物理学专业招收调剂
已经有12人回复
-
VASP 的一组 GPU / CPU 基准测试记录
已经有0人回复
-
美国顶级物理期刊《应用物理快报》（APL）的编辑欺骗和歧视作者及AIP的官僚主义傲慢
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab非线性拟合求助
已经有7人回复
- MATLAB求解非线性方程组
已经有5人回复
- 求Matlab解非线性方程代码
已经有12人回复
- 求助matlab解一个非线性方程组
已经有8人回复
- matlab解非线性方程组解不开
已经有3人回复
- 拟牛顿法求解复杂（带根号的）非线性方程组
已经有17人回复
- 求助matlab 非线性方程组全部解！
已经有6人回复
- matlab求解非线性方程组，错误提示怎么解决
已经有5人回复
- 求助 MATLAB解方程组-fslove
已经有7人回复
- 求助matlab---fsolve解非线性方程组
已经有6人回复
- 用不动点迭代求非线性方程组，求得的不是想要的解，怎么办
已经有12人回复
- matlab求解非线性方程组，求助！
已经有6人回复
- 求助！matlab用fsolve函数求解非线性方程组的问题！
已经有19人回复
- matlab求解非线性方程组
已经有16人回复
- 求助有关Matlab有约束非线性最优化问题
已经有9人回复
- 求matlab求解一个方程组的问题
已经有6人回复
- 求高人指点用matlab求解非线性方程组，解决了追加100金币；
已经有11人回复
- matlab的fsove 命令求解非线性方程组
已经有6人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】用mathematica 5.0求解一个非线性方程组失败，特发帖求助！
已经有5人回复
- 【求助】求教matlab解非线性方程组
已经有9人回复
thermalE
木虫 (正式写手)
- 应助: 24 (小学生)
- 金币: 1999.2
- 红花: 8
- 帖子: 864
- 在线: 357.2小时
- 虫号: 2057299
- 注册: 2012-10-13
- 性别: GG
- 专业: 传递过程
2楼2013-05-22 14:56:33
小木虫zb
木虫 (正式写手)
- 应助: 11 (小学生)
- 金币: 1560.9
- 红花: 5
- 帖子: 345
- 在线: 142.1小时
- 虫号: 1227944
- 注册: 2011-03-10
- 性别: GG
- 专业: 分离过程
3楼2013-05-23 19:40:56
5