全部版块 我的主页
论坛 提问 悬赏 求职 新闻 读书 功能一区 经管百科 爱问频道
2035 0
2016-04-01
悬赏 200 个论坛币 未解决

[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]


二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群