sum_test.m

来自「MATSNL is a package of MATLAB M-files fo」· M 代码 · 共 16 行

clear;
Q=5;
sum_sim=0;
p=0.9;
for i=1:Q
%     sum_sim=sum_sim+i*p*(1-p)^(i-1);
    sum_sim=sum_sim+i*p^(i);
    sum_sim_log(i)=sum_sim;
end

sum_thy=(p/(1-p))*(1-Q*p^Q+(p/(1-p))*(1-p^(Q-1)));

disp(['sum_sim: ',num2str(sum_sim)]);
disp(['sum_thy: ',num2str(sum_thy)]);

plot(sum_sim_log);