P = 30; % 发射功率 sigma_b = 1; % Bob端AWGN的方差 sigma_e = 1; % Eve端AWGN的方差 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的信道参数 Nsamp = 1e6; x = randi([0 1],Nsamp,1) d_AB = 3^0.5; % Alice到Bob的距离 d_AE = 2; % Alice到Eve的距离 eta_AB = 2; % Alice到Bob的衰落系数 eta_AE = 2; % Alice到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 )^2) / (theta * P * d_AE^(-eta_AE) * norm(h_AE )^2 + sigma_e^2); % Eve的信噪比给以上内容添加matlab代码，得到gamma_b和gamma_e随theta变化的曲线图

好的，以下是添加的代码：

P = 30;   % 发射功率
sigma_b = 1;   % Bob端AWGN的方差
sigma_e = 1;   % Eve端AWGN的方差
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的信道参数
Nsamp = 1e6;
x = randi([0 1],Nsamp,1);
d_AB = 3^0.5;   % Alice到Bob的距离
d_AE = 2;   % Alice到Eve的距离
eta_AB = 2;   % Alice到Bob的衰落系数
eta_AE = 2;   % Alice到Eve的衰落系数

theta = linspace(0, 1, 101);
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 )^2) ./ (theta * P * d_AE^(-eta_AE) * norm(h_AE )^2 + sigma_e^2);   % Eve的信噪比

plot(theta, gamma_b, theta, gamma_e);
title('Gamma_B and Gamma_E vs. Theta');
xlabel('Theta');
ylabel('Gamma_B and Gamma_E');
legend('Gamma_B', 'Gamma_E');

这段代码会生成一个图像，横轴是 `theta`，纵轴是 `gamma_b` 和 `gamma_e`，分别表示 Bob 和 Eve 的信噪比随功率分配因子 `theta` 变化的曲线。