Imaginary parts of complex X and/or Y arguments ignored

Question

Toghrul Nasirli 2020-4-4

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/515408-imaginary-parts-of-complex-x-and-or-y-arguments-ignored

编辑： dpb 2020-4-7

在 MATLAB Online 中打开

Hi, I wrote this code for my homework and it works but I get this error.

Warning: Imaginary parts of complex X and/or Y arguments ignored

> In EthaneSatcurve (line 36)

Warning: Imaginary parts of complex X and/or Y arguments ignored

> In EthaneSatcurve (line 38)

I cannot solve this problem, could you please help me?

function EthaneSatcurve()
Pc = 706.5;% psi
Tc = 549.92;% R
Vc = 0.0783;% ft3/lbmole
alpha = [ 1.124 ];
R = 10.732; % psi ft3/(lbmole-R)
a = (0.42727*(R.^(2)).*(Tc.^(2)))./Pc;
b = (0.08664 * R * Tc) / Pc;
V = linspace(b*1.2,100,1000);   
T= -10;
T = T + 460; % temperature in R
Psatguess=600;
Psat = fzero(@SRKLoopAreaDif, Psatguess);
fSRKEOSp = @(T,Vx,alpha)((R*T)./(Vx-b)) - (a*alpha./(Vx.*(Vx+b)));
P = fSRKEOSp(T,V,alpha);
plot(V,P); 
xlim([0 Vc*100]);
ylim([0 500]);
xlabel('Volume, ft^3');
ylabel('Pressure, psi');
hold all;
plot(Vsat(1),Psat,'bo',Vsat(2),Psat,'bo',Vsat(3),Psat,'bo')
hold all;
plot([Vsat(1) Vsat(3)],[Psat Psat]);
function Z = SRKLoopAreaDif(Px)
SRKp = [Px -R*T +(a*alpha-b*R*T-b^(2)*Px) -(a*b*alpha)];
Vsat = roots(SRKp);
Vsat = sort(Vsat);
if any(imag(Vsat))
	Z = 100;
	return
end
	
fSRKEOSpvInt = @(Vx)((a.*log(b+Vx)-a.*log(Vx)+b.*R*T.^(3/2)*log(Vx-b))/(b.*T.^(1/2)));
i1 = fSRKEOSpvInt(Vsat(1));
i2 = fSRKEOSpvInt(Vsat(2));
i3 = fSRKEOSpvInt(Vsat(3));
j1 = (Vsat(2) - Vsat(1)) * Px;
j2 = (Vsat(3) - Vsat(2)) * Px;
AreaLeft = j1 - (i2 - i1);
AreaRight = (i3 - i2) - j2;
Z = AreaLeft - AreaRight;
end
end

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

dpb 2020-4-4

编辑：dpb 2020-4-7

在 MATLAB Online 中打开

We can't run your code because

>> EthaneSatCurve
Error:File: EthaneSatCurve.m Line: 26 Column: 6
Identifier 'Vsat' is not a function or a shared variable. To share 'Vsat' with nested function,
initialize it in the current scope. For more
information, see Sharing Variables Between Parent and Nested Functions. 
>> 

and we don't know for certain which lines are 36 and 38 but

plot(Vsat(1),Psat,'bo',Vsat(2),Psat,'bo',Vsat(3),Psat,'bo')
hold all;
plot([Vsat(1) Vsat(3)],[Psat Psat]);

look like likely candidates; particularly when you have specific code

function Z = SRKLoopAreaDif(Px)
SRKp = [Px -R*T +(a*alpha-b*R*T-b^(2)*Px) -(a*b*alpha)];
Vsat = roots(SRKp);
Vsat = sort(Vsat);
if any(imag(Vsat))
	Z = 100;
	return
end

to fix up a problem that your solution to Vsat returns complex results. That can't be right besides the problem of the return value from the function that calculates Vsat doesn't return that result.

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Harshendra Shah 2020-4-7

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/515408-imaginary-parts-of-complex-x-and-or-y-arguments-ignored#answer_424497

Hi Toghrul,

This issue might be occurring because you are trying to plot complex numbers using the "plot" function.

To debug this issue, you can do the following:

Stop MATLAB on the line where the warning occurs.
Verify that the two vectors or matrices you pass to the PLOT function, are not complex.
If you want to plot the real part of a vector versus the complex part, pass the vector as a single complex vector to the PLOT function.
If you want to plot the magnitude of the elements, use the ABS function on the vector before passing it to the PLOT function.

You can refere to the following MATLAB Answers link for the same:

https://www.mathworks.com/matlabcentral/answers/96197-why-do-i-get-warning-messages-about-imaginary-parts-of-complex-x-and-or-y-arguments-ignored