SUBS problem

[Jan Houska]

Dobry den,

skoro se mi zda, jako kdyby ten vyraz ztratil presnost (presneji receno 
rozliseni) okolo cisla 1. V pripade MATLABu R2006a je to myslim 
matematicky spravne:

>> subs(y,x,1-eps)
ans =
    -1
>> subs(y,x,1)
ans =
   NaN
>> subs(y,x,1+eps)
ans =
     0

zatimco MATLAB R2006b ztratil rozliseni:

>> subs(y,x,1-9*eps)
ans =
   NaN
>> subs(y,x,1)
ans =
   NaN
>> subs(y,x,1+10*eps)
ans =
   NaN

tedy cisla od 1-9*eps do 1+10*eps bere jako 1, a teprve pri vetsich 
odchylkach od 1 vyhodnoti vyraz spravne. To povazuji za chybu v MATLABu 
R2006b a nahlasim to, diky za upozorneni.

Zdravi

Jan Houska



> zjistili jsme nasledujici vypocetni nesrovnalosti mezi verzi Matlabu
> 2006a a 2006b.
> Po zadani (potrebny Symbolic Math Toolbox):
> 
> syms x
> y = heaviside(x)-2*heaviside(x-1/2)+heaviside(x-1);
> subs(y,x,1+2*eps)
> 
> 
> %% ve verzi 2006a dostaneme
> 
> ans =
> 
> 0
> 
> 
> %% ve verzi 2006b dostaneme
> 
> ans =
> 
> NaN
> 
> 
> proc? neni to chyba? nebo byla to chyba?
> 
> Jedine co jsme byli schopni o problemu zjistit je ze funkce SUBS byla z
> 2006a do 2006b zmenena. Soubor $MATLAB\toolbox\symbolic\@sym\subs.m ma
> nasledujici verze
> 
> 2006a %   $Revision: 1.39.4.4 $  $Date: 2005/03/07 17:42:10 $
> 2006b %   $Revision: 1.39.4.5 $  $Date: 2006/06/20 20:52:42 $
> 
> Dekujem
> -- 
> Ales PAVELKA
> 
> ICQ:171-964-427
> http://alespavelka.wz.cz/
> http://uprt.vscht.cz/pavelka
> http://dsp.vscht.cz/pavelka
> http://www.volny.cz/ales.pavelka
> http://www.volny.cz/pavelka.ales
> 
> Tomorrow's the day My bride's gonna come
> Oh, oh, are we gonna fly Down in the easy chair!
> 
> 
> 
> -- 
> Pokerová mánie - 6 freerollových turnajů o 5000 USD!
> http://im.impact.as/sign/betway/
> 
> 
> _______________________________________________
> CsMUG mailing list
> CsMUG at vscht.cz
> http://listserv.vscht.cz/mailman/listinfo/csmug


--
Jan Houska                                           HUMUSOFT s.r.o.
houska at humusoft.com                                  Pobrezni 20
http://www.humusoft.com                              186 00 Praha 8
tel: ++ 420 284 011 730                              Czech Republic
fax: ++ 420 284 011 740