什么阻止我的代码生成正确的结果
5视图(30天)
显示旧的评论
你好我在matlab编写以下方程:
情节应该类似这样:
我的代码是打击,我uncluded代码绘制仿真结果(下限),所以我们可以使用它作为参考
%定义参数
M_array = 1:50;%的BS天线
SNR0_dB = 0;%的信噪比期望的问题(dB)
SNR0 = 10 ^ (SNR0_dB / 10);%的信噪比期望的问题(在线性范围内)
beta_bar_dB = -10;% Inter-cell干扰强度(dB)
beta_bar = 10 ^ (beta_bar_dB / 10);% Inter-cell干扰强度(在线性范围内)
dH = 1/2;%为齿龈半波间距
π= 180;
num_realizations = 5000;
%计算SE_NLOS_LOWER_BOUND
SE_NLOS_LOWER_BOUND = log2 (1 + ((M_array-1)。/ (beta_bar + (1 / SNR0))));
%计算每个SE_LOS M和问题的角度
SE_LOS = 0(1、长度(M_array));
为i = 1:长度(M_array)
G = 0 (1、num_realizations);
为j = 1: num_realizations
Q = 2 *π* randn ();
W = 2 *π* randn ();
如果信德(Q) ~ =信德(W)
G (j) =(信德(π* dH * M_array(我)*(信德(Q)信德(W)))) ^ 2 / (M_array(我)*(信德(π* dH *(信德(Q)信德(W)))) ^ 2);
其他的
G (j) = M_array(我);
结束
结束
SE_LOS (i) =意味着(log2 (1 + M_array(我)/ (beta_bar *意味着(G) + 1 / SNR0)));
结束
%绘制结果
网格在;持有在;
情节(M_array SE_LOS,“红色”);
情节(M_array SE_NLOS_LOWER_BOUND,“蓝”);
%,仿真结果
包含(“废话天线(M)的数量”);
ylabel (的平均SE(比特/秒/ Hz) ');
传奇(“洛杉矶”,仿真结果下界的,“位置”,“西北”);
阴谋的结果是:
我不理解为什么《曲线下降NLOS_LOWER_BOUND下曲线,图中所示的im试图复制,洛杉矶曲线应该大于NLOS_LOWER_BOUND贯穿,也在《crukidness曲线M值更高但我相信它可以消除通过增加的数量实现(我只使用了5000,因为更高的值花太多时间来计算),我也有M_array = 1:50代替1:10 0来节约时间,因为需要太长时间来计算。
我想知道有什么错误在我的代码,当然非常感谢所有的帮助
