word恒星距离的计算方法摘要光速不变原理是相对论的基础，但为什么不变？假设光子为了保持速度恒定能够自动衰减它的能量，利用这个假设推导出了红移的计算方法：多普勒红移zD=v/(c-v),引力红移zg=GM/(rc2-GM)，距离红移zd= exp(L/4669)-1，并指出了多普勒公式中的错误和根据广义相对论推导出的引力红移公式的瑕疵利用普森公式、距离与红移的关系及两个Ia型超新星参数，求出了关联系数a及星际消光系数b用两种方法计算出了155个恒星的距离并比较了它们的差别证明了光子在星际中传播时能量衰减很小，但它却是引起红移的主要原因；星际的绝对消光也很小，但忽略它会引起很大的误差最后，解释了引起计算误差的原因，提出了测量关联系数的方法，说明了哈勃常数是距离的函数，分析了类星体的发光原因,并估算了类星体的距离和半径关键词：距离红移，关联系数，消光系数，哈勃常数，史瓦西半径1 引 言 / 光子在真空中传播时，其速度是常数，它是相对论的基石，并已被多次证明是正确的为什么光速是常数？光子是如何实现的？假设光子为了保持速度的恒定能够自动衰减（或增加）自身的能量1.1 根据假设，如果光源以速度v离开观察者，则它的动量为P = mv，光子为了保证速度的恒定，也必须降低能量克服光子的初始动能（f为发射源的原始发射频率，f’为接收到的频率，h为普朗克常数，c为光速，m为光子运动质量）即：h*f - P*c=h*f’ （1）把P=m*v和m=h*f/c2 代入上式得 ：h*f - h*f*v/c=h*f’， f’=(1-v/c)*f。
zD=[c/(c-v)]-1=v/(c-v)，（2）或：vD=c*z/(1+z) （3）根据多普勒原理，观察者和发射源彼此远离时的频率关系为(v0为观察者移动速度，vs为发射源移动速度，V为介质速度）：f’=[(V-vo)/(V+vs)]*f （4）当v0=0，vs=0.6V时，f’=f*而当v0=0.6V，vs=0时，f’=f*彼此靠近时的关系为：f’=[(V+vo)/(V-vs)]*f （5）当v0=0，vs=0.6V时，f’=2.5*f 而当v0=0.6V，vs=0时，f’=1.6*f这明显是个错误，在现有的实验条件下很容易验证，宇宙学红移就是根据多普勒原理推导出了z=v/c，当z>1时必须用洛伦兹变换，但在同一坐标系下为什么变换呢？多普勒公式应改为：彼此远离时：f’=[(V-vo)/V]*[(V-vs)/V]*f （6）彼此靠近时：f’=[(V+vo)/V]*[(V+vs)/V]*f （7）红移和速度的关系也应改为：z=f/f’-1=[V/(V-v)]-1=v/(V-v) （8）或：v=V*z/(1+z) （9）1.2 根据前面的假设，光子在一个引力场中，为了保证速度的恒定，必须衰减能量用于克服引力，则下式成立（g为重力加速度，r为星球的半径，）：h*f - m*g*r =h*f’ (10)把m=h*f/c2，g=G*M/r2，（G为引力常数，M为星球质量）代入上式得：h*f -h*f*G*M/r*c2=h*f’，f’=（1-G*M/r*c2)*f ,z=f/f’-1=G*M/（r*c2-G*M） (11)根据广义相对论推得引力红移公式为[1]：zg=[1- 2G*M/(r*c2)] -1/2-1，简化后的公式为[2]：zg=GM/rc2用本文的能量衰减论所得的等式为：zg=GM/(rc2-GM)。
为了比较哪个是合理的，分别把地球（M=5.965*1024kg,r=6.371*106m)、太阳(M=1.9887*1030 kg, r=6.955*108m)和当r=2GM/c2时的数据分别代入，计算得到结果如表1所示．表1 不同理论计算的引力红移量比较Table 1. parison of gravitational redshift in different theoretical calculationsGravity redshiftearthsunr=2GM/c2[1- (2G*M/r*c2)] -1/2-16.9526895352*10-102.12339151*10-6∞GM/（rc2-GM）6.9526896728*10-102.12338925*10-61GM/rc26.9526896680*10-102.12338474*10-6由表1中可以看出：当r=2G*M/c2时，根据广义相对论推导的公式红移为无穷大，也就是说，在史瓦西半径上光子是无法逃逸的如果光子的运动方向是史瓦西半径球面的切线方向，则光子永远不能逃离这个球面，但如果光子的运动方向是垂直于球面向外，光子为了保持速度的恒定，就会自动降低能量克服引力。
设光子的频率为f，这时光子的能量为：E=h*f ，史瓦西半径球面的势能为：Es=m*g*r=(h*f/c2)*(G*M/r2)*r= (h*f/c2)*G*M/r = (h*f/c2)*G*M/(2*G*M/c2) = h*f/2可见，光子完全有能力克服引力逃逸出来，这说明根据广义相对论推导出的红移公式是有瑕疵的，从史瓦西半径球面垂直逃逸所产生的红移为z=(E/Es)-1=1光子的逃逸半径为：Re=G*M/c2 ，通过上面的分析，说明假设是正确的 2 计 算物体的传播是需要能量的，光子也不例外根据上面的假设，设光子每传播百万秒差距的衰减量为k，则（EL为光子传播L百万秒差距后所具有的能量，E0为光子的初始能量，L为光子传播的距离）：EL=(1-k)L*E0 （12）光子在传播L百万秒差距后的距离红移为： zd=f/fn-1=（E0/EL）-1=[1/(1-k)L]-1=（1-k)-L-1 （13） 对zd+1=1/(1-k)L两边取自然对数可得：Ln(1+zd)= - L*Ln(1-k)，即：L=[ -1/Ln(1-k)]*Ln(1+zd)，令 红移-距离的关联系数 a= -1/Ln(1-k) （14）得：L=a*Ln(1+zd) （15）如果单位时间和单位面积，在理想状态下从恒星接收的光子数为1，设光子每传播百万秒差距所减少的数为n，则实际所能接收的光子数为：Y=(1-n) L（16） 根据普森公式，星际间和波长无关的绝对消光为[3]：Aa= -2.5* lg(Y)= - 2.5L*lg(1-n) （17）令 和波长无关的消光系数 b= - 2.5*lg(1-n) （18） 则：Aa=b*L （19）L=10^[(m-M-25-Aa)/5] (m为视星等、M为绝对星等)， 或：5*lg(L)=m-M-25-b*L （20）对于同一颗恒星，两种计算方法所得的距离应该是相同的，把(15)代入(20)得：5*lg[a*Ln(z+1)]=m-M-25-b*a*Ln(1+z)取两个恒星的参数代入上式可得方程组：（SN1999fv的参数为：m=24.5，M=-19.27，z=1.19。
SN2007uh的参数为：m=22.21，M=-19.8，z=0.53，也可以取其它中、高红移的Ia超新星 )5*lg[a*Ln(1.19+1)]=24.5+19.27 - 25 - b*a*Ln(1+1.19)5*lg[a*Ln(0.53+1)]=22.21+19.8 - 25 - b*a*Ln(1+0.53)用上面的两个参数代入到下面的等式可以得到两种计算恒星距离的方法：DZ=4669*Ln(z+1) （21）5*lg(DL)=m-M-25-0.0002622*DL（21）计算的结果如表2所示(从开放的超新星目录中随机抽取了155个Ia超新星，但为了节省篇幅，只列出了10个，其它的显示在图1中)：表2 恒星的距离及两种计算方法的百分差Table 2. the distance of stars and the percentage deviation between two methods of calculationNamezmMDz(Mpc)DL(Mpc)%v(km/s)UDS10Wil-1.96%6108SN150G-1.42%4160SN2003ak0.90%-2550SN2002hp-0.35%866SCP-06F120.13%-297SNLS-04Dlow0.50%-980SNLS-04D4hu0.49%-788SN1997cd-19-0.97%1179SN2009jo-19-1.36%1083SN2011ct-19-2.94%892为了更直观地显示出计算结果，把表1中的数据表示在图1中，其上面的曲线为红移与距离（是DZ的，DL的曲线几乎与之重合）的关系，下面的曲线为用两种方法计算所得结果的百分差。
当z值小于0.5时，DZ和DL之间出现偏差的几率较大，当z值大于1.6时，误差有增大的趋势，但由于观测到恒星数目较少，没有分析出原因，表2中的最后一列为恒星的径向速度5000 L（Mpc）400030002000100010 (1-DZ/DL)% z -10图 1 .红移与距离的关系及两种计算结果的百分差Figure 1. The relationship between the redshift and the distance， and the percentage deviationof the two calculated results3讨论3.1各参数的物理意义：由a= -1/Ln(1-k) ， a=4669 ，得k=1-exp(-1/4669)=0.000214156，它表示每百万秒差距能量的衰减为百分之0.0214，光子传播一年的能量损失率为：△E%=1-exp[Ln(1-k)/3261600]=6.57*10-11，所产生的红移为：z=[1/(1-△E%)]-1=6.57*10-11，相当于速度为v=0.02m/s的恒星所产生的多普勒红移，这就是科学界一直忽略不计的原因。
计算绝对星等时，已经考虑到了消光，但只是根据一系列统计化的色余方法，对和波长无关的消光考虑的少由b= - 2.5*lg(1-n) ，b=0.000262177 ，得n = 0.00024 ，它表示每百万秒差距光子数的损失为百分之0.024，也就是说：光子传播一年的损失为：△N=1-10^[lg(1-n)/3261600]=7.4*10-11个（每年每100亿个光子损失不到1个），这也是科学界一直忽略不计的原因恒星一直在运动，但我们只能观察到视向运动，当z值小于0.5时，DZ和DL之间出现偏差的几率较大，这是由于多普勒红移与距离红移具有可比性，在表2中，运动速度最大的是UDS10Wil，其值为6108千米/秒，方向是离我们而去，所产生的多普勒红移为：6108/(300000-6108)=0.02083当z值大于0.5时，速度的影。