I am having trouble constructing a 'for loop'. Can you suggest me on how to modify it?

2 次查看(过去 30 天)
I have written the Matlab code based on the results I obtained. However, the code is not producing the expected result, and I am having trouble pinpointing my mistake. Can you please help me identify where I went wrong and assist me in correcting it?
I want to do a summation.
I211_sum=0;
I22_sum =0;
% Loop for I211
for k= 1: L
for i = 1:k
for m = 1: i-1
I211_sum = I211_sum ...
+ ( nchoosek(L,k) * ( (-1)^k ) * (exp(-H1*k) ) ...
* (exp(B)) * (P^-k) * nchoosek(k,i) * ((B/P)^i)...
* ( (exp(-D) * ( (factorial(m-1)/factorial(i-1)) * ((-A)^(i-1-m)) * ((E- (B/P))^(-m)) ) )...
+ ( ((-A)^(i-1)) * (1/factorial(i-1)) * ((exp(-(B*A)/ P))) * (expint(A * (E- (B/P)))) ) ) ) ;
end
end
end
I211(r) = I211_sum;
% Loop for I22
for k= 1: L
I22_sum = I22_sum + ( nchoosek(L,k) * ((-1)^k) * (P1^k) * (1-( exp(-F) )) * ((exp(-H*k))/(G^k)) ) ;
end
I22(r) = I22_sum;
  1 个评论
Voss
Voss 2024-8-29
Tip: You can use the nchoosek function to avoid having to write three separate factorial function calls.
Example:
L = 5;
k = 3;
factorial(L) / (factorial(k)*factorial(L-k) )
ans = 10
nchoosek(L,k)
ans = 10

请先登录,再进行评论。

采纳的回答

Voss
Voss 2024-8-29
编辑:Voss 2024-8-29
One thing I notice is that
(-A)^i-1
should be
(-A)^(i-1)
EDIT: that's now corrected in the question.
  3 个评论
Voss
Voss 2024-8-29
编辑:Voss 2024-8-29
And the code appears to have an extra *A here
* (exp(B)) * A * (P^-k)
% ^^^
that doesn't seem to be in the image.
EDIT: that's now corrected in the question.

请先登录,再进行评论。

更多回答(0 个)

类别

Help CenterFile Exchange 中查找有关 Loops and Conditional Statements 的更多信息

标签

产品


版本

R2021b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by