[Scilab-users] non linear optimization why produce Nan

[Stéphane Mottelet]

You have to remove the first line of Value, because of the header "x,y"  
in your sheet:

--> Value
ans  =

    Nan    Nan   Nan
    0.     97.   Nan
    20.    74.   Nan
    40.    83.   Nan
    60.    84.   Nan
    80.    84.   Nan
    100.   85.   Nan
    120.   86.   Nan
    140.   84.   Nan
    160.   84.   Nan
    (...)
    2820.   585.    Nan
    2840.   303.    Nan
    2860.   28.     Nan
    2880.  -132.    Nan
    2900.  -205.    Nan
    2920.  -214.    Nan
    2940.  -170.    Nan
    2960.  -107.    Nan
    2980.  -62.     Nan
    3000.  -32.     Nan
    3020.  -2.      Nan
    3040.   22.     Nan
    3060.   40.     Nan
    3080.   54.     Nan

Just insert

Value(1,:) = []

after calling read_xls

S.


Le 02/04/2019 à 10:48, Stéphane Mottelet a écrit :
> Hello,
>
> Le 02/04/2019 à 10:44, karoli a écrit :
>> Hello mottelet
>> I tested it, residual results is ok
> No it is not ok
>> as bellow shown;
>> -->dyke_error(p0,Value, i, y)
>>   ans  =
>>        Nan
> you have a "Not a Number" in first component
>>      - 332.01567
>>    - 312.19667
>>    - 324.46484
>>    - 328.8238
>>    - 332.27738
>>    - 336.82964
>>    - 341.48485
>>    - 343.24755
>>    - 347.12254
>>    - 353.11492
>>    - 355.2301
>>    - 359.47381
>>    - 360.85216 ........ continue
>>
>> but once I run leastsq function again, I got same results
>>
>> -->[f,p,g]=leastsq(dyke_error,p0)
>>   g  =
>>        Nan
>>      Nan
>>      Nan
>>      Nan
>>      Nan
>>   p  =
>>        53.8
>>      1527.
>>      43.7
>>      55.7
>>      3424.7
>>   f  =
>>        Nan
>>
>> I still don't understand what's wrong with my code.
>>
>> regards
>> karoli
>>
>>
>>
>> -- 
>> Sent from: 
>> https://antispam.utc.fr/proxy/2/c3RlcGhhbmUubW90dGVsZXRAdXRjLmZy/antispam.utc.fr/proxy/1/c3RlcGhhbmUubW90dGVsZXRAdXRjLmZy/mailinglists.scilab.org/Scilab-users-Mailing-Lists-Archives-f2602246.html
>> _______________________________________________
>> users mailing list
>> users at lists.scilab.org
>> https://antispam.utc.fr/proxy/2/c3RlcGhhbmUubW90dGVsZXRAdXRjLmZy/antispam.utc.fr/proxy/1/c3RlcGhhbmUubW90dGVsZXRAdXRjLmZy/lists.scilab.org/mailman/listinfo/users 
>>
>
>

-- 
Stéphane Mottelet
Ingénieur de recherche
EA 4297 Transformations Intégrées de la Matière Renouvelable
Département Génie des Procédés Industriels
Sorbonne Universités - Université de Technologie de Compiègne
CS 60319, 60203 Compiègne cedex
Tel : +33(0)344234688
http://www.utc.fr/~mottelet