求高手解答用Matlab 解复杂方程组的解
syms x y z a b c
[x,y,z]=solve('72.8*(1+cos(a))=2*(sqrt(x*21.8)+sqrt(25.5*y)+sqrt(25.5*z))','50.8*(1+cos(b))=2*sqrt(x*50.8)','48*(1+cos(c))=2*(sqrt(x*29)+sqrt(47*y)+sqrt(1.92*z))','x,y,z')
x =
12.700000000000000000000000000000+25.400000000000000000000000000000*cos(b)+12.700000000000000000000000000000*cos(b)^2
y =
-.98416931170728906540241988580218*cos(c)+9.7643439805890752790240685157559*cos(a)*cos(b)-23.462398366292790782350185322144*cos(c)*cos(b)+7.1488996766620738069752563394849*cos(b)^2+3.3341639162497969191291895638814*cos(a)^2+.40958164398297044880227538559881*cos(a)+.59974496762043211062439587258181*cos(b)-16.023090169105698668480172257920*cos(c)*cos(a)+19.250659611845600192713968847131*cos(c)^2+.12578649948056747012125686189223e-1
z =
.15194565438596826916006412565125e-193*(.58477278352492431981436322974990e98+.58477278352492431981436322974990e98*cos(a)-.26731041191058757560106827855831e98*((cos(b)+1.)^2)^(1/2)-1.*(-.64771139107889760381210034820223e194*cos(c)+.64262081202967154457022801098613e195*cos(a)*cos(b)-.15441309237244940736312726425476e196*cos(c)*cos(b)+.47049056490306928474554525841306e195*cos(b)^2+.21943134403702282352778341593461e195*cos(a)^2+.82783874266682206962299704487537e192+.26955798481907364090410143645319e194*cos(a)+.39471018111316040430045885264655e194*cos(b)-.10545277016218098275353846992100e196*cos(c)*cos(a)+.12669437431192070703927236534687e196*cos(c)^2)^(1/2))^2
______________________________
就是代入求值就行了，例如：（接着上面的求X）
x=subs(x,{a,b,c},{77.7,33.3,48.5})
x =
127/10+127/5*cos(333/10)+127/10*cos(333/10)^2
__________________________________
其实在一开始就给定a,b,c的值也是可以的，如果不要求一定要求解析解的话，楼下说的fsolve倒也是一种选择
用matlab解复杂方程组
function hhh
To=400
Tf=0.75*To+0.25*302
V=0.6*4.6*0.3/sqrt(1+1)*sqrt(2*9.8*2.5*(Tf-302)/302)
Nu=4.9+0.0606*(V/4.6/0.04*0.04/0.000016*0.701*2*0.04/2.5)^1.2/(1+0.0909*(V/4.6/0.04*0.04/0.000016*0.701*2*0.04/2.5)^0.7*0.701^0.17)
h=0.0267*Nu/2.5
q=1.128*V*1005*(Tf-302)/0.75/4.6/2.2
eq1=sym('(Tg-302)/(1/23.26+0.01/1.1)+h*(Tg-Tf)+hwg*(Tg-Tw)=0.7*600');
eq2=sym('h*(Tg-Tf)+h*(Tw-Tf)=q');
eq5=sym('h*(Tw-Tf)+1/(1/8.72+0.2/0.49)*(Tw-299)+hwg*(Tw-Tg)=0');
%eq7=sym('5.67*10^8*(Tg^2+Tw^2)*(Tw+Tg)/(1/0.8+1/0.94-1)=hwg');
[Tg,Tw,hwg]=solve(eq1,eq2,eq5,'Tg,Tw,hwg');
Tg=subs(Tg),Tw=subs(Tw),hwg=subs(hwg)
hwg1=5.67*10^8*(Tg^2+Tw^2)*(Tw+Tg)/(1/0.8+1/0.94-1)
x0=[Tg h Tf hwg Tw q To Nu V]
%x0=[11.2698 0.9291 455.5052 -1.8440 -11.2698 113.7074 487.3698 97.0550 0.0031]
fsolve(@myfuns,x0)
function eqs=myfuns(x)
Tg=x(1);h=x(2);Tf=x(3);hwg=x(4);Tw=x(5);q=x(6);To=x(7);Nu=x(8);V=x(9);
eq1=(Tg-302)/(1/23.26+0.01/1.1)+h*(Tg-Tf)+hwg*(Tg-Tw)-0.7*600;
eq2=h*(Tg-Tf)+h*(Tw-Tf)-q;
eq3=1.128*V*1005*(Tf-302)/0.75/4.6/2.2-q;
eq4=0.75*To+0.25*302-Tf;
eq5=h*(Tw-Tf)+1/(1/8.72+0.2/0.49)*(Tw-299)+hwg*(Tw-Tg);
eq6=0.6*4.6*0.3/sqrt(1+1)*sqrt(2*9.8*2.5*(Tf-302)/302)-V;
eq7=5.67*10^8*(Tg^2+Tw^2)*(Tw+Tg)/(1/0.8+1/0.94-1)-hwg;
eq8=0.0267*Nu/2.5-h;
eq9=4.9+0.0606*(V/4.6/0.04*0.04/0.000016*0.701*2*0.04/2.5)^1.2/(1+0.0909*(V/4.6/0.04*0.04/0.000016*0.701*2*0.04/2.5)^0.7*0.701^0.17)-Nu;
eqs=[eq1;eq2;eq3;eq4;eq5;eq6;eq7;eq8;eq9];
结果：
Optimizer appears to be converging to a point which is not a root.
Norm of relative change in X is less than max(options.TolX^2,eps) but
sum-of-squares of function values is greater than or equal to sqrt(options.TolFun)
Try again with a new starting guess.
ans =
1.0e+004 *
-3.6696 0.0000 0.0376 -0.0002 3.6696 2.2196 0.0400 0.0020 0.0002
求一个win7 能用的maple 8软件 的下载
首先，请尽量卸载之前安装的office2010版本，使用下面工具卸载(卸载更干净)：
并且，卸载家庭版和学生版（阉割版）和试用版，然后用以下方案。
以下为几个OfFice个人解决方案，都是最合理的，可以自由选择：
方案一 （安装office 2010免激活版）
Office Professional Plus 2010免激活版（包含32位和64位）：
由VOL版本封装，集成了“MAK”密钥，跟方案三很相似。只是做了封装。
方案二 （安装Office 2010企业版+KMS激活——推荐）：
Microsoft Office 2010 Pro VOL简体中文正式版+迷你KMS激活:
KMS 激活方法：
方案三（安装Office 2010专业版+密钥激活）：
Microsoft Office Pro 2010 简体中文正式版下载：
最新Microsoft Office Pro 2010 简体中文正式版激活密钥(密钥激活过程)：
目前网上大多数密钥已封，这是最新的密钥。看运气吧，有的人能直接激活，有的人要电话激活，有的则不行。
附教程：
《Microsoft Office 2010实用技巧宝典》微软官方光盘下载：
因为密钥的激活次数有限，个人推荐方案二（完全按方案二），希望对你有帮助
