定日镜场余弦效率python代码

from math import *
import math
import matplotlib.pyplot as plt
import numpy as np
from tqdm import tqdm
import pandas as pd
from 太阳高度角和方位角 import cal_gamma,cal_alpha
def cal_Lambda(x,y,tower_height):
    #参数：定日镜中心坐标，吸收塔高度
    #计算定日镜中心与集热器中心连线与数值方向夹角，弧度制
    return atan(sqrt((x**2+y**2))/tower_height)

def cal_Thita(x,y):
    #参数：定日镜中心坐标
    #计算定日镜中心的方位角，弧度制，从正北顺时针旋转
    a=0
    if(x>0 and y != 0):
        a=atan(x/y)
    elif(x<0 and y != 0):
        a=atan(x/y)+pi
    elif(x>0 and y == 0):
        a=pi/2
    elif(x<0 and y == 0):
        a=pi *3/ 2
    elif (x ==0 and y > 0):
        a = 0
    elif (x == 0 and y < 0):
        a = pi *2
    return a
def cal_cosLoss(position_data,ST,N,tower_height,day):
#参数：position_data:所有定日镜的坐标,读取文件得来镜中心
#time:一天的某个时刻
#N:定日镜数量
#tower_height:吸收塔高度
#day:距离春分的天数
#功能：计算某天整个