【自我挑戰63】

Joseflin

求a b 值：   ( 註：两圆相切)

WJQ6080

15.21

19.08

对不对

Joseflin

请参考：

asong83

求解？？？？？？？？？？？？？？？？？？？？？？？

Joseflin

请高手出招。

qjchen

不懂几何求解，只能先给个结果，是利用python编程的，不过还是要寻求几何解，不过似乎很难，还请高手指教

python代码（activepython+scipy）

2. import os
3. import math
4. from scipy import *
5. from scipy.optimize import fsolve
6. def dis(p0,p1):
7.     p0x=p0[0]
8.     p0y=p0[1]
9.     p1x=p1[0]
10.     p1y=p1[1]
11.     return math.sqrt((p0x-p1x)**2+(p0y-p1y)**2)
12. def p2l(p1,p2,p3):
13.     return [dis(p1,p2),dis(p2,p3),dis(p3,p1)]
14. def in_rb(p1,p2,p3):
15.     l=p2l(p1,p2,p3)
16.     a=l[0]
17.     b=l[1]
18.     c=l[2]
19.     p=0.5*(a+b+c)
20.     return [math.sqrt((p-a)*(p-b)*(p-c)/p),p-c]
21. def func(x):
22.     point1=[0,0]
23.     point2=[x[0],0]
24.     point3=[x[1],45]
25.     point11=[100,0]
26.     res1=in_rb(point1,point2,point3)
27.     res2=in_rb(point11,point2,point3)
28.     out=[res1[0]-2*res2[0]]
29.     out.append(res1[1]-res2[1])
30.     return out
31. x02 = fsolve(func, [50, 50])
32. print x02

复制代码