clear all; clc;
close all
radeg = 180/pi;
derad=1/radeg;
twpi = 2*pi;
Melm= 7;
kelm = 6; %
dd = 0.5; %
d=0:dd:(Melm-1)*dd; %
iwave = 2; % number of DOA
theta = [40 50]; % DOA
n = 1024 %
A=exp(-j*twpi*d.'*sin(theta*derad));%%%% direction matrix
%构造相干信号源
S0=randn(iwave-1,n);
S=[S0(1,:);S0];
X0=A*S;
X=awgn(X0,10,'measured');
Rxxm=X*X'/n
issp=2;
%空间平滑
if issp==1
Rxx=ssp(Rxxm,kelm);
elseif issp==2
Rxx=mssp(Rxxm,kelm);
else
Rxx=Rxxm;
kelm=Melm;
end
[EV,D]=eig(Rxx);
EVA=diag(D)'; [EVA,I]=sort(EVA);
EVA=fliplr(EVA); EV=fliplr(EV(:,I));
for iang = 1:361
angle(iang)=(iang-181)/2;
phim=derad*angle(iang);
a=exp(-j*twpi*d(1:kelm)*sin(phim)).';
L=iwave;
En=EV(:,L+1:kelm);
SP(iang)=(a'*a)/(a'*En*En'*a);
end
SP=abs(SP);
SPmax=max(SP);
SP=10*log10(SP/SPmax);
figure
h=plot(angle,SP);
set(h,'Linewidth',2)
xlabel('angle (degree)')
ylabel('magnitude (dB)')
axis([-90 90 -60 0])
set(gca, 'XTick',[-90:10:90],'YTick',[-60:10:0])
grid on
hold on
legend('平滑MUSIC')
function crs = mssp(cr, K)
% modified spatial smoothing
[M,MM]=size(cr);
N=M-K+1;
J = fliplr(eye(M));
crfb = (cr + J*cr.'*J)/2;
crs = zeros(K,K);
for in =1:N
crs = crs + crfb(in:in+K-1,in:in+K-1)
end
crs = crs / N;
end
%End mssp.m
function crs = ssp(cr, K)
% spatial smoothing
[M,MM]=size(cr);
N=M-K+1;
crs = zeros(K,K);
for in =1:N
crs = crs + cr(in:in+K-1,in:in+K-1)
end
crs = crs / N;
end
DOA.zip_DOA_DOA改进_capon_平滑延迟_改进DOA
版权申诉
114 浏览量
2022-07-13
22:23:28
上传
评论
收藏 4KB ZIP 举报
![avatar](https://profile-avatar.csdnimg.cn/2416af5c19524431b870352d943af459_weixin_42659196.jpg!1)
周楷雯
- 粉丝: 80
- 资源: 1万+
最新资源
- 一段Scratch语言的计算斐波那契数列的步骤流程
- 一段VisualBasic语言的计算斐波那契数列的代码
- RUBY-DEFAULT-CONTAINS
- 基于Go语言写的一段计算斐波那契数列的代码!
- vue-material v1.0.0-beta-10.2
- TbReactions(BaseModel)
- Blumentals Rapid PHP 2025 v18.1.0.264 是一款功能强大的PHP编辑器和开发工具
- do not edit .js files directly - edit src/index.jst
- request-validator.cpython-37
- pygments.lexers.javascript ~~~~~~~~~~~~~~~~~~~~~~~~~~
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈
![feedback](https://img-home.csdnimg.cn/images/20220527035711.png)
![feedback](https://img-home.csdnimg.cn/images/20220527035711.png)
![feedback-tip](https://img-home.csdnimg.cn/images/20220527035111.png)
评论0