[交流] 求助：多变量曲线拟合，没有经验公式已有5人参与

3楼: Originally posted by dyhseu at 2015-03-22 16:47:42
您好！我是一名大二的学生，用1stopt 1.5有一个问题想请教您，因为刚刚注册不能发帖和密语，所以找到您，麻烦您帮忙解答下谢谢！~代码如下：
Constant     W=25000,b=0.3,D=1,R=0.5;                             ...

function solve_nonleqs

clear all;clc
format long
global i

ispan=0.1:0.05:0.9;
m=length(ispan);
for j=1:m
    i=ispan(j);
x0=0.5;
options=optimset('TolFun',1e-6);
[o(j) fval(j)] =fsolve(@eqs,x0,options);


W=25000;b=0.3;D=1;R=0.5;                                                      
c=5000;phi=(30/180*pi);n=0.8;kc=10000;kphi=1000000;k=0.015;  %土壤参数

                                                            %滑转率
                                                    %接触角，挂钩牵引力

sigma=(kc/b+kphi)*(0.5*R*o(j)^2)^n*(1-n/3);                    %驱动轮正应力
tau=c+sigma*tan(phi);                                                     %驱动轮剪应力
eta=sigma/tau;
zeta=i*R/k;


FDP(j)  =b*R*tau*(eta*(cos(o(j))-1)+sin(o(j))-1/(1+zeta^2)*(zeta*cos(o(j))+sin(o(j))-zeta*exp(-zeta*o(j)))) ;    %由o求FDP
end
o=o'  
FDP=FDP',
cancha=fval'

figure(1),plot(ispan,FDP,'MarkerEdgeColor',[1 0 0],'MarkerSize',12,'Marker','o',...
    'LineWidth',3);


% Create xlabel
xlabel('i','FontSize',24);

% Create ylabel
ylabel('FDP','FontSize',24);


function f=eqs(x)
global i

o=x;

W=25000;b=0.3;D=1;R=0.5;                                                      
c=5000;phi=(30/180*pi);n=0.8;kc=10000;kphi=1000000;k=0.015;  %土壤参数

                                                            %滑转率
                                                    %接触角，挂钩牵引力

sigma=(kc/b+kphi)*(0.5*R*o^2)^n*(1-n/3);                    %驱动轮正应力
tau=c+sigma*tan(phi);                                                     %驱动轮剪应力
eta=sigma/tau;
zeta=i*R/k;



f= 0.4*W- ( b*R*tau*((1-cos(o))+eta*sin(o)-1/(1+zeta^2)*(zeta*sin(o)-cos(o)+   exp(-zeta*o))) );       %已知W求o