主权项 |
一种循环平稳泊松信号快速模拟方法,其特征在于:采用高斯模型拟合一类循环平稳泊松信号的标准累积轮廓,所述高斯拟合方法生成泊松信号的仿真信号的方法就是生成在周期内具有高斯分布特性的光子序列,产生光子序列仿真信号的流程如下:1).利用RXTE探测器的实测数据,获取信号流量强度,计算单个周期内的信号平均光子流量强度λ<sub>s</sub>和噪声平均流量强度λ<sub>b</sub>;2).对泊松信号的标准累积轮廓进行高斯拟合,通过调整成分数使得高斯拟合曲线与标准轮廓的拟合程度最好,即使拟合得到的轮廓与标准累积轮廓的差异在均方意义下最小,记录高斯拟合成分的参数,包括成分数C,幅度系数α<sub>i</sub>、均值μ<sub>i</sub>和方差δ<sub>i</sub>(1≤i≤C);3).在一个周期内分别对高斯成分曲线积分,求出曲线包含的面积a<sub>i</sub>,以及计算总面积<img file="FDA0000976927300000011.GIF" wi="266" he="87" />4).求出各子高斯成分曲线包含的面积在总面积中所占的比例p<sub>i</sub>=a<sub>i</sub>/A,并由此计算出各子高斯成分的平均光子流量强度f<sub>i</sub>=λ<sub>s</sub>·p<sub>i</sub>;5).根据f<sub>i</sub>计算各子高斯成分在一个周期内的信号光子数N<sub>i</sub>=poissrnd(f<sub>i</sub>),其中poissrnd(·)是泊松随机数生成函数,并根据输入的流量强度返回光子数;6).对于每个高斯成分f(φ;{α<sub>i</sub>,μ<sub>i</sub>,δ<sub>i</sub>}),根据其均值μ<sub>i</sub>和方差δ<sub>i</sub>产生一个高斯随机数ph=normrnd(μ<sub>i</sub>,δ<sub>i</sub>),其中ph看作是到达探测器的一个信号光子,其数值表示该信号光子在轮廓中的位置,循环N<sub>i</sub>次,产生高斯成分f(φ;{α<sub>i</sub>,μ<sub>i</sub>,δ<sub>i</sub>})在一个周期内的N<sub>i</sub>个信号光子,对这N<sub>i</sub>个信号光子时间从小到大排序得到该周期内的光子时间序列,根据下式:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>h</mi><mrow><mo>(</mo><mi>φ</mi><mo>;</mo><mo>{</mo><msub><mi>α</mi><mi>i</mi></msub><mo>,</mo><msub><mi>μ</mi><mi>i</mi></msub><mo>,</mo><msub><mi>δ</mi><mi>i</mi></msub><mo>}</mo><mo>)</mo></mrow><mo>=</mo><munderover><mo>Σ</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>C</mi></munderover><msub><mi>f</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>φ</mi><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000976927300000021.GIF" wi="590" he="127" /></maths>合成子高斯轮廓光子序列,最终生成拟合轮廓的光子序列。 |