⛄一、获取代码方式
获取代码方式1:
完整代码已上传我的资源:【MIMO通信】基于matlab MIMO-WINNER模型【含Matlab源码 4948期】
点击上面蓝色字体,直接付费下载,即可。
获取代码方式2:
付费专栏信号处理(Matlab)
备注:
点击上面蓝色字体付费专栏图像处理(Matlab),扫描上面二维码,付费299.9元订阅海神之光博客付费专栏,凭支付凭证,私信博主,可免费获得5份本博客上传CSDN资源代码(有效期为订阅日起,三天内有效);
点击CSDN资源下载链接:5份本博客上传CSDN资源代码
⛄二、部分源代码
%输入: scene: 场景选择 1 为 UMi, 2 为 UMa
% d2D: 小区半径
% U: 移动端天线数
% S: 基站端天线数
%输出: H: 信道矩阵
function H = channel_3DMIMO(scene,d2D,U,S)
clusterE = pi/2+(2rand-1)45/180pi; % phi_zod
clusterD = pi/2+(2rand-1)60/180pi; % phi_aod
[sigma_ASD,sigma_ZSD,sigma_DS,sigma_SF,sigma_ASA,sigma_ZSA,N_cluster,N_ray,…
c_ASD,m_ZSD,Dsp,Pcs,m_ZSD_offset,c_ZSA,c_ASA,m_xpr,s_xpr]=generate_para(scene,d2D);
[pha_aod_n_m,pha_zod_n_m,pha_aoa_n_m,pha_zoa_n_m,p]=angles(clusterE,clusterD,…
sigma_ASD,sigma_ZSD,sigma_ASA,sigma_ZSA,sigma_DS,N_cluster,N_ray,c_ASD,c_ASA,c_ZSA,m_ZSD,Dsp,Pcs,m_ZSD_offset);
% 产生天线交叉极化比
m_XPR = m_xpr;
s_XPR = s_xpr;
XPR = 10.^(normrnd(m_XPR,s_XPR,U,S)/10);
alpha = pi/4;
% 产生四种极化相位
phi_VV = pi.*rand(N_ray,N_cluster);
phi_VH = pi.*rand(N_ray,N_cluster);
phi_HV = pi.*rand(N_ray,N_cluster);
phi_HH = pi.*rand(N_ray,N_cluster);
% 产生信道系数
t = [1:100];
CarrierFrequency = 2e9;
wavelength = 3e8/CarrierFrequency;
ds = 0.5wavelength;
du = 0.5wavelength;
k_CONST = 2pi/wavelength; % 波矢
v = 10; % m/s
theta_v = pirand;
% 信道矩阵初始化
H = zeros(U,S,length(t));
h_m(t) = zeros(size(t));
h_n(t) = zeros(size(t));
h(t) = zeros(size(t));
for u=1:U
for s=1:S
for n=1:N_cluster
for m=1:N_ray
h_m(t) = sqrt(p(n)/N_ray)[cos(alpha)sin(pha_zoa_n_m(m,n))+sin(alpha)sin(pha_zoa_n_m(m,n))cos(pha_aoa_n_m(m,n)),cos(alpha)cos(pha_zoa_n_m(m,n))] …
[exp(1iphi_VV(m,n)),sqrt(XPR(u,s))exp(1iphi_VH(m,n));sqrt(XPR(u,s))exp(1iphi_HV(m,n)),exp(1iphi_HH(m,n))]…
[cos(alpha)sin(pha_zod_n_m(m,n))+sin(alpha)sin(pha_zod_n_m(m,n))cos(pha_aod_n_m(m,n));cos(alpha)cos(pha_zod_n_m(m,n))]…
exp(1i(k_CONSTdssin(pha_aod_n_m(m,n))sin(pha_zod_n_m(m,n))))…
exp(1i(k_CONSTdusin(pha_aoa_n_m(m,n))sin(pha_zoa_n_m(m,n)))).…
exp(1i*k_CONST)vsin(pha_zoa_n_m(m,n)*cos(pha_aoa_n_m(m,n)-theta_v).*t);
h_n(t) = h_n(t)+h_m(t);
end
h(t) = h(t)+h_n(t);
end
H(u,s,:) = h(t);
end
end
end
⛄三、运行结果
⛄四、matlab版本及参考文献
1 matlab版本
2014a
2 参考文献
[1] 沈再阳.精通MATLAB信号处理[M].清华大学出版社,2015.
[2]高宝建,彭进业,王琳,潘建寿.信号与系统——使用MATLAB分析与实现[M].清华大学出版社,2020.
[3]王文光,魏少明,任欣.信号处理与系统分析的MATLAB实现[M].电子工业出版社,2018.
3 备注
简介此部分摘自互联网,仅供参考,若侵权,联系删除