三正弦叠加高斯白噪声信号源进行频率估计 （Pisarenko_music_esprit）

分别用Pisarenko谐波估计法、MUSIC算法和ESPRIT算法对三正弦叠加高斯白噪声信号源进行频率估计，附件中有代码说明，编译环境为matlab。

% ESPRIT
clc
clear all
close all 

Nr = 7;
ldata = 128; % 数据样本长度32/64/128
w = 2 * pi * [0.05 0.40 0.42]; % 频率
N_Var = 0.32; 
Z = zeros（NrNr）;
Z（2:Nr1:Nr-1） = blkdiag（111111）;%1111）;%1111）;
I = blkdiag（1111111）;%1111）;%1111）;
valuediag = zeros（203）;
frq_mean = zeros（N_Var3）;
frq_var = zeros（N_Var3）;
k = 1
    for i  = 1:20
    smp = fix（100*rand）;%随机初始相位
    Nise = randn（1ldata） * sqrt（0.32+（k-1）*.01）;
    n = smp:smp+ldata-1;
    X = 2 * cos（ w（1）*n ）+3 * cos（ w（2）*n ）  + 1.2 * cos（ w（3）*n ）;
    Y = （X + Nise）‘;

    Rline = xcorr（Y‘）/ldata;
    Rxx = toeplitz（Rline（ldata:ldata+Nr-1））;

    Rxy1 = toeplitz（Rline（ldata:ldata+Nr   ））;
    Rxy = Rxy1（1:Nr2:Nr+1）;

    [eigvcteigval] = eig（Rxx）;
    lamdamin = eigval（11）;
    C

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1458  2011-01-19 01:22  代码和说明\Esprit.m

     文件       2220  2011-01-19 10:16  代码和说明\Esprit_eig.m

     文件       2026  2011-01-19 10:38  代码和说明\ESprit_tls.m

     文件       1340  2011-01-18 23:58  代码和说明\music1.m

     文件       2164  2011-01-19 09:55  代码和说明\music_evaluate.m

     文件        876  2011-01-18 21:41  代码和说明\pisarenko.m

     文件       1642  2011-01-19 01:47  代码和说明\pisko_evaluate.m

     文件       1431  2011-01-19 10:50  代码和说明\Tls_Esprit.m

     文件        684  2011-09-30 14:27  代码和说明\代码说明.txt

     目录          0  2011-09-30 14:25  代码和说明

----------- ---------  ---------- -----  ----

                13841                    10