Raising a matrix to a power
42 views (last 30 days)
Show older comments
I want to raise the Tr value to the power of 1/2. I keep getting the same error.
0 Comments
Answers (2)
Guillaume
on 7 Sep 2019
The problem is not the .^0.5, it's the ^2, which probably should be .^2. Note that x.^0.5 is the same as sqrt(x).
However, If Tr is a matrix, then it's unlikely that alpha(Tr) is a valid target for assignment.
4 Comments
madhan ravi
on 9 Sep 2019
Also note don’t name variable in the name of MATLAB’s inbuilt function. alpha() is an inbuilt function.
Guillaume
on 9 Sep 2019
Yes, as said, the current error is trivially fixed by replacing ^2 by .^2. As for the alpha(Tr) =, it's not clear what you intended to do with that.
Bruno Luong
on 9 Sep 2019
% This script will generate a plot of pressure P (in bar) as a function of molar
% volume V (in cm^3/mol) for various values of temperature T (in K) for
% ethylene as described by the Peng-Robinso n equation of state.
T=[260;270;280;282.4;290;300;310];
V=[50:400];
R= 83.14; % cm^3bar/mol K
Tc= 282.4; %K
Tr= T./Tc;
a= 5.001*10^6;% barcm^6/mol^2
b= 36.24; % cm^3/mol
k= 0.5098;
alpha = (1+k*(1-((Tr).^0.5)).^2);
P= (((R*T)./(V-b))- (a*alpha)./((V.^2)+((2*b*V)-b.^2 ))) ;
surf(V,alpha,P)
7 Comments
madhan ravi
on 10 Sep 2019
Edited: madhan ravi
on 10 Sep 2019
The comment wasn’t meant to you , it was meant for the OP who reads it, if you took it that way then I apologise. As I said earlier it is not a problem for you but for others it might.
Walter Roberson
on 10 Sep 2019
I see posts from people having trouble with alpha roughly once a year.
See Also
Categories
Find more on Graphics Performance in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!