3.1.4　程序设计

3.1.4.1　计算流程

图3.1　功率谱计算流程

3.1.4.2　程序代码

（1）功率谱分析子程序。程序输入信度r0 （0.10、0.05、0.025、0.01 或0.005）、原始序列数据X（n）、序列长度N，即可分析计算序列的周期。子程序没有设计通过变量输出，如果需要可以将原始谱、平滑谱、平均谱和否定域上界等中间结果及计算的序列周期通过子程序变量输出。计算功率谱分析计算的VisualBasic6.0程序代码如下：

（3）GAMMQ函数。函数根据不完全Γ函数的余函数的收敛快慢，调用子程序GSER或GCF来计算余函数近似值。代码如下：

（4）GSER子程序。该子程序用级数展开式计算不完全Γ函数近似值。程序输入Γ函数的参数A、自变量X1，输出用级数展开式计算不完全Γ函数近似值GAMSER、Γ函数的对数值GLN。其程序代码如下：

（5）GCF子程序。该子程序用于连分式展开计算不完全Γ函数的余函数近似值。程序输入参数A、X1、GLN，输出参数GAMMCF。程序代码如下：