how can i write attached equation in MATLAB code
1 Ansicht (letzte 30 Tage)
Ältere Kommentare anzeigen
Akzeptierte Antwort
Walter Roberson
am 17 Nov. 2021
syms Delta_P_lo Delta_x phi_lo Delta_P_f x LB UB
Delta_P_f = Delta_P_lo * (1/Delta_x .* int(phi_lo.^2, x, LB, UB, 'hold', true))
The 'hold', true is there to prevent int() from immediately doing the integral. is independent of x so int() would treat it as being constant in x and would immediate do the integration, leaving you with something that did not include the integration sign.
If is dependent on x then the code would need to be changed a little to reflect that.
I had to introduce LB and UB as the bounds of integration because it happens that the 'hold' feature would be misinterpreted as being bounds if you were to try to do indefinite integration while 'hold' is on.
3 Kommentare
Walter Roberson
am 21 Nov. 2021
Bearbeitet: Walter Roberson
am 21 Nov. 2021
syms d g mu_1 mu_v rho_1 rho_v sigma chi
GAMMA = (rho_1/rho_v)^sym(1/2) / (mu_v / mu_1)^sym(1/8)
N_conf = 1/d * (sigma/ (g*(rho_1 - rho_v)))^sym(1/2)
phi_lo__2 = 1 + (sym(4.3) * GAMMA^2 - 1) * (N_conf * chi^sym(7/8) * (1-chi)^sym(7/8) + chi^sym(1/8))
syms Delta_chi Delta_P_f e_lo G L R
f_lo = (sym(0.79) * R * e_lo - sym(1.64))^(-2)
Delta_P_lo = 2*f_lo * G^2 * L / (rho_1 * d)
syms LB UB
Delta_P_f = Delta_P_lo * (1/Delta_chi .* int(phi_lo__2, chi, LB, UB, 'hold', true))
dPf = simplify(release(Delta_P_f))
Note that the Γ that shows up in the result is not the same as the Γ that you define at the bottom left of the equations. The one in the results is the Gamma function -- the extension of factorials.
Weitere Antworten (0)
Siehe auch
Kategorien
Mehr zu Calculus finden Sie in Help Center und File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!