prince00840 发表于 2009-6-6 11:40

求助二次规划优化

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

xiezhh 发表于 2009-6-6 12:27

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

prince00840 发表于 2009-6-6 12:30

回复 沙发 xiezhh 的帖子

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

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

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

prince00840 发表于 2009-6-11 21:19

%先建立目标函数文件，并保存为myfun.m：function f = myfun(x)
x=';
eta=0.05;gamma=10;
A=[1 -1 1 000000 0 -1 0
   0 1 0 -1 -1 -1 -1 -11 0 1 -1
   0 0 -1 00000 -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 04.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=18.0803;x=55.568

Aeq=';

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

prince00840 发表于 2009-6-27 17:13

就是转化成标准格式的时候不会:@(
页: [1]
查看完整版本: 求助二次规划优化