求高手翻译一下这个Matlab程序的意思,要详细一点,谢谢了 10

clearallcloseallderad=pi/180;radeg=180/pi;twpi=2*pi;kelm=20;%dd=0.5;%d=-(kelm-1)/2*d... clear all
close all
derad = pi/180;
radeg = 180/pi;
twpi = 2*pi;
kelm = 20; %
dd = 0.5; % 
d=-(kelm-1)/2*dd:dd:(kelm-1)/2*dd; % 
iwave = 3; % number of DOA
theta1 = [10 20 30];
theta2 = [15 25 35];% DOA
snr=10; % input SNR (dB)
n=100; %
A0=exp(j*twpi*d.'*(sin(theta1*derad).*cos(theta2*derad)))/sqrt(kelm);
A1=exp(j*twpi*d.'*(sin(theta1*derad).*sin(theta2*derad)))/sqrt(kelm);%%%% direction matrix
S=randn(iwave,n);
X0=[];
for im=1:kelm
X0=[X0;A0*diag(A1(im,:))*S];
end
X=awgn(X0,snr,'measured');
L=iwave;
J1=eye(kelm-1,kelm);
J2=flipud(fliplr(J1));
Q=qq(kelm);
Y=kron(Q',Q')*X;
Q0=qq(kelm-1);
K1=real(Q0'*J2*Q);
K2=imag(Q0'*J2*Q);
I=eye(kelm);
Ku1=kron(I,K1);
Ku2=kron(I,K2);
Kv1=kron(K1,I);
Kv2=kron(K2,I);
E=[real(Y),imag(Y)];
Ey=E*E'/n;
[V,D]=eig(Ey);
EVAs =diag(D).';
[EVAs,I0] = sort(EVAs);
EVAs=fliplr(EVAs);
EVs=fliplr(V(:,I0));
Es=EVs(:,1:L);
fiu=pinv(Ku1*Es)*Ku2*Es;
fiv=pinv(Kv1*Es)*Kv2*Es;
F=fiu+j*fiv;
[VV,DD]=eig(F);
EVA = diag(DD).';
u=2*atan(real(EVA))/pi;
v=2*atan(imag(EVA))/pi;
theta10=asin(sqrt(u.^2+v.^2))*radeg
theta20=atan(v./u)*radeg
展开
 我来答
静默澜世
2013-06-07
知道答主
回答量:10
采纳率:0%
帮助的人:6.2万
展开全部
它是一种高效的工程计算语言。
已赞过 已踩过<
你对这个回答的评价是?
评论 收起
navylq
2013-06-07 · TA获得超过266个赞
知道答主
回答量:25
采纳率:0%
帮助的人:24.9万
展开全部
你懒得可以!!!
已赞过 已踩过<
你对这个回答的评价是?
评论 收起
推荐律师服务: 若未解决您的问题,请您详细描述您的问题,通过百度律临进行免费专业咨询

为你推荐:

下载百度知道APP,抢鲜体验
使用百度知道APP,立即抢鲜体验。你的手机镜头里或许有别人想知道的答案。
扫描二维码下载
×

类别

我们会通过消息、邮箱等方式尽快将举报结果通知您。

说明

0/200

提交
取消

辅 助

模 式