Which documentation and example can I use for iterative solution.

2 次查看（过去 30 天）

Tadele 2023-1-12

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1892795-which-documentation-and-example-can-i-use-for-iterative-solution

编辑： Tadele 2023-1-12

Dear expert,

I was struglling with my matlab model. I have a cascade of functions which each set of function recieves the outputs of the other function as an input arguments. To run all these cascaded functions, the first set of functions needs guess values. The outputs of this functions set should be used as inputs for the next set of functions and continue like this untill the final part. Then the outputs of the last set of functions should be the same with guessed values in the begining using iterative approach. This will be the solution of my model. However, I don't know how to make an iterative solution for this problem which will stop iterating when the guess value is equal to the final calculated value (the difference between them converges to zero).

Can someone help me with examples or which documentation can I use, please? I thank you a lot in advance.

N.B I can share codes but it is too much for help.

12 个评论
显示 10更早的评论隐藏 10更早的评论

Tadele 2023-1-12

在 MATLAB Online 中打开

That is very clear.

And of course "fsolve" solves the equation defined by "find_CP_MSNP" shown below.

My last question is:

Can we say that the solver is using ''Cp_guess_arb(1:3)'' as an initial value for the cascaded functions and trying to satisfy the equations defined in "find_CP_MSNP" that is set to converge the condition defined by "sol"? (Is my approach correct?)

Equation defined by "find_CP_MSNP" looks as shown below.

function [sol] = find_CP_MSNP( Cp_guess, phi_start, Jv, Kc, DP,phi, Z, z_pol, X_mem)
Cp_guess_MS = Cp_guess; 
[y0, yp0]  = decic(@MS_AA_imp, 0, [phi_start 0], [1 1 0 0 0], [0 0 0 0 0], [],[],Cp_guess_MS, Jv, Z);
stepsMS = linspace(0 , z_pol, 500) ;
[zsol_MS, ysol_MS] = ode15i(@(z, x, dx_dz)MS_AA_imp(z, x, dx_dz ,Cp_guess_MS, Jv, Z), stepsMS, y0, yp0);
C_surf = ysol_MS(end,:);
%-------------------------------------------------------------------------------
% define anonymous function 'partionAnnonymous' which accepts
% C_surf(out put from the top) as an input for the next function
zz=[0.00001 0.00001 0.00001];
partitonAnonymous = @(z)partition(z,C_surf(1:3));
C_entrance = fsolve(partitonAnonymous,zz);
%-------------------------------------------------------------------------------
% define anonymous function 'nernstPlanckAnonymous'  which accept
% C_entrance as an input arguement for the next function
nernstPlanckAnonymous = @(z)NernstPlanck_Vol(z,C_entrance);
C_end = fsolve(nernstPlanckAnonymous,zz);
%-------------------------------------------------------------------------------
% define anonymous function 'secondpartitionAnonymous'  which accept
% C_end as an input arguement and computes the Cp_calc
secondpartitionAnonymous = @(z)partitionsecond(z,C_end(1:3)); 
Cp_calc = fsolve(secondpartitionAnonymous,zz);
%-------------------------------------------------------------------------------
%% the last par that I want the solver [Cp_calc_arb, fval, exitflag] solves
sol = ((y(1:3) - Cp_guess(1:3)).^2);
end