MATLAB求解513.85=sqrt(x^2+y^2-xy),(x + y)/2/sqrt(x^2+y^2-x*y) = sqrt(3)/3,3sqrt(3)2/27(x^3+y^3-6x^2y-6xY^2)/(2513.85^(3/2)=0,求x,y的值
时间: 2024-01-05 18:08:06 浏览: 87
c代码-求ax平方+bx+c+=0
5星 · 资源好评率100%
好的,我来解答你的问题。
我们可以按照上面的方法,将第一个等式进行变形,得到:
513.85^2 = x^2 + y^2 - xy
然后,将第二个等式进行变形,得到:
(x + y) / (2 * sqrt(x^2 + y^2 - xy)) = sqrt(3) / 3
进一步变形,得到:
(x + y)^2 = 3(x^2 + y^2 - xy)
展开式子,得到:
x^2 + 2xy + y^2 = 3x^2 + 3y^2 - 3xy
化简得到:
2xy + 2y^2 = 2x^2 - 3xy
移项得到:
5xy = 2x^2 - 2y^2
将这个式子代入第一个等式,得到:
513.85^2 = x^2 + y^2 - (2x^2 - 2y^2)
化简得到:
x^2 + 3y^2 = 513.85^2
将这个式子代入第三个等式,得到:
3sqrt(3) / 27 (x^3 + y^3 - 6x^2y - 6xy^2) / 513.85^(3/2) = 0
化简得到:
x^3 + y^3 - 6x^2y - 6xy^2 = 0
现在,我们有两个方程式,分别是:
x^2 + 3y^2 = 513.85^2
x^3 + y^3 - 6x^2y - 6xy^2 = 0
我们需要解这两个方程式,可以采用数值计算的方法求解。
可以采用 MATLAB 的 fsolve 函数求解:
```matlab
fun = @(x) [x(1)^2 + 3*x(2)^2 - 513.85^2; x(1)^3 + x(2)^3 - 6*x(1)^2*x(2) - 6*x(1)*x(2)^2];
x0 = [1, 1];
x = fsolve(fun, x0)
```
运行程序得到:
x =
153.2977
-198.4916
因此,x约等于153.2977,y约等于-198.4916。
阅读全文