How can i solve these equations in Matlab?

  • Topic: MATLAB 
  • Thread starter koii123
  • Start date
  • Tags Tags
    Matlab
4 replies · 2K views · last reply
koii123
Messages
8
Reaction score
0
Equation 1:
Equation: eta_R = eta_0 times e raised to the power alpha* times p, with eta_R and eta_0 as symbols, and alpha* and p in the exponent.

Equation 2:
Typeset equation: alpha* p = [ln eta_0 + 9.67] times a curly-brace expression containing ((theta - 138)/(theta_0 - 138)) raised to the power -S_0, multiplied by (1 + 5.1 x 10^(-9) p)^Z minus 1.

Equation 3:
Equation Z = α divided by 5.1 × 10^-9 times [ln η0 + 9.67].

Equation 4:
Equation S0 = β(θ0 − 138) divided by (ln η0 + 9.67), with S0 on the left side of the equals sign.

Equation 5:
Equation: β = [ln η0 + 9.67] [1 + 5.1 × 10^-9 p]^Z [S0 / (θ0 − 138)]


I would like to solve
upload_2015-5-17_22-9-25.png
, (not the one with a star)
and i want to define values for p, η0, ηr, θ0, θr.
i have tried to solve them, but didn't work,
can someone please help and show me step by step?
 
Last edited:
What exactly are you trying to do? I count 10 unknowns (α, α*, p, η0, ηR, θ, θ0, S0, Z, β), and only 5 equations. So the system is underdetermined and there are an infinite number of possible solutions. You either need to fix some of the values or find more equations until the number of equations equals the number of unknowns. What exactly have you tried so far?
 
phyzguy said:
What exactly are you trying to do? I count 10 unknowns (α, α*, p, η0, ηR, θ, θ0, S0, Z, β), and only 5 equations. So the system is underdetermined and there are an infinite number of possible solutions. You either need to fix some of the values or find more equations until the number of equations equals the number of unknowns. What exactly have you tried so far?
Actually, i have values for p, η0, ηr, θ0, θr,
for example:
p=101325,
η0=4.43,
ηr=0.181,
θ0=313.15,
θr=373.15,
also, as you can see, it is possible to solve α , because these 5 equations can be combined into 1
 
OK, that makes more sense. So can you do the following?

(1) Solve equation 4 for the ratio β/S0.
(2) Divide equation 5 by S0 to get β/S0 on the LHS. Then plug the result from step 1 into this and solve this equation for Z.
(3) Take the result from step 2 and plug it into equation 3, then solve this equation for α.

Edit:

Actually, I see now that it is even simpler than that. You can combine equations 4 and 5 to give:

(1+5.1E-9*p)^Z = 1

Solve this for Z and plug it into equation 3 to give α.
 
Last edited:
Looking more at your equations. Equations 4 and 5 together imply that:

(1+5.1E-9*p)^Z = 1

Since p = 101325, (1+5.1E-9*p) = 1.0005, so this equation is only satisfied for Z = 0. Equation 3 then implies that α=0. Are you sure this is all correct, or am I missing something?