a=[ ];%建立a的一个数组
x=[650 340 400 800 300 600 720 480 440 540 750];
y=[2.48 4.45 4.52 1.38 4.65 2.96 2.18 4.04 4.2 3.1 1.5];
q=length(x);
for i=1:q-1
if x(i)>500
a(i)=1;
else
a(i)=0;
end
end
%%
得出
a =
1 0 0 1 0 1 1 0 0 1 1
利用1stopt进行拟合,输入:
Parameters b1,b2,b3;
Variable x,a,y;
Function y=b1+b2*x+b3*(x-500)*a;
data;
6.5 1 2.48
340 0 4.45
400 0 4.52
800 1 1.38
300 0 4.65
600 1 2.96
720 1 2.18
480 0 4.04
440 0 4.2
540 1 3.1
750 1 1.5
结果
迭代数: 20
计算用时(时:分:秒:微秒): 00:00:00:281
优化算法: 麦夸特法(Levenberg-Marquardt) + 通用全局优化法
计算结束原因: 达到收敛判断标准
均方差(RM