main.m

有约束2a+b^2+tan(c)=15 目标函数a*exp(b/x)+c=y 的最小二乘解

function [ x ] = main(  )
%MAIN Summary of this function goes here
%   Detailed explanation goes here
global data1
global data2
global theda
data1=[0.2 1 2 3 5 7 11 15]
data2=[5.05 8.88 11.63 12.93 14.15 14.73 15.30 15.60]
x=[11.362697524217458  -1.077276645684700   5.000409076117723]
P=1
theda=1
c=10
k=0
while theda*P>0.001
k=k+1
alpha=0.01
beta=10
d=[1 1 1]
X=1
f_i=1
while (X'*f_i')'*(X'*f_i')>0.01
    X=A(x)
    f_i=f(x)
    F_1=SumF(x)
    d=-(X'*X+alpha*eye)^(-1)*X'*f_i'
    x=x+d'
    F_2=SumF(x)
    if F_1>F_2
        alpha=alpha/beta
    else
        alpha=alpha*beta
    end
end
P=P_k(x)
theda=theda*c
end