資源簡介
信號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
- 上一篇:anysend發包工具
- 下一篇:計算機考研方向
評論
共有 條評論