【急！！！】求解代码哪里出问题了，附题
已知某一400mm普通机床，电机功率为7.5kw，主轴的内径d-=48mm，作用在主轴外伸端处的载荷F=15000N，许用挠度=0.124mm，需用切应力=240MPa，允许转角=0.01rad。主轴材料的密度rou=7.8x10^-6kg/mm^3，主轴材料的弹性模量为E=2x10^5MPa,剪切模量为G=8x10^5MPa。主轴的计算转速为50r/min，经过计算主轴的最大输出扭矩为T=1432Nm。根据简历的主轴优化模型，结合Matlab优化工具箱对主轴结构进行优化设计。
下面是代码：
(1）编写目标函数的m文件：objfun.m,返回x处的函数值f
function f=objfun(x)
f=6.126*(10e-006)*(x(2)+x(3))(x(1)^2-48^2);
(2)因设计约束含非线性约束，需编写一个描述非线性约束的nonlinconstr.m文件：
function =nonlinconstr(x)
g(1)=0.5093*x(2)^2*(x(2)+x(3))/(x(1)^4-49^4)-0.124;
g(2)=7162.5/(x(1)^3-48^3)-240;
g(3)=0.5093*x(1)*x(2)/(x(1)^4-48^4)-0.01;
ceq=[];
(3)给定变量的初值，并调用优化函数：
x0=;
A=
b=[-200;650;-50;180;-80;150]
Ib=
options=optimset('display','iter','latgescale'.'off');
=fmincon('objfun',x0,A,b,[],[],Ib,[],'nonlinconstr',options)
但是输入进去有错误，求高手指出错误，十万火急啊！！！！！！！
页:
[1]