[Octave-bug-tracker] [bug #37508] fminsearch misses in deconvolving 2 gaussian peaks

[anonymous]

URL:
  <http://savannah.gnu.org/bugs/?37508>

                 Summary: fminsearch misses in deconvolving 2 gaussian peaks
                 Project: GNU Octave
            Submitted by: None
            Submitted on: fre  5 okt 2012 14.55.14
                Category: None
                Severity: 3 - Normal
                Priority: 5 - Normal
              Item Group: Matlab Compatibility
                  Status: None
             Assigned to: None
         Originator Name: CFB
        Originator Email: address@hidden
             Open/Closed: Open
         Discussion Lock: Any
                 Release: 3.6.1
        Operating System: GNU/Linux

    _______________________________________________________

Details:

1st report from my part. I hope to post in the right place.

I am trying to write code to deconvolve superimposed peaks.
I start by setting up a simulation of 2 gaussian peaks and superimposing them
(instead of importing real chromatographic data - best to use "clean" model
first, right?). Then I try to deconvolve by minimizing residuals.

Main program is exp2.m, which uses functions gaussian.m and fitexp2.m (see
attachments).

Using Matlab at work, it works very nicely. I get the exact same values back
that I put in (at 4 decimals for all 6 parameters, 3 from each peak). In GNU
Octave, it's a mess. 

I haven't hacked "Matlab" code in approx 10 years, so maybe there are better
approaches to solve the problem (please enlighten me with any good ideas).
However, I would expect them to produce the same results...

Again, I hope I'm doing this the right way. Sorry if it's the wrong place or
any other error. I'm relatively new to Linux and supernew to Octave. I have
checked for similar post as best I could.

Thanks,
CFB



    _______________________________________________________

File Attachments:


-------------------------------------------------------
Date: fre  5 okt 2012 14.55.14  Name: exp2.m  Size: 644 B   By: None

<http://savannah.gnu.org/bugs/download.php?file_id=26703>
-------------------------------------------------------
Date: fre  5 okt 2012 14.55.14  Name: gaussian.m  Size: 160 B   By: None

<http://savannah.gnu.org/bugs/download.php?file_id=26704>
-------------------------------------------------------
Date: fre  5 okt 2012 14.55.14  Name: fitexp2.m  Size: 191 B   By: None

<http://savannah.gnu.org/bugs/download.php?file_id=26705>

    _______________________________________________________

Reply to this item at:

  <http://savannah.gnu.org/bugs/?37508>

_______________________________________________
  Meddelandet skickades via/av Savannah
  http://savannah.gnu.org/