matlab 中awgn函数怎么用clcclear allformat longN=200;%%快拍数doa=[20 40]/180*pi; %%信号到达角,w=[pi/4 pi/3]';%%信号频率M=8;%%阵元数P=length(w); %%信号个数,也可以用特征分解的大特征值来决定lambda=150;%波长d=lam
来源：学生作业帮助网 编辑：六六作业网 时间：2024/12/22 15:33:17
matlab 中awgn函数怎么用clcclear allformat longN=200;%%快拍数doa=[20 40]/180*pi; %%信号到达角,w=[pi/4 pi/3]';%%信号频率M=8;%%阵元数P=length(w); %%信号个数,也可以用特征分解的大特征值来决定lambda=150;%波长d=lam
matlab 中awgn函数怎么用
clc
clear all
format long
N=200;%%快拍数
doa=[20 40]/180*pi; %%信号到达角,
w=[pi/4 pi/3]';%%信号频率
M=8;%%阵元数
P=length(w); %%信号个数,也可以用特征分解的大特征值来决定
lambda=150;%波长
d=lambda/2;%阵元间距
snr=15;%%信噪比
%%%导向向量
B=zeros(P,M);
for k=1:P
B(k,:)=exp(-j*2*pi*d*sin(doa(k))/lambda*[0:M-1]);
end
B=B';
%%%导向向量
xx=2*exp(j*(w*[1:N]));
x=B*xx;
%%%%噪声平均 因为matlab产生的噪声不太好,不是严格意义上的白噪声,
%%%平均后结果更好
[pp,ppp]=size(x);
xx=zeros(pp,ppp);
cycle=200;
for kkk=1:cycle
xx=xx+awgn(x,snr);
end
x=xx/cycle;
%%%%噪声平均结束
R=x*x'; %数据协方差矩阵
%针对相干源的时候进行平衡
J=fliplr(eye(M));
R=R+J*conj(R)*J;
%%以下是ESPRIT程序
Rxx=R(1:M-1,1:M-1);%%%M-1维的自相关函数
Rxy=R(1:M-1,2:M);%%%M-1维的互相关函数
b=[zeros(1,M-2);eye(M-2)];
b=[b zeros(M-1,1)];
Cxx=Rxx-min(eig(Rxx))*eye(M-1);
Cxy=Rxy-min(eig(Rxx))*b;
a=eig(Cxx,Cxy);
%找出最接近1的a值其对应的角度即为φ
a1=abs(abs(a)-1);
for i=1:P
[c,d]=min(a1);
a1(d)=1000;
bb(i)=a(d);
a(d)=1000;
end
if P>1
disp('The angles of signals are')
else
disp('The angle of signal is')
end
DOA=asin(angle(bb)/pi)/pi*180
其中的
for kkk=1:cycle
xx=xx+awgn(x,snr);
end
有什么错误吗
matlab 中awgn函数怎么用clcclear allformat longN=200;%%快拍数doa=[20 40]/180*pi; %%信号到达角,w=[pi/4 pi/3]';%%信号频率M=8;%%阵元数P=length(w); %%信号个数,也可以用特征分解的大特征值来决定lambda=150;%波长d=lam
AWGN：在某一信号中加入高斯白噪声
y = awgn(x,SNR) 在信号x中加入高斯白噪声.信噪比SNR以dB为单位.x的强度假定为0dBW.如果x是复数,就加入复噪声.
y = awgn(x,SNR,SIGPOWER) 如果SIGPOWER是数值,则其代表以dBW为单位的信号强度；如果SIGPOWER为'measured',则函数将在加入噪声之前测定信号强度.
y = awgn(x,SNR,SIGPOWER,STATE) 重置RANDN的状态.
y = awgn(…,POWERTYPE) 指定SNR和SIGPOWER的单位.POWERTYPE可以是'dB'或'linear'.如果POWERTYPE是'dB',那么SNR以dB为单位,而SIGPOWER以dBW为单位.如果POWERTYPE是'linear',那么SNR作为比值来度量,而SIGPOWER以瓦特为单位.