maxiaoweitj 发表于 2009-5-29 14:52

大侠们 求助啊 我的毕业设计遇到的问题！！！

=solve('a*(x0-x1)+b*(y0-y1)+c*(z0-z1)=0','(x0-x1)^2+(y0-y1)^2+(z0-z1)^2=(x0-x2)^2+(y0-y2)^2+(z0-z2)^2','(x0-x1)^2+(y0-y1)^2+(z0-z1)^2=(x0-x3)^2+(y0-y3)^2+(z0-z3)^2')
这是求空间一圆柱的截面上圆的中心中标
a,b,c是已经求出的法向量加上所测得的圆上3个点的坐标
三个方程三个未知数应该能求出来的呀
为什么matlab 算出的结果乱七八糟的、、、、、、、...
x0 =

[ (a*x3+(-2*b^2*y0*y1-2*c^2*z0*z1+b^2*y0^2+b^2*y1^2+c^2*z0^2+c^2*z1^2+2*b*y0*c*z0-2*b*y0*c*z1-2*b*y1*c*z0+2*b*y1*c*z1-2*y0*y1*a^2+y1^2*a^2-2*z0*z1*a^2+z1^2*a^2+2*y0*y3*a^2-y3^2*a^2+2*z0*z3*a^2-z3^2*a^2)^(1/2))/a]
[ (a*x3+(-2*b^2*y0*y1-2*c^2*z0*z1+b^2*y0^2+b^2*y1^2+c^2*z0^2+c^2*z1^2+2*b*y0*c*z0-2*b*y0*c*z1-2*b*y1*c*z0+2*b*y1*c*z1-2*y0*y1*a^2+y1^2*a^2-2*z0*z1*a^2+z1^2*a^2+2*y0*y3*a^2-y3^2*a^2+2*z0*z3*a^2-z3^2*a^2)^(1/2))/a]
[ (a*x3-(-2*b^2*y0*y1-2*c^2*z0*z1+b^2*y0^2+b^2*y1^2+c^2*z0^2+c^2*z1^2+2*b*y0*c*z0-2*b*y0*c*z1-2*b*y1*c*z0+2*b*y1*c*z1-2*y0*y1*a^2+y1^2*a^2-2*z0*z1*a^2+z1^2*a^2+2*y0*y3*a^2-y3^2*a^2+2*z0*z3*a^2-z3^2*a^2)^(1/2))/a]
[ (a*x3-(-2*b^2*y0*y1-2*c^2*z0*z1+b^2*y0^2+b^2*y1^2+c^2*z0^2+c^2*z1^2+2*b*y0*c*z0-2*b*y0*c*z1-2*b*y1*c*z0+2*b*y1*c*z1-2*y0*y1*a^2+y1^2*a^2-2*z0*z1*a^2+z1^2*a^2+2*y0*y3*a^2-y3^2*a^2+2*z0*z3*a^2-z3^2*a^2)^(1/2))/a]


y0 =

[ -(-a*x3-(-2*b^2*y0*y1-2*c^2*z0*z1+b^2*y0^2+b^2*y1^2+c^2*z0^2+c^2*z1^2+2*b*y0*c*z0-2*b*y0*c*z1-2*b*y1*c*z0+2*b*y1*c*z1-2*y0*y1*a^2+y1^2*a^2-2*z0*z1*a^2+z1^2*a^2+2*y0*y3*a^2-y3^2*a^2+2*z0*z3*a^2-z3^2*a^2)^(1/2)-b*y0+b*y1-c*z0+c*z1)/a]
[ -(-a*x3-(-2*b^2*y0*y1-2*c^2*z0*z1+b^2*y0^2+b^2*y1^2+c^2*z0^2+c^2*z1^2+2*b*y0*c*z0-2*b*y0*c*z1-2*b*y1*c*z0+2*b*y1*c*z1-2*y0*y1*a^2+y1^2*a^2-2*z0*z1*a^2+z1^2*a^2+2*y0*y3*a^2-y3^2*a^2+2*z0*z3*a^2-z3^2*a^2)^(1/2)-b*y0+b*y1-c*z0+c*z1)/a]
[ -(-a*x3+(-2*b^2*y0*y1-2*c^2*z0*z1+b^2*y0^2+b^2*y1^2+c^2*z0^2+c^2*z1^2+2*b*y0*c*z0-2*b*y0*c*z1-2*b*y1*c*z0+2*b*y1*c*z1-2*y0*y1*a^2+y1^2*a^2-2*z0*z1*a^2+z1^2*a^2+2*y0*y3*a^2-y3^2*a^2+2*z0*z3*a^2-z3^2*a^2)^(1/2)-b*y0+b*y1-c*z0+c*z1)/a]
[ -(-a*x3+(-2*b^2*y0*y1-2*c^2*z0*z1+b^2*y0^2+b^2*y1^2+c^2*z0^2+c^2*z1^2+2*b*y0*c*z0-2*b*y0*c*z1-2*b*y1*c*z0+2*b*y1*c*z1-2*y0*y1*a^2+y1^2*a^2-2*z0*z1*a^2+z1^2*a^2+2*y0*y3*a^2-y3^2*a^2+2*z0*z3*a^2-z3^2*a^2)^(1/2)-b*y0+b*y1-c*z0+c*z1)/a]


z0 =

[ 1/2/a*(2*a*x3+2*(-2*b^2*y0*y1-2*c^2*z0*z1+b^2*y0^2+b^2*y1^2+c^2*z0^2+c^2*z1^2+2*b*y0*c*z0-2*b*y0*c*z1-2*b*y1*c*z0+2*b*y1*c*z1-2*y0*y1*a^2+y1^2*a^2-2*z0*z1*a^2+z1^2*a^2+2*y0*y3*a^2-y3^2*a^2+2*z0*z3*a^2-z3^2*a^2)^(1/2)+2*(-2*y0*y1*a^2+2*z0*z2*a^2-2*b^2*y0*y1-2*z0*z1*a^2+2*y0*y2*a^2-2*c^2*z0*z1+y1^2*a^2+z1^2*a^2-y2^2*a^2-z2^2*a^2+b^2*y0^2+b^2*y1^2+c^2*z0^2+c^2*z1^2+2*b*y0*c*z0-2*b*y0*c*z1-2*b*y1*c*z0+2*b*y1*c*z1)^(1/2))]
[ 1/2/a*(2*a*x3+2*(-2*b^2*y0*y1-2*c^2*z0*z1+b^2*y0^2+b^2*y1^2+c^2*z0^2+c^2*z1^2+2*b*y0*c*z0-2*b*y0*c*z1-2*b*y1*c*z0+2*b*y1*c*z1-2*y0*y1*a^2+y1^2*a^2-2*z0*z1*a^2+z1^2*a^2+2*y0*y3*a^2-y3^2*a^2+2*z0*z3*a^2-z3^2*a^2)^(1/2)-2*(-2*y0*y1*a^2+2*z0*z2*a^2-2*b^2*y0*y1-2*z0*z1*a^2+2*y0*y2*a^2-2*c^2*z0*z1+y1^2*a^2+z1^2*a^2-y2^2*a^2-z2^2*a^2+b^2*y0^2+b^2*y1^2+c^2*z0^2+c^2*z1^2+2*b*y0*c*z0-2*b*y0*c*z1-2*b*y1*c*z0+2*b*y1*c*z1)^(1/2))]
[ 1/2/a*(2*a*x3-2*(-2*b^2*y0*y1-2*c^2*z0*z1+b^2*y0^2+b^2*y1^2+c^2*z0^2+c^2*z1^2+2*b*y0*c*z0-2*b*y0*c*z1-2*b*y1*c*z0+2*b*y1*c*z1-2*y0*y1*a^2+y1^2*a^2-2*z0*z1*a^2+z1^2*a^2+2*y0*y3*a^2-y3^2*a^2+2*z0*z3*a^2-z3^2*a^2)^(1/2)+2*(-2*y0*y1*a^2+2*z0*z2*a^2-2*b^2*y0*y1-2*z0*z1*a^2+2*y0*y2*a^2-2*c^2*z0*z1+y1^2*a^2+z1^2*a^2-y2^2*a^2-z2^2*a^2+b^2*y0^2+b^2*y1^2+c^2*z0^2+c^2*z1^2+2*b*y0*c*z0-2*b*y0*c*z1-2*b*y1*c*z0+2*b*y1*c*z1)^(1/2))]
[ 1/2/a*(2*a*x3-2*(-2*b^2*y0*y1-2*c^2*z0*z1+b^2*y0^2+b^2*y1^2+c^2*z0^2+c^2*z1^2+2*b*y0*c*z0-2*b*y0*c*z1-2*b*y1*c*z0+2*b*y1*c*z1-2*y0*y1*a^2+y1^2*a^2-2*z0*z1*a^2+z1^2*a^2+2*y0*y3*a^2-y3^2*a^2+2*z0*z3*a^2-z3^2*a^2)^(1/2)-2*(-2*y0*y1*a^2+2*z0*z2*a^2-2*b^2*y0*y1-2*z0*z1*a^2+2*y0*y2*a^2-2*c^2*z0*z1+y1^2*a^2+z1^2*a^2-y2^2*a^2-z2^2*a^2+b^2*y0^2+b^2*y1^2+c^2*z0^2+c^2*z1^2+2*b*y0*c*z0-2*b*y0*c*z1-2*b*y1*c*z0+2*b*y1*c*z1)^(1/2))]

>> 。。。
哪位知道哪里出问题了 、、、？？？？？

无水1324 发表于 2009-5-29 16:24

考虑一下数值解法吧，非线性的方程可能没有准确的解表达式

maxiaoweitj 发表于 2009-5-31 10:36

回复 沙发 无水1324 的帖子

LS能否指导一下？？多谢！！

friendchj 发表于 2009-5-31 11:43

help fsolve
页: [1]
查看完整版本: 大侠们 求助啊 我的毕业设计遇到的问题！！！