solve求方程组的问题 。
显示 更早的评论
syms C1 C2 C3 x y;
c1=[2.0810;3.6763];c2=[1.7604;4.9588];
xc=2;%yc=4;
yc=(xc-c1(1))/(c2(1)-c1(1))*(c2(2)-c1(2))+c1(2);
f=C1*x^2+C2*x*y+C3*y^2-(2*C1*xc+C2*yc)*x-(2*C3*yc+C2*xc)*y+(C1*xc^2+C2*xc*yc+C3*yc^2-1);
y1=7.3080*x;
t1=subs(f,y,y1);
d1=coeffs(collect(vpa(t1)),x);
deta1=d1(2)^2-4*d1(3)*d1(1)
y2=0.9744*x;
t2=subs(f,y,y2);
d2=coeffs(collect(vpa(t2)),x);
deta2=d2(2)^2-4*d2(3)*d2(1)
y3=-0.6943*x+4.1660;
t3=subs(f,y,y3);
d3=coeffs(collect(vpa(t3)),x);
deta3=d3(2)^2-4*d3(3)*d3(1)
[C1,C2,C3]=solve(deta1,deta2,deta3)
为什么这样求出来的解是错误的，而如果把yc＝4，就能求出正解，但是yc=(xc-c1(1))/(c2(1)-c1(1))*(c2(2)-c1(2))+c1(2)也是等于4的呀。
采纳的回答
votatev
2022-11-25
0 个投票
计算机是有舍入误差的噢，鉴于的计算对数值是敏感的。如果计算方式都是一模一样的，并排除了程序存在问题的可能性，那么，就可以考虑存在由于舍入误差引入计算结果差异的可能性。
更多回答（0 个）
类别
在 帮助中心 和 File Exchange 中查找有关 Symbolic Math Toolbox 的更多信息
