Итерационный метод - Matlab - PullRequest
0 голосов
/ 27 мая 2020

Пожалуйста, помогите мне с этим.

У меня есть такой псевдокод:

-

Итерационный метод с 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, как указано выше?

Большое спасибо!

Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...