技术领域
[0001] 本
发明属于医疗
电子仪器的测试计量技术领域，涉及一种低波峰因数多频率正弦信号的合成方法。
背景技术
[0002]
生物电阻抗(Bioimpedance)是生物对象对所施加的交流
电流的流动所呈现出的宏观抑制特性，其反应了生物系统的微观电特性。相对其他传统的测量方法，生物电阻抗测量法具有低成本、无创性、实用性和快捷性，从而成为生物医学中提取生理、病理信息的重要手段，对研究
生物组织病变、人体成分评估、
疾病早期监测等方面都有着广泛的运用前景。
[0003] 早期生物电阻抗测量主要采用分时测量来获得阻抗谱信息，如直接采用数字频率合成器(DDS)产生
激励信号，其特点是在一个测量时刻，激励信号只有一种频率，因此需要分时扫频才能完成所有激励频率的测量。扫频法虽然简单易行，其拥有较高的
信噪比(SNR)，但由于如心跳、呼吸、血液流动等动态变化使基于扫频的分时测量法不能及时响应生物系统的变化从而无法得到生物的瞬时特性，失去了宝贵的诊断信息。
[0004] 为了更方便的进行动态生物电阻抗测量，采用多频率
同步信号成为电阻抗测量的趋势。多频率同步信号可以在很短的时间内同时对被测阻抗进行不同频率的激励，从而可以同时获得多个频率下的频率响应。如拥有较强的抗干扰性的伪随机信号m序列可以在连续的
频谱分布下得到系统的阻抗信息，但一般而言，阻抗的测量只需要在一定
频率范围(通常为5kHz至1MHz)选择有限个稀疏、宽频的频率分布点。因此拥有连续频谱的m序列会因为
能量有限而降低各个频率下的幅值，从而降低对应的SNR。
[0005] 多频率正弦(multi-sine)信号是一种通过包含有限个不同频率、幅值、
相位的
正弦波从而使其具有任意频率分布与幅值分布的MFS。由于可以自定义频率分布和幅值分布，其很高的能量利用率和灵活性使其在工业测量中得到了广泛的运用。但没有进行优化过的multi-sine信号拥有较大的波峰因数(CF)，会破坏测量系统的线性特征，因此如何合成具有低CF的multi-sine信号是当今研究的主要方向。一般而言CF的优化方法主要分为以下两类：
解析法和
迭代法。解析法可以直接通过公式计算出各个谐波的优化相位，但主要针对于等幅，频率均匀分布的频谱，并且CF值并没有达到最优。而迭代法一般都可以得到比解析法更优的CF值，代价是
算法复杂且计算量很大。迭代法的优点是通过大量的迭代运算，CF值可以降到很低，但其算法一般很复杂且运算时间过长。因此如何综合两种方法的优点，设计出一种针对任意幅值、频率分布的multi-sine信号，都具有高效、通用的低波峰因数合成方法对于实际的运用非常关键。
发明内容
[0006] 本发明的目的是提供一种低波峰因数多频率正弦信号的合成方法，解决了
现有技术中存在的低波峰因数多频率正弦信号的合成方法复杂、效率低的技术问题。
[0007] 本发明所采用的技术方案是，一种低波峰因数多频率正弦信号的合成方法，包括以下步骤：
[0008] 步骤1：将各次谐波频率fk输入计算机，k代表谐波的次数；
[0009] 步骤2：获取各次谐波频率对应的将功率pk的和归一化后的ak，并将其输入计算机；
[0010] 步骤3：对各次谐波的初相位 进行Schroeder编码，并将其输入计算机；
[0011] 步骤4：根据步骤1、步骤2和步骤3得到的参数fk、ak、 通过IDFT生成初始时域信号 N为一个周期内的码元数；
[0012] 步骤5：通过迭代过程优化初相位，最终获取最优相位 结合步骤1、步骤2中的各次谐波的fk、ak，通过IDFT生成最终的时域信号，即低波峰因数的多频率正弦信号；
[0013] 步骤6：将最终得到低波峰因数的多频率正弦信号输出。
[0014] 本发明的特点还在于，
[0015] 步骤1中fk的表达式如下：
[0016] fk＝k×f1 (1)；
[0017] 其中，k＝1,...,K，K为谐波的总次数，f1为基波频率。
[0018] 步骤2中ak的获取方法为：
[0019] 将各次谐波对应的功率pk的和归一化为1，即 令 因为则
[0020]
[0021] 步骤3中各次谐波的初相位 为：
[0022]
[0023] 其中， 为基波初相位，取0～2π的随机数。
[0024] 步骤4中的IDFT过程中码元产生的频率fs不小于信号谐波最高频率fK的64倍。
[0025] 步骤5中生成最优相位 得到最终时域信号的具体方法为：
[0026] 5.1，首先构造时域信号钳位函数：
[0027] 时域信号钳位函数采用
原型为一个上下限给定的对数函数y＝loga(i-b)，其中，i为当前迭代次数，i＝0,....,X，X为总的迭代次数，钳位值y的下限设为0.7，上限设为1，时域信号钳位函数中参数a和b的计算方法为：
[0028] 分别将(0，0.7)，(X，1)代入对数函数y＝loga(i-b)，通过Matlab中的solve()函数求解得到a、b的值，得到时域信号的钳位函数；
[0029] 5.2，采用步骤5.1中构造的钳位函数y＝loga(i-b)，对当前时域信号进行钳位，得到钳位后的信号，其中i＝0,1,2,....,X，当i＝0时，当前时域信号即为步骤4中生成的初始时域信号，具体操作为：
[0030] 计算当前
阈值H(i)， 当 大于H(i)时， 当 小于(i)
等于H 时， 保持不变；
[0031] 5.3，对钳位后的信号进行FFT，得到新的各次谐波对应的相位 结合步骤1、步骤2中的各次谐波的fk、ak通过IDFT生成新的时域信号，返回步骤5.2进行下一次迭代；
[0032] 重复步骤5.2～5.3，直到X次迭代结束，得到最优相位 下的最终时域信号，即低波峰因数的多频率正弦信号。
[0033] 本发明的有益效果是，本发明采用Schroeder编码和V.D.O迭代算法的思想得到对于生物电阻抗测量非常理想的低波峰因数多频率正弦信号，合成过程简单、易实现。本发明的低波峰因数多频率正弦信号合成方法与构成信号的基波频率无关，且当合成信号的各次谐波幅值相等时，得到的最终波峰因数与幅值的选取也无关，对合成低波峰因数的不同谐波次数的多频率正弦信号具有通用性。
附图说明
[0035] 图2是本发明具体实施方式中使用的钳位函数曲线；
[0036] 图3是采用随机相位合成的31次谐波多频率正弦信号
波形；
[0037] 图4是采用Schroeder相位合成的31次谐波多频率正弦信号波形；
[0038] 图5是采用本发明方法合成的31次谐波多频率正弦信号波形；
[0039] 图6是采用随机相位合成的63次谐波多频率正弦信号波形；
[0040] 图7是采用Schroeder相位合成的63次谐波多频率正弦信号波形；
[0041] 图8是采用本发明方法合成的63次谐波多频率正弦信号波形；
[0042] 图9是采用随机相位合成的127次谐波多频率正弦信号波形；
[0043] 图10是采用Schroeder相位合成的127次谐波多频率正弦信号波形；
[0044] 图11是采用本发明方法合成的127次谐波多频率正弦信号波形；
[0045] 图12是采用随机相位合成的255次谐波多频率正弦信号波形；
[0046] 图13是采用Schroeder相位合成的255次谐波多频率正弦信号波形；
[0047] 图14是采用本发明方法合成的255次谐波多频率正弦信号波形。
具体实施方式
[0048] 下面结合附图和具体实施方式对本发明进行详细说明。
[0049] 本发明的一种低波峰因数多频率正弦信号的合成方法，包括以下步骤：
[0050] 步骤1：将各次谐波频率fk输入计算机，其中，fk表达式如下：
[0051] fk＝k×f1 (1)；
[0052] 其中，k＝1,...,K，k为谐波的次数,K为谐波的总次数，f1为基波频率，本发明该具体实施方式中取1Hz；
[0053] 步骤2：获取各次谐波频率对应的将功率pk的和归一化后的ak，并将其输入计算机，具体过程为：
[0054] 将各次谐波对应的功率pk的和归一化为1，即 令 因为则
[0055]
[0056] 步骤3：对各次谐波的初相位 进行Schroeder编码，并将编码后的 输入计算机，具体过程如下：
[0057]
[0058] 其中， 为基波初相
角，取0～2π的随机数；
[0059] 通过Schroeder相位可以优化初始相位，减少后续的迭代次数。
[0060] 步骤4：生成初始波形，具体方法如下：
[0061] 根据步骤1、步骤2和步骤3得到的参数fk、ak、 通过IDFT生成初始时域信号N为一个周期内的码元数；
[0062] 对于离散序列xn，若一个周期的
采样点数N不够多，则得到的CF值会与真实值相差甚远，实验表明，若IDFT过程中码元产生频率fs不小于信号最高频率fK的64倍，可以保证得到的CF值拥有4位
精度，所以，本发明的IDFT过程中采用的码元产生频率fs不小于信号最高频率fK的64倍。
[0063] 步骤5：通过迭代过程获取最优相位 得到最终的时域信号，具体过程如下：
[0064] 5.1，首先构造时域信号钳位函数：
[0065] 本发明中采用时域信号钳位函数的原型为一个上下限给定的对数函数y＝loga(i-b)，其中，i为当前迭代次数，i＝0～1000，钳位值y的下限设为0.7，上限设为1，时域信号钳位函数中参数a和b的计算方法为：
[0066] 当i＝0，y取0.7即：
[0067] 0.7＝loga-b (4)；
[0068] 当i＝1000，y取1，即：
[0069] 1＝loga1000-b (5)；
[0070] 联立(4)、(5)得：
[0071] a-a0.7-1000＝0 (6)；
[0072] 建立函数y＝a-a0.7-1000，通过Matlab中的solve()函数可以求解得：(i+137.80)
a≈1137.80，b≈-137.80，得到时域信号钳位函数为y＝log1137.80 ，如图2所示，横轴为迭代次数i，纵轴为钳位函数值y；
[0073] 5.2，采用步骤5.1中构造的钳位函数：y＝log1137.80(i+137.80)对当前时域信号进行钳位，得到钳位后的信号，其中i＝0,1,2,...,1000，当i＝0时，当前时域信号即为步骤4中生成的初始时域信号，具体操作为：
[0074] 计算当前阈值H(i)， 当 大于H(i)时， 当 小于等(i)
于H 时， 保持不变；
[0075] 5.3，对钳位后的信号进行FFT，得到新的各次谐波对应的相位 结合步骤1、步骤2中的各次谐波的fk、ak通过IDFT生成新的时域信号，返回步骤5.2进行下一次迭代；
[0076] 重复步骤5.2～5.3，直到1000次迭代结束，得到最优相位 下的最终时域信号，即低波峰因数的多频率正弦信号。
[0077] 为使下一次迭代时，时域信号超过对应迭代次数的阈值时钳位到该阈值，若钳位值一直是一个常量且取值过大，如取当前时域信号峰值的0.999倍，则会使算法收敛过慢，若取值过小，如取当前时域信号峰值的0.7倍，会在迭代的后期产生震荡，而得不到最优值，所以将时域信号的钳位函数值y的下限设为0.7，上限设为1。
[0078] 步骤6：将最终得到的低波峰因数的多频率正弦信号输出。
[0079] 其中，步骤2中将各次谐波pk的和归一化为1的原理如下：
[0080] 时域连续信号x(t)的波峰因数CF：
[0081]
[0082] 其中，xpeak代表信号峰值，xRMS代表信号在有效频带内的均方根值(RMS)。
[0083] 通过Parseval's Theorem和公式(7)，离散时域信号的CF值可以表示为[0084]
[0085] 将pk归一化后得到的multi-sine信号的CF则可以简化为：
[0086]
[0087] 因此，对于将各次谐波频率对应的功率pk的和归一化后的multi-sine频谱，其CF的计算简化为寻找离散multi-sine序列xn的峰值，可以看出对于任意给定谐
波数量K时，归一化
功率谱的multi-sine信号的波峰因数只与该信号的峰值有关而与信号的有效值无关。
[0088] 综上所述，从将各次谐波频率对应的功率pk的和归一化后的时域信号可以方便的看出该信号下的CF值，也就是说最终得到的时域信号的峰值，即为该信号下的CF值。
[0089] 图3到图14的纵坐标归一化后的最大函数值即为最终得到的时域信号的峰值，也就是多频率正弦信号的CF值，从图3到图14中可以看出不同谐波次数的多频率正弦信号的CF值分别为2.6674、1.7811、1.4164、3.0700、1.6945、1.4009、3.2844、1.6727、1.4020、3.5379、1.6632、1.3914，因此，采用本发明方法合成的多频率正弦信号的CF值均低于采用随机相位合成的多频率正弦信号的CF值及仅采用Schroeder相位而没有钳位迭代的多频率正弦信号的CF值，证明了本发明的方法可以获得低波峰因数的多频率正弦信号。
