In spite of being real, I get the error: Warning: Imaginary parts of complex X and/or Y arguments ignored.

Question

raha ahmadi 2022-5-6

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1713390-in-spite-of-being-real-i-get-the-error-warning-imaginary-parts-of-complex-x-and-or-y-arguments-ig

评论： raha ahmadi 2022-5-7

Hi every one. I have a system of ODE. My solution is real (isreal (A_ODE)=1), but I get this error: Warning: Imaginary parts of complex X and/or Y arguments ignored. Actually I expected the steady- state soIution but I got unstable system. I dont know, is this error massage affected on the answer or not I attached the code and really appreciate any help.

thanks in advance

LaserODE

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

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Walter Roberson 2022-5-6

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1713390-in-spite-of-being-real-i-get-the-error-warning-imaginary-parts-of-complex-x-and-or-y-arguments-ig#answer_958830

在 MATLAB Online 中打开

dNis = Nis - Nth;
Ais = sqrt((AInj^2 - gammaN/gammaP*dNis)/(1 + g*dNis/gammaP));

dNis is negative.

That makes the numerator for the next line a value minus a negative value, which is positive, so the numerator is okay.

But with dNis being sufficiently negative, the denominator is negative, so overall you are taking square root of a negative value.

phi_Tr = zeros(size(A_Tr));

You do not assign any other value to phi_Tr, so phi_Tr is all zero at the time you semilogy(), which means you are trying to plot the log of 0.

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

raha ahmadi 2022-5-7

编辑：raha ahmadi 2022-5-7

I know from the steady state solution and physics of the problem that phi must be limited on the set:

. It seems I must impose this constraint on the solution

raha ahmadi 2022-5-7

Dear @Walter Roberson

The problem is solved. these lasers after injection experience non stable situations. I tried to impose constraints on phi, but I couldnt. In steady state solution after calculation we choose between the answers and choose those satisfy the physical conditions. I try to repeat this strategy hear but nothing happens. I changed one key parameter (rInj) and reach the answer.

again very thanks for your help

with the best wishes

请先登录，再进行评论。