asufit.c

ASUFIT－Matlab－全局拟合程序

/*
 * MATLAB Compiler: 2.0.1
 * Date: Tue May 08 21:28:20 2001
 * Arguments: "-B" "sgl" "-m" "-W" "mainhg" "-L" "C" "asufit.m" "absfun.m"
 * "absfunfree.m" "absfunshift.m" "absfunwidth.m" "dispfit.m" "dispfitdisp.m"
 * "dispfitfree.m" "dispfitwidth.m" "seqmodfree.m" "spcfun.m" "spcfun1.m"
 * "atamult.m" "aprecon.m" 
 */
#include "Asufit.h"
#include "asufit.h"
#include "cla.h"
#include "colorbar.h"
#include "inputdlg.h"
#include "listdlg.h"
#include "lsqnonlin.h"
#include "shading.h"
#include "surf.h"
#include "title.h"
#include "xlabel.h"
#include "xlim.h"
#include "ylabel.h"
#include "ylim.h"
#include "zlabel.h"
#include "zlim.h"

extern mxArray * OPTIONS;
extern mxArray * OPT_STOP;
extern mxArray * S;
extern mxArray * U;
extern mxArray * V;
extern mxArray * abs3dAxes;
extern mxArray * abs3dFig;
extern mxArray * abs3dView;
extern mxArray * absAbsorbBox;
extern mxArray * absAbsorbText;
extern mxArray * absAllTaoBox;
extern mxArray * absAllTaoText;
extern mxArray * absCurrentTimeBox;
extern mxArray * absCurrentTimeText;
extern mxArray * absCurrentWlBox;
extern mxArray * absCurrentWlText;
extern mxArray * absDataGrid1;
extern mxArray * absDataGrid2;
extern mxArray * absDisper1Hypertext;
extern mxArray * absDisper2Hypertext;
extern mxArray * absDisper3Hypertext;
extern mxArray * absDisperCoeffsBox;
extern mxArray * absDisperCoeffsText;
extern mxArray * absDisperText;
extern mxArray * absDispersion;
extern mxArray * absErrorBox;
extern mxArray * absErrorText;
extern mxArray * absFilename;
extern mxArray * absInd1;
extern mxArray * absInd2;
extern mxArray * absIntPlot1;
extern mxArray * absIntPlot2;
extern mxArray * absIntSetAxes1;
extern mxArray * absIntSetAxes2;
extern mxArray * absLeftEnd1;
extern mxArray * absLeftEnd2;
extern mxArray * absMaxItBox;
extern mxArray * absMaxItText;
extern mxArray * absNumAllTaoBox;
extern mxArray * absNumAllTaoText;
extern mxArray * absNumTaoBox;
extern mxArray * absNumTaoFixedBox;
extern mxArray * absNumTaoFixedText;
extern mxArray * absNumTaoText;
extern mxArray * absPlot1;
extern mxArray * absPlot2;
extern mxArray * absPlotKinSpec;
extern mxArray * absPulseWidthBox;
extern mxArray * absPulseWidthText;
extern mxArray * absResidGrid1;
extern mxArray * absResidGrid2;
extern mxArray * absResidPlot1;
extern mxArray * absResidPlot2;
extern mxArray * absRightEnd1;
extern mxArray * absRightEnd2;
extern mxArray * absShiftGuessBox;
extern mxArray * absShiftGuessText;
extern mxArray * absTaoAnsBox;
extern mxArray * absTaoAnsText;
extern mxArray * absTaoBox;
extern mxArray * absTaoFixedBox;
extern mxArray * absTaoFixedText;
extern mxArray * absTaoText;
extern mxArray * absTog1;
extern mxArray * absTog2;
extern mxArray * absYTicks;
extern mxArray * abscheckbox1;
extern mxArray * abscheckbox2;
extern mxArray * abscheckbox3;
extern mxArray * abscheckboxtext;
extern mxArray * abstimeaxis;
extern mxArray * abstimelin;
extern mxArray * action1button;
extern mxArray * action2button;
extern mxArray * action3button;
extern mxArray * ampClosebutton;
extern mxArray * ampDownbutton;
extern mxArray * ampGrid;
extern mxArray * ampInd;
extern mxArray * ampPlotAll;
extern mxArray * ampSpecFile;
extern mxArray * ampSpecPlot;
extern mxArray * ampSpecWindow;
extern mxArray * ampSwitchbutton;
extern mxArray * ampTextBox;
extern mxArray * ampUpbutton;
extern mxArray * amps;
extern mxArray * analysisType;
extern mxArray * average;
extern mxArray * backupTypes;
extern mxArray * chisq;
extern mxArray * components;
extern mxArray * data;
extern mxArray * dataBackup;
extern mxArray * dataFit;
extern mxArray * dataIntervalStatus;
extern mxArray * dataMax;
extern mxArray * dataMin;
extern mxArray * date;
extern mxArray * disperCoeffs;
extern mxArray * dispersionStatus;
extern mxArray * dispersionVector;
extern mxArray * dispvecguess;
extern mxArray * emit;
extern mxArray * excite;
extern mxArray * expFile;
extern mxArray * fTol;
extern mxArray * fit;
extern mxArray * fitError;
extern mxArray * fitconst;
extern mxArray * fitmodel;
extern mxArray * glob;
extern mxArray * horizRotBox;
extern mxArray * irf;
extern mxArray * irfFilenames;
extern mxArray * kineticsFile;
extern mxArray * leftInt;
extern mxArray * mainWindow;
extern mxArray * maxIt;
extern mxArray * name;
extern mxArray * numTao;
extern mxArray * numTaoFixed;
extern mxArray * numTimePoints1;
extern mxArray * numTimePoints2;
extern mxArray * numVecs;
extern mxArray * optionsMenu;
extern mxArray * optionsbutton;
extern mxArray * paramsAns;
extern mxArray * pathname;
extern mxArray * perform;
extern mxArray * plot3D;
extern mxArray * plotMenu;
extern mxArray * plotbutton;
extern mxArray * pulseWidth;
extern mxArray * restore;
extern mxArray * restoreName;
extern mxArray * revert;
extern mxArray * rightInt;
extern mxArray * shift;
extern mxArray * shiftGuess;
extern mxArray * spcAllTaoBox;
extern mxArray * spcAllTaoText;
extern mxArray * spcAmpAnsBox;
extern mxArray * spcAmpAnsText;
extern mxArray * spcDataGrid;
extern mxArray * spcDateBox;
extern mxArray * spcDateText;
extern mxArray * spcEmitBox;
extern mxArray * spcEmitText;
extern mxArray * spcErrorBox;
extern mxArray * spcErrorText;
extern mxArray * spcExciteBox;
extern mxArray * spcExciteText;
extern mxArray * spcFilenames;
extern mxArray * spcFitSave;
extern mxArray * spcGlobal;
extern mxArray * spcGlobalErrorBox;
extern mxArray * spcGlobalErrorText;
extern mxArray * spcInd;
extern mxArray * spcIntPlot;
extern mxArray * spcIntSetAxes;
extern mxArray * spcLeftEnd;
extern mxArray * spcMaxItBox;
extern mxArray * spcMaxItText;
extern mxArray * spcNumAllTaoBox;
extern mxArray * spcNumAllTaoText;
extern mxArray * spcNumTaoBox;
extern mxArray * spcNumTaoFixedBox;
extern mxArray * spcNumTaoFixedText;
extern mxArray * spcNumTaoText;
extern mxArray * spcPlot;
extern mxArray * spcPlotMode;
extern mxArray * spcPromptBox;
extern mxArray * spcPromptText;
extern mxArray * spcResidGrid;
extern mxArray * spcResidPlot;
extern mxArray * spcRightEnd;
extern mxArray * spcShiftNumBox;
extern mxArray * spcShiftNumText;
extern mxArray * spcSpectrumBox;
extern mxArray * spcSpectrumText;
extern mxArray * spcTaoAnsBox;
extern mxArray * spcTaoAnsText;
extern mxArray * spcTaoBox;
extern mxArray * spcTaoFixedBox;
extern mxArray * spcTaoFixedText;
extern mxArray * spcTaoText;
extern mxArray * spcTog;
extern mxArray * spcshift;
extern mxArray * spcshifttext;
extern mxArray * spctpchBox;
extern mxArray * spctpchText;
extern mxArray * spec;
extern mxArray * spectrumFile;
extern mxArray * subtract;
extern mxArray * svdAnal;
extern mxArray * svdErrorBox;
extern mxArray * svdErrorText;
extern mxArray * svdGlobalTogHypertext;
extern mxArray * svdGlobalToggle;
extern mxArray * svdGrid1;
extern mxArray * svdGrid2;
extern mxArray * svdGrid3;
extern mxArray * svdNumVecsBox;
extern mxArray * svdNumVecsText;
extern mxArray * svdSurfNew;
extern mxArray * svdSurfOld;
extern mxArray * svdTimeBox;
extern mxArray * svdTimeText;
extern mxArray * svdWavelengthBox;
extern mxArray * svdWavelengthText;
extern mxArray * svdWeightBox;
extern mxArray * svdWeightText;
extern mxArray * svdbutton;
extern mxArray * t;
extern mxArray * tao;
extern mxArray * taoAns;
extern mxArray * taoFixed;
extern mxArray * textBox;
extern mxArray * tpch;
extern mxArray * tzero;
extern mxArray * vertRotBox;
extern mxArray * wavelength;
extern mxArray * weight;
extern mxArray * xTol;

static double __Array0_r[4] = { .05, .05, .9, .9 };

static double __Array1_r[3] = { 0.8, 0.8, 0.8 };

static double __Array2_r[4] = { 0.06, .3, 0.88, 0.55 };

static double __Array3_r[4] = { 0.06, .2, 0.88, 0.07 };

static double __Array4_r[4] = { 0.05, .37, 0.425, 0.47 };

static double __Array5_r[4] = { 0.525, .37, 0.425, 0.47 };

static double __Array6_r[4] = { 0.05, .12, 0.425, 0.2 };

static double __Array7_r[4] = { 0.525, .12, 0.425, 0.2 };

static double __Array8_r[4] = { .07, .45, 0.4, 0.42 };

static double __Array9_r[4] = { .53, .45, 0.4, 0.42 };

static double __Array10_r[4] = { 0.05, .47, 0.35, 0.37 };

static double __Array11_r[4] = { 0.45, .47, 0.35, 0.37 };

static double __Array12_r[4] = { 0.45, .05, 0.35, 0.37 };

static double __Array13_r[4] = { .07, .45, 0.86, 0.42 };

static double __Array14_r[3] = { 0.8, 0.8, 0.8 };

static double __Array15_r[3] = { 0.0, 0.0, 0.0 };

static double __Array16_r[4] = { .03, .92, .57, .05 };

static double __Array17_r[3] = { 0.8, 0.8, 0.8 };

static double __Array18_r[3] = { 0.0, 0.0, 0.0 };

static double __Array19_r[4] = { .6, .94, .1, .03 };

static double __Array20_r[4] = { .7, .94, .1, .03 };

static double __Array21_r[4] = { .8, .94, .1, .03 };

static double __Array22_r[4] = { .9, .94, .1, .03 };

static double __Array23_r[3] = { 0.8, 0.8, 0.8 };

static double __Array24_r[3] = { 0.0, 0.0, 0.0 };

static double __Array25_r[4] = { .6, .91, .1, .03 };

static double __Array26_r[4] = { .7, .91, .1, .03 };

static double __Array27_r[4] = { .8, .91, .1, .03 };

static double __Array28_r[4] = { .9, .91, .1, .03 };

static double __Array29_r[3] = { 0.8, 0.8, 0.8 };

static double __Array30_r[3] = { 0.0, 0.0, 0.0 };

static double __Array31_r[4] = { .05, .14, .1, .03 };

static double __Array32_r[3] = { 0.8, 0.8, 0.8 };

static double __Array33_r[3] = { 0.0, 0.0, 0.0 };

static double __Array34_r[4] = { .05, .11, .1, .03 };

static double __Array35_r[3] = { 0.8, 0.8, 0.8 };

static double __Array36_r[3] = { 0.0, 0.0, 0.0 };

static double __Array37_r[4] = { .15, .14, .1, .03 };

static double __Array38_r[3] = { 0.8, 0.8, 0.8 };

static double __Array39_r[3] = { 0.0, 0.0, 0.0 };

static double __Array40_r[4] = { .15, .11, .1, .03 };

static double __Array41_r[3] = { 0.8, 0.8, 0.8 };

static double __Array42_r[3] = { 0.0, 0.0, 0.0 };

static double __Array43_r[4] = { .25, .14, .1, .03 };

static double __Array44_r[3] = { 0.8, 0.8, 0.8 };

static double __Array45_r[3] = { 0.0, 0.0, 0.0 };

static double __Array46_r[4] = { .25, .11, .1, .03 };

static double __Array47_r[3] = { 0.8, 0.8, 0.8 };

static double __Array48_r[3] = { 0.0, 0.0, 0.0 };

static double __Array49_r[4] = { .35, .14, .1, .03 };

static double __Array50_r[3] = { 0.8, 0.8, 0.8 };

static double __Array51_r[3] = { 0.0, 0.0, 0.0 };

static double __Array52_r[4] = { .35, .11, .1, .03 };

static double __Array53_r[3] = { 0.8, 0.8, 0.8 };

static double __Array54_r[3] = { 0.0, 0.0, 0.0 };

static double __Array55_r[4] = { .45, .14, .1, .03 };

static double __Array56_r[3] = { 0.8, 0.8, 0.8 };

static double __Array57_r[3] = { 0.0, 0.0, 0.0 };

static double __Array58_r[4] = { .45, .11, .1, .03 };

static double __Array59_r[3] = { 0.8, 0.8, 0.8 };

static double __Array60_r[3] = { 0.0, 0.0, 0.0 };

static double __Array61_r[4] = { .55, .14, .1, .03 };

static double __Array62_r[3] = { 0.8, 0.8, 0.8 };

static double __Array63_r[3] = { 0.0, 0.0, 0.0 };

static double __Array64_r[4] = { .55, .11, .1, .03 };

static double __Array65_r[3] = { 0.8, 0.8, 0.8 };

static double __Array66_r[3] = { 0.0, 0.0, 0.0 };

static double __Array67_r[4] = { .65, .14, .1, .03 };

static double __Array68_r[3] = { 0.8, 0.8, 0.8 };

static double __Array69_r[3] = { 0.0, 0.0, 0.0 };

static double __Array70_r[4] = { .65, .11, .1, .03 };

static double __Array71_r[3] = { 0.8, 0.8, 0.8 };

static double __Array72_r[3] = { 0.0, 0.0, 0.0 };

static double __Array73_r[4] = { .75, .14, .1, .03 };

static double __Array74_r[3] = { 0.8, 0.8, 0.8 };

static double __Array75_r[3] = { 0.0, 0.0, 0.0 };

static double __Array76_r[4] = { .75, .11, .1, .03 };

static double __Array77_r[3] = { 0.8, 0.8, 0.8 };

static double __Array78_r[3] = { 0.0, 0.0, 0.0 };

static double __Array79_r[4] = { .85, .14, .1, .03 };

static double __Array80_r[3] = { 0.8, 0.8, 0.8 };

static double __Array81_r[3] = { 0.0, 0.0, 0.0 };

static double __Array82_r[4] = { .85, .11, .1, .03 };

static double __Array83_r[3] = { 0.8, 0.8, 0.8 };

static double __Array84_r[3] = { 0.0, 0.0, 0.0 };

static double __Array85_r[4] = { .05, .08, .25, .03 };

static double __Array86_r[3] = { 0.8, 0.8, 0.8 };

static double __Array87_r[3] = { 0.0, 0.0, 0.0 };

static double __Array88_r[4] = { .3, .08, .85, .03 };

static double __Array89_r[3] = { 0.8, 0.8, 0.8 };

static double __Array90_r[3] = { 0.0, 0.0, 0.0 };

static double __Array91_r[4] = { .05, .05, .25, .03 };

static double __Array92_r[3] = { 0.8, 0.8, 0.8 };