QQ登录

只需要一步,快速开始

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

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

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

2

主题

0

听众

4

积分

升级  80%

该用户从未签到

自我介绍
学理的青年一枚~
跳转到指定楼层
1#
发表于 2013-6-23 14:00 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
想要优化一个函数,程序如下:. ^9 J6 J* R0 P7 t# I  u
%给定初值- @7 |8 P  [& `! I# W
cm=3;& c; b5 O& X. g) L
cp=5;) d# o- G. t2 o3 p- C( i9 B! b
l=20;4 @6 I6 G# ]- |7 M/ Z5 \$ L. N3 Y+ J* C
h=0.05;" r! A! ?8 _) i% |
p0=1/h;
8 b7 r$ I. F: N: j4 Dq0=2*l/h;
/ L: B2 g3 U; n4 d8 |: inn=0;
4 t+ t" ^$ x( |. acc=cm*l*l
& \; p, X/ v: {5 z: G7 s% X" Q%主程序
5 x) H* |3 x( B6 }for p=1:h:p0* K* C) k( d  K2 [' }
    for q=1:h:q0
: y. k! p% V# w! K        a=p*h;* O4 X/ F6 U4 \" r, D2 G% U
        b=q*h;
% \) Q& _9 Q0 f8 ^. a        for n=1:12*l/b) %行向量/ n! g8 ~5 w1 o
            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;
8 {3 f, _5 r$ v" r4 _7 k            if c(2)>=c(1)' m: ^% p& N4 F) p0 G' v) }  J4 A
                nn=1;
/ \0 Y" R) z' P1 P+ C            elseif c(2*l/b)-c(2*l/b-1)<=0! E* e1 ~2 @. k$ D! E2 q
                nn=2*l/b;
1 H2 H& e4 C) `8 V            else
+ w' \' z' l( T6 P                for k=1:12*l/b)1 i+ t. X4 G5 S9 Z
                    if (c(k)-c(k-1)<=0 && c(k+1)-c(k)>=0)1 ?' C) @) s0 f& j2 z* s* \; G+ c
                        nn=k;& d5 U0 c! c) _: ^( G( X. m# Q! ?
                    end" D! _- X7 n7 |5 o/ t' s
                end
& o2 o2 Z7 Q6 g- A- A1 S' w            end
' O5 i& |0 n! I- C4 L) J; b9 p        end* H! c& @: N2 Q( e$ E5 N* m
        if(c(nn)<cc)) ~2 D' ^0 O& h) C* e0 x
            n=nn;5 [/ E, \% k6 \, T1 J! N
            aa=a;* k" l5 L9 R7 R/ U5 ~
            bb=b;2 z% A, N8 U. h6 b1 k
           cc=c(nn)
/ E) |+ T  ^' M. k  y+ B" y        end9 U2 b7 t: z1 z
    end1 K9 U: p! v8 S, j( c: t( C( @' Q2 f% E
end* e! c: p- }  y
fprinft('%d',n)' `9 X+ @) Q. X: A* O5 {
fprintf('%d',cc)" E9 v. }7 p7 c) u% O5 ?7 p8 P
fprintf('%d',aa)
9 `% b* i6 f& o; C1 sfprintf('%d',bb)
) G8 ?+ l7 c5 k6 n6 V/ |+ |  ?8 ^' b. a6 W
% y, b1 Q* g- p! c
运行不出来,错误提示:, X  O2 N! z8 E, a
cc =
4 U' e2 Y' H8 @0 a0 e9 }
4 \  [9 k- X/ g) W' z        1200/ T+ t5 u8 C3 I
+ @) d( V* l+ e, t. X  S4 p3 ~
??? Attempted to access c(2); index out of bounds because numel(c)=1.
- t  p5 `% p7 C1 J- ^1 o2 F: Q2 k7 M4 [3 N2 O

" X& }4 Y& \0 [这是什么问题呀......
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-6-9 08:43 , Processed in 0.317949 second(s), 60 queries .

回顶部