fminbnd的优化问题

羽逐云

function f=H(n,a1)  
c=23.43/180*pi;
& y$ M( V/ |9 Y( x7 s( va0=32.05*pi/180;
h=zeros(n,1);
%a1=zeros(1,n);
    t=0;
) J& t9 i$ S9 U) h# R    for k = 1: n
        b=asin(sin(c)*cos(pi*k/n));
) n3 {. O- [( X0 }9 P. I' c; w        A=sqrt(1-(sec(a0))^2*(sin(b))^2);
8 b+ B' |2 Q" t6 r" r3 _; F        B=sin(b)*acos(tan(a0)*tan(b));
        e=atan(B/A);
4 D. [: i& ^5 P        h(k)=sqrt(A^2+B^2)*cos(a0+e-a1);
! \' q( y* S2 n5 ^        %T(k)=k/n;
# G9 `4 T$ N3 v% W6 [5 {        %a1(k)=a1(k)*180/pi;
" q3 Y" F1 r4 e        t=t+h(k);
    end
  I3 j0 p5 A* L5 C    f=189.8325-t;

$ I/ q8 y4 P" X" v4 C) Q4 P然后用下面程序调用，怎么老出错呢，希望高手能解答一下

. w+ G2 R4 \) J& ?n=182;
/ b  }  }$ L! \' Ya1=fminbnd(@H,0,pi/2,[],n)
) A" U; ]2 e( w/ y% R) [8 T8 M& j不知道哪里出错了，谢谢帮助了

吖信

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