SparseDOA.zip_DOA L1-SVD_DOA估计_L1——SVD DOA_l1 svd doa_l1_svd

T=300; M=8; % …………………产生信号……………… nsource=3; % s=zeros(nsource,T); % amp=[3 3 3]; % for i=1:nsource % s(i,:)=rpiid(T,'exp') * amp(i); % end % s=abs(s); % s(2,:)= s(1,:); % s(3,:)= s(1,:); t=0:1/T:(1-1/T); s1=abs(cos(2*pi*15*t+pi/3)); s2=abs(cos(2*pi*25*t+pi/17)); s3=abs(cos(2*pi*20*t+pi/11)); % % s4=abs(cos(2*pi*45*t+pi/17)); % % s5=abs(cos(2*pi*55*t+pi/8)); % % s6=abs(cos(2*pi*65*t+pi/13)); % % s7=abs(cos(2*pi*75*t+pi/19)); % % % s8=abs(cos(2*pi*90*t+pi/23)); s=3*[s1;s2;s3]; % % % sigmae=10^((10*log10(9)-10)/20); % % noise=sigmae*randn(8,5000); % % noise=noise-mean(noise,2)*ones(1,5000); % % s=2*ones(nsource,T); % % s=awgn(s,5); % % s=s-mean(s,2)*ones(1,T); % ……………………产生传输矩阵A…………………… bearing=[10 30 60]; bearing = bearing * pi / 180; % degrees to radians A=exp(sqrt(-1)*(0:M-1)'*pi*sin(bearing)); % % ……………………产生阵列观测矩阵…………………… y=A*s; % y=awgn(y,10); y=awgn(y,50,'measured',1234,'dB'); % y=y+noise(:,1000:999+T); % y1=y.^2; % ……………………产生过完备Ac…………………… % t=48:0.1:52; % o=-89:1:90; o=1:90; N=length(o); theta=o*pi/180; m=0:(M-1); Ac=exp(sqrt(-1)*m'*pi*sin(theta)); Ac_r=real(Ac); % Ac_i=imag(Ac); % y_r=real(y); % y_i=imag(y); % y1_r=real(y1); % y1_i=imag(y1); AA=Ac; Ac_n=abs(conj(AA')*Ac).^2; z=conj(AA')*y; znew=abs(z).^2; %………………求解稀疏源…………………… % s_re=pinv(Ac_r)*mean(y1_r,2); p=0.1; r=10^(-3); b=zeros(1,N); H=zeros(N,N); s_re=ones(N,1); C=2*conj(Ac_r')*Ac_r; s0=zeros(N,1); % S=zeros(N,100); % k=1; %--------------Lp------------------------ % while (norm(s_re(:,1)-s0)/norm(s_re(:,1)))>0.0001 % s0=s_re(:,1); % for i=1:N % b(i)=p/(((abs(s_re(i,1)))^2+r)^(1-p/2)); % end % B=diag(b); % H=C+B*1; % s_re(:,1)=pinv(H)*2*conj(Ac_r')*real(y); % S(:,k)=s_re; % k=k+1; % end % plot(10*log10(abs(s_re(1:90))/max(abs(s_re(1:90))))); %---------------SVD---------------------- [S V D]=svd(y); % yy=S*V*[ones(1,3) zeros(1,297)]'; % cvx_begin % variable sr(90); % minimize( norm(Ac*sr-yy)+1*norm(sr,1) ); % cvx_end cvx_begin variable sr(90); minimize( norm(sr,1)); subject to norm(Ac_n*sr-mean(znew,2))<=0.7; cvx_end