Question

I m试图做一些适合装饰。我有这样的职能:

function F = cdf_3p_model(a,data)
F=1-((1-a(5)-a(6)).*(exp(-abs(data)./a(1)))+((1-a(4)-a(6)).*(exp(-abs(data)./a(2))))+((1-a(4)-a(5)).*(exp(-abs(data)./a(3)))));

以及

function [a residual] = cdf_fit_3p(x,y)
a0 = [10 1 0.1 0.3 0.3 0.3];
lb = [0 0 0 0 0 0];
ub = [];
curvefitoptions = optimset( Display , final , MaxFunEvals ,100000, MaxIter ,50000);
[a, residual] = fmincon(@cdf_3p_model,a0,x,y,lb,ub,curvefitoptions);
end

我确定了初步参数、泡沫、灯泡,但我也如何声明:

a(1) > a(2) > a(3)
a(5) + a(6) +a(7) = 1

Answer 1

我认为,你更有机会使用诸如fmincon 。允许你具体说明否则可能无法做到的制约因素。你们可以通过把模型和数据差异的L2-规范化纳入最小的季度。

Answer 2

通常,我要说,“你的职能中的错误条款,如果这些条件得不到满足,就给真正的可怕的分数”。然而,你的条件使各种可允许的参数范围如此之广,我认为,如果你这样做,就会使电离板变得永远无法调和。我要说的是,“安全”不是你们的正确解决办法。

Answer 3

你们必须从更适于马塔雷克的一套参数中计算你们的参数。

例如,你可以改写仪式。

a(1) > a(2) > a(3)
a(5) + a(6) + a(7) = 1

页: 1

a(3) = p(1)
a(2) = p(1) + p(2)
a(1) = p(1) + p(2) + p(3)
a(4) = p(4)
a(5) = p(5)
a(6) = p(6)
a(7) = 1 - p(5) - p(6)

with

lb = [0 0 0 0 0 0]
ub = [Inf Inf Inf Inf 1 1]

Well, it s not perfect, because it allows a(7) 页: 1 low 页: 1 -1 instead of 0. But it includes your other constraints.

友情链接