加入 Gitee
与超过 1200万 开发者一起发现、参与优秀开源项目,私有仓库也完全免费 :)
免费加入
文件
该仓库未声明开源许可证文件(LICENSE),使用请关注具体项目描述及其代码上游依赖。
克隆/下载
data_sim_test.m 1.22 KB
一键复制 编辑 原始数据 按行查看 历史
Xixi22 提交于 2021-03-09 15:56 . 新版本
% 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;
Loading...
马建仓 AI 助手
尝试更多
代码解读
代码找茬
代码优化