QQ登录

只需要一步,快速开始

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

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

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

2

主题

0

听众

4

积分

升级  80%

该用户从未签到

自我介绍
学理的青年一枚~
跳转到指定楼层
1#
发表于 2013-6-23 14:00 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
想要优化一个函数,程序如下:; b- w! A/ W* e9 S; Y, Y
%给定初值$ S, z1 \) f) ^
cm=3;
( k, G# P: T/ L+ h3 y4 mcp=5;. @, e8 x' H  u3 n0 m
l=20;
2 W! q* }5 H5 N6 _/ p& k3 j8 sh=0.05;  \' f* Z. C* Q+ v; ~4 X  F9 V
p0=1/h;( @' ?' f# u: Y
q0=2*l/h;
* C! N. p8 B: n9 v; ]; Wnn=0;7 e' G+ @( ]. ?0 b& d6 P# Q3 h
cc=cm*l*l
- A1 j. [% ?; E1 n0 U& J+ K%主程序! R' l' \$ `5 Y
for p=1:h:p0
& K! g0 R8 H3 K6 m& w9 L% c  j    for q=1:h:q0. B' ^, g: w3 A+ h4 v: Y$ p
        a=p*h;
6 E) r2 ]( H% V  b0 t2 k        b=q*h;5 W* a6 Q" J# L6 P. D4 Y' `- I
        for n=1:12*l/b) %行向量- i4 i) k6 b- U* _  u) L3 {2 c  U
            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;) Y; u( e; \- r8 V4 O5 l, ~
            if c(2)>=c(1)2 a/ ]0 ^' ?2 j0 z7 F
                nn=1;
  k5 `. d+ L; }# w# W1 K# q, x            elseif c(2*l/b)-c(2*l/b-1)<=0& U* d% |6 z/ G; X" I5 Q5 U
                nn=2*l/b;+ \' ^" E2 v3 K8 u8 p6 n# o! H  b
            else
! P+ O, H1 |) I4 @% X                for k=1:12*l/b)! y2 A2 R0 a% ?
                    if (c(k)-c(k-1)<=0 && c(k+1)-c(k)>=0)
8 C$ r9 U1 D+ Y' n: ?                        nn=k;
* S, p* o, p. s2 [+ _5 e( y                    end
5 l$ ]9 N+ t  I# D5 z* R5 l                end
$ F' q6 Y' y( Y1 V: z7 ^            end2 w3 ]. P; Z1 O2 x* i5 T) W
        end
, ?3 ^' h* I) g        if(c(nn)<cc)" E3 `6 m; T) {- z1 s" D! G' v
            n=nn;/ a8 ~* |. x% N! ^- \. D8 d
            aa=a;1 Z; Q- ~. N. r4 V9 B* Z9 G
            bb=b;
  }" I. S! R/ {, @: Z; ]1 I           cc=c(nn)
/ M2 z# F( [, L        end
6 h. W* w" n  o' @/ a3 f+ j% D    end6 b( L$ p9 [0 i/ v8 i
end9 ]9 G- Y1 p2 L  n$ D
fprinft('%d',n)/ Z& d, m) |$ c2 L% o5 z) o
fprintf('%d',cc)+ K. C$ N; i. |( a4 G) s) M
fprintf('%d',aa)0 a4 _$ g; Z- _
fprintf('%d',bb)
2 C" U3 o. r$ n) W0 e( v- X* \8 q6 _' R. f( Q
2 D; Z* a+ E0 s  f2 `
运行不出来,错误提示:
; `! Z) o8 M: m6 H9 Mcc =/ A6 ~8 Z' z, a/ W. k! v1 v  }
0 P/ W2 h$ ~4 u6 P* c  a
        1200
2 [8 }: c4 X/ K- G$ x& w7 g& ~& Q1 h' d8 Z# @% h
??? Attempted to access c(2); index out of bounds because numel(c)=1.7 _* ]/ r6 q6 e  U* W9 D

$ D, B* H* a* e; ?, W9 v. W1 m
% F2 D& h/ a, A/ }1 `. s这是什么问题呀......
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-7-6 16:26 , Processed in 0.446737 second(s), 59 queries .

回顶部