代码拉取完成,页面将自动刷新
% Data Simulation
clear all;
derad = pi/180; %角度->弧度
M = 4; % 阵元个数
N = 3;
K = 512; % 载波数
path = 2; % 信源数目
theta1 = [150 90]; % 待估计角度aoa
theta2 = [30 90];
c = 3e8;
tof = [6 2.5]/c;
dcps = pi/3;
snr = 10; % 信噪比
ex1 = 0:M-1;
ex2 = 0:N-1;
ex3 = 0:K-1;
d = 0.26; % 阵元间距
d1 = d*ex1;
d2 = d*ex2;
delta_f = 312.5e3;
f = delta_f*ex3;
A1 = exp(-1i*2*pi*d1.'*cos(theta1*derad)); %方向矢量aoa
A2 = exp(-1i*2*pi*d2.'*cos(theta2*derad)); %方向矢量aod
A3 = exp(-1i*2*pi*f.'* tof); % 时延矢量tof
A4 = exp(-1i*dcps);
%信道状态矩阵aoa_aod
H_aoad = A1 * A2';
H_aoad = awgn(H_aoad,snr,'measured'); %将白色高斯噪声添加到信号中
%信道状态矩阵aoa_tof
H_aoatof = A1 * A3';
H_aoatof = awgn(H_aoatof,snr,'measured'); %将白色高斯噪声添加到信号中
%信道状态矩阵aoa_aod_tof
A11 = permute((repmat(A1,1,1,N)),[1,3,2]);
A22 = permute((repmat(A2,1,1,M)),[3,1,2]);
H_adtof = reshape(reshape(A11 .* A22,M*N,path) * A3',M,N,K);
H_adtof = awgn(H_adtof,snr,'measured'); %将白色高斯噪声添加到信号中
%信道状态矩阵aoa_aod_tof_dcps
H_adtp = H_adtof * dcps;
此处可能存在不合适展示的内容,页面不予展示。您可通过相关编辑功能自查并修改。
如您确认内容无涉及 不当用语 / 纯广告导流 / 暴力 / 低俗色情 / 侵权 / 盗版 / 虚假 / 无价值内容或违法国家有关法律法规的内容,可点击提交进行申诉,我们将尽快为您处理。