QQ登录

只需要一步,快速开始

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

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

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

2

主题

0

听众

4

积分

升级  80%

该用户从未签到

自我介绍
学理的青年一枚~
跳转到指定楼层
1#
发表于 2013-6-23 14:00 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
想要优化一个函数,程序如下:7 ^( i/ }* \7 f% Y$ T  S* U
%给定初值
1 ?5 o" O) g+ l) ^, s' rcm=3;
: ?' P5 l1 _$ ^5 V  _cp=5;
3 f" O# ~4 b7 l) s' m. C  {l=20;5 g, t2 b" c  [5 q7 ^6 P
h=0.05;
3 B+ B0 B% ?# a1 {( L- r# U; Dp0=1/h;
- I/ _9 Q( B& ~" i# E) C) n1 mq0=2*l/h;
8 K& W6 r! T! u. V6 L& Qnn=0;
6 E/ Z  H4 a- m4 B+ M& c' Occ=cm*l*l' h1 U# n5 i' N2 j4 r3 W
%主程序4 H: V1 j! K3 i
for p=1:h:p00 A, `: y1 L1 Y& c/ p+ i3 ?0 X' c4 I
    for q=1:h:q0
' y( V2 Y' T4 c: w, F        a=p*h;5 ^1 E) W$ y3 m; d
        b=q*h;8 b2 i' C( K* O; J2 V" N: S
        for n=1:12*l/b) %行向量
& V" a9 o$ B" A) R! D. o: A6 w            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 z) p& x5 ?1 E" q5 \            if c(2)>=c(1)
/ Z! `) i5 F5 _* A                nn=1;9 X6 E3 |- F7 H* o$ v+ k' L8 {
            elseif c(2*l/b)-c(2*l/b-1)<=0  j/ ?! C* F6 ^. W* m& }; n; o, U
                nn=2*l/b;
  o& k! Z- E" W5 N$ D! @' R2 ^1 a            else. |; ~  {8 N/ A) r& p) b3 R/ Z9 E; J
                for k=1:12*l/b)
; e6 u, k1 T' A/ X5 }- T$ C" @                    if (c(k)-c(k-1)<=0 && c(k+1)-c(k)>=0)
+ m5 ?' ^. c& \2 m' Z9 t" l! C                        nn=k;4 [& R. p* A& u" q& ^8 C
                    end
, I" _) k0 i4 s1 Z1 Q) T                end
' ^8 N! K( y8 l5 \/ b- @( U9 \            end+ \) z3 L- L% s, Q- q
        end) L# S- C3 \$ Z1 P7 p, z4 p
        if(c(nn)<cc)
6 z$ h$ t9 i( R9 @3 g( P# ^' k            n=nn;4 C3 F% [" F" j; j! T3 A
            aa=a;
: J5 e& i9 A: h2 P+ f            bb=b;
4 ^8 B  A7 y% _1 @; V; ~9 P           cc=c(nn)3 f6 |- T1 m! b  p$ H7 V/ e0 t
        end+ n$ F8 ?! E0 Z$ Y4 B( d
    end
6 I4 E: I5 G! U$ G# }end
  `: u; d4 F' S/ e8 _/ f0 Kfprinft('%d',n)
+ d" _( T! e9 M& Cfprintf('%d',cc)
& `- }9 M4 B* `: _; cfprintf('%d',aa)
% ^  R+ T* y$ c* ~fprintf('%d',bb)7 J) Z  G( i' h  ~

2 a# ^( B& N9 r/ }" o6 e" q. E6 w( t$ _  D4 g8 j% i
运行不出来,错误提示:
" |  M! n% d8 R6 {1 ecc =
% k# l' s: R; f0 P. O! S  R- }7 [, s7 A
        1200
- ]( Z3 a3 r" C7 k: e$ \' V! H: W5 u5 K9 `& F. h2 g4 s
??? Attempted to access c(2); index out of bounds because numel(c)=1.
9 ?6 `, V8 T( g' I0 n" b/ @! s% P6 }( Q4 I( k& U; ~- @" y8 G
& T( z% M1 l' s$ T1 F; b
这是什么问题呀......
zan
转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信
2#
无效楼层,该帖已经被删除
您需要登录后才可以回帖 登录 | 注册地址

qq
收缩
  • 电话咨询

  • 04714969085
fastpost

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

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

蒙公网安备 15010502000194号

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

GMT+8, 2025-8-17 02:28 , Processed in 1.743806 second(s), 64 queries .

回顶部