求助二次规划优化

prince00840

怎么用序列二次规划，我查资料说要化成标准形式
help quadprog
不知道怎么把图中的式（9）化成二次规划的标准形式，能有高手明白吗？请您赐教:@)
就是问matlab编程时怎么变

xiezhh

求最大似然估计可以试一试mle函数。

prince00840

是把极大似然估计转化成了序列二次规划来求:lol

给出一组数据，怎样检测她的粗差（大的误差）的位置呢

[ 本帖最后由 ChaChing 于 2009-6-27 19:28 编辑 ]

prince00840

%先建立目标函数文件，并保存为myfun.m：function f = myfun(x)
x=[190.26 174.596 3.139 32.7 33.47 7.25 0.316 100.24 28.629 23.8]';
eta=0.05;gamma=10;
A=[1 -1 1 0  0  0  0  0  0 0 -1 0
   0 1 0 -1 -1 -1 -1 -1  1 0 1 -1
   0 0 -1 0  0  0  0  0 -1 -1 0 1];
sigma=[
    6.1549 0 0 0 0 0 0 0 0 0
    0 0.4529 0 0 0 0 0 0 0 0
    0 0 0.8814 0 0 0 0 0 0 0
    0 0 0 0.1305 0 0 0 0 0 0
    0 0 0 0 0.2283 0 0 0 0 0
    0 0 0 0 0 0.0326 0 0 0 0
    0 0 0 0 0 0 0.0015 0 0 0
    0 0 0 0 0 0 0  4.318 0 0
    0 0 0 0 0 0 0 0 1.9852 0
    0 0 0 0 0 0 0 0 0 0.1262];
f=ln((1-eta)*exp(-0.5*((x(i)-y(i)).^2)/(sigma^2))+eta/gamma*exp(-0.5*((x(i)-y(i))/(gamma*sigma)).^2));
y[0]=18.0803;x[0]=55.568

Aeq=[1,-1,-1,-1,-1,-1,-1,]';

------------------------------------
[x,fval,exitflag] = quadprog(HH,f,A,b,Aeq,beq,lb,ub,[],options)
:@(

prince00840

就是转化成标准格式的时候不会:@(