请教高手MATLAB约束非线性规划/一M文件的调试M文件如下,怎么调试呀,老有错误,% 两级斜齿轮传动中心距优化设计% 1---减速器中心距优化设计主程序% 设计变量的初始值x0=[2;4;18;20;6.4;10];% 设计变

来源:学生作业帮助网 编辑:六六作业网 时间:2024/12/23 09:11:26
请教高手MATLAB约束非线性规划/一M文件的调试M文件如下,怎么调试呀,老有错误,%两级斜齿轮传动中心距优化设计%1---减速器中心距优化设计主程序%设计变量的初始值x0=[2;4;18;20;6.

请教高手MATLAB约束非线性规划/一M文件的调试M文件如下,怎么调试呀,老有错误,% 两级斜齿轮传动中心距优化设计% 1---减速器中心距优化设计主程序% 设计变量的初始值x0=[2;4;18;20;6.4;10];% 设计变
请教高手MATLAB约束非线性规划/一M文件的调试
M文件如下,怎么调试呀,老有错误,
% 两级斜齿轮传动中心距优化设计
% 1---减速器中心距优化设计主程序
% 设计变量的初始值
x0=[2;4;18;20;6.4;10];
% 设计变量的下界与上界
lb=[2;3.5;14;16;5.8;8];
ub=[5;6;22;22;7;15];
% 使用多维约束优化命令fmincon(调用目标函数jsqyh_g和非线性约束函数jsqyh_g)
% 不定义线性不等式约束中设计变量的系数矩阵a=[]和常数项向量b=[]
% 没有等式约束,侧参数:系数矩阵Aeq=[]和常数项向量beq=[]
[x,fn]=fmincon(@jsqyh_f,x0,[],[],[],[],lb,ub,@jsqyh_g);
disp ' '
fprintf(1,' 高速级齿轮副模数 Mn1=%3.4f mm\n',x(1))
fprintf(1,' 低速级齿轮副模数 Mn2=%3.4f mm\n',x(2))
fprintf(1,' 高速级小齿轮齿数 z1=%3.4f \n',x(3))
fprintf(1,' 低速级小齿轮齿数 z3=%3.4f \n',x(4))
fprintf(1,' 高速级齿轮副传动比 i1=%3.4f \n',x(5))
fprintf(1,' 齿轮副螺旋角 beta=%3.4f 度\n',x(6))
fprintf(1,' 减速器中心距 a12=%3.4f mm\n',fn)
% 调用多维约束优化非线性约束函数(jsqyh_g)计算最优点x*的性能约束函数值g=jsqyh_g(x);
disp ' '
fprintf(1,' 高速级齿轮副解除疲劳强度约束函数值 g1=%3.4f \n',g(1))
fprintf(1,' 低速级齿轮副解除疲劳强度约束函数值 g2=%3.4f \n',g(2))
fprintf(1,' 高速级大齿轮齿根弯曲强度约束函数值 g3=%3.4f \n',g(3))
fprintf(1,' 低速级大齿轮齿根弯曲强度约束函数值 g4=%3.4f \n',g(4))
fprintf(1,' 大齿轮齿顶与轴不干涉几何约束函数值 g5=%3.4f \n',g(5))
>> % 2----两级斜齿轮减速器总中心距的目标函数(jsqyh_f)
function f=jsqyh_f(x);
hd=pi/180;
a1=x(1)*x(3)*(1+x(5));
a2=x(2)*x(4)*(1+31.5/x(5));
cb=2*cos(x(6)*hd);
f=(a1+a2)/cb;
% 3----两级斜齿轮减速器优化设计的非线性不等式约束函数(jsqyh_g)
function [g,ceq]=jsqyh_g(x);
hd=pi/180;
g(1)=cos(x(6)*hd)^3-3.079e-6*x(1)^3*x(3)^3*x(5);
g(2)=x(5)^2cos(x(6)*hd)^3-1.701e-4*x(2)^3*x(4)^3;
g(3)=cos(x(6)*hd)^2-9.939e-5*(1+x(5))*x(1)^3*x(3)^2;
g(4)=x(5)^2*cos(x(6)*hd)^2-1.076e-4*(31.5+x(5))*x(2)^3*x(4)^2;
g(5)=x(5)*(2*(x(1)+50)*cos(x(6)*hd)+x(1)*x(2)*x(3))-x(2)*x(4)*(31.5+x(5));

请教高手MATLAB约束非线性规划/一M文件的调试M文件如下,怎么调试呀,老有错误,% 两级斜齿轮传动中心距优化设计% 1---减速器中心距优化设计主程序% 设计变量的初始值x0=[2;4;18;20;6.4;10];% 设计变
改了几处,通了.
function hahah
% 两级斜齿轮传动中心距优化设计
% 1---减速器中心距优化设计主程序
% 设计变量的初始值
x0=[2;4;18;20;6.4;10];
% 设计变量的下界与上界
lb=[2;3.5;14;16;5.8;8];
ub=[5;6;22;22;7;15];
% 使用多维约束优化命令fmincon(调用目标函数jsqyh_g和非线性约束函数jsqyh_g)
% 不定义线性不等式约束中设计变量的系数矩阵a=[]和常数项向量b=[]
% 没有等式约束,侧参数:系数矩阵Aeq=[]和常数项向量beq=[]
[x,fn,exitflag]=fmincon(@jsqyh_f,x0,[],[],[],[],lb,ub,@jsqyh_g)
disp '
高速级齿轮副模数 Mn1=2.0461 mm
低速级齿轮副模数 Mn2=3.6059 mm
高速级小齿轮齿数 z1=18.5158
低速级小齿轮齿数 z3=16.0000
高速级齿轮副传动比 i1=5.8000
齿轮副螺旋角 beta=8.0000 度
减速器中心距 a12=317.4186 mm
g =
0.0000 -0.0000 -1.0042 -15.1854 -761.8045
高速级齿轮副解除疲劳强度约束函数值 g1=0.0000
低速级齿轮副解除疲劳强度约束函数值 g2=-0.0000
高速级大齿轮齿根弯曲强度约束函数值 g3=-1.0042
低速级大齿轮齿根弯曲强度约束函数值 g4=-15.1854
大齿轮齿顶与轴不干涉几何约束函数值 g5=-761.8045

请教高手MATLAB约束非线性规划/一M文件的调试%1 目标函数function f=abfun(x);f=pi*130*(x(1)-0.75*x(3))^2/4+pi*[80^2-(x(1)-0.75*x(3))^2]*x(2)/4;%2 非线性约束funtion [c,ceq]=abfung(x);c(1)=40000-[0.75*pi*(x(1)-0.75*x(3))*x(2)]*10; 请教高手MATLAB约束非线性规划/一M文件的调试M文件如下,怎么调试呀,老有错误,% 两级斜齿轮传动中心距优化设计% 1---减速器中心距优化设计主程序% 设计变量的初始值x0=[2;4;18;20;6.4;10];% 设计变 matlab 非线性规划问题 求解约束非线性规划中使用MATLAB时为什么要设定初始值X0? MATLAB基础问题麻烦高手帮个忙,请问解有约束非线性规划问题,用MATLAB中的哪个函数?再麻烦各位把函数的调用格式也写清楚吧, 如何利用MATLAB解非线性规划 有约束的非线性规划怎么求解 运筹学 matlab一道有约束的非线性规划问题max gs.t.(a-sqrt(a^2-4*g^2*x^2))/(2*b*g*x) 请教关于用Matlab在非线性约束条件下的最优解在以下约束条件下cos(5*x)+cos(5*y)+cos(5*z)+cos(5*m)=0;cos(7*x)+cos(7*y)+cos(7*z)+cos(7*m)=0;cos(11*x)+cos(11*y)+cos(11*z)+cos(11*m)=0;使M=cos(x)+cos(y)+cos(z)+cos(m)最大的求 求高手给下面非线性规划的lingo的代码 matlab fmincon函数用法高手帮我看下这段有什么问题,为什么运行不出来目标函数和非线性约束不能放在一个M文件里面么?function f = myfun(x)f = 0.192457*1e-4*(x(2)+2)*x(1)^2*x(3);function [c,ceq] = mycon(x)c(1)=350 请教matlab高手,线性拟合 非线性约束规划 我编制了 目标函数和约束函数 为什么在matlab中还是出错 求大神看看这不给力的function f=myfun(x)f=0.192457*le-4*(x(2)+2)*x(1)^2*x(3);function [c,cep]=mycon(x)c(1)=350-163*x(1)^(-2.86)*x(3)^0.86;c(2)=1 如何用matlab求解非线性规划问题 最好给几个例子 机会约束规划模型 一道关于Matlab非线性约束问题的求解,需要有m文件的代码具体的我会上传附件,我主要是忘记怎么编了,我星期一要去讲的 matlab 函数 fmincon求解非线性规划,求这个函数的最小值 function f = myfun(x)f = (pi*(x(3)+2)*x(2)*x(1)^2)/4;运行后返回(约束在最后面有写)Warning:Trust-region-reflective method does not currently solve this type of prob matlab 非线性约束求最小值求目标函数f 的最小值,约束条件为 x>0,y>0,de