Incorrect values when calculating the optimum value using while loop
1 vue (au cours des 30 derniers jours)
Afficher commentaires plus anciens
bassant tolba
le 20 Nov 2022
Commenté : bassant tolba
le 21 Nov 2022
Dear all,
Please I'm trying to get the optimum power allocation Pk using while loop.. but it seems that I wrongly writing this algorithm and all values are equal to 1 without making alot of iterations... please can any one help me to correct the code.. I appreciate your help.
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/1200813/image.png)
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/1200863/image.png)
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/1200818/image.png)
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/1200823/image.png)
This is my code:
%parameters
TN=10
lamda=[5.3899,2.8511,1.0846,0.4003,2.8364,0.0988,0.0470,0.3417,0.5460,0.5992];
segma_squared=10^-9; %variance
tau=0.1;
Pt=1;
beta= [0.3414,0.1707,0.1138,0.0854,0.0683,0.0569,0.0488,0.0427,0.0379,0.0341];
epsilon=0.001
%%%%%%%%%%%starting of algorithm 1%%%%%%%%%%%%%%%
max_iterations=1000;
power_allocation_initial=[5.3899,4.435,1.1024,0.4004,2.854,0.109988,0.0470,0.3417,0.3450,0.6592];
for i=2:max_iterations
for k=1:TN
gamma_infinity(i,k)=(lamda(k)*(1-tau^2)*power_allocation_initial(k)*Pt)/(segma_squared); % calculation initial value of gamma_inf(0)
w_k(k)=log2(power_allocation_initial(k));
while (gamma_infinity(i,k)-gamma_infinity(i-1,k) >= epsilon)
i=i+1;
gamma_infinity(i,k)=gamma_infinity(i-1,k);
k_(k)=(gamma_infinity(i,k))/(1+gamma_infinity(i,k));
v_k(k)=log2(1+gamma_infinity(i,k))-((gamma_infinity(i,k))/(1+gamma_infinity(i,k)))*log2(gamma_infinity(i,k));
eita_wk(i,k)=-k_(k) *w_k(k)-k_(k)*log2(lamda(k)*(1-tau^2))*Pt + 2*k_(k)*log2(sqrt(segma_squared))-v_k(k);
w(i,k)=beta(k)*eita_wk(i,k);
ww(i)=sum(w(i,k));
[p_opt(k),fval] = linprog(ww(i),[],[],[],[],0,1);
end
end
end
6 commentaires
Torsten
le 20 Nov 2022
Modifié(e) : Torsten
le 20 Nov 2022
Don't you see from the description of the algorithm that the optimization problem has to be solved within the while-loop and not outside ?
And now that you search for a vector p of length 10, your upper and lower bounds also have to be vectors of length 10. But you set them to 0 and 1 - so just one value.
And what is the for-loop good for ? Isn't it a superfluous duplicate of the while-loop ?
Réponse acceptée
Torsten
le 20 Nov 2022
Modifié(e) : Torsten
le 21 Nov 2022
With regards to for loop, I used it to be able to do j-1 as a previous step
But the steps (variable i) are done in the while loop - so no further for loop is needed.
I don't understand the underlying problem. Here is a suggestion for code as I understand the algorithm (although problem (25) is quite strange since the extra terms (2nd and 3rd) in the definition of Gamma don't depend on the optimization variable omega):
TN=10;
lamda=[5.3899,2.8511,1.0846,0.4003,2.8364,0.0988,0.0470,0.3417,0.5460,0.5992];
segma_squared=10^-9; %variance
tau=0.1;
Pt=1;
beta= [0.3414,0.1707,0.1138,0.0854,0.0683,0.0569,0.0488,0.0427,0.0379,0.0341];
epsilon=0.001;
power_allocation_initial=[5.3899,4.435,1.1024,0.4004,2.854,0.109988,0.0470,0.3417,0.3450,0.6592];
gamma_infinity_im1=(lamda.*(1-tau^2).*power_allocation_initial*Pt)/(segma_squared); % calculation initial value of gamma_inf(0)
imax = 20;
i = 0;
error = 1;
while error > epsilon && i < imax
i = i+1;
kappa = gamma_infinity_im1./(1+gamma_infinity_im1);
v = log2(1+gamma_infinity_im1)-kappa.*log2(gamma_infinity_im1);
f = -beta.*kappa; % Don't know if there are other terms in gamma(omega) that depend on p except the first
omega = linprog(f,[],[],[],[],-Inf(10,1),log2(Pt)*ones(10,1));
P = 2.^omega
gamma_infinity_i = lamda.*(1-tau^2).*P.'*Pt/segma_squared;
error = sum(abs(gamma_infinity_im1-gamma_infinity_i))
gamma_infinity_im1 = gamma_infinity_i;
end
3 commentaires
Torsten
le 21 Nov 2022
I'm quite sure we both still don't understand the model correctly.
So I suggest you let the MATLAB code pause for a while and read the background article from which the algorithm was developed again carefully.
Do you have a link to the article in question ?
Plus de réponses (0)
Voir également
Catégories
En savoir plus sur Matrix Indexing dans Help Center et File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!