声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1010|回复: 0

[编程技巧] 【急!!!】求解代码哪里出问题了,附题

[复制链接]
发表于 2012-5-21 12:18 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
已知某一400mm普通机床,电机功率为7.5kw,主轴的内径d-=48mm,作用在主轴外伸端处的载荷F=15000N,许用挠度[y]=0.124mm,需用切应力[t]=240MPa,允许转角[sita]=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 [g,ceq] =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=[140,130,550];
A=[0,0,-1;0,0,1;-1,0,0;1,0,0;0,-1,0;0,1,0]
b=[-200;650;-50;180;-80;150]
Ib=[0,0,0]
options=optimset('display','iter','latgescale'.'off');
[x,fval,exitflag,ouput]=fmincon('objfun',x0,A,b,[],[],Ib,[],'nonlinconstr',options)

但是输入进去有错误,求高手指出错误,十万火急啊!!!!!!!
回复
分享到:

使用道具 举报

您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-29 05:50 , Processed in 0.058813 second(s), 18 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表