while loop and for loop get different result

2 次查看（过去 30 天）
LIN kai
LIN kai 2018-8-2
编辑： Stephen23 2018-8-3
This is a simple problem.I think I had coded correct.But it's still have some problem. The problem is why these three plot are different?

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

回答（3 个）

Masoud
Masoud 2018-8-2
I saw no difference in figures, but your loops may be different because you wrote : "while t<40" which does not include the 40, it self. But "for" contains the upper bound, which is 40. In other word, your first loop is done 4000 times, however, the "for" loop runs 4001 times. This may be the main result of difference that you are looking for.
Bob Thompson
Bob Thompson 2018-8-2
The while loop doesn't produce the same plot because the values of t are not exactly the same as the values desired. If you look at the final value of t in long format you get 40.000000000000611. Although this is effectively 40, when you subtract 40 from it you do not get 0, so the dirac() command does not result in 0.
Unfortunately, I do not know why the while loop produces this error, but it is the reason you aren't getting the same plot.
  2 个评论
Stephen23
Stephen23 2018-8-2
编辑：Stephen23 2018-8-3
"Unfortunately, I do not know why the while loop produces this error..."
Because of this line:
t=t+0.01;
The value 0.1 cannot be represented using binary floating point numbers, so that addition is actually
t = t+0.1+someError
This addition (with error) is repeated lots of times, thus introducing more and more error into the calculated values. In short, using a loop is a bad way to define those values. MATLAB offers much more numerically robust methods via linspace and the colon operators.
Learn about the practical limits of floating point numbers:
http://www.mathworks.com/help/matlab/matlab_prog/floating-point-numbers.html
http://www.mathworks.com/matlabcentral/answers/69-why-does-1-2-3-1-3-not-equal-zero
http://www.mathworks.com/matlabcentral/answers/57444-faq-why-is-0-3-0-2-0-1-not-equal-to-zero
http://www.mathworks.com/company/newsletters/news_notes/pdf/Fall96Cleve.pdf
http://www.mathworks.com/matlabcentral/answers/41536-why-is-2-24-10-22-4-not-equal-to-0
http://www.mathworks.com/company/newsletters/news_notes/pdf/Fall96Cleve.pdf
http://blogs.mathworks.com/loren/2006/08/23/a-glimpse-into-floating-point-accuracy/
And some external links on this topic:
http://www.lahey.com/float.htm
http://matlabgeeks.com/tips-tutorials/floating-point-comparisons-in-matlab/
https://docs.oracle.com/cd/E19957-01/806-3568/ncg_goldberg.html
Bob Thompson
Bob Thompson 2018-8-2
Huh, I suppose the binary makes sense. I had been able to track it down to that line as well, I just didn't have a good reason of why there was error.

请先登录，再进行评论。

LIN kai
LIN kai 2018-8-3
Thanks for everyone's answer. I'm already understand.

类别

在 Help Center 和 File Exchange 中查找有关 Loops and Conditional Statements 的更多信息

标签

产品

Community Treasure Hunt

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

Start Hunting!

Translated by