把它用matlab表示出来

好的，那么我将给出一些Matlab代码，用于实现上述内容中的一些公式和计算。请注意，以下代码仅供参考，实现细节和具体参数需要根据具体情况进行调整和修改。

首先，我们可以定义一些常量和参数，如发射功率P、功率分配因子ϑ、AWGN的方差σ2等：

P = 1;   % 发射功率
theta = 0.5;   % 功率分配因子
sigma_b = 0.1;   % Bob端AWGN的方差
sigma_e = 0.2;   % Eve端AWGN的方差

然后，我们可以生成一些随机信号和噪声，如期望信号预编码向量u、有用信号s、人工噪声预编码矩阵W和人工噪声向量z：

N_A = 4;   % Alice的天线数
N_B = 1;   % Bob的天线数
N_E = 1;   % Eve的天线数
h_AB = randn(N_A, N_B);   % Alice到Bob的信道参数
h_AE = randn(N_A, N_E);   % Alice到Eve的信道参数
u = randn(N_A, 1);   % 期望信号预编码向量
s = sqrt((1 - theta) * P) * randn(N_A, 1);   % 有用信号
W = randn(N_A, N_A);   % 人工噪声预编码矩阵
z = sqrt(theta * P) * randn(N_A, 1);   % 人工噪声向量

接下来，我们可以根据公式计算Bob和Eve接收到的信号y_b和y_e，以及对应的信噪比γ_b和γ_e：

d_AB = 1;   % Alice到Bob的距离
d_AE = 2;   % Alice到Eve的距离
eta_AB = 3;   % Alice到Bob的衰落系数
eta_AE = 4;   % Alice到Eve的衰落系数
y_b = sqrt((1 - theta) * P * d_AB^(-eta_AB)) * h_AB * (u + W * z) + sigma_b * randn(1);   % Bob的接收信号
y_e = sqrt(theta * P * d_AE^(-eta_AE)) * h_AE * (u + W * z) + sigma_e * randn(1);   % Eve的接收信号
gamma_b = ((1 - theta) * P * d_AB^(-eta_AB) * norm(h_AB)^2) / sigma_b^2;   % Bob的信噪比
gamma_e = (theta * P * d_AE^(-eta_AE) * norm(h_AE * u)^2) / (theta * P * d_AE^(-eta_AE) * norm(h_AE * W)^2 + sigma_e^2);   % Eve的信噪比

以上代码实现了该通信系统的一些基本计算，但是具体的实现还需要根据实际情况进行调整和修改。