Empty sym: 0-by-1
14 vues (au cours des 30 derniers jours)
Afficher commentaires plus anciens
When I run this code, my answer for A comes out as 'Empty sym: 0-by-1'
clc;clear;format compact;
syms V
P=[0.98 1.97 4.93 9.86 49.36 98.69]; %atm
T=[573 573 573 573 573 573]; %K
R=461.5; %Pa*m3/mol*K
R=convpres(R,'Pa','atm');
R=R*1000; %L*atm/mol*K
GI=(P*V)-(R*T);
A=solve(GI,V)
0 commentaires
Réponse acceptée
Torsten
le 12 Oct 2023
Use
A = double(arrayfun(@(i)solve(GI(i),V),1:numel(P)))
instead of
A = solve(GI,V)
1 commentaire
Dyuman Joshi
le 12 Oct 2023
which can be simplified to
out = double(arrayfun(@(x) solve(x,V), GI))
Plus de réponses (1)
Steven Lord
le 12 Oct 2023
Let's look at the equations you're trying to solve.
syms V
P=[0.98 1.97 4.93 9.86 49.36 98.69]; %atm
T=[573 573 573 573 573 573]; %K
R=461.5; %Pa*m3/mol*K
R=convpres(R,'Pa','atm');
R=R*1000; %L*atm/mol*K
GI=(P*V)-(R*T)
The first element of GI has a solution, as does the second element of GI. Is the solution for those two elements exactly the same? To ask that question another way, is there a value of V that satisfies all the elements in GI simultaneously?
A1=solve(GI(1), V)
A2 = solve(GI(2), V)
isAlways(A1 == A2)
No, there isn't. So your original solve call that's looking for that one value of V that satisfies all the elements in GI simultaneously is correct to return an empty sym.
Now if you wanted a vector of values of V that satisfy each corresponding element in GI, that's possible to compute.
A = zeros(size(GI), 'sym');
for whichTerm = 1:numel(GI)
A(whichTerm) = solve(GI(whichTerm), V);
end
A
Do they satisfy their equations?
result = zeros(size(GI), 'sym');
for whichEquation = 1:numel(GI)
result(whichEquation) = subs(GI(whichEquation), V, A(whichEquation));
end
result
Yes. Are they close to the same value?
vpa(A, 5)
Nope. Not even same order of magnitude.
1 commentaire
Walter Roberson
le 12 Oct 2023
This is an important point: solve is always looking for simultaneous solution of all the inputs. It is designed for simultaneous equation solving, and does not have any option for solving each equation individually.
Voir également
Catégories
En savoir plus sur Number Theory dans Help Center et File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!