Python计算离轴三反系统的初始结构

输入面间距d1,d2,d3和焦距f的值，可计算出遮
拦系数和放大系数，进而代入半径公式得出半径

http://www.optzmx.com/thread-23103-1-1.html

import sympy as sp


d1 = -1200

d2 = 1200

d3 = -1300

f = -2000

a1, a2, b1, b2 = sp.symbols('a1,a2,b1,b2')
eq = [
    sp.Eq(f * (1 - a1) / b1 / b2, d1),
    sp.Eq(f * a1 * (1 - a2) / b2, d2),
    sp.Eq(a1 * a2 * f, d3),
    sp.Eq(b1*b2-b2*(1+b1)/a1+(1+b2)/a1/a2, 0)
]

solution = sp.solve(eq, (a1, a2, b1, b2))

print(solution)

r1 = 2*f/b1*b2
r2 = 2*a1*f/b2/(1+b1)
r3 = 2*a1*a2*f/(1+b1)