<html>
<head>
<meta content="text/html; charset=UTF-8" http-equiv="Content-Type">
</head>
<body text="#000000" bgcolor="#FFFFFF">
<div class="moz-cite-prefix">Hello,<br>
<br>
The fit works with the following:<br>
<pre style="font-family:Monospaced;font-style:normal;font-size:12.0;"><span style="color:rgb(100,174,100);font-style:italic;">//xm;ym</span>
<span style="color:rgb(0,0,0);">d</span> <span style="color:rgb(92,92,92);">=</span> <span style="color:rgb(74,85,219);">[</span>
<span style="color:rgb(188,143,143);">2.0</span> <span style="color:rgb(188,143,143);">99.9449173761</span>
<span style="color:rgb(188,143,143);">1.0</span> <span style="color:rgb(188,143,143);">99.8097145719</span>
<span style="color:rgb(188,143,143);">0.5</span> <span style="color:rgb(188,143,143);">97.9769654482</span>
<span style="color:rgb(188,143,143);">0.25</span> <span style="color:rgb(188,143,143);">36.4046069104</span>
<span style="color:rgb(188,143,143);">0.1</span> <span style="color:rgb(188,143,143);">1.4872308463</span>
<span style="color:rgb(188,143,143);">0.071</span> <span style="color:rgb(188,143,143);">0.5207811718</span>
<span style="color:rgb(188,143,143);">0.063</span> <span style="color:rgb(188,143,143);">0.3705558338</span>
<span style="color:rgb(74,85,219);">]</span><span style="color:rgb(0,0,0);">;</span>
<span style="color:rgb(0,0,0);">xm</span> <span style="color:rgb(92,92,92);">=</span> <span style="color:rgb(0,0,0);">d</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(255,170,0);">:</span><span style="color:rgb(0,0,0);">,</span><span style="color:rgb(188,143,143);">1</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(0,0,0);">;</span>
<span style="color:rgb(0,0,0);">ym</span> <span style="color:rgb(92,92,92);">=</span> <span style="color:rgb(0,0,0);">d</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(255,170,0);">:</span><span style="color:rgb(0,0,0);">,</span><span style="color:rgb(188,143,143);">2</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(0,0,0);">;</span>
<span style="color:rgb(100,174,100);font-style:italic;">//Script used to fit data</span>
<span style="color:rgb(0,0,0);">miny</span> <span style="color:rgb(92,92,92);">=</span> <span style="color:rgb(50,185,185);">min</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(0,0,0);">ym</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(0,0,0);">;</span>
<span style="color:rgb(0,0,0);">maxy</span> <span style="color:rgb(92,92,92);">=</span> <span style="color:rgb(50,185,185);">max</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(0,0,0);">ym</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(0,0,0);">;</span>
<span style="color:rgb(0,0,0);">minx</span> <span style="color:rgb(92,92,92);">=</span> <span style="color:rgb(50,185,185);">min</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(0,0,0);">xm</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(0,0,0);">;</span>
<span style="color:rgb(0,0,0);">maxx</span> <span style="color:rgb(92,92,92);">=</span> <span style="color:rgb(50,185,185);">max</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(0,0,0);">xm</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(0,0,0);">;</span>
<span style="color:rgb(100,174,100);font-style:italic;">//Function describing problem (I want to find best fit by finding n()).</span>
<span style="color:rgb(176,24,19);">function</span> <span style="color:rgb(131,67,16);font-weight:bold;">g1</span><span style="color:rgb(92,92,92);">=</span><span style="color:rgb(0,0,0);">f1</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(131,67,16);font-weight:bold;">x</span><span style="color:rgb(0,0,0);">, </span><span style="color:rgb(131,67,16);font-weight:bold;">n</span><span style="color:rgb(74,85,219);">)</span>
<span style="color:rgb(131,67,16);font-weight:bold;">g1</span> <span style="color:rgb(92,92,92);">=</span> <span style="color:rgb(0,0,0);">miny</span> <span style="color:rgb(92,92,92);">+</span> <span style="color:rgb(74,85,219);">(</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(0,0,0);">maxy</span><span style="color:rgb(92,92,92);">-</span><span style="color:rgb(0,0,0);">miny</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(92,92,92);">./</span> <span style="color:rgb(255,170,0);">..</span>
<span style="color:rgb(74,85,219);">(</span> <span style="color:rgb(74,85,219);">(</span><span style="color:rgb(188,143,143);">1</span> <span style="color:rgb(92,92,92);">+</span> <span style="color:rgb(74,85,219);">(</span><span style="color:rgb(131,67,16);font-weight:bold;">n</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(188,143,143);">1</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(92,92,92);">*</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(131,67,16);font-weight:bold;">x</span><span style="color:rgb(92,92,92);">/</span><span style="color:rgb(131,67,16);font-weight:bold;">n</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(188,143,143);">2</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(92,92,92);">.^</span><span style="color:rgb(131,67,16);font-weight:bold;">n</span><span style="color:rgb(74,85,219);">(</
span><span style="color:rgb(188,143,143);">3</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(74,85,219);">)</span> <span style="color:rgb(92,92,92);">+</span> <span style="color:rgb(74,85,219);">(</span><span style="color:rgb(188,143,143);">1</span><span style="color:rgb(92,92,92);">-</span><span style="color:rgb(131,67,16);font-weight:bold;">n</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(188,143,143);">1</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(92,92,92);">*</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(131,67,16);font-weight:bold;">x</span><span style="color:rgb(92,92,92);">/</span><span style="color:rgb(131,67,16);font-weight:bold;">n</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(188,143,143);">2</span><span style="color:rgb(74,85,219);">)</span><span style=
"color:rgb(74,85,219);">)</span><span style="color:rgb(92,92,92);">.^</span><span style="color:rgb(131,67,16);font-weight:bold;">n</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(188,143,143);">4</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(92,92,92);">.^</span><span style="color:rgb(131,67,16);font-weight:bold;">n</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(188,143,143);">5</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(74,85,219);">)</span>
<span style="color:rgb(176,24,19);">endfunction</span>
<span style="color:rgb(176,24,19);">function</span> <span style="color:rgb(131,67,16);font-weight:bold;">e</span><span style="color:rgb(92,92,92);">=</span><span style="color:rgb(0,0,0);">G</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(131,67,16);font-weight:bold;">n</span><span style="color:rgb(0,0,0);">, </span><span style="color:rgb(131,67,16);font-weight:bold;">z</span><span style="color:rgb(74,85,219);">)</span>
<span style="color:rgb(131,67,16);font-weight:bold;">e</span> <span style="color:rgb(92,92,92);">=</span> <span style="color:rgb(50,185,185);">abs</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(131,67,16);font-weight:bold;">z</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(188,143,143);">2</span><span style="color:rgb(74,85,219);">)</span> <span style="color:rgb(92,92,92);">-</span> <span style="color:rgb(0,0,0);">f1</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(131,67,16);font-weight:bold;">z</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(188,143,143);">1</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(0,0,0);">,</span><span style="color:rgb(131,67,16);font-weight:bold;">n</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(74,85,219);">)</span> // Note here that we take the modulus
<span style="color:rgb(176,24,19);">endfunction</span>
<span style="color:rgb(0,0,0);">n0</span> <span style="color:rgb(92,92,92);">=</span> <span style="color:rgb(74,85,219);">[</span><span style="color:rgb(188,143,143);">1</span> <span style="color:rgb(0,0,0);">;</span> <span style="color:rgb(188,143,143);">1</span><span style="color:rgb(0,0,0);">;</span> <span style="color:rgb(188,143,143);">1</span><span style="color:rgb(0,0,0);">;</span> <span style="color:rgb(188,143,143);">1</span><span style="color:rgb(0,0,0);">;</span> <span style="color:rgb(188,143,143);">1</span><span style="color:rgb(74,85,219);">]</span><span style="color:rgb(0,0,0);">;</span>
<span style="color:rgb(0,0,0);">Zz</span> <span style="color:rgb(92,92,92);">=</span> <span style="color:rgb(74,85,219);">[</span><span style="color:rgb(0,0,0);">xm</span><span style="color:rgb(92,92,92);">.'</span> <span style="color:rgb(0,0,0);">;</span> <span style="color:rgb(0,0,0);">ym</span><span style="color:rgb(92,92,92);">.'</span><span style="color:rgb(74,85,219);">]</span><span style="color:rgb(0,0,0);">;</span>
<span style="color:rgb(74,85,219);">[</span><span style="color:rgb(0,0,0);">n0_opt</span><span style="color:rgb(0,0,0);">,</span><span style="color:rgb(0,0,0);">err</span><span style="color:rgb(74,85,219);">]</span> <span style="color:rgb(92,92,92);">=</span> <span style="color:rgb(174,92,176);">datafit</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(0,0,0);">G</span><span style="color:rgb(0,0,0);">,</span> <span style="color:rgb(0,0,0);">Zz</span><span style="color:rgb(0,0,0);">,</span> <span style="color:rgb(0,0,0);">n0</span><span style="color:rgb(74,85,219);">)
</span>
<span style="color:rgb(174,92,176);">clf</span>
<span style="color:rgb(50,185,185);">plot2d</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(188,143,143);">"</span><span style="color:rgb(188,143,143);">ll</span><span style="color:rgb(188,143,143);">"</span><span style="color:rgb(0,0,0);">,</span> <span style="color:rgb(0,0,0);">xm</span><span style="color:rgb(0,0,0);">,</span><span style="color:rgb(0,0,0);">ym</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(0,0,0);">;</span>
<span style="color:rgb(0,0,0);">c</span> <span style="color:rgb(92,92,92);">=</span> <span style="color:rgb(0,0,0);">e</span><span style="color:rgb(92,92,92);">.</span><span style="color:rgb(170,170,170);">children</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(188,143,143);">1</span><span style="color:rgb(74,85,219);">)</span>
<span style="color:rgb(0,0,0);">c</span><span style="color:rgb(92,92,92);">.</span><span style="color:rgb(170,170,170);">line_mode</span> <span style="color:rgb(92,92,92);">=</span> <span style="color:rgb(188,143,143);">"</span><span style="color:rgb(188,143,143);">off</span><span style="color:rgb(188,143,143);">"</span><span style="color:rgb(0,0,0);">;</span>
<span style="color:rgb(0,0,0);">c</span><span style="color:rgb(92,92,92);">.</span><span style="color:rgb(170,170,170);">mark_mode</span> <span style="color:rgb(92,92,92);">=</span> <span style="color:rgb(188,143,143);">"</span><span style="color:rgb(188,143,143);">on</span><span style="color:rgb(188,143,143);">"</span><span style="color:rgb(0,0,0);">;</span>
<span style="color:rgb(0,0,0);">c</span><span style="color:rgb(92,92,92);">.</span><span style="color:rgb(170,170,170);">mark_size_unit</span> <span style="color:rgb(92,92,92);">=</span> <span style="color:rgb(188,143,143);">"</span><span style="color:rgb(188,143,143);">point</span><span style="color:rgb(188,143,143);">"</span><span style="color:rgb(0,0,0);">;</span>
<span style="color:rgb(0,0,0);">c</span><span style="color:rgb(92,92,92);">.</span><span style="color:rgb(170,170,170);">mark_size</span> <span style="color:rgb(92,92,92);">=</span> <span style="color:rgb(188,143,143);">4</span><span style="color:rgb(0,0,0);">;</span>
<span style="color:rgb(50,185,185);">plot2d</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(0,0,0);">xm</span><span style="color:rgb(0,0,0);">,</span><span style="color:rgb(0,0,0);">f1</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(0,0,0);">xm</span><span style="color:rgb(0,0,0);">,</span><span style="color:rgb(0,0,0);">n0_opt</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(74,85,219);">)</span>
</pre>
We get:<br>
-->[n0_opt,err] = datafit(G, Zz, n0) <br>
err =<br>
0.0485009 <br>
n0_opt =<br>
0.9641563 <br>
0.3183242 <br>
- 7.833407 <br>
- 2.3794352 <br>
0.5056993 <br>
<br>
// + Figure<br>
<img src="cid:part1.09030200.09040204@free.fr" alt="" width="290"
height="249"><br>
<br>
Not so bad.<br>
<br>
Regards<br>
Samuel<br>
<br>
Le 17/12/2014 21:14, Maihem a écrit :<br>
</div>
<blockquote cite="mid:1418847245075-4031549.post@n3.nabble.com"
type="cite">
<pre wrap="">Thanks for answer but it seems that implementation of this only allowed to
avoid error cause calculated parameters still don't fit to data.
Parameters calculated by implementation real on e gives fit parameters:
n(1)= 1.264264
n(2)= - 0.1891816
n(3)= 1.6592188
n(4)= 0.9977104
n(5)= 1.0404815
but this parameters are wrong.
For example I fitted this equation using Gnuplot and I recived:
n(1)= 0.960002
n(2)= 0.331008
n(3)= -8.26945
n(4)= -1.9777
n(5)= 0.423854
and this results from gnuplot fits this function to data
I'm little confused why I can't receive similar solution from scilab
Łukasz
--
View this message in context: <a class="moz-txt-link-freetext" href="http://mailinglists.scilab.org/Problem-with-function-fitting-tp4031546p4031549.html">http://mailinglists.scilab.org/Problem-with-function-fitting-tp4031546p4031549.html</a>
Sent from the Scilab users - Mailing Lists Archives mailing list archive at Nabble.com.
_______________________________________________
users mailing list
<a class="moz-txt-link-abbreviated" href="mailto:users@lists.scilab.org">users@lists.scilab.org</a>
<a class="moz-txt-link-freetext" href="http://lists.scilab.org/mailman/listinfo/users">http://lists.scilab.org/mailman/listinfo/users</a>
</pre>
</blockquote>
<br>
</body>
</html>