小弟有一题请教！！！

网中剑

各位大仙：
     小弟有一题请教：
非线性模型的参数估计问题 .
方程为:  
y=a*Exp(-b/x1)*(1-x2)^c+d*Exp(-f/x1)*[(1-x2)^g]*x2^h
其中 a b c d f g h 为未知参数.
观测数据如下:
     y                X1              X2       
0.41        448.2        0.176       
0.9205        454.2        4.41e-002       
1.799        458.2        9.69e-002       
2.159        461.2        0.1322       
2.611        643.2        0.1806       
3.347        466.2        0.2555       
3.582        467.2        0.2996       
4.008        468.2        0.348       
4.10        470.2        0.3921       
4.276        471.2        0.4405       
       
初始值：a=    1.121379e+014     b =181534       
            c=0.9                          d= 5.011866e+016
           f=171632                     g=0.8                                   
           h=         1.2       







[URL=http://http://m156.163.com/cgi/ldmsapp/非线性模型的参数估计问题.doc?sid=PAUasBFWUUQAUvfg&mid=MX6kKExF71000-_UpFEk&part=2&funcid=readpart&msglevel=0&filename=非线性模型的参数估计问题.doc





[em05][em05][em05][em05][em05]

[此贴子已经被作者于2003-8-16 17:42:59编辑过]

xfgasa2

I don't know

dcyu

非线性的，挺难的，没有可以参考的数学模型，只能编程搜索，目标函数就是最小二乘值，用退火算法或遗传算法可以解决，不过程序很不好编。

golden

先建立个符号变量,y=a*exp(-b/x1)*(1-x2)^c+d*exp(-f/x1)*[(1-x2)^g]*x2^h
再用nlinfit(x,y,'表达式 ' ,beat());beat是初始系数.看帮助

chen017

yong ssps软件编

golden

y=[ .41 .9205 1.799 2.159 2.611 3.347 3.582 4.008 4.10 4.276]';
x=[448.2 0.176
    454.2 4.41e-002
    458.2 9.69e-002
    461.2 0.1322
    643.2 0.1806
    466.2 0.2555
    467.2 0.2996
    468.2 0.348
    470.2 0.3921
     471.2 0.4405 ];   
beat0=[    1.121379e+014  181534   0.9   5.011866e+016  171632  0.8   1.2 ];
beatfit=nlinfit(x,y,'mmm',beat0)
ｍｍｍ.ｍ子程序如下
function yy=model(beat0,x)
        a=beat0(1);
        b=beat0(2);
        c=beat0(3);
        d=beat0(4);
        f=beat0(5);
        g=beat0(6);
        h=beat0(7);      
        x1=x(:,1);
        x2=x(:,2);         
        yy=a*exp(-b/x1)*(1-x2).^c+d*exp(-f/x1)*[(1-x2).^g]*x2.^h;
答案
1.0e+016 *

    0.0112
    0.0000
    0.0000
    5.0119
    0.0000
    0.0000
    0.0000
不知道对不对，有什么高见告诉我呀！

golden

我也知道答案有问题，但没有找出哪儿出错了，哪位高手帮我指点一下。

网中剑

非常感谢!! golden
让我们继续努力！！