声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1403|回复: 2

[工具箱] 一道优化的题,过早的退出,

[复制链接]
发表于 2008-8-27 20:55 | 显示全部楼层 |阅读模式

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

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

x
function maxtrans = transmissibility_zk(zk)
m1=260.8;
m2=14.5;
k1=334259;
z1 =0.03;
c1=2*sqrt(m1*k1)*z1;
z2=zk(1);
k2=zk(2);
c2=2*sqrt(k2*m2)*z2;
w1=0.5*2*pi;
w2=10*2*pi;
w = w1:(w2-w1)/500:w2;
T =( (c1*j*w + k1) .* (m2*(j*w).*2 + c2*j*w + k2) ) ./...
    ( (m1*(j*w).^2 +( c1 +c2 )*j*w +k1 + k2).*(m2*(j*w).^2 + c2*j*w +k2) - ( c2*j*w + k2).^2 );
% figure
% plot(w,abs(T))
% figure
% plot(w,angle(T)*180/pi)
maxtrans=max(abs(T));
这个是目标函数
以下是调用语句

[x,fval,exitflag,output,lambda,grad,hessian]=fmincon(@transmissibility_zk,[0.1,13865],[ ],[ ],[ ],[ ],[0.1,10000],[0.7,100000])
结果没有得到要求的目标值,
只迭代了6次就退出了
Optimization terminated: magnitude of directional derivative in search
direction less than 2*options.TolFun and maximum constraint violation
  is less than options.TolCon.
No active inequalities.

x =
  1.0e+004 *
   0.000036464769068   1.386498505477171

fval =
  11.592576914326367

exitflag =
     5

output =
       iterations: 6
        funcCount: 23
     lssteplength: 1
         stepsize: 0.001372954668234
        algorithm: 'medium-scale: SQP, Quasi-Newton, line-search'
    firstorderopt: 9.591579437255859e-004
          message: [1x172 char]

lambda =
         lower: [2x1 double]
         upper: [2x1 double]
         eqlin: [0x1 double]
      eqnonlin: [0x1 double]
       ineqlin: [0x1 double]
    ineqnonlin: [0x1 double]

grad =
  1.0e-003 *
   0.959157943725586
   0.537485388391948

hessian =
  18.401867691914799  -0.004367311374674
  -0.004367311374674   0.999885314836522
回复
分享到:

使用道具 举报

发表于 2008-8-27 21:01 | 显示全部楼层
在现有的优化限制条件下没有解,改一下限制条件吧
发表于 2008-8-29 11:00 | 显示全部楼层
如果用几个参数做fmincon优化,应该将这些优化求解的参数的数量级用新变量替换为一致,
在进行迭代的函数中还原计算。否则这个优化函数将得不到正确的优化解,你只会看到
一个变量有变化,另一个几乎不变。
用fminsearch优化函数就没有这个数量级的限制,但是不能提供上下界约束。

评分

1

查看全部评分

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

本版积分规则

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

GMT+8, 2024-12-26 01:51 , Processed in 0.179810 second(s), 21 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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