[发明专利]奇异谱分析截断频谱信号重建方法及其应用

摘要

本方法关键是建立信号谱奇异谱分析数学模型,从采集到的已有的截断频谱数据中抽取模型参数(奇异点和奇异谱分析的加权系数),再根据奇异谱分析数学模型,推出截断频谱所隐含的高频分量,解决截断频谱信号重建中的伪影问题。不论含噪声的还是无噪声的、对实际的还是对计算机模拟的截断频谱数据的实验结果,都表明了它是一种高精度的高频信号重建方法,重建精度和速度大大优于现有重建方法。它适用于磁共振截断频谱图象重建和其它截断频谱信号重建。

主权项

1、奇异谱分析截断频谱信号重建方法其特征在于采用下列操作步骤为：(其中，F=[f^(0),f^(1),...,f^(n-1)]T]]>为截断频谱，对应的数字信号为f＝[f(0)，f(1)，...，f(N-1)]T，n＜N，)(1)运用下式计算各尺度下的小波系数：cj,k=1NΣu=0n-1f^(u)ψ^j,k(u)‾,]]>k＝0，1，...，N-1，其中尺度j取0，1，2，3；(2)检索极值将随着尺度的增加而增加的极点，即确定奇异点bi， i＝1，2，...，Q。(3)用伪逆矩阵法求解以下方程组，求加权系数y＝{yb1，yb2，...，ybQ}T：(4)用下式计算截断频谱的各高频分量：f^(u)=Σi=1QybiWbi^(u),]]>u＝n、n+1，...，N-1；(5)对频谱矢量F=[f^(0),f^(1),...,f^(n-1),f^(n),...,f^(N-1)]T]]>进行Fourier 反变换，得重建信号。