bugGNU Octave - Bugs: bug #37508, fminsearch misses in deconvolving...

 
 

bug #37508: fminsearch misses in deconvolving 2 gaussian peaks

Submitter:  None
Submitted:  Fri 05 Oct 2012 02:55:14 PM UTC
   
 
Category:  Octave Package Severity:  3 - Normal
Priority:  5 - Normal Item Group:  Matlab Compatibility
Status:  Fixed Assigned to:  None
Originator Name:  CFB Originator Email:  -email is unavailable-
Open/Closed:  * Closed Release:  * 3.6.1
Operating System:  * GNU/Linux Fixed Release:  None
Planned Release:  None
* Mandatory Fields

Add a New Comment Rich Markup
   

Fri 14 Feb 2014 02:52:30 AM UTC, comment #2: 

The fminsearch is now a function in Octave itself. I've run the attached exp2 script in Octave 3.8.0 and the results (a1m, m1m, s1m) match the input parameters (a1, m1, s1), so I think this one is fixed now.

If I've misunderstood the problem, or if you still have problems with fminsearch in Octave 3.8.0, please comment or submit a new bug report.

Mike Miller <mtmiller>
Group Member
Fri 05 Oct 2012 03:03:18 PM UTC, comment #1: 

This looks like a problem with an Octave Forge package, not a bug
in Octave.  The correct place to report this problem is the
Octave Forge mailing list, octave-dev@lists.sourceforge.net.

A copy of this update to the bug report is being sent to that
list so that they are aware of it.  Please follow up there.

Jordi GutiƩrrez Hermoso <jordigh>
Group Member
Fri 05 Oct 2012 02:55:14 PM UTC, original submission:  

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

Anonymous

 

(Note: upload size limit is set to 16384 kB, after insertion of the required escape characters.)

Attach Files:
   
   
Comment:
   

Attached Files
file #26703:  exp2.m added by None (644B - text/x-objcsrc)
file #26704:  gaussian.m added by None (160B - text/x-objcsrc)
file #26705:  fitexp2.m added by None (191B - text/x-objcsrc)

 

Depends on the following items: None found

Items that depend on this one: None found

 

Carbon-Copy List
  • -email is unavailable- added by mtmiller (Posted a comment)
  • -email is unavailable- added by jwe (Updated the item)
  • -email is unavailable- added by rik5 (Updated the item)
  • -email is unavailable- added by jordigh (Posted a comment)
  • -email is unavailable- added by jordigh
  • -email is unavailable- added by None (Submitted the item)
  •  

    There are 0 votes so far. Votes easily highlight which items people would like to see resolved in priority, independently of the priority of the item set by tracker managers.

    Only group members can vote.

     

    Follow 9 latest changes.

    Date Changed by Updated Field Previous Value => Replaced by
    2014-02-14 mtmiller StatusNone Fixed
        Open/ClosedOpen Closed
    2013-05-16 jwe CategoryNone Octave Package
        StatusOctave Forge None
    2012-10-12 rik5 StatusNone Octave Forge
    2012-10-05 jordigh Carbon-Copy- Added "octave forge list" <octave-dev@lists.sourceforge.net>
    2012-10-05 None Attached File- Added exp2.m, #26703
        Attached File- Added gaussian.m, #26704
        Attached File- Added fitexp2.m, #26705

    Back to the top

    Powered by Savane 3.13-758e.
    Corresponding source code