声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1712|回复: 7

[综合讨论] solve解方程的问题

[复制链接]
发表于 2009-10-29 18:53 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
我用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 编辑 ]
回复
分享到:

使用道具 举报

发表于 2009-10-29 21:39 | 显示全部楼层
help double
 楼主| 发表于 2009-10-30 18:06 | 显示全部楼层

回复 沙发 ChaChing 的帖子

问题解决了,谢谢版主!这个论坛真是太好了!继续学习!
 楼主| 发表于 2009-10-31 12:28 | 显示全部楼层

回复 楼主 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 编辑 ]
发表于 2009-10-31 12:49 | 显示全部楼层

回复 地板 yuanping09 的帖子

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 编辑 ]
 楼主| 发表于 2009-10-31 14:10 | 显示全部楼层

回复 5楼 ChaChing 的帖子

ChaChing真的是很厉害啊,谢谢你!
发表于 2009-10-31 14:33 | 显示全部楼层

回复 6楼 yuanping09 的帖子

客气! 个人并不厉害, 仅是分享我会的!
亦希望楼主能与会员们一起分享
发表于 2009-11-9 20:37 | 显示全部楼层
俺也遇到了这个问题 谢谢chaching
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-26 07:42 , Processed in 0.091396 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表