solve解方程的问题

yuanping09

我用solve函数解方程为什么返回值不能给出小数表示的结果呢？
如：给定a=3,b=4,c=1,d=1,
clear; syms x; syms y; syms z; syms w;
a=input('please input a= '); b=input('please input b= ');
c=input('please input c= '); d=input('please input d= ');
eq1=(x+y+z)*w/(x+y+z+w)-a; eq2=(x+y+w)*z/(x+y+z+w)-b;
eq3=(x+w+z)*y/(x+y+z+w)-c; eq4=(w+y+z)*x/(x+y+z+w)-d;
[solvx,solvy,solvz,solvw] = solve(eq1,eq2,eq3,eq4,x,y,z,w)
这样一个简单的程序解得的结果是这样表示的：
solvx =

11/4+1/4*17^(1/2)
11/4-1/4*17^(1/2)
              -1/2
              -1/2
怎么转换成一个小数来表示最终的结果啊，请大家帮帮忙啊？

[ 本帖最后由 ChaChing 于 2009-10-30 21:22 编辑 ]

ChaChing

help double

yuanping09

问题解决了，谢谢版主！这个论坛真是太好了！继续学习！

yuanping09

大家好，还是这个程序，在这里还想请教两个问题：
1.实际中只要x,y,z,w都为正数的那组解，怎么在显示的结果中只保留x,y,z,w同时都为正数的那组解啊？
   因为要用这组解来计算e=(x+z)*(y+w)/(x+y+z+w); f=(x+w)*(y+z)/(x+y+z+w);
2.而且直接计算e,f会提示错误??? Error using ==> mtimes
Inner matrix dimensions must agree.
Error in ==> abcd1 at 19
e=(x+z)*(y+w)/(x+y+z+w);
应该怎么修改啊？
谢谢！

[ 本帖最后由 ChaChing 于 2009-10-31 12:51 编辑 ]

ChaChing

a=3;b=4;c=1;d=1;
syms x y z w;
eq1=(x+y+z)*w/(x+y+z+w)-a; eq2=(x+y+w)*z/(x+y+z+w)-b;
eq3=(x+w+z)*y/(x+y+z+w)-c; eq4=(w+y+z)*x/(x+y+z+w)-d;
[solvx,solvy,solvz,solvw] = solve(eq1,eq2,eq3,eq4,x,y,z,w);
x1=double(solvx); y1=double(solvy); z1=double(solvz); w1=double(solvw);
Sol=[x1,y1,z1,w1]; ii=all(Sol>0,2);
x1=x1(ii); y1=y1(ii); z1=z1(ii); w1=w1(ii);
e=(x1+z1)*(y1+w1)/(x1+y1+z1+w1)
f=(x1+w1)*(y1+z1)/(x1+y1+z1+w1)

[ 本帖最后由 ChaChing 于 2009-10-31 12:52 编辑 ]

yuanping09

ChaChing真的是很厉害啊，谢谢你！

ChaChing

客气! 个人并不厉害, 仅是分享我会的!
亦希望楼主能与会员们一起分享

zhaozhao0726

俺也遇到了这个问题 谢谢chaching