用波尔共振仪研究受迫振动

#相位差理论值python公式

import math

n=1

while True:

    if n==1:

        rou=eval(input('阻尼系数 rou='))

    w0=eval(input(f"{n}. w0="))

    w=eval(input('    w='))

    if w<w0:

        print(f"结果：{180*(math.atan(2*rou*w/(w0**2-w**2)))/math.pi}\n")

    elif w>w0:

        print(f"结果：{180*(math.atan(2*rou*w/(w0**2-w**2))+math.pi)/math.pi}\n")

    elif w==w0:

        print(f"结果：{180*(math.pi/2)/math.pi}\n")

    n+=1

    if n==8:

        n=1