信号DCT字典稀疏表示

信号DCT字典稀疏表示，匹配追踪算法（MP）搜索稀疏系数，重构误差达到精度停止循环。

%稀疏表示方法分类相控阵的故障
%产生测量数据
clear;close all;
N=6;
B=8;
lambda=0.2;
d=0.51*lambda;
r0=fix（2*（N-1）^2*d^2/lambda*10）/10;
r=5;%每一类的样本数，训练样本
r1=3;%测试样本
[y_temptheta_b]=nearfield_test_u1（Nr0Bzeros（1N）00）;
Y=zeros（length（y_temp）r*（2^N-1））;
T=zeros（1r*（2^N-1））;
Y_t=zeros（length（y_temp）r1*（2^N-1））;
T_t=zeros（1r1*（2^N-1））;
vec=zeros（N2^N-1）;%产生损坏情况的（01）向量
%产生训练样本
for i=1:2^N-1
    for k=1:r
        vec（:i）=de2bi（2^N-iN）‘;
        [Y（:（i-1）*r+k）theta]=nearfield_test_u1（Nr0Bvec（:i）0.010.01）;
        T（1（i-1）*r+k）=i;
    end
end
%产生测试样本
for i=1:2^N-1
    for k=1:r1
        vec（:i）=de2bi（2^N-iN）‘;
        [Y_t（:（i-1）*r1+k）theta]=nearfield_test_u1（Nr0Bvec（:i）0.010.01）;
        T_t（1（i-1）*r1+k）=i;
    end
end

%DCT字典
K=256;
sigma = 25; 
C = 1.15;
Pn=ceil（sqrt（K））;
bb=ceil（sqrt（size（Y1）））;
DCT=zeros（bbPn）;
for k=0:1:Pn-1
    V=cos（[0:1:bb-1]‘*k*pi/Pn）;
    if k>0
        V=V-mean（V）;
    end
    DCT（:k+1）=V/norm（V）;
end
DCT=kron（DCTDCT）;
Y=[Y;zeros（size（DCT1）-size（Y1）size（Y2））];
Coefs = OMPerr（DCTY0.01）;

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

     文件       1368  2013-08-01 16:20  nearfield_test_u1.m

     文件       1154  2013-08-06 10:49  nearfield_sparse.m

     文件       1083  2006-12-11 15:08  OMPerr.m

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

                 3605                    3