Question

怎么用matlab或C语言编写程序求解下列方程组 急求啊 谢谢

((xn+sn/2)^2+(yn-50)^2)^0.5=60-sn
(yn-50)/(xn-sn/2)=25/(sn/2-(((160-s1)^2-50^2)^0.5+s1)/2
其中s1到sn为已知量，n取1—10 求 xn yn
最好用for循环

Answer 1

clear
syms x y
s(10)=0; % 已知量

xs(2,10)=0;
ys(2,10)=0;
for n=1:10
[sx,sy]=solve( ((x+s(n)/2)^2+(y-50)^2)^0.5==60-s(n),(y-50)/(x-s(n)/2)==25/(s(n)/2-(((160-s(1))^2-50^2)^0.5+s(1)))/2 );

xs(:,n)=sx;ys(:,n)=sy;
end
xs,ys

追问

??? Error using ==> horzcat
Conversion to char from logical is not possible.
Error in ==> C:\MATLAB6p5p1\toolbox\symbolic\solve.m
On line 77 ==> eqns = [eqns ',' v];
您能把算出来的答案直接发给我吗？谢谢

追答

安装个matlab 2013，matlab 2014. 吧

追问

额 可是我现在急需答案啊。

追答

给 s1到sn的值给我吧。

追问

s1=10.9659s2=15.0499s3=17.8956s4=20.0375s5=21.6853s6=22.9456s7=23.8799s8=24.5255s9=24.9048s10=25.03

追答

clear
syms x y
s=[10.9659,15.0499,17.8956,20.0375,21.6853,22.9456,23.8799,24.5255,24.9048,25.03];
xs(2,10)=0;
ys(2,10)=0;
for n=1:10
[sx,sy]=solve( ((x+s(n)/2)^2+(y-50)^2)^0.5==60-s(n),(y-50)/(x-s(n)/2)==25/(s(n)/2-(((160-s(1))^2-50^2)^0.5+s(1)))/2 );

xs(:,n)=sx;ys(:,n)=sy;
end
xs,ys

xs =

Columns 1 through 7

-54.2492 -52.1746 -50.7252 -49.6315 -48.7884 -48.1423 -47.6626
43.4431 37.3504 33.1032 29.9050 27.4436 25.5602 24.1635

Columns 8 through 10

-47.3307 -47.1356 -47.0711
23.1981 22.6308 22.4436

ys =

Columns 1 through 7

55.1183 55.1881 55.2376 55.2753 55.3045 55.3270 55.3437
46.7473 47.4081 47.8798 48.2413 48.5233 48.7412 48.9041

Columns 8 through 10

55.3553 55.3621 55.3644
49.0173 49.0841 49.1062

追问

谢谢。