clear?all;
close?all;
clc;
lambda=60;
dd=lambda/2;??????%陣元間距離，取為入射波長的一半
Kp=80;????????????%采樣快拍數
Nnum=9;?%%陣列1陣元數量
fs=20*10^3;

theta=[626];
AA=[233];
SignalNum=length（theta）;
Fc=[2*10^35*10^32*10^3];?????%入射信號頻率
SNR=10;??
Aratio=sqrt（10^（SNR/10））;????%信號幅度與噪聲幅度比值，并假設信號幅度為1
%?thetatest=（0*pi/180:1*pi/180:90*pi/180）;???%theta角度搜索范圍
thetatest=[-60:1:60]*pi/180;
thetanum=length（thetatest）;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
SignalVector=zeros（SignalNumKp）;
Xt=zeros（NnumKp）;
????M=Nnum;

????T_Vector=（1:Kp）/fs;
d=dd*（0:M-1）;
w=2*pi/lambda*sin（thetatest）;
A=exp（-1j*d‘*w）;

%?%?%??for?k2=1:thetanum
%?%?%?????for?k1=1:M
%?%?%?????????Att（k1）=exp（-j*（k1-1）*2*pi*dd*sin（thetatest（k2））/lambda）;
%?%?%?????????A（k1k2）?=Att（k1）;
%?%?%?????end
%?%?%???end

%%構造觀測矩陣Y
for?k2=1:SignalNum
????for?k1=1:Nnum??????????????????????????
????????At（k1）=exp（-j*（k1-1）*2*pi*dd*sin（theta（k2）*pi/180）/lambda）;
????????Atemp（k1k2）=At（k1）;???
????end?????
end
%%%構造信號矩陣和噪聲矩陣
for?k1=1:SignalNum
????SignalVector（k1:）=exp（j*2*pi*Fc（k1）.*T_Vector）;??%信號
end
Xtt=Atemp*SignalVector;

%NoiseVector=sqrt（0.5）*（randn（NnumK）+j*randn（NnumK））;
%--------------------加噪聲方法---1
for?kk=1:Nnum
???Xt（kk:）=awgn（Xtt（kk:）SNR‘measured‘）;
end
%-------------------加噪聲方法---2
%?????for?k1=1:Nnum
%?????????NoiseVector（k1:）=SignalNum/Aratio.*noisecg（Kp）‘;
%?????????Xt（k1:）=Xtt（k1:）+NoiseVector（k1:）;
%?????end?

?Rx=（Xt*Xt‘）./Kp;
?

%?%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
[V?D]?=?eig（Rx）;????????????????%???X*V?=?V*D?
DD?????=?diag（D）;???????????????%?對角陣變矢量
sigm_n=min（DD）;????????%最小特征值^作為?的估計

I=eye（M）;??
W12=sqrt（Kp）*（kron（（Rx^（-0.5））.‘Rx^（-0.5）））;
Y=Rx-sigm_n*I;
y=W12*Y（:）;

Fai=W12*kron（IA）;????

II=ones（1thetanum）;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%???
belta=sqrt（chi2inv（0.99999M*M））

%?%?%?cvx_begin
%?%?%?????variable?B（thetanumM）?
%?%?%??????variable?b（thetanum1）;????
%?%?%???????minimize（norm（b1））
%?%?%?????%?minimize?（norm（sqrt（sum（（B.^2）‘））1））
%?%?%????subject?to????
%?%?%????for?kk=1:1:thetanum
%?%?%??????????norm（B（kk:））<=b（kk1）;
%?%?%????end???
%?%?%??????temp0=Rx-A*B-sigm_n*I;
%?%?%??????temp=W12*vec（temp0）;
%?%?%???????norm（temp）<=belta?;
%?%?%?%????norm（y-Fai*B（:））<=23?;
%?%?%??cvx_end
%?%?%?
%?%?%??figure（）
%?%?%??plot（0:1:90b）;
%%%%%%%%%%%%%%%%%%%%%%----BW*vec（R-AB-I）---%%%%%%%%%%%%%%%%%%%%??
?belta=sqrt（chi2inv（0.9999M*M））;
??AP=[];
????for?i=1:1:M
????????AP=blkdiag（APA）;
????end
????
%?I2=eye（thetanum）;
%??APP=[];
%?cvx_begin
%?????variable?B（thetanum*M）?complex;
%??????variable?b（thetanum）;
%????
%???????minimize（norm（b1））
%?????%?minimize?（norm（sqrt（sum（（B.^2）‘））1））
%????subject?to??
%????????for?i=1:1:thetanum
%?????????for?j=1:1:M
%?????????????APP=blkdiag（APPI2（i:））;
%?????????end
%?????????b（i）>=norm（APP*B）;
%?????????APP=[];
%?????end
%?????
%?????????temp0=

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????4277??2020-09-14?14:43??L1_SRACV_Sparse_CVX.m