多重信号分类MUSIC 算法

基于子空间的多重信号分类算法MUSIC，可估计信号的到达角、离开角

function y=awgn（varargin）
%AWGN Add white Gaussian noise to a signal.
%   Y = AWGN（XSNR） adds white Gaussian noise to X.  The SNR is in dB.
%   The power of X is assumed to be 0 dBW.  If X is complex then 
%   AWGN adds complex noise.
%
%   Y = AWGN（XSNRSIGPOWER） when SIGPOWER is numeric it represents 
%   the signal power in dBW. When SIGPOWER is ‘measured‘ AWGN measures
%   the signal power before adding noise.
%
%   Y = AWGN（XSNRSIGPOWERS） uses S which is a random stream handle to
%   generate random noise samples with RANDN. If S is an integer then
%   resets the state of RANDN to S. The latter usage is obsoleted and may
%   be removed in a future release.  If you want to generate repeateable
%   noise samples then provide the handle of a random stream or use reset
%   method on the default random stream.  Type ‘help RandStream‘ for more
%   information.
%
%   Y = AWGN（XSNRSIGPOWERSTATE） resets the state of RANDN to STATE.
%   This usage is deprecated and may be removed in a future release.
%
%   Y = AWGN（... POWERTYPE） specifies the units of SNR and SIGPOWER.
%   POWERTYPE can be ‘db‘ or ‘linear‘.  If POWERTYPE is ‘db‘ then SNR
%   is measured in dB and SIGPOWER is measured in dBW.  If POWERTYPE is
%   ‘linear‘ then SNR is measured as a ratio and SIGPOWER is measured
%   in Watts.
%
%   Example 1: 
%        % To specify the power of X to be 0 dBW and add noise to produce
%        % an SNR of 10dB use:
%        X = sqrt（2）*sin（0:pi/8:6*pi）;
%        Y = awgn（X100）;
%
%   Example 2: 
%        % To specify the power of X to be 3 Watts and add noise to
%        % produce a linear SNR of 4 use:
%        X = sqrt（2）*sin（0:pi/8:6*pi）;
%        Y = awgn（X43‘linear‘）;
%
%   Example 3: 
%        % To cause AWGN to measure the power of X and add noise to
%        % produce a linear SNR of 4 use:
%        X = sqrt（2）*sin（0:pi/8:6*pi）;
%        Y = awgn（X4‘measured‘‘linear‘）;
%
%   Example 4: 
%        % To specify the power of X to be 0 dBW add noise to produce
%        % an SNR of 10dB and utilize a local random stream use:
%        S = RandStream（‘mt19937ar‘‘seed‘5489）;
%        X = sqrt（2）*sin（0:pi/8:6*pi）;
%        Y = awgn（X100S）;
%
%   Example 5: 
%        % To specify the power of X to be 0 dBW add noise to produce
%        % an SNR of 10dB and produce reproducible results use:
%        reset（RandStream.getDefaultStream）
%        X = sqrt（2）*sin（0:pi/8:6*pi）;
%        Y = awgn（X100S）;
%
%
%   See also WGN RANDN RandStream/RANDN and BSC.

%   Copyright 1996-2008 The MathWorks Inc.
%   $Revision: 1.9.4.7 $  $Date: 2009/01/05 17:45:01 $ 

% --- Initial checks
error（nargchk（25nargin‘struct‘））;

% --- Value set indicators （used for the string flags）
pModeSet    = 0;
measModeSet = 0;

% --- Set default values
sigPower = 0;
pMode    = ‘db‘;
measMode = ‘specify‘;
state    = [];

% --- Placeholder for the signature string
sigSt

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        1152  2013-09-03 21:23  doa_music.m
     文件        6522  2013-09-03 16:19  wgn.m
     文件        8109  2013-09-03 08:17  awgn.m