信号相位匹配原理的正弦信号参数总体最小二乘估计方法

合集下载

一种提取噪声中正弦信号的总体最小二乘法_喻胜

∑i =1 ∑j =1 |X ij -Y ij |
( E)
m
n
/2 2 1
( 15)
文献[ 4] 的定理二指出 : 在以上( 15) 式定义的 Frobenius 量最小的意义下 , 与( 14) 式实际观测矩阵 H 最接近的秩为 2K 的增强数据矩阵 H 可按如下的迭代方法求出( 算法 2) : ① 将 Hankel 数据矩阵 H 作奇异值分解展开并取前 2K 个最大奇异值而将其余奇异值置0 : H(E) = ② 将( 16) 式所得 H
2k
-k )=0 n =2K +1 , …, N
( 4) ( 5) ( 6) ( 7)
式中的线性预测系数{ ak } 具有如下的中心元素对称性 : a 2K -k =a k k = 1 , … , K 且 a 0 =a 2K = 1 若能由观测序列求出{ ak } 的估计值 , 则以下 z 的多项式 Χ ( z )=
2 参数。这里假定 e ( n) 为具有零均值、方差为 σ 的宽带噪声。将纯谐波过程
本文于 1999 年 2 月收到。喻胜 : 博士研究生 ; 闫波 : 讲师 ; 陈光
: 教授 , 博导。
第 2 期一种提取噪声中正弦信号的总体最小二乘法 · 7 · y( n)= 写作 y( n) =
k
ω kn + k ) ∑k =1 Ak cos(
n
k
( 2) ( 3)
= 1 ( bkzk ∑ k
+bk zk ) n = 1, …, N
＊
＊
式中 bk = ( 1/2) Ak exp( j k) ; zk =exp( jω k) 。很明显 , 当被测过程的各频率量 ω k 求出后 , 以上 N 个方程组成的方程组( 3) 是复振幅 bk 的线性方程 , 不难由传统的线性最小二乘法求出诸 { bk } 。所以 , 本文的讨论仅限于由观测序列{ x( n) } 估计各正弦信号频率。我们知道 , 对任意频率的 K 个纯净实正弦信号( 2) 可以用一个 AR( 2K ) 线性预测模型来描述 : n ∑k =0 aky(

现代信号课件第4章最小二乘滤波

归一化均方误差性能评估
NMSE越小，说明滤波器的性能越好，信号处理的效果越接近原始信号。
归一化均方误差（NMSE）是另一种衡量滤波器性能的指标，它表示信号经过滤波器处理后的误差相对于原始信号的均方误差的比例。
NMSE的计算公式为：$NMSE = frac{MSE}{MSE_{total}}$，其中$MSE_{total}$为原始信号的均方误差。
加权最小二乘滤波
加权最小二乘滤波是在线性最小二乘滤波的基础上引入了权重因子，以调整误差的权重。
通过给不同的误差项赋予不同的权重，加权最小二乘滤波能够更好地适应不同的噪声分布和信号特性。
加权最小二乘滤波在处理具有不同特性的信号和噪声时能够获得更好的滤波效果。
03
最小二乘滤波的算法实现
递归最小二乘滤波
04
在控制系统中，最小二乘滤波用于系统辨识和参数估计等。
02
最小二乘滤波的数学模型
线性最小二乘滤波
线性最小二乘滤波是一种常用的信号处理方法，通过最小化误差的平方和来估计信号中的未知参
数。
它假设信号和噪声之间存在线性关系，通过解线性方程组来得到
最优估计值。
线性最小二乘滤波具有简单、稳定和快速收敛等优点，适用于多
THANKS FOR WATCHING
感谢您的观看
信噪比性能评估
信噪比（SNR）是衡量滤波器在噪声干扰下性能的重要指标，它表示信号与噪声的功率比值。
SNR越大，说明滤波器对噪声的抑制能力越强，信号处理的效果越好。
SNR的计算公式为：$SNR = 10log_{10}frac{P_s}{P_n}$，其中$P_s$为信号功率，$P_n$为噪声功率。
自适应滤波算法优化

第四章线性系统参数估计的最小二乘法

下面讨论更为一般的情况。假设在t1, t2, …, tm时刻对Y及X的观测值序列已经被我们获得，并且用
y(i), x1(i), x2(i), x3(i), … i = 1,2, …, m 来表示这些观测数据。显然，可以用 m 个方程组来表示量测数据与估计值之间的关系
⎧ y(1) = θ1x1(1) +θ 2 x2 (1) +L+θ n xn (1)
从图中可看到，前两条线都仅能满足两个点的要求，而对其它点的误差都很大，其 6 个点的误差平方累计分别为 0.49 和 0.42。第三条线能满足三个点的要求，但误差平方累计更大，为 1.58。显然我们需要找到一条更为理想的直线来取得较小的误差。例如图中的红色短划线，它的方程为 y=1.697 + 0.294x，误差平方累计为 0.25。这条线是怎样得到的呢？它是用最小二乘法得到的。
z
−2
，在其输入端加入 M 序列输入后
所得到的输出输入数据见下表，请利用这些数据辨识出系统的传递函数的系数。
k
1
2
3
4
5
6
7
8
9
10
输入 u
1
0
1
1
0
0
1
1
1
0
输出 y -0.45 -0.01
1.15
2.56
1.92
-0.30 -0.80 0.91 2.92 2.40
解：已知系统阶数 n=2，有 4 个未知数。将式（4.4）展开 y(k) = −a1 y(k −1) − a2 y(k − 2) + b0u(k) + b1u(k −1) 根据要求，观测次数 N>2n+1，取 N 为 6，k=3

最小二乘参数辨识方法及原理

v ( k ) 是均值为 0 的随机噪声。
2.2 一般最小二乘法原理及算法
z (k ) a i y (k i) bi u (k i) v (k )
i 1 i 1 n n
如果定义
h ( k ) [ y ( k 1), y ( k 2 ), , y ( k n ), u ( k 1), u ( k 2 ), , u ( k n )]
1 1 1
1 1 1
1
1
1
z1 1 1 ( z 1 z 2 ) 2 z2
r 1 0 0 1 1 4 r 1 1 1 1
2、最小二乘辨识方法的基本概念
通过试验确定热敏电阻阻值和温度间的关系
t (C ) R ( )
t1 R1
t2 R2

tN
1
tN RN
RN
1
R a bt
• 当测量没有任何误差时，仅需2个测量值。 • 每次测量总是存在随机误差。
y i R i v i 或 y i a bt v i
v i y i R i 或 v i＝ y i a bt i
常见做法：
太复杂使
max | y i R i |
1 i N
N
最小 /* minimax problem */ 不可导，求解困难
使 |y
i 1
i
Ri |
最小
最小
使 |y
i 1
m
i
Ri |
H
2
1 1
r R 0
0 4r

基于最小二乘法的信号处理与参数估计

基于最小二乘法的信号处理与参数估计随着科技的不断发展，我们对信号处理与参数估计的需求越来越高。

在这个背景下，基于最小二乘法的信号处理与参数估计成为了一种非常重要的方法。

本文将深入探讨基于最小二乘法的信号处理与参数估计，从数学角度分析其原理与优劣，以及在实际应用中的具体方法与注意事项。

一、基本原理最小二乘法是一种数学优化方法，它的基本思想是通过最小化误差的平方和来求解参数估计值。

在信号处理与参数估计中，最小二乘法的应用十分广泛，比如回归分析、滤波器设计、谱估计等。

以回归分析为例，最小二乘法可以用来计算线性回归模型中的系数，使得观测数据与模型预测值之间的误差最小。

具体来讲，假设我们有$n$个样本$(x_i,y_i)$，其中$x_i$为自变量，$y_i$为因变量。

在最小二乘法中，我们要求解的是模型参数$\beta$，使得残差平方和最小，即：$$\min_{\beta}\sum_{i=1}^{n}(y_i-\beta^Tx_i)^2$$这里$\beta$是一个列向量，$x_i$是一个行向量，$^T$表示向量的转置。

根据最小二乘法的原理，我们可以推导出$\beta$的解析表达式：$$\beta=(X^TX)^{-1}X^Ty$$其中，$X$为$n\times p$的设计矩阵，每行表示一个样本，每列表示一个自变量。

$y$为$n\times 1$的观测向量，每个元素表示一个因变量。

二、优劣分析最小二乘法具有以下优点：1. 理论基础扎实。

最小二乘法是一种数学优化方法，其理论基础非常扎实，在数学上可以比较轻松地进行推导和分析。

2. 稳定性好。

最小二乘法的解析表达式通常具有良好的数学性质，比如对称性、非负性等，因此可以保证算法的稳定性。

3. 可扩展性强。

最小二乘法可以用于各种不同的信号处理与参数估计问题，比如回归分析、滤波器设计、谱估计等。

但是，最小二乘法也存在一些缺点：1. 对异常值敏感。

在实际应用中，如果数据中存在异常值，那么最小二乘法会受到很大的干扰，导致估计结果不准确。

线性系统参数估计的最小二乘方法

y X e
un (1) un (2) un ( N )
从而 ˆ ( X T X )1 X T y 如何推广到 MIMO 系统？
8
2. 非线性系统的辨识
y(k ) a0 a1u (k ) a2u 2 (k ) e(k )
选取
x1 ( k ) u (k )
得新解
ˆ(k 1) [ X T (k 1) X (k 1)]1 X T (k 1) y (k 1)
12
其中
X (k ) y (k ) , y (k 1) , X (k 1) T x (k 1) y (k 1)
20
辅助变量法
y (k ) ai y (k i ) biu (k i ) e(k )
i 1 i 1 n n
传统做法是由
X T X X T y X T e
得到
ˆ ( X T X )1 X T y ( X T X )1 X T e ，
ˆ。如果 X 与 e 不相关，则 E ( )
i 1 i 1 i 1 n n n
令 A(q 1 ) 1 aq i ， B(q ) bi q i ，则系统写成
i 1 i 1
n
n
A(q 1 ) y (k ) B(q 1 )u (k ) e(k )
其中
n
e(k ) A(q 1 ) (k )
引入变白滤波器
e( N ) 将进入ˆ 而使它偏离
2. 如 N 2n ，此时观测方程个数大于参数个数，这是一个解超定方程问题，可取
ˆ(i )) 2 J ( N ) ( y (i ) T (i )

信号检测与估计理论(8)第八章最小二乘估计

−1
T T −1 T
)( ( = x ( I − H( H H) H ) x
)
（8-11）
上式推导中，利用了 I − H(H T H ) H T 是幂等矩阵的性质，即 A2 = A 。
−1
J min 也可写成另一种形式
Jmin = xT x − xT H HT H HT x
ˆ = xT x − Hθ
w[n] = E (w[n]) + w' [n]
式中 w ' [n] 是零均值噪声，这时观测数据则为
x[n] = A + E (w[n]) + w ' [n]
显然，样本均值估计实际上是 A + E (w[n]) 的估计，这时的 LSE 是有偏估计，不是 MVU 估计。也就是说，必须假设观测数据是由确定信号和零均值噪声组成，在信号参量的正确选择下，误差 e[n] = x[n] − s[n] 在平均意义下趋于零，则（8-1）式的最小是合理的。如果假设的直流量信号模型不正确，例如描述观测数据的模型是 x[n] = A + Bn + w[n] ，则这个模型误差也将引起 LSE 是有偏的。
对于信号模型是待估计参量的线性函数的最小二乘情况我们称为线性最小二乘问题，如例 8-1。否则，这个问题称为非线性最小二乘问题，如例 8-2。求解非线性 LS 问题是通过栅格寻优或迭代最小化方法（前面已介绍）。应该注意的是线性最小二乘问题是指信号是待估计参量的线性函数，信号本身不必是线性的。
p ×1
维的，信号
s = [s[0] s[1] " s[ N − 1]]T 是待估计参量的线性函数，假设
s = Hθ
（8-8）

基于最小二乘的变采样率超低频正弦扫频信号幅值识别方法

基于最小二乘的变采样率超低频正弦扫频信号幅值识别方法朱景振;蒋瑜;陶俊勇【摘要】Aiming at the problem of long signal period, nonstationary and real time amplitude identification of ultra low frequency (0.1~1 Hz) sine sweep signal of shaking table, variety rate sampling method was used to transform the sine sweep signal into standard sinusoidal signal approximately in a very short time. Then use the least square(LS) fitting method to identify the amplitude of the signal in real time. Simulation results show that the method is effective.%针对振动台超低频(0.1~1 Hz)正弦扫频振动信号周期长、非平稳、幅值实时识别困难的问题,采用变采样率的方法将极短时间内的正弦扫频信号近似变成频率恒定的标准正弦信号,结合最小二乘拟合方法,可以有效识别非整周期正弦扫频信号的幅值.仿真结果表明所研究的方法是有效的.【期刊名称】《强度与环境》【年(卷),期】2017(044)004【总页数】5页(P60-64)【关键词】超低频正弦扫频;变采样率;最小二乘;幅值识别【作者】朱景振;蒋瑜;陶俊勇【作者单位】国防科技大学机电工程与自动化学院,长沙 410073;国防科技大学机电工程与自动化学院,长沙 410073;国防科技大学机电工程与自动化学院,长沙410073【正文语种】中文【中图分类】TB534+.2正弦扫频振动测试是一类基本的环境与可靠性试验类型，对测试产品的动力学特性和保障产品的环境适应性及可靠性水平具有重要作用。

第二十四讲：最小二乘估计、波形估计-课件

z(i)H (i)X 0w (i)
z k z(1 ) z(2 ) ... z(k )T
H k H ( 1 ) H (2 ) ... H (k )T
zk HkX0wk
X ˆ0 (k ) [(H k)T H k] 1 (H k)T z k
批处理算法，运算量太大。
递推算法：
X ˆ 0 ( k ) X ˆ 0 ( k 1 ) K ( k ) ( z ( k ) H ( k ) X ˆ 0 ( k 1 ) )
{ z(k), k= n0, n0+1,...,nf }对区间内的某一个时刻 n(n0<n<nf)的信号进行估计，内插也称为平滑。
数据
n0
n
nf
sˆ ( n )
波形估计宜采用可建立递推算法的线性最小均方估计或最小二乘估计。
z ( n ) s ( n ) v ( n ) n 0 ,1 ,...,N 1
A
n0
Aˆ
1
N 1
z(n)
N n0
例：正弦信号频率的估计
s(n)cos2f0n
N1
J(f0) (z(n)cos2f0n)2 n0
最小化难以得到闭合性形式的解，原因是信号与未知参数f0之间存在高度的非线性关系。
zHθv
zz1,z2,...,zNT
θ1,2,...,MT
vv1,v2,...,vNT
θ ˆlsw(H TW H )1H TW z
讨论：
（1）当观测噪声的均值为零时，最小二乘与加权最小二乘是无偏估计。
E[θˆls ] (HT H)1HT E[z] (HT H)1HT E[Hθ v] (HT H)1HT Hθ θ
（2）估计的方差阵
V a r ( θ l s ) E { [ θ θ ˆ l s ] [ θ θ ˆ l s ] T } ( H T H ) 1 H T R H ( H T H ) 1

最小二乘估计

r
缺点：对所有测量数据同等看待
最小二乘估计——线性最小二乘估计

•
•
加权最小二乘法对不同置信度的测量值采用加权的办法分别对待强调可靠数据的重要性
LS改为：J ( ) ( x H )
T
W ( x H )
N×N维加权矩阵
最小二乘估计——线性最小二乘估计

加权最小二乘法
ˆ 此时的LSE为： ( H T WH ) 1 H T Wx
ˆ 一般的LSE为： ( H T H ) 1 H T x
此时估计的均方误差为： ( H T WH ) 1 H TWRWH ( H TWH ) 1 一般的均方误差为：
( H T H ) 1 H T RH ( H T H ) 1
最小二乘估计——线性最小二乘估计

加权最小二乘法
W diag[w(0), w(1), , w( N 1)]
改写形式
ˆ[ N ] N A[ N 1] 1 x[ N ] ˆ A N 1 N 1
1 ˆ ˆ ˆ A[ N ] A[ N 1] ( x[ N ] A[ N 1]) N 1
老的估计
修正项
最小二乘估计——序贯最小二乘估计

若对观测数据加权，w( n)
1
2 n
最小二乘估计——线性最小二乘估计
当系统测量噪声V是均值为0，方差为R时
性质1. 最小二乘估计即是无偏估计，有： ( ) E ˆ
ˆ) E[ ( H T H ) 1 H T x] E ( E[( H T H ) 1 ( H T H ) ( H T H ) 1 H T x] ( H H ) H E ( H x)
n 1 n 1 n 1

信号相位匹配原理的正弦信号参数的最小二乘估计

信号相位匹配原理的正弦信号参数的最小二乘估计
孙进才;朱维杰;孙轶源;刘理
【期刊名称】《自然科学进展》
【年(卷),期】2004(014)010
【摘要】根据信号相位匹配原理,提出了一种正弦信号频率、振幅和相位参数的最小二乘估计方法.推导了利用单传感器接收信号的参数估计的最小二乘估计计算公式,给出了已知频带内的未知频率,振幅和相位的信号参数估计的搜索算法.分析了最小二乘法使用的方程数、信噪比和采样频率确定后FFT的序列长度对参数估计精度的影响.理论和仿真结果说明,该方法不仅能降低估计频率带宽以外的噪声,而且还可降低被估计频率信号带宽内的噪声,提高了低信噪比时的信号参数估计精度.该算法简单、快速,具有工程应用前景.
【总页数】5页(P1204-1208)
【作者】孙进才;朱维杰;孙轶源;刘理
【作者单位】西北工业大学航海工程学院,西安,710072;西北工业大学航海工程学院,西安,710072;西北工业大学航海工程学院,西安,710072;西北工业大学航海工程学院,西安,710072
【正文语种】中文
【中图分类】N941
【相关文献】
1.利用信号相位匹配原理的正弦信号参数估计 [J], 孙轶源;朱维杰;孙进才;李志舜
2.正弦信号参数估计的最小二乘算法 [J], 孙进才;朱维杰;孙轶源;刘理
3.一种基于相位匹配的高精度估计正弦信号频率新方法 [J], 杜文超;孟小芬;王国宏;张治杰
4.信号相位匹配原理的正弦信号参数总体最小二乘估计方法 [J], 刘理;朱维杰;孙进才
5.自相关域信号相位匹配原理的正弦信号参数估计 [J], 陈志菲;孙进才
因版权原因，仅展示原文概要，查看原文内容请购买。

最小残差法的方向估计

最小残差法的方向估计
陈志菲;孙进才;侯宏
【期刊名称】《西北工业大学学报》
【年(卷),期】2011(029)003
【摘要】针对基于相位匹配原理的最小二乘法(LSMSPM)无法估计具有相同中心频率的窄带入射源的波达方向问题,提出了最小残差法(MRM)的方向估计方法.首先将阵列接收信号转换到频域,在中心频点处将相位匹配后的数据模型转换为估计频谱参数的线性方程组,采用加权最小二乘法解该方程组,以残差平方和的倒数作为MRM方法的方向估计算子.MRM方法提高了LSMSPM方法的方向估计性能,相同中心频率的等功率窄带双源入射时其角度分辨能力优于常规波束形成方法.
【总页数】4页(P347-350)
【作者】陈志菲;孙进才;侯宏
【作者单位】西北工业大学航海学院,陕西西安710072;西北工业大学航海学院,陕西西安710072;西北工业大学航海学院,陕西西安710072
【正文语种】中文
【中图分类】TN911.7
【相关文献】
1.最小残差法加速局部加权LSSVM求解及其应用 [J], 林超;王华杰
2.窄带信号波达方向估计仿真实验 [J], 廖昌俊;潘晔;刘继芝
3.基于LoRa信号的行走距离和方向估计 [J], 常兆心;张扶桑;牛凯;张大庆
4.基于奇异值分解的虚拟阵列波达方向估计算法 [J], 徐朋豪;高春林;董华
5.一种高精度的神经网络波达方向估计算法 [J], 李国彬;朱清超;宋晓鸥
因版权原因，仅展示原文概要，查看原文内容请购买。

1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。

１引言．
文献〔提出了利用三元阵的信号相位匹配原理的相干干扰抵消原理，能降低噪声和相１］干干扰的影响。文献「基于信号相位匹配原理进行正弦信号参数估计，频率和幅度估计精２］度较高，三个序列中但当有两个序列的噪声“ 相” 接近 “ 相” 在该搜索频点同或同时，将出现奇异峰，在低信噪比条件下尤为严重。因文献［增加了此，［３］方程数，利用最小二乘（）（ＬＳ
一２－１８
Ｘ）凡（＋（）ｍ＝ｋ戈ｋ（ｋ）
根据文献〔［证明结果，４５Ｊ」对于信号有以下关系
（２）
Ｓ（＝， ’）Ｓ（ｍＳｋ一＋二）ｋ（嵘（，）） ‘ ｋ
当Ｎ＞＞Ｍ时，有
＜３）
１）（鲁ＳＩ‘ ；／，（Ｉ｝ｋＩ‘ Ｓ
ｅｓ
， ‘
ｅ
，ＡＩ｛ＭＮ ’Ａ一
自月口万．．．．．．．．．卜．
２正弦信号参数总体最小二乘估计．
设接收信号为数字信号ｘｎ，传感器的（）信号与噪声是线性可加的，依次取Ｍ段长度为Ｎ
的信号得
ｘ（＝ｒ（１＋［（１．Ｓ＋ｍ一）ｎ＋ｍ一）Ｚｉ）］ｉ］
ＦＴ得Ｆ
（）１
其中ｉ，．．－１．．ｓ为单弦＝１．．ｍ，．．（频正信号，（斯白声。０，．Ｎ１２．Ｍ，２．）＝．ｉ）ｎ）ｉ为高噪对式（）１作
关键词：信号相位匹配原理，总体最小二乘，奇异值分解，信号增强
Ｍｅｏｏｔａｌｓｓｕｒｏｐｒｍｔｅｉａｏｆｓｕｏａｓｎｌｕｉｔｄｏｌｔａｅｆａｅｒｍｔｎｉｓｉｌａｓｇｈｆｅｑｓａｔａｅｓｉｏｎｄｉｓｔｒｇｎ
ｏｓｎｌｃｅ－ａ．ｅｅｄｅｔｉｄｎｉｔｓｅｔａｃｒｙｈｐｒｅｒｆａｍｔｄｈｓＳｍａｎｍｎｓｅｈｐｐｒｈｃｕｃｏｔａｍｔｓｉｇａｈｐｅｏｍｓｏｎａ，ｉｅａｆａｅｅｃｎｉｐｖｄｔｒｕｉｔｍｔｄｔａｌｓｓａｓＡｅｔａｎｔｏｓｖｄａｂｍｒｅｆｈｂｓｇｅｏｏｏｌｔｒ．ｒｔｇｂｅｅｅｏｕｅｙｎｈｒｅｈｆｅｑｅｆｒｉｈｔａｕｔｅｅｒｄｔｕｉｔｔｈｉｅｓｎｌａｃｅｔｔｅｅｗｌｂｉｐｏｅ．ｓｕｔｎａｓｇｅｎｕｏｉａｅｈｎｅｎｈｆｃｉｅｒｄＴｅｌｉａｎｈｃｑｆｅｇｎｍ，ｅｔｌｍｖｈｉａｏｍｒｕｓｗｔｍｔｄｓｔｉｔｐｐｒａｔｐｒｒａｃｔｎｍｔｄｅｌｓｏｔｔｅｏｐｅｅｎａｅｈｓｅｅｅｏｎｅｔｅｏｏｓｔｈｈｈａｅｈｒｎｄｈｅｅａｂｔｒｆｍｈｈａｅｈｆ
方法求解，可部分消除噪声对被估信号的影响，从而提高了低信噪比情况下的正弦信号参数估计精度。如果考虑到观测数据的误差，那么它的解是不准确的。因此，本文在此基础上进
行了改进，采用总体最小二乘（Ｓ（）Ｔ估计，减小Ｌ它能误差的影响。另外，信噪比较低时，利
用信号增强技术对信号进行预处理后再进行估计，也可改善估计性能。
ｌｓｓｕｒ．ｅｔａｓａｑｅＫｙｒｓｃｌｏｓｎｌｃｅ－ａ，ｌｔａｓＳｇｌｖｌｄｃｍｏｉｎｅｗｄＰｉｉｅｉａｍｔｄｈｓＴｔｌｓｓｒ，ｕｒｅｏｐｓｉ，ｏｒｐｆｎｇａｈｐｅｏｅｑｅｉａａｅａａｕｎｕｔｏＳｇａｅｈｎｅｅｔｉｌａｃｍｎｎｎ
令（）Ｘ（）Ｎ＇，Ｎ（）Ｎ（）Ｎ－ｏ－，将（式边凡ｋ＝ｍ（）Ａ，，ｏ益（）１，４右的ｏｋＷｍｋ０－０ｋ＝ｋ＋１叫ｍｋｏｋ）））
Ｓｋ移左然求平得（简省括中ｏ，ｏ至边，后模的方为便，去号的ｋ）（）
‘
．
一
１
ｅｓ
Ｊ
一２－１９
而采总最二（Ｓ６法解。Ｌ方主想是：到Ｓ最解应当用“ 体小乘”Ｔ）方求ＴＳ法的要法找的优Ｌ１］
使得以一卜
（ＥＳ＝Ｐ￡Ａ＋）一
（３）１
方中扰矩Ｅ：总误｝－ｚ刁这｝为阵Ｆｂｉ范，程的动阵和的体差ＥＩ最、里．矩的ｒｎ数将ＩＦ，｛，｝；ｏｕｅｓ
１１２１．
一
只
，
户
ｌｗｅ一
门
厂ｔ十
门１，ｓｅｗｅｓｅｓｅｓｅｗｅ
．
君
－－
一
Ｉ２一ＮＮ ’ＡＡ｝｝，Ｉ３一ＮＮ ’ＡＡＩ｝２
１，
』
ｅｓ
自，
Ｌ
氏
一
凡
刁
Ｌ＝－Ｔｕｖ（３３３
将式（）写为矩阵形式９
Ｐ＝Ａ・Ｓ＋６
（０１）
（１Ｉ）
厂
Ａ・Ｓ尸一￡
厂一
Ｒ（只ｅ几一）
！．
门
其中Ａ
一一
．．！．．．Ｌ
Ｒ（凡）ｅ几一
Ｉ（只ｍ几一）ｌ（几）ｍ凡一
１，ｅ．，ｗｅ．．ｌｌｌｓｅ
（ＬｄＣｄ０Ｃ＋）＝Ｌ＝（５１）
由１，式（）采用ＴＳ５Ｌ方法求解期望信号Ｓ最优解问的题等价于：找到一个最小范数扰动
矩使：ＣＬ非秩阵；外，Ｃ的空量ｄ最一元应阵Ｌ得Ｃ二＋为满矩此解出：零间向的后个素
ＡｓａｔａｃｒｙｐｒｅｒｍｔｎｓｕｏａｓｎｌｃｎｐｏｅｉｔｃｓｏｂｔｃＴｅｕｃｏａｍｔｅｉａｏｆｉｓｉｌａａｂｒｖｄｈａｅｒｈｃａｆａｅｓｉｏｎｄｉｓｅｔｒｇｎｅｆｌｅｒｉｏｓｎｌｏｅｕｉｔｍｔｄｌｓｓｕｒ，ｃｉｂｓｏｔｐｎｉｅｏｒｏｉａｔｎｉｂｓｇｅｏｏｅｔａｓｗｉｓｅｎｒｃｌｗａｆｏｓｙｎｈｔｇｅｈｆｑｅｈｈａｄｈｉｐａｅ
信号相位匹配原理的正弦信号参数总体最小二乘估计方法
刘理朱维杰孙进才
（西北工业大学航海工程学院，西安，７０７）１０２
摘要利用信号相位匹配原理进行正弦信号最小二乘估计可提高低信噪比时的估计精度，在此基础上进行改进，应用总体最小二乘方法对正弦信号进行参数估计，能进一步提高精度。利用信号增强技术对观测序列进行预处理，能改善估计效果。仿真实验表明，本文方法性能优于最小二乘方法。
Ｓ
一一
！！ｌｅｗｅＬ
门
Ｒ（，ｓｅｅｅ）ＳｗｅＩ（，ＪＭＳ）
１气 ‘
ｌ
１１１ ‘ ．一且ｒ
Ｒ（Ｐ＿ｌ（一１ｅＭ一Ｍｍ氏凡＿Ｐ１））
．．Ｊ
．
尸
一一
Ｉ几Ｐ２｝一ＩＰ３
那么，期望信号Ｓ的最小二乘解为
Ｓ（＇一Ｔ＝十二Ａ）ＡＡＰＡ，
文献［也正是采用了这种方法。［３］
（２１）
但是，传统的最小二乘方法问题中，系数矩阵Ａ是没有噪声干扰的，而在本文问题中，
Ａ和尸均受到观测误差的影响。所以这里求解期望信号Ｓ传统最小二乘方法并不准确，时，
若知号率ｋ时根信相匹原［ｌ用叫，）式（两己信频为ｏ，据号位配理ｌ，ｔｌｔ２一乘３边，ｋ）
可
（
６
得（
ｔｌ
Ｘ（）Ｎ１＝，ｏＮ（）Ｎ－ｋ叫ｍ’ ｍ（）Ｓｋ＋，ｏｍ（）一（ｋＷｍｋ（）ｋ＋ｌ，０－，ｏ ‘ ）
ｐｉｃｌｏｓｎｌｃｅ－ｈｓｒｉｅｉａｍｔｄｐａｅｎｐｆｇａｈ
ＬｕＨＷｅｉＳｎｃｉｉＬｉＵｉｕＪａＺｊｅｉｎ
（ｏｅｅｒｅｉｅｎＮｒｗｓｒＰｙｃｎａＵｉｓｙｉ７０２，ｔｅｅｌｅｈｉｌｖｒｔＸ＇，０７）ＣｌｇｏＭａｎＥｇｅｒｇｏｈｔｏｔｃｎｅｉ，ｎ１ｌｆｉｎｎｉｎａ
Ｓ和Ｓ众之间的相位差满足
（）４
Ａ二、１Ｏ竺、
其中ｍ１，，＝，二Ｍａ２－，
Ｎ
（５）
一般取ＦＴ样本长Ｆ的度Ｎ远大于程组，方数Ｍ忽略Ｓ（）ｍ对估计信号Ｓ（）影响。ｋ，ｋ的
这样方程（）可改写为２
Ｘ，ＳｋＮ－＋Ｘ “，Ｎ－ｋ，），（）Ｎ一＋ｍ（（＝（Ｗｍｋｋ）１， ‘ １）
当为非零值以便于利用ｄ的末元素归一化求出Ｓｏ以上最优化问题可以由矩阵Ｃ的奇异值分解求出。将Ｃ写成奇异值分解向量积展开形