[size=13.3333px]已经建立了如下三个方程组,想解未知数:
syms a b c g12 g23 g21 g32 g13 t12 t23 t13 t31 t21 x y z
eq1=sym('a/(x+y*g12+z*g13)+(b*t12^k)/(x*g12+y+z*g23)+c*t12^k/(x*t13+y*g23+z)=(a*t12^k)/(x+y*g12+z*g13)+b/(x*g12+y+z*g23)+(c*t23^k)/(x*g13+y*g23+z)')
eq2=sym('a/(x+y*g12+z*g13)+b*t12^k/(x*g12+y+z*g23)+c*t13^k/(x*g13+y*g23+z)=a*t13^k/(x*g13+y*g23+z)+b*t23^k/(x*g13+y*g23+z)+c/(x*g12+y*g23+z)')
eq3=sym('z=1-x-y')
[x,y,z]=solve(eq1,eq2,eq3,'x','y','z')
[size=13.3333px]但实际运行以后,貌似没有解析解(Explicit solution could not be found[size=13.3333px]),只能数值模拟看变量之间的关系。求助内容:想要的结果有三个:
[size=13.3333px](1)x、y的关系[size=13.3333px](如果无解析解就画出相关关系图)。[size=13.3333px]数值模拟时,赋值可考虑为:a=b=c=0.33 g12=[size=13.3333px]g21=1.2247 g23=g32=g13=2.8284 t12=[size=13.3333px]t21=1.5 t23=t13=t31=2 k=1.5 (注,从经济含义上:x、y、z均大于等于0,小于等于1,a、b、c均大于0小于1,g、t均大于1,k大于等于1小于等于2)
[size=13.3333px](2)x/y、g13的关系。[size=13.3333px]数值模拟时,赋值可考虑为:a=b=c=0.33 g12=[size=13.3333px]g21=1.2247 t12=[size=13.3333px]t21=1.5 t23=t13=t31=2 k=1.5 g23=g32=g13[size=13.3333px](注,从经济含义上:x、y、z均大于等于0,小于等于1,a、b、c均大于0小于1,g、t均大于1,k大于等于1小于等于2)
[size=13.3333px](3)x与a的关系。[size=13.3333px]数值模拟时,赋值可考虑为:b=c=a/2 g12=[size=13.3333px]g21=1.2247 t12=[size=13.3333px]t21=1.5 t23=t13=t31=2 k=1.5 g23=g32=g13=2.8284[size=13.3333px](注,从经济含义上:x、y、z均大于等于0,小于等于1,a、b、c均大于0小于1,g、t均大于1,k大于等于1小于等于2)
[size=13.3333px]由于matlab不熟悉,也没有进行过变量数值模拟,如果上述各种情形赋值如果不合理可调整。需要请提供实现上述结果的matlab命令源代码。感激不尽!!
[size=13.3333px]