fminbnd的优化问题

羽逐云

function f=H(n,a1)  
# C1 d: k* s' K* a  i: T& ic=23.43/180*pi;
a0=32.05*pi/180;
8 x1 \. Z$ h/ b5 p, j  Qh=zeros(n,1);
%a1=zeros(1,n);
- W6 ^" W6 z' v, G4 n+ y& c4 d    t=0;
; `6 l9 y/ i/ t  F/ Q7 e# \    for k = 1: n
        b=asin(sin(c)*cos(pi*k/n));
6 {# h6 s. F: B1 H        A=sqrt(1-(sec(a0))^2*(sin(b))^2);
9 P% }$ t- r' x2 Q) J        B=sin(b)*acos(tan(a0)*tan(b));
" [* S3 A) j9 A        e=atan(B/A);
4 {' Y  _. _9 S0 f' i% b        h(k)=sqrt(A^2+B^2)*cos(a0+e-a1);
# C% Z- O! p& R- B        %T(k)=k/n;
        %a1(k)=a1(k)*180/pi;
        t=t+h(k);
    end
    f=189.8325-t;
- K; d$ Z- l" C) r# |
然后用下面程序调用，怎么老出错呢，希望高手能解答一下

n=182;
* o1 O5 r1 y+ X' C: p: }a1=fminbnd(@H,0,pi/2,[],n)
4 ^9 ^8 x" B8 {% @' C# J不知道哪里出错了，谢谢帮助了

吖信

报错？还是结果错了？我刚试试没报错啊~