charlotte.w 发表于 2007-4-3 15:39

验证方程的解怎么回代呢

如题, 要验证方程的解怎么回代呢? 谢谢

[ 本帖最后由 eight 于 2007-4-3 15:44 编辑 ]

xjzuo 发表于 2007-4-3 16:11

定义为函数后,代入即可.
如还不清楚,请贴出方程.

charlotte.w 发表于 2007-4-3 16:20

解关于B的方程f=[((20.*sin(A)-20.*cot(B).*(1-cos(A))-L.*sin(e+B)+tan(B).*(L*cos(e+B)-950))./(cot(B)+tan(B))+20.*(1-cos(A))).^2+(-20.*cot(B).*(1-cos(A))-cot(B).*((20.*sin(A)-20.*cot(B).*(1-cos(A))-L.*sin(e+B)+tan(B).*(L.*cos(e+B)-950))./(cot(B)+tan(B)))).^2-70^2]

得到的f=   0.00197951734026,不为零。 这样得到的B值还可信么

要让f为0 要怎么做呢

charlotte.w 发表于 2007-4-3 16:24

我就是要解这个方程, 但是得到的B值和真实差太多. 我怀疑应该是因为解方程时函数值不为零造成的?

跟B的初值有关?

[ 本帖最后由 charlotte.w 于 2007-4-3 16:26 编辑 ]

gerry 发表于 2007-4-3 17:32

数值解法很少有完全等于零,达到你需要的精度就可以了,不知道你是用什么方法解的,但Matlab内置的函数通常可以在输入时指定精度

charlotte.w 发表于 2007-4-4 10:25

gerry: 我用的fsolve
options = optimset('Display','iter')
fsolve(@myfun,B0,options,AA)
我要解的精度到10e-5, 要怎么设置呢?谢谢
页: [1]
查看完整版本: 验证方程的解怎么回代呢