Python数学方程

Question

我试图用python得到数学方程的解，如下所示：

# Open netCDF file
fd1 = nc.Dataset(nc_file1, 'r')

# Read variables from the netCDF file
rng1 = fd1.variables['range'][:]
tim1 = fd1.variables['time'][:]
pwr1 = fd1.variables['pwr'][:]
dpl1 = fd1.variables['dpl'][:]
nfft1 = fd1.variables['nfft'][0]
pn1 = fd1.variables['pnoise'][:]

# Close netCDF file

fd1.close()

# Specify beam 
ibeam1 = 0

# Time convertion from seconds to hours
tim1 = tim1/3600.0

# Select data and transpose
p_plot1 = pwr1[ibeam1]
for it1 in range(len(tim1)):
   p_plot1[it1] = p_plot1[it1] - pn1[ibeam1][it1] - 10.*np.log10(nfft1)
p_plot1 = p_plot1.transpose() 

**#Determine Height**
A=6378.1370
Tetha=math.cos((37.5 * math.pi) / 180)
for j in range(len(tim1)):
   for i in range(len(rng1)):
       D[i]=sqrt(A**2+(rng1[i]**2)+(2*A*(rng1[i])*Tetha)
       height1(i,j)=(D[i]-A)

从数学方程式（确定高度）我得到了错误。它可能是python中表达式的语法。实际上，我关注的是三个数据变量，即tim1、rng1和pwr1。从rng1数据中，我将得到高度1，只需使用该方程将rng1转换为高度1。 提前谢谢你。在

Answer 1

问题是争吵

D[i]=sqrt(A**2+(rng1[i]**2)+(2*A*(rng1[i])*Tetha)

括号不匹配，就这样在末尾加一个就行了。在

^{pr2}$