matlab程序没报错,但是出来的图片是一条直线

Question

clc;
BW = 1e6;
Pcir =188*1e-6;
E0=1e-4;
pe=1e-3;
G0 =-70; G0 = 10^(G0/10);
N0 = -140; N0= (10^(N0/10))*1e-3;
a=0.4;c1=0.2;c2=1.6;c3=1;c4=1;
d1=1000;d2=5000;
h1 = G0*d1.^(-3);
h2 = G0*(d2-d1).^(-3);
G1 = h1/N0;
G2 = h2/N0;
gt = 1/sqrt(G1)+1/sqrt(G2);

M=[5:5:50];
for i=1:length(M)
gamma = [(2.^(c3*M-c4))/c2]*log(c1/pe);
EE(M)=(a*BW*M)/((gt.^2)*gamma+4*a*Pcir+a*BW*E0*(2.^M));
plot(M,EE(M));
end

Answer 1

clc;

clear,close all

BW = 1e6; 

Pcir =188*1e-6; 

E0=1e-4;

pe=1e-3;

G0 =-70;  G0 = 10^(G0/10);

N0 = -140;  N0= (10^(N0/10))*1e-3;

a=0.4;c1=0.2;c2=1.6;c3=1;c4=1;

d1=1000;d2=5000;

h1 = G0*d1.^(-3);

h2 = G0*(d2-d1).^(-3);

G1 = h1/N0;

G2 = h2/N0;

gt = 1/sqrt(G1)+1/sqrt(G2);

M=[5:5:50];

for i=1:length(M)

    gamma = ((2.^(c3*M(i)-c4))/c2)*log(c1/pe);

EE(i)=(a*BW*M)/((gt.^2)*gamma+4*a*Pcir+a*BW*E0*(2.^M));

% plot(M,EE(M)),hold on;

end

figure,plot(M,EE)

追问

EE(i)=(a*BW*M（i）)/((gt.^2)*gamma+4*a*Pcir+a*BW*E0*(2.^M（i）))是否要这样修改

Answer 2

M=[5:5:50];  

for i=1:length(M)

    gamma = [(2.^(c3.*M(i)-c4))/c2]*log(c1/pe);

    EE=(a*BW.*M(i))/((gt^2).*gamma+4*a*Pcir+a*BW*E0.*(2.^M(i)));

    EEdata(i) = EE;

end

plot(M,EEdata);

最后一段代码还是修改一下吧。。我出的图

也不知道对不对，最后一个小建议，在你进行循环计算的时候，最好建立一个空的矩阵，将你每次循环算的数据存到新的矩阵里面，这样可能会好一些

追问

理论分析应该是一个抛物线的形式，有一个最优点能达到最大值