假设引力场的速度为光速C，那么，星系之间的作用一般都是要通过几十年几百年以上的光速距离才能达到，这个时候星系之间的运动平衡状态早破坏了。这种无序运动导致星系之间是一片混乱，事实上，宇宙是和谐的，是有序的，假设它们之间的作用场为引力场，自然会存在一种超光速的（瞬时速度）引力作用场维持它们之间的运动平衡。
由此可以推测，引力场速度是超光速的一个稳定速度。

而所有的光速实验证明，光速都存在着电磁场的参与，而电磁场的理论速度也是光速C。那么，引力场的速度就不一定是光速。

在目前的情况下，可以假设引力场速度为光速参与下的速度，所以，这个引力场达朗贝尔公式还是可以有存在的意义。

电磁场波动方程为
▽²A－εμ ∂ ²A/∂t² =0
▽²φ－εμ ∂²φ/∂t² =0
电磁波方程的通解一般为
φ(r，t ) = C +∮E.dr ＋ c².∮B.dt ＋ φ₀ .exp[i ω.(t±r/c)]
A(r，t ) = C +∮B.d r＋∮E.dt ＋A₀ .exp[i ω.(t±r/c)]
其中电场E 、磁B、C都是积分常矢。
假设光子能量场W(r，t ) 的波动方程为
▽²W － (1/c²)∂²W/∂t² = 0
则其解为
W(r，t ) = C + ∮F·dr + ∮P·dt＋W₀·exp[i ω·(t±r/c)]
其中力F 、功率P、C都是积分常矢
不妨设一个矢量场T(r，t ) 为
T(r，t ) ＝M(r，t ) +φ(r，t )
假设它们的波速都为光速C，则复合矢量场T(r，t ) 的波动方程为
▽²T－(1/c²) ∂²T/∂²t= 0
那上式的解为
T(r，t ) = C +∮E.dr ＋ c².∮B.dt + ∮F·dr + ∮P·dt＋ T₀ .exp[i ω.(t±r/c)]
不妨令C=0，则有
T(r，t ) = ∮E.dr ＋ c².∮B.dt + ∮F·dr + ∮P·dt＋ T₀ .exp[i ω.(t±r/c)]
可知电磁场或光波在空间中沿着任意直线传播。
在一个周期内，上式可表示为
T(r，t ) = φ＋Ac+ p.c＋h.ω＋ T₀ .exp[i ω.(t±r/c)]
当一个光子穿过一个杨氏缝时，屏幕上就会显示一个量子点，当许多个光子穿过一个杨氏缝时，各个光子在穿过这点时由于在周期性的路径位置的不同，叠加时就会表现出一个干涉的波动纹路。
光子是一个量子波，而不是粒子。光子波本身含有电场，磁场，动量，角动量，能量的物理性质。

由于时空都是量子化的，令dr=ndr₀，dt=ndt₀，L=nh则
u=dr/dt=ndr₀/ndt₀=dr₀/dt₀
假设核外运动电子的牛顿总动能量E为
E=m₀.u²/2－L²/2m₀.r
那么，上式有
E=m₀.u²/2－h²/2m₀.r₀²
所以，牛顿动能量也是量子化的，它们是一个最小的量子化值。
对于核外运动电子，根据向心力相等的原则，有
m₀.u²/r = ke²/r²
得
r=L²/m₀. ke²
令r=n.r₀，L=nh，则得
r₀ =nh²/m₀. ke²
代入到牛顿动能量公式中，有
E=m₀.u²/2－m₀. (ke²)²/2n².h²
如果处于外磁场B中，对波长最小长路径r₀积分，则核外运动电子的牛顿总能量有
E=m₀.u²/2 + eB.h/2m₀－m₀. (ke²)²/2n².h²

目前的薛定谔方程，最终的描述仍然是牛顿物理量。

很显然，量子力学的理论需要重建。

设复合波函数y(x，t )的波动方程为
d²y/dx²－(1/c²)d²y /dt²= 0
不妨把上式分解为
.d²y/dx² + k². φ = 0
d²y /dt² + ω².φ = 0
其中，ω² = k².c² 。
当 复合波函数 y(r，t ) 等于 φ(r，t )时，则上式化为
.d²φ/dx² + k². φ = 0
d²φ/dt² + ω².φ = 0
这就是常见的波动方程，它的解很显然为
φ(r，t ) = φ₀ .exp[i ω.(t±r/c )]
但时，当 y(r，t ) 不等于 φ(r，t )时，其解为
y(x，t )=Bx＋bt＋C + φ₀ .xp[i ω.(t±x/c )]
设复合波函数y(x，t )的达朗贝尔方程为
d²y/dx²－(1/c²)d²y /dt² = f(x)
可以移项为下面的形式
【 d²y/dx² - f(x) 】－(1/c²)d²y /dt² = 0
不妨令，【 d²y/dx² - f(x) 】=d²φ/dx² ，则
d²φ/dx²－(1/c²)d²y /dt² = 0
上式不妨分解为
d²φ/dx² + k². g = 0
d²y/dt² + ω².g = 0
即
d²y/dx² - f(x) + k². g = 0
d²y/dt² + ω².g = 0
当y不等于g时，很显然，上面第二式的解
y(t ) = bt＋C + g₀.exp[i ω.(t±x/c )]
很显然，则第一式的解为
y(x ) = －∮∮ f(x) ds＋Bx＋C +g₀ .exp[i ω.(t±r/c)]
当y等于g时，
d²y/dx² - f(x) + k². y = 0
d²y/dt² + ω².y = 0
第二项显然为
y(t ) = y₀.exp[i ω.(t±x/c )]
很显然，上面第一式的解就有点复杂了。

设复合波函数y(x，t )的波动方程为
d²y/dx² － (1/c²)d² w /dt²= 0
当 复合波函数 y(r，t ) = w(r，t )时，则有
d²y/dx² － (1/c²)d² y /dt²= 0
令A = k². g = ω²/c². g ，不妨把上式分解为
.d²y/dx² + k². g = 0
d² y /dt² + ω².g = 0
其中，ω² = k².c² 。
当 复合波函数 y(r，t ) = g(r，t )时，它的解很显然为
y(x，t ) =g(r，t ) = g₀ .exp[i ω.(t±r/c )]
当 y(r，t ) 不等于 g(r，t )时，其解为
y(x，t ) = Bx＋bt＋C + g₀.xp[i ω.(t±x/c )]
当复合波函数 y(r，t ) 不等于 w(r，t )时，则有
.d²y/dx² + k². g = 0
d² w /dt² + ω².g = 0
不妨令
d²y/dx² =d²φ/dx² - f(x)
d²w/dt² =d²φ/dt² - f(t)
则有
d²φ/dx² - f(x) + k². g = 0
d²φ/dt² - f(t) + ω².g = 0
当复合波函数φ不等于g时，很显然，上式的解为
φ(x ) = ∫ ∫f(x) ds＋Bx＋C + φ₀ .exp[i k.r]
φ(t ) = ∫ ∫f(t) dtdt＋bt＋c + φ₀ .exp[i ω.t]
当复合波函数φ等于g时，上面的方程化为如下常见的形式
d²φ/dx² - f(x) + k². φ = 0
d²φ/dt² - f(t) + ω².φ = 0
上面其实就是二阶常系数非线性微分方程，很容易求解。

▽²φ ＋ ∂ B/∂t = －ρ/ε
它又可以表示为
▽²φ ＋ ∂ (▽×A )/∂t = －ρ/ε
上面这个方程在电磁场理论中有相似的描述，它的表达式为
▽²φ ＋ ∂ (▽. A)/∂t = －ρ/ε
很显然，在这个方程的产生过程中，它们的数学物理意义都不同。

对波尔量子理论的一点个人看法
波尔理论三个假设：
第一假设是“原子定态假设”
玻尔原子理论认为，原子只能处于一系列不连续的能量状态中。而且在这些状态中，电子虽做变速运动，但并不向外辐射电磁波（即处于能量守恒的稳态），这样的相对稳定的状态称为定态（补充一下，对于原子来说，电子所处的定态是有很多个的）。
第二假设是“电子跃迁假设”
玻尔原子理论认为，电子绕核转动处于定态时不辐射电磁波，但电子在两个不同定态间发生跃迁时，却要辐射（吸收）电磁波（光子）；其需要的能量及外界光的频率（或者向外界释放的光的频率）由两个定态的能量差值决定hv=E2-E1。
第三假设是“轨道量子化假设”
由于能量状态的不连续，因此电子绕核转动的轨道半径也不能任意取值，玻尔原子理论认为在某个半径下的运动，原子的能量（和半径的大小）只能取到固定的若干个数值。
上面这个假设似乎没有涉及到时空的本质，我在波尔假设的基础上，提出时空的量子化假设，然后来推导一遍。
设时空是量子化的，则有
dr=ndr₀
dt=ndt₀
u=dr/dt=ndr₀/ndt₀=dr₀/dt₀
角动量L=nh
.........................................................................................1式
假设质量为m₀的自由粒子的运动来自于能量场的作用，那么，它所表现出来的作用力F为
F＝m₀.(a +u²/r )
它的运动轨迹是一条曲线
核外电子动体的加速度a取决于粒子速度u的大小，所以
m₀a=m₀.u²/r - ke²/r²
对上式积分，得
m₀.u²/2 = - L²/2m₀.r² + ke²/r
因此，核外运动电子的总能量E为
E=m₀.u²/2 + L²/2m₀.r² - ke²/r
.........................................................................................2式
在外磁场B中，其总能量E为
E=m₀.u²/2 + L²/2m₀.r² - ke²/r ± eBL/m₀
.........................................................................................3式
把1式代入3式中，eBL/m₀=eBh/m₀得
E=m₀.u²/2 + h²/2m₀.r₀² - ke²/r ± eB.h/m₀
.........................................................................................4式
如果粒子处于定态，则 m₀a=0，有
m₀.u²/r - ke²/r²=0
所以
r =L²/m₀.ke²
代入时空量子化条件，有
r =n²h²/m₀.ke²
即
r₀ =nh²/m₀.ke²
.........................................................................................5式
把5式代入4式中，整理后得
E=m₀.u²/2－m₀.(ke²)² /2n².h² ± eB.h/m₀
.........................................................................................6式
上式第一项表示核外运动电子的动能量m₀.u²/2是一个不变值，处于波尔理论的定态假设之中。第二项表示核外运动电子的量子化向心势能量值，每一次向心势能量的变化都是量子化的，表示波尔理论的能量跃迁假设，第三项表示核外运动电子处于外磁场B中的约束能量值，在经典量子理论中解释为运动电子的自旋能量。
令u=ke² /h，则有
E=【m₀.(ke²)² /2h²】 （1 - 1/n²） ± eB.h/m₀
.........................................................................................7式
这就是氢原子核外运动单电子的牛顿能量系统的表达式。

对于简谐振动，作用振子上的作用力f为
f=K.X
它的动能E为
E＝∫fdX＝K.X²/2
假设振动谐振子动能势能E为
E＝m.u²/2＝m.ω²X²/2
两相比较，显然，有
K＝m.ω²
很显然，能量E对X的微分表达式为
d²E/dx² = K
又在简谐振动中，振动位移X可以表示为
X(r，t ) = X₀ [sin (ω.t±ω.r/c )＋cos(ω.t±ω.r/c )]
其中r≤x，c为波速。
由于时间轴与X轴相差九十度角，假设时空轴是对称的，所以，可设
X(r，t ) = X₀ [sin (ω.t)＋cos(k.r)]
或
X(r，t ) = X₀ [sin (k.r)＋cos(ω.t)]
其中ω²= c².k²。
u=dX/dt = ωX
a=du/dt=ω²X = u²/X
谐振子的向心势能量V为
V＝∫fdX＝m∫(u²/X)dX＝- L²/2m.X² ＝－m.ω²X²/2
假设谐振子不仅存在动能，还存在势能，所以，可以假设
d²E/dx² － K = d²H/dx²
把 d²H/dt² ＝c². d²H/dx² 代入上式，得
d²E/dx² －（1/c²）d²H/dt² = K
可以把上式分解为
d²E/dx² + k².g = K
d²H/dt² + ω²g = 0
其中ω²= c².k²。
因此，上面方程的解为
E(r) = K.X²/2＋∮Fdr＋g₀ .exp[±i k.r]
H(t ) = ∮Pdt＋ g₀ .exp[±i ω.t]
我们知道， K.X²/2是一个最小量子值，跟量子数无关。可设
K.X²/2 =h.ω/2
∮FdX= n h.ω
因此，E(r) 的能量子本征值为
E(r) = ( n＋ 1/2)h.ω＋g₀ .exp[±i 2nπ]
E(t) = n.h.ω＋g₀ .exp[±i 2nπ]
这就是谐振子的能量本征值。

设复合波函数的波动方程
d²y/dx² － (1/c²)d² w /dt²= 0
当 复合波函数 y(r，t ) = w(r，t )时，则
d²y/dx² － (1/c²)d² y /dt²= 0
或
d²w/dx² － (1/c²)d²w /dt²= 0
令A = -k². g = -ω²/c². g ，不妨把上式分解为
.d²y/dx² + k². g = 0
d² w /dt² + ω².g = 0
其中，ω² = k².c² 。
当 复合波函数 y(r，t ) =w(r，t ) = g(r，t )时，很显然，它的解为
y(x，t )=w(r，t ) =g(r，t ) = g₀ .exp[i ω.(t±r/c )]
当 y(r，t )或 w(r，t )不等于 g(r，t )时，其解为
y(x) = Bx＋C + g₀.xp[i ω.(t±x/c )]
w(t ) = bt＋C + g₀.xp[i ω.(t±x/c )]
当复合波函数 y(r，t ) 不等于 w(r，t )时，则有
.d²y/dx² + k². g = 0
d² w /dt² + ω².g = 0
不妨令
d²y/dx² =d²φ/dx² - f(x)
d²w/dt² =d²Ψ/dt² - f(t)
则有
d²φ/dx² - f(x) + k². g = 0
d²Ψ/dt² - f(t) + ω².g = 0
当复合波函数φ不等于g时，很显然，上式的解为
φ(x ) = ∫ ∫f(x) ds＋Bx＋C + g₀ .exp[i k.r]
Ψ(t ) = ∫ ∫f(t) dtdt＋bt＋c + g₀ .exp[i ω.t]
当复合波函数φ等于g时，上面的方程化为如下常见的形式
d²φ/dx² - f(x) + k². φ = 0
d²Ψ/dt² - f(t) + ω².Ψ = 0
上面其实就是二阶常系数非线性微分方程，参考教科书，很容易求解。
实际上，这就是广义的波动方程以及广义的波函数。

如果不能自己创造出一套数学方法，根本就不应该去研究理论物理。
因为人家走过的路，不会再有路留给你走。

引力场方程
不妨假设自由粒子静止质量为m₀的能量场M(r，t )的自由运动方程为
▽² M－ (1/c2)∂² M /∂t²= － g
.................................................................................10式
g是自由粒子m₀的质量场强度，所以，上式也可理解为质量场引力场方程。它的一般数学解为
M(r，t ) = －∮ g.ds＋f × r＋∮F.dr＋∮P.dt＋W₀ .exp[i ω.(t±r/c)]
其中F、f、P为积分常矢。从上式可以看出，能量场可以分两部分，一部分是跟静止质量有关的能量场，另一部分是跟波动质量有关的量子化能量场。不管是机械波，光波，声波，其波函数所表示的物理量都可以是量子化的，由此可以看到，物质波函数的物理意义也并不是唯一确定的。复合粒子能量场M(r，t ) 等动体具有波粒二象性。
因此，自由粒子m₀的运动速度u为
u =(1－∮g.ds/|M|)c=(1－m0.c²/|M|)c
当∂² M/∂t²= 0时，上式化为
▽²M = － g
上式可理解为牛顿引力场方程。如果考虑到能量场函数M(r，t )的复合性，令M(r) =M₀(r)＋W(r)上式可以分解为
▽² M₀ = － g
▽× (▽× W) = 0
考虑到静止质量场M₀的独立性，上面第一式的特解为
M₀ (r ) = －∮g.ds = m₀.c²
第二式的特解为
W (r) =f × ∮dr
上式与力矩的物理作用相同，其中f为积分常矢，所以
M (r) = m₀.c²＋f × ∮dr
上式第二项是旋量表达式，对于静止质量m₀而言，它的周围可以存在一个带旋度的能量场，任何处于力场f中的质量体都能受到这个力f的作用。在原子核内，力场f 对其它粒子作用表现的最明显。
当g=0时，10式化为
▽² M － (1/c2)∂² M /∂t²= 0
上式可理解为光的能量场波动方场，其一般特解为
M(r，t ) = ∮F.dr＋∮P.dt＋W₀ .exp[i ω.(t±r/c)]
其中F、P都是积分常矢。在一个波长和周期的积分上下限中，有
∮F.dr=∮(p×ω).dr=p.c
∮P.dt=∮(d²h/dt²).dt=dh/dt=h.ω
因此，在r轴上，光波函数M(r) 的特解为
M(r)= p.c＋W₀.exp(±i k.r)
在t轴上，光波函数M(t)的特解为
M(t)= h.ω＋W₀.exp(±iω.t)
由上式的方程解可知，在欧氏空间中，光沿着任意直线传播。在r —t轴上，光波具有偏振的物理特性。
从光波动方程的解可以看出来，光是一种量子化的波，它不具有粒子性。
从这个引力场方程可以看到，黑洞的产生跟引力场强度无关，只跟物质场的波函数特性有关。
物质波函数（波动的能量场）同时也反映了暗物质能量场或引力波的存在，它导致了处于引力波能量场中的物质产生了量子化物理特性。