| 查看: 687 | 回复: 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 ] |
回复此楼» 猜你喜欢
德国亥姆霍兹Hereon中心招收两位医用镁合金腐蚀与LPSO相变方向2026公派博士生
已经有0人回复
-
推荐一款可以AI辅助写作的Latex编辑器SmartLatexEditor，超级好用，推荐试试
已经有11人回复
-
物理学I论文润色/翻译怎么收费?
已经有94人回复
-
2026-CJ开始申报了
已经有1人回复
-
西安电子科学大学杭州研究院刘丽香教授招收智能多模态传感器和微型储能器件方向博士
已经有13人回复
-
重庆交通大学光子学微结构与器件课题组招收2026年硕士研究生信息
已经有0人回复
-
一志愿郑大材料学硕298分，求调剂
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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
30