NaN and Infs popping in my matrices after a few iterations of EM algorithm

2 次查看（过去 30 天）
Jérôme
Jérôme 2012-11-20
Hello ;
So, I'm implementing the EM algorithm in Matlab, but my matrices quickly end up contaminated by NaN and Inf values. I think it might be caused by matrix inversions, but I'm not sure it's the only reason.
Here is the code : em_algo.m http://pastebin.com/vJGqQPw7 iterates until convergence (which never happens due to my issue) and calls smoother.m http://pastebin.com/yD9Je37i at each iteration. The NaNs and Infs start popping around the ~8th iteration, usually.
I guess in there somewhere I'm doing something unholy with my matrices, but I really have no clue about what's wrong. I trust your expertise.
Thanks in advance for the help!

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

回答（3 个）

José-Luis
José-Luis 2012-11-20
It sounds like you are running into numerical stability issues. Look at the following example, that might help you understand what is happening:
val = 1;
while (val ~= Inf)
val = val / 10^-100
end
The NaN's that pop up can be a result of
0/0
Jérôme
Jérôme 2012-11-22
Thanks José-Luis. Indeed, something like that might be the source of my problem. But I have a hard time finding out WHERE it happens in the code. :/
Another issue (or most likely another aspect of the same one) is the fact that some values seem to grow out of bonds which probably leads to Inf values. But again, I have no idea why this is happening.
Matt J
Matt J 2012-11-23
编辑：Matt J 2012-11-23
Set a conditional breakpoint at line 45 (making it conditional on i==8 or whenever the NaNs start to pop up). When the code stops there, see whether inv(A) has sensible values.

类别

在 Help Center 和 File Exchange 中查找有关 Logical 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by