QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 3836|回复: 0
打印 上一主题 下一主题

[问题求助] matlab问题 求大家指点

[复制链接]
字体大小: 正常 放大
菠萝饭        

2

主题

0

听众

4

积分

升级  80%

该用户从未签到

自我介绍
学理的青年一枚~
跳转到指定楼层
1#
发表于 2013-6-23 14:00 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
想要优化一个函数,程序如下:
3 ~* c+ a" e& V, ~4 g7 ~0 H%给定初值9 ~  S/ ]5 E. U  R& y
cm=3;
2 d; f0 E0 j1 P7 [8 _  Dcp=5;
9 @; ?" \: f$ bl=20;- {3 v0 O2 ]( P
h=0.05;
# i: v8 L( Z1 N8 _p0=1/h;- `) W( u4 ^5 p/ J- n' `1 O& o
q0=2*l/h;
; `' q. ~) K5 y, y! U" ?nn=0;# u' D  H# X# e! X, _4 w# a$ B
cc=cm*l*l4 N9 v* y9 C6 P0 c+ u' T. F
%主程序$ b1 V6 ~" ~$ t3 U% Q
for p=1:h:p0# D, R& _- K$ L
    for q=1:h:q0, s/ d, f  ]3 g0 S. ?. X6 G3 K4 ^
        a=p*h;
: z8 d1 W1 ]; y+ ]/ w! c        b=q*h;
' m# X. W- m% \6 u        for n=1:12*l/b) %行向量
5 a2 z  U! n3 I* `" v$ T$ a9 S: r            c(n)=cm*b*b/4*((1-a^n)/(1-a)+2*n/(1-a)-2*(1-a^n)/(1-2*a+a*a))+cm*(a^n*(l-n*b/2)*(l-n*b/2)+(1-a^n)*b*(l-n*b/2)/(1-a))+n*cp;
2 c6 \/ K& U. ^; X, [& V            if c(2)>=c(1)  P# E  f2 ~5 k8 o* J0 X$ T4 b
                nn=1;
) ?3 M7 t9 }( l) \7 z4 w/ v            elseif c(2*l/b)-c(2*l/b-1)<=0( R' X% |6 F* F7 d3 {9 B1 l
                nn=2*l/b;( y3 M' A+ l: c3 T
            else! y$ g; Q  L1 p! _2 J
                for k=1:12*l/b)
' h, g" y* ^+ }  ?9 o                    if (c(k)-c(k-1)<=0 && c(k+1)-c(k)>=0)' [+ F5 P% v/ V
                        nn=k;
" w; N4 A$ g" y* w9 W" p  z6 R                    end
  c! K4 W+ v# }6 h                end$ h% L% |4 F9 x, H2 d
            end% ]1 N: e/ n! b! R) F% V
        end* F3 J4 Y8 _6 `+ i" A3 u
        if(c(nn)<cc)4 E  L2 }/ h7 {
            n=nn;
; k7 U6 i5 z0 Q6 R$ k' K) s; f            aa=a;+ s% g- B# |' l% `/ {6 [3 e9 z- O, j
            bb=b;
3 a4 j( ~! O( @7 T           cc=c(nn)0 t5 g- a  Q  ?2 V
        end
% h" q7 s, O: s    end) S% R5 b7 x( O* U
end
) j7 I+ D2 `3 |* qfprinft('%d',n)- i: J  L$ m. |5 b
fprintf('%d',cc)
$ K/ i! E; h( m. a; K' r( Zfprintf('%d',aa)
; ]6 u) C. w* w) e) Y: i! t+ jfprintf('%d',bb)5 A+ L7 l( A( q8 w( B

* x, ]) h  }" o5 e9 d1 f, M1 \5 O! ?3 n' x, g) q
运行不出来,错误提示:8 Y) R9 G0 E8 p- D# m, i& _1 E
cc =
0 {, q; m; E0 Q7 t9 x' E% ^
; C; b. Q* C1 {/ }0 v( H! m        1200. K$ O. E  b1 X5 ]: X
. U: W2 D9 O2 r1 u
??? Attempted to access c(2); index out of bounds because numel(c)=1.
0 I2 V% u9 P* G5 b- _! w- E3 e- H8 C: ~
8 q5 j, M$ g' H
这是什么问题呀......
zan
转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信
您需要登录后才可以回帖 登录 | 注册地址

qq
收缩
  • 电话咨询

  • 04714969085
fastpost

关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

手机版|Archiver| |繁體中文 手机客户端  

蒙公网安备 15010502000194号

Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

GMT+8, 2026-8-23 06:54 , Processed in 0.314128 second(s), 58 queries .

回顶部