作业帮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/11/05 11:56:19
![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](/uploads/image/z/8325536-32-6.jpg?t=matlab+%E4%B8%ADawgn%E5%87%BD%E6%95%B0%E6%80%8E%E4%B9%88%E7%94%A8clcclear+allformat+longN%3D200%3B%25%25%E5%BF%AB%E6%8B%8D%E6%95%B0doa%3D%5B20+40%5D%2F180%2Api%3B+%25%25%E4%BF%A1%E5%8F%B7%E5%88%B0%E8%BE%BE%E8%A7%92%2Cw%3D%5Bpi%2F4+pi%2F3%5D%27%3B%25%25%E4%BF%A1%E5%8F%B7%E9%A2%91%E7%8E%87M%3D8%3B%25%25%E9%98%B5%E5%85%83%E6%95%B0P%3Dlength%28w%29%3B+%25%25%E4%BF%A1%E5%8F%B7%E4%B8%AA%E6%95%B0%2C%E4%B9%9F%E5%8F%AF%E4%BB%A5%E7%94%A8%E7%89%B9%E5%BE%81%E5%88%86%E8%A7%A3%E7%9A%84%E5%A4%A7%E7%89%B9%E5%BE%81%E5%80%BC%E6%9D%A5%E5%86%B3%E5%AE%9Alambda%3D150%3B%25%E6%B3%A2%E9%95%BFd%3Dlam)
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以瓦特为单位.