有关MATLAB求解三元二次方程组的问题?
有关MATLAB求解三元二次方程组的问题?
>> syms x y z
>> eq1=(18*(2*x - 30))/(5*((y - z)^2 + (x - 15)^2)^(1/2)) + (36*x)/(5*((y - 5)^2 + x^2)^(1/2));eq2=(18*(2*y - 10))/(5*((y - 5)^2 + x^2)^(1/2)) + (18*(2*y - 2*z))/(5*((y - z)^2 + (x - 15)^2)^(1/2)) + 36/5;eq3=(141*(2*z - 16))/(10*((z - 8)^2 + 25)^(1/2)) - (18*(2*y - 2*z))/(5*((y - z)^2 + (x - 15)^2)^(1/2));[x,y,z]=solve(eq1,eq2,eq3,'x','y','z')
请求大神看一下哪里有错误,我求不出结果
>> syms x y z
>> eq1=(18*(2*x - 30))/(5*((y - z)^2 + (x - 15)^2)^(1/2)) + (36*x)/(5*((y - 5)^2 + x^2)^(1/2));eq2=(18*(2*y - 10))/(5*((y - 5)^2 + x^2)^(1/2)) + (18*(2*y - 2*z))/(5*((y - z)^2 + (x - 15)^2)^(1/2)) + 36/5;eq3=(141*(2*z - 16))/(10*((z - 8)^2 + 25)^(1/2)) - (18*(2*y - 2*z))/(5*((y - z)^2 + (x - 15)^2)^(1/2));[x,y,z]=solve(eq1,eq2,eq3,'x','y','z')
请求大神看一下哪里有错误,我求不出结果
赞 爱恨一瞬间 幼苗
共回答了24个问题采纳率：87.5% 举报
syms x y z
>> eq1=(18*(2*x - 30))/(5*((y - z)^2 + (x - 15)^2)^(1/2)) + (36*x)/(5*((y - 5)^2 + x^2)^(1/2));
>> eq2=(18*(2*y - 10))/(5*((y - 5)^2 + x^2)^(1/2)) + (18*(2*y - 2*z))/(5*((y - z)^2 + (x - 15)^2)^(1/2)) + 36/5;
>> eq3=(141*(2*z - 16))/(10*((z - 8)^2 + 25)^(1/2)) - (18*(2*y - 2*z))/(5*((y - z)^2 + (x - 15)^2)^(1/2));
>> [x,y,z]=solve(eq1,eq2,eq3,'x','y','z')
经过在matlab7.0中尝试,发现没有问题,
x =
1861137696/586981625+2173/1080500*(-3+5*3^(1/2)+30/2173*2173^(1/2))^3-26751/54025*3^(1/2)-160506/117396325*2173^(1/2)+36791/540250*(-3+5*3^(1/2)+30/2173*2173^(1/2))^2
y =
13/2-5/2*3^(1/2)-15/2173*2173^(1/2)
z =
-6519/108050*(-3+5*3^(1/2)+30/2173*2173^(1/2))^2+451757/54025-2173/324150*(-3+5*3^(1/2)+30/2173*2173^(1/2))^3+17834/10805*3^(1/2)+107004/23479265*2173^(1/2)
这里得到的是符号解,转化为数值解；
>> x=subs(x),y=subs(y),z=subs(z)
x =
5.4593
y =
1.8481
z =
7.3564
>> eq1=(18*(2*x - 30))/(5*((y - z)^2 + (x - 15)^2)^(1/2)) + (36*x)/(5*((y - 5)^2 + x^2)^(1/2));
>> eq2=(18*(2*y - 10))/(5*((y - 5)^2 + x^2)^(1/2)) + (18*(2*y - 2*z))/(5*((y - z)^2 + (x - 15)^2)^(1/2)) + 36/5;
>> eq3=(141*(2*z - 16))/(10*((z - 8)^2 + 25)^(1/2)) - (18*(2*y - 2*z))/(5*((y - z)^2 + (x - 15)^2)^(1/2));
>> [x,y,z]=solve(eq1,eq2,eq3,'x','y','z')
经过在matlab7.0中尝试,发现没有问题,
x =
1861137696/586981625+2173/1080500*(-3+5*3^(1/2)+30/2173*2173^(1/2))^3-26751/54025*3^(1/2)-160506/117396325*2173^(1/2)+36791/540250*(-3+5*3^(1/2)+30/2173*2173^(1/2))^2
y =
13/2-5/2*3^(1/2)-15/2173*2173^(1/2)
z =
-6519/108050*(-3+5*3^(1/2)+30/2173*2173^(1/2))^2+451757/54025-2173/324150*(-3+5*3^(1/2)+30/2173*2173^(1/2))^3+17834/10805*3^(1/2)+107004/23479265*2173^(1/2)
这里得到的是符号解,转化为数值解；
>> x=subs(x),y=subs(y),z=subs(z)
x =
5.4593
y =
1.8481
z =
7.3564
1年前 追问
7
可是我这边显示BUSY 就是不出答案 怎么办？
可能相似的问题
-
1年前2个回答
-
求这个三元二次方程组解法a+b=-4ab+c=2ac=1лл
1年前2个回答
-
1年前1个回答
-
1年前2个回答
你能帮帮他们吗