| 查看: 1372 | 回复: 0 | ||
[求助]
MATLAB求解多元非线性方程组
|
|
我想求解一个非常复杂的方程组，原理很简单，就是解三个方程求出三个未知数a、b、c，但每一个方程都包含exp、log、sqrt函数，不知道能否用MATLAB求出解来，我试着写了程序，但是运行不出来，我第一次用MATLAB，不知道写的对不对，请各位高手帮忙解出该方程组来，我有急用，谢谢了~ syms a b c I=3*a+4*b+3*c; f=-4*0.377/1.2*I*log(1+1.2*sqrt(I)); f1=-2*0.377*(sqrt(I)/(1+1.2*sqrt(I))+2/1.2*log(1+1.2*sqrt(I))); B1=-0.03942+0.97225/2/I*(1-(1+2*sqrt(I))*exp(-2*sqrt(I))); B11=0.97225/2/I/I*(-1+(1+2*sqrt(I)+2*I)*exp(-2*sqrt(I))); C1=0.0179/sqrt(8); B2=0.23825+2.9605/0.98/I*(1-(1+1.4*sqrt(I))*exp(-1.4*sqrt(I)))-30.905/72/I*(1-(1+12*sqrt(I))*exp(-12*sqrt(I))); B21=2.9605/0.98/I/I*(-1+(1+1.4*sqrt(I)+0.98*I)*exp(-1.4*sqrt(I)))-30.905/72/I/I*(-1+(1+12*sqrt(I)+72*I)*exp(-12*sqrt(I))); C2=0.003; B3=0.00488+0.4911/2/I*(1-(1+2*sqrt(I))*exp(-2*sqrt(I))); B31=0.4911/2/I/I*(-1+(1+2*sqrt(I)+2*I)*exp(-2*sqrt(I))); C3=-0.00116/sqrt(8); g=(1-(1+2*sqrt(I))*exp(-2*sqrt(I)))/2/I; g1=(-1+(1+2*sqrt(I)+2*I)*exp(-2*sqrt(I)))/2/I/I; w=1+(I*f1-f+4*a*(a+b+c)*(B1+I*B11+4*(a+b+c)*C1)+2*b*(a+b+c)*(B2+I*B21+4*(a+b+c)*C2)+4*c*(a+b+c)*(B3+I*B31+4*(a+b+c)*C3))/(3*a+2*b+3*c); r1=exp(f1+(4/3*a+4/3*(a+b+c))*(B1+2*(a+b+c)*C1)+4/3*c*(B3+2*(a+b+c)*C3)+2/3*b*(B2+2*(a+b+c)*C2)+2*a*(a+b+c)*(2*B11+4/3*C1)+2*c*(a+b+c)*(2*B31+4/3*C3)+b*(a+b+c)*(2*B21+4/3*C2)+8/3*c*(0.125167+1.078184*g)+4/3*b*(2.521475-15.811274*g)+(4/3*c*(a+b+c)+4/6*a*c)*(-0.0080597)+(2/3*b*(a+b+c)+1/3*a*b)*(-1.016358)+1/3*c*b*0.0672253+4*a*c*1.078184*g1+2*a*b*(-15.811274)*g1+2*c*b*10.163009*g1); r2=exp(2*f1+(b+(a+b+c))*(B2+2*(a+b+c)*C2)+2*a*(B1+2*(a+b+c)*C1)+2*c*(B3+2*(a+b+c)*C3)+2*a*(a+b+c)*(4*B11+2*C1)+b*(a+b+c)*(4*B21+2*C2)+2*c*(a+b+c)*(4*B31+2*C3)+2*a*(2.521475-15.811274*g)+2*c*(-0.0692694+10.163009*g)+(a*(a+b+c)+1/2*a*b)*(-1.016358)+(c*(a+b+c)+1/2*b*c)*0.0672253+a*c*(-0.0080597)+4*a*b*(-15.811274)*g1+4*c*b*10.163009*g1+8*a*c*1.078184*g1); r3=exp(f1+(4/3*c+4/3*(a+b+c))*(B3+2*(a+b+c)*C3)+4/3*a*(B1+2*(a+b+c)*C1)+2/3*b*(B2+2*(a+b+c)*C2)+2*a*(a+b+c)*(2*B11+4/3*C1)+2*c*(a+b+c)*(2*B31+4/3*C3)+b*(a+b+c)*(2*B21+4/3*C2)+8/3*a*(0.125167+1.078184*g)+4/3*b*(-0.0692694+10.163009*g)+(4/3*a*(a+b+c)+4/6*a*c)*(-0.0080597)+(2/3*b*(a+b+c)+1/3*c*b)*0.0672253+1/3*a*b*(-1.016358)+4*a*c*1.078184*g1+2*a*b*(-15.811274)*g1+2*c*b*10.163009*g1); eq1=b*r2*(exp(-(3*a+2*b+3*c)*w/55.51))^7-0.076849; eq2=b*c*r2*r3*(exp(-(3*a+2*b+3*c)*w/55.51))^6-0.012198; eq3=a*r1*(exp(-(3*a+2*b+3*c)*w/55.51))^10-0.092085; [a,b,c]=solve(eq1,eq2,eq3,'a,b,c') |
回复此楼» 猜你喜欢
316求调剂
已经有6人回复
-
346求调剂[0856]
已经有7人回复
-
一志愿山大07化学 332分 四六级已过 本科山东双非 求调剂！
已经有3人回复
-
310求调剂
已经有3人回复
-
机械专硕299求调剂至材料
已经有4人回复
-
070300化学319求调剂
已经有7人回复
-
08工科 320总分 求调剂
已经有6人回复
-
一志愿天津大学化学工艺专业（081702）315分求调剂
已经有12人回复
-
307求调剂
已经有3人回复
-
265求调剂
已经有3人回复
找到一些相关的精华帖子，希望有用哦~
matlab解非线性方程组解不开
已经有3人回复
- matlab求解非线性方程组
已经有3人回复
- matlab求解非线性方程组，错误提示怎么解决
已经有5人回复
- 求助 MATLAB解方程组-fslove
已经有7人回复
- 请教用matlab求解一个非线性偏微分方程组的数值解
已经有8人回复
- 求助matlab---fsolve解非线性方程组
已经有6人回复
- matlab非线性方程组求解，急！！！
已经有5人回复
- 用不动点迭代求非线性方程组，求得的不是想要的解，怎么办
已经有12人回复
- 初学者，请教大家关于matlab解多元二次方程组
已经有4人回复
- matlab求解非线性方程组，求助！
已经有6人回复
- 如何使用matlab求解非线性方程组的所有整数解？
已经有9人回复
- Matlab求解偏微分方程组
已经有22人回复
- matlab求解非线性方程组
已经有16人回复
- matlab数值求解边界条件微分方程组
已经有7人回复
- 求matlab求解一个方程组的问题
已经有6人回复
- 求高人指点用matlab求解非线性方程组，解决了追加100金币；
已经有11人回复
- 请教Matlab高手解多元非线性方程组中的参数
已经有6人回复
- matlab的fsove 命令求解非线性方程组
已经有6人回复
- MATLAB非线性方程组求解~~~反应釜模拟
已经有9人回复
- 【求助】多元方程组求解及拟合
已经有9人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 非线性方程组的迭代法（数值计算高手请进）
已经有7人回复
- 【求助】求教matlab解非线性方程组
已经有9人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
