Question

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
有什么错误吗

Answer 1

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以瓦特为单位。

追问

我的有什么问题吗
我觉得没错啊