Labfans是一个针对大学生、工程师和科研工作者的技术社区。 论坛首页 | 联系我们(Contact Us)
MATLAB爱好者论坛-LabFans.com
返回   MATLAB爱好者论坛-LabFans.com > 基础科学 > 物理学
物理学 A discussion board for Physics.
回复
 
主题工具 显示模式
旧 2009-05-14, 19:07   #1
kaka768
初级会员
 
注册日期: 2009-01-23
年龄: 35
帖子: 2
声望力: 0
kaka768 正向着好的方向发展
问题 求助!matlab解符号方程出错!

代码如下:
syms Ry Rb Rc Ra Rd Rf X E Iy Ib
>> [Ry Iy Ib]=solve('Iy*Ry+X*Rf=Ib*Rb','(Iy-X)*Rc=X*Rf+(X+Ib)*Ra','(Iy+Ib)*Rd+Iy*Ry+(Iy-X)*Rc=E','Ry','Iy','Ib')

Ry =

(-Rd*Rc*X+Rc*E-Rd*X*Rf-Rd*Ra*X-Rc*X*Ra)/(Rb*Rc+Rd*Ra+Rd*Rc+Rc*Ra)



Iy =

(X*Rf*Ra+Rb*Rc*X+Rb*X*Rf+Rb*Ra*X+Rd*Rc*X+Rd*X*Rf+Rd*Ra*X+Rc*X*Ra+E*Ra)/(Rb*Rc+Rd*Ra+Rd*Rc+Rc*Ra)



Ib =

(-Rb*X*Rf*Rd-X*Rf*Ra*Rd-X*Rf*Ra*Rc+Rb*Rc*E-Rb*Ra*X*Rd-Rb*Ra*X*Rc-X*Rf*Rb*Rc-X*Rf*Rd*Rc-Rb*Rc*Rd*X)/(X*Rf*Ra+Rb*Rc*X+Rb*X*Rf+Rb*Ra*X+Rd*Rc*X+Rd*X*Rf+Rd*Ra*X+Rc*X*Ra+E*Ra)
其中I,R表示电流、电阻,X也为电流,E为电势!

电阻Ry =(-Rd*Rc*X+Rc*E-Rd*X*Rf-Rd*Ra*X-Rc*X*Ra)/(Rb*Rc+Rd*Ra+Rd*Rc+Rc*Ra)
其量纲为电流!而电流Ib =(-Rb*X*Rf*Rd-X*Rf*Ra*Rd-X*Rf*Ra*Rc+Rb*Rc*E-Rb*Ra*X*Rd-Rb*Ra*X*Rc-X*Rf*Rb*Rc-X*Rf*Rd*Rc-Rb*Rc*Rd*X)/(X*Rf*Ra+Rb*Rc*X+Rb*X*Rf+Rb*Ra*X+Rd*Rc*X+Rd*X*Rf+Rd*Ra*X+Rc*X*Ra+E*Ra)
其量纲为电阻!!!!
这是为什么?哪里出错了!高手指点!!!
kaka768 当前离线   回复时引用此帖
旧 2010-07-18, 22:45   #2
yisanming
初级会员
 
yisanming 的头像
 
注册日期: 2010-07-18
年龄: 37
帖子: 10
声望力: 15
yisanming 正向着好的方向发展
默认 回复: 求助!matlab解符号方程出错!

主要是matlab中的solve函数求解出结果后会按照变量的字母顺序列出,最后相当于对Iy、Ib、Ry进行了一次重新赋值。[Ry Iy Ib]=solve('Iy*Ry+X*Rf=Ib*Rb','(Iy-X)*Rc=X*Rf+(X+Ib)*Ra','(Iy+Ib)*Rd+Iy*Ry+(Iy-X)*Rc=E','Ry','Iy','Ib'),解出来之后Ry=ans(Ib)、Iy=ans(Iy)、Ib=ans(Ry)。
yisanming 当前离线   回复时引用此帖
回复

主题工具
显示模式

发帖规则
不可以发表新主题
不可以发表回复
不可以上传附件
不可以编辑自己的帖子

启用 BB 代码
论坛启用 表情符号
论坛启用 [IMG] 代码
论坛禁用 HTML 代码



所有时间均为北京时间。现在的时间是 05:17


Powered by vBulletin
版权所有 ©2000 - 2025,Jelsoft Enterprises Ltd.