MATLAB 曲线拟合（优化）

function f=fun(x)
global K L x thetamax alpha
for i=1:61
f=0; %函数值初始化
betae=atan(tan(alpha(i))/(1-(K/L)*tan(alpha(i)))); %计算转角期望值
A(i)=2*x(1).^2*sin(x(1)+alpha(i)); %计算实际转角
B(i)=2*K*x(1)-2*x(1).^2*cos(x(2)+alpha(i));
C(i)=2*x(1).^2-4*x(1).^2*(cos(x(2)).^2+4*K*x(1)*cos(x(2))-2*K*x(1)*cos(x(2)+alpha(i)));
theta3=2*acot((A(i)+sqrt(A(i).^2+B(i).^2-C(i).^2))/B(i)+C(i));
beta(i)=x(2)+theta(i)-pi;
if alpha(i)<=pi/18 %计算目标函数值
f(i)=1.5*abs(beta(k)-betae(k));
elseif alpha(i)<=pi/180;f(i)=abs(beta(k)-betae(k));
else f(i)=0.5*abs(beta(k)-betae(k));
end
end

global K L thetamax alpha
K=input( '输入两主销中心线间距(单位mm)K='); %提示输入结构参数
L=input('输入轴距(单位mm)L=');
thetamax=input( '输入外转向轮最大转角(单位为度)');
x0(1)=input( '输入初始点的第1个分量(臂长、mm)');
x0(2)=input('输入初始点的第2个分

matlab曲线拟合 怎样用matlab进行曲线拟合？ 求教matlab曲线拟合 [求助] 用MATLAB怎么做曲线拟合（高手留步） 如何用MATLAB实现曲线拟合？ MATLAB的曲线拟合具体实例 用matlab曲线拟合时遇到的问题 matlab线性优化 matlab工具箱优化问题 matlab求解优化问题