solving system of multiple equations
1 visualizzazione (ultimi 30 giorni)
Mostra commenti meno recenti
My code is:
syms BF1 BF2 BF3 BF4 CO COreq
eqn1 = BF3(i)==0.012*CO;
eqn2 = BF1(i)==0.92*CO;
eqn3 = BF4(i)==MaxBFvsk*PctBF*1000;
eqn4 = BF2(i)==(100*MuVO2)/(0.8*210);
eqn5 = COreq==BF1(i)+BF2(i)+BF3(i)+BF4(i);
eqn6 = CO==min(COmax,COreq);
[BF1(i),BF2(i),BF3(i),BF4(i),CO,COreq]=solve(eqn1,eqn2,eqn3,eqn4,eqn5,eqn6,BF1(i),BF2(i),BF3(i),BF4(i),CO,COreq);
Where MaxBFvsk, PctBF, MuVO2, COmax are all constants.
And it doesn't work because of equation 6, the error says:
Error using sym/min (line 98) Input arguments must be convertible to floating-point numbers.
Error in test (line 148) eqn6 = CO==min(COmax,COreq); %ml/min CO=HRR*SV is another method
I am not sure how to correct this problem.
0 Commenti
Risposta accettata
Star Strider
il 7 Ott 2015
‘Where MaxBFvsk, PctBF, MuVO2, COmax are all constants.’
You have to define them after the syms call. If ‘COmax’ and ‘COreq’ contain any symbolic variables (I did not run your code), that line will also throw an error.
You might also want to use the vpasolve function instead of solve.
5 Commenti
Star Strider
il 8 Ott 2015
What are the approximate values for: MaxBFvsk, PctBF, MuVO2, COmax ?
What are the equations you’re starting with? It looks like some sort of metabolic or circulatory model, but I don’t recognise it. It’s fine with me if you post a PDF of the paper or references you’re working from.
Since you already Accepted my Answer, I’ll work with you as much as I can to get a solution that works, if one is possible.
Più risposte (1)
John D'Errico
il 7 Ott 2015
Modificato: John D'Errico
il 7 Ott 2015
Look at equation 6, and then think about the error message.
eqn6 = CO==min(COmax,COreq);
Try this:
syms x y
min(x,y)
Error using sym/min (line 98)
Input arguments must be convertible to floating-point numbers.
The point is, min is invalid there in your equation. Min requires numeric arguments, not symbolic ones.
0 Commenti
Vedere anche
Categorie
Scopri di più su Equation Solving in Help Center e File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!