Пожалуйста, помогите мне с этим.
У меня есть такой псевдокод:
Итерационный метод с 3 циклами упомянутый в этом псевдокоде, могу ли я использовать итерацию с фиксированной точкой, как показано ниже:
for n = 1:500
Ex(n) = Ex(1);
maxIter = 3;
%--------Computation using Fixed point Iteration------%
Ex(nz) = Ex(n);
for i = 1 : maxIter
Phi_Kx(nz) = Ex(nz).^2;
P_LDx(nz) = (4/(2+(omega^2*dt^2)))*P_LDx(nz) - P_LDx(nz) + (2*omega^2*dt^2/(2+(omega^2*dt^2)))*Ex(nz);
Xi_Rx(nz) = ((4-2*dt^2)/(dt*BW+2))*Xi_Rx(nz) +((dt*BW-2)/(dt*BW+2))*Xi_Rx(nz) +((2*dt^2*BW)/(dt*BW+2))* Ex(nz);
N(nz) =((2*tau_c-dt)/(2*tau_c+dt))* N(nz) +((tau_c*dt)/(2*tau_c-dt))*(10/(4*hbar*omega))*Ex(nz).^4;
A2 = epsilon_2.*S_Kx(nz) + Xi_Rx(nz) - (3*N(nz)-N(nz-1))+(c*n0*dt/4)*(alpha_loss+alpha_fca) + (n0^2*c^2.*epsilon*beta_tpa/4).*S_Kx(nz);
D_NLx = Dx - P_LDx;
Ex = (A1./A2).*Ex(nz)- (D_NLx(nz)-D_NLx(nz-1))./(epsilon.*A2);
end
end
В этом l oop, Phi_Kx, P_LDx, Xi_Rx, N являются функциями E (n), которые также вычисляются before
Нужно ли мне включать их в l oop, как указано выше?
Большое спасибо!