matrix_x.m

来自「此程序是在MATLAB上实现用于短波定位的MUSIC算法仿真」· M 代码 · 共 34 行

function [ x ] = matrix_X( M , D , r , t , w0 , s0 )

%  &*&*&*&*&*&*&*&*&*&*&*&*&*&*&*&*&*&*&*&*&*&*&*&*&*
%   M--天线根数（均匀圆阵）
%   D--信号个数
%   r--圆阵半径
%   t--采样时刻


c=3*10^8;             %光速

st = 0 : 2*pi/D : 2*pi-2*pi/D; %方位角360度平分成D个

fi = 0 : pi/2/D : pi/2-pi/2/D; %仰角90度平分成D个

tao = zeros( M , D );          

k = 1 : M;

for i = 1 : D
    tao(:,i) = r/c * cos( 2*pi*(k-1)'/M - st(i) ) * cos( fi(i) );% 生成一个M*D的时延矢量矩阵
end

A=exp( (-tao)*w0*1j );        % 生成一个M*D的方向矢量矩阵

s = s0 : 10000 : s0+10000*(D-1);  %产生D个基带信号

s = cos(s*t)*exp(w0*t*1j);    %产生D个基带信号的调制信号

nis = randn(M,1)*0.1;

nis = hilbert(nis)/2^0.5;     %产生M个噪声

x = A * s' + nis;             %模拟产生一个时刻的M*D矩阵的信号