MUSIC算法详解

接收信号模型

假设目标1在点A处信号模型为：$$x_1[t]=e^{j2\pi ft}\text{\hspace{1em}(1)}$$
因为：$\lambda =\frac{c}{f}$，所以经过路径L后相位变化为：$\phi =2\pi f\frac{L}{c}$
则经过距离L后到达天线1的信号模型为：$$y[t]=A.x_1[t].e^{j2\pi f\frac{L}{c}}\text{\hspace{1em}(2)}$$
假设阵列有 n 个天线，假设他们之间间隔固定为d，经过时间延迟后，到达不同天线的距离为：
$$\{\begin{array}{l}{L}_0=L\\ L_1=L+d.cos(\theta )\\ ⋮\\ L_n=L+(n)\ast d.cos(\theta )\end{array}\text{\hspace{1em}(3)}$$

对应天线在不同时刻接收到的信号模型为：
$$\{\begin{array}{l}{y}_1[\tau ]=A_1.x_1[t].e^{j2\pi f\frac{L_1}{c}}\\ y_2[\tau ]=A_2.x_1[t].e^{j2\pi f\frac{L_2}{c}}\\ ⋮\\ y_n[\tau ]=A_n.x_1[t].e^{j2\pi f\frac{L_n}{c}}\end{array}\text{\hspace{1em}(4)}$$

A表示在不同的路径传输引起接收支路之间的幅度不平衡因素，在这里考虑到的是天线方向图在相同方位处的增益不一致，信号在内部通道之间传输的幅度响应不一致。

不同的距离引起的相位差增益在不同的天线处表示为：
{ ϕ 0 = 1 ϕ 1 = e j 2 π f d × c o s ( θ ) c ϕ 2 = e j 2 π f 2 × d × c o s ( θ ) c = { e j 2 π f d × c o s ( θ ) c } 2 ⋮ ϕ n = { e j 2 π f d × c o s ( θ ) c } n (5) \begin{cases} \phi_0 = 1\\ \phi_1 = e^{j 2 \pi f \frac{d×cos(\theta)}{c}}\\ \phi_2 = e^{j 2 \pi f \frac{2×d×cos(\theta)}{c}}=\{e^{j 2 \pi f \frac{d×cos(\theta)}{c}}\}^2 \\ \vdots\\ \phi_n = \{e^{j 2 \pi f \frac{d×cos(\theta)}{c}}\}^n \\ \end{cases} \tag5 ⎩ ⎨ ⎧​ϕ0​=1ϕ1​=ej2πfcd×cos(θ)​ϕ2​=ej2πfc2×d×cos(θ)​={ej2πfcd×cos(θ)​}2⋮ϕn​={ej2πfcd×cos(θ)​}n​(5)

假设：${\gamma }_i=A_i\ast e^{j2\pi f\frac{L}{c}}$，我们假设，每个接收支路的天线在相同方向的增益一致，内部信号传输通道幅度被校准到一致，则$A_0\to A_n$的值都是一致。

将以上公式进行符号简化:
$$\{\begin{array}{l}{y}_0[\tau ]={\gamma }_1.x_1[t].{\varphi }_0\\ y_1[\tau ]={\gamma }_1.x_1[t].{\varphi }_1\\ y_2[\tau ]={\gamma }_1.x_1[t].{{\varphi }_1}^2\\ ⋮\\ y_n[\tau ]={\gamma }_1.x_1[t].{{\varphi }_1}^n\end{array}\text{\hspace{1em}(6)}$$

Music算法模型

多信号情况下的信号模型

假设信号个数为 $j$
$$\{\begin{array}{l}{y}_0[\tau ]={\gamma }_1.x_1[t]+{\gamma }_2.x_2[t]...+{\gamma }_j.x_j[t]\\ y_1[\tau ]={\gamma }_1.x_1[t].{\varphi }_1+{\gamma }_2.x_2[t].{\varphi }_2...+{\gamma }_j.x_j[t].{\varphi }_j\\ y_2[\tau ]={\gamma }_1.x_1[t].{{\varphi }_1}^2+{\gamma }_2.x_2[t].{\varphi }_2^2...+{\gamma }_j.x_j[t].{\varphi }_j^2\\ ⋮\\ y_n[\tau ]={\gamma }_1.x_1[t].{{\varphi }_1}^n+{\gamma }_2.x_2[t].{\varphi }_2^n...+{\gamma }_j.x_j[t].{\varphi }_j^n\end{array}\text{\hspace{1em}(7)}$$

假设对多天线进行了m次的信号采样，则采样矩阵如下表示
则可以转换成矩阵的方式进行表示为：

$$\left[\begin{array}{ccc}y[0{]}_0& y[1{]}_0& y[2{]}_0...y[m{]}_0\\ y[0{]}_1& y[1{]}_1& y[2{]}_1...y[m{]}_1\\ ⋮& & \\ y[0{]}_j& y[1{]}_j& y[2{]}_j...y[m{]}_j\end{array}\right]=\left[\begin{array}{cccc}{\varphi }_1^0& {\varphi }_2^0& ...& {\varphi }_j^0\\ {\varphi }_1^1& {\varphi }_2^1& ...& {\varphi }_j^1\\ ⋮& & & \\ {\varphi }_1^n& {\varphi }_2^n& ...& {\varphi }_j^n\end{array}\right]\times \left[\begin{array}{cccc}{\gamma }_1.x_1[1]& {\gamma }_1.x_1[2]& ...& {\gamma }_1.x_1[m]\\ {\gamma }_2.x_2[1]& {\gamma }_2.x_2[2]& ...& {\gamma }_2.x_2[m]\\ ⋮& & & \\ {\gamma }_j.x_j[1]& {\gamma }_j.x_j[2]& ...& {\gamma }_j.x_j[m]\end{array}\right]\text{\hspace{1em}(8)}$$

将以上矩阵进行简化：
$$Y=\Phi X\text{\hspace{1em}(9)}$$
$Y$矩阵是已知的采样点，最终的目的是求 $\Phi$;
正常情况下 $\Phi$=$Y{X}^{-1}$在这里的最困难的点是$X$矩阵的值是未知的。

对信号的预先处理

对于信号的处理，可以从两方面进行考虑：

1）Amphtude
2）Phase

假设对一个信号的幅度和相位进行调整，综合为：$C=a\times e^{j\varphi }$
假设对信号进行调整并叠加后
$$y_1\times C_1+y_2\times C_2...y_j\times C_j=0\text{\hspace{1em}(10)}$$
将式子10 ，代入式子7中，并简化，可得：
{ γ 1 C 1 + γ 1 ϕ 1 C 2 + γ 1 ϕ 1 2 C 3 . . . + γ 1 ϕ 1 j − 1 C j } x 1 + { γ 2 C 1 + γ 2 ϕ 2 C 2 + γ 1 ϕ 2 2 C 3 . . . + γ 2 ϕ 2 j − 1 C j } x 2 + ⋮ { γ m C 1 + γ m ϕ m C 2 + γ m ϕ m 2 C 3 . . . + γ m ϕ m j − 1 C j } x j + = 0 (11) \begin{matrix} \{\gamma_1 C_1 + \gamma_1 \phi_1 C_2 + \gamma_1 \phi_1^2 C_3 ... +\gamma_1 \phi_1^{j - 1} C_j \} x_1 + \\ \{ \gamma_2 C_1 + \gamma_2 \phi_2 C_2 + \gamma_1 \phi_2^2 C_3 ... +\gamma_2 \phi_2^{j - 1} C_j \} x_2 +\\ \vdots\\ \{ \gamma_m C_1 + \gamma_m \phi_m C_2 + \gamma_m \phi_m^2 C_3 ... +\gamma_m \phi_m^{j - 1} C_j \} x_j + \\ \end{matrix} = 0 \tag{11} {γ1​C1​+γ1​ϕ1​C2​+γ1​ϕ12​C3​...+γ1​ϕ1j−1​Cj​}x1​+{γ2​C1​+γ2​ϕ2​C2​+γ1​ϕ22​C3​...+γ2​ϕ2j−1​Cj​}x2​+⋮{γm​C1​+γm​ϕm​C2​+γm​ϕm2​C3​...+γm​ϕmj−1​Cj​}xj​+​=0(11)
假设信号都不相关的情况下，则不存在一组非0数 { K 1 , K 2 , K 3... K j } \{K_1,K2,K3 ... Kj\} {K1​,K2,K3...Kj},使得$K_1\ast x_1+K2\ast x_2+K_3\ast x_3=0$;所以在信号不相关的情况下，则该组数只能为： { 0 , 0 , 0 , . . . 0 } \{0,0,0,...0\} {0,0,0,...0}。

意即：
$$\left[\begin{array}{c}{\gamma }_1{C}_1+{\gamma }_1{\varphi }_1{C}_2+{\gamma }_1{\varphi }_1^2{C}_3...+{\gamma }_1{\varphi }_1^{j-1}{C}_j\\ {\gamma }_2{C}_1+{\gamma }_2{\varphi }_2{C}_2+{\gamma }_2{\varphi }_2^2{C}_3...+{\gamma }_2{\varphi }_2^{j-1}{C}_j\\ ⋮\\ {\gamma }_m{C}_1+{\gamma }_m{\varphi }_m{C}_2+{\gamma }_m{\varphi }_m^2{C}_3...+{\gamma }_m{\varphi }_m^{j-1}{C}_j\end{array}\right]=0\text{\hspace{1em}(12)}$$
划简公式为：
$$\left[\begin{array}{c}{C}_1+{\varphi }_1{C}_2+{\varphi }_1^2{C}_3...+{\varphi }_1^{j-1}{C}_j\\ C_1+{\varphi }_2{C}_2+{\varphi }_2^2{C}_3...+{\varphi }_2^{j-1}{C}_j\\ ⋮\\ C_1+{\varphi }_m{C}_2+{\varphi }_m^2{C}_3...+{\varphi }_m^{j-1}{C}_j\end{array}\right]=0\text{\hspace{1em}(13)}$$
如果已经找到了关于C的系数，则根据以上式子可以解关于$\varphi$的值。但是有一个前提条件是接收信号的个数必须小于天线的个数

最终将问题转换为：
$$\left[\begin{array}{cccc}{y}_1[1]& y_2[1]& ...& y_n[1]\\ y_1[2]& y_2[2]& ...& y_n[2]\\ ⋮& & & \\ y_1[m]& y_2[m]& ...& y_n[m]\end{array}\right]\left[\begin{array}{c}{C}_1\\ C_2\\ ⋮\\ C_n\end{array}\right]=\left[\begin{array}{c}0\\ 0\\ ⋮\\ 0\end{array}\right]\text{\hspace{1em}(12)}$$
这和一般的解方程组的方式不一样的地方是这个不是方阵，未知数小于方程个数。
是一个典型的解超定方程组系统（overdetermind system）。

超定方程组求解

TODO ！！！