fminbnd的优化问题

羽逐云

function f=H(n,a1)  
c=23.43/180*pi;
a0=32.05*pi/180;
h=zeros(n,1);
; z# {8 l& ?$ L' b4 N8 H' E%a1=zeros(1,n);
    t=0;
    for k = 1: n
        b=asin(sin(c)*cos(pi*k/n));
0 s0 t& a' Q/ ^        A=sqrt(1-(sec(a0))^2*(sin(b))^2);
2 j7 Q4 e2 r1 z+ g5 P* ]6 i0 t        B=sin(b)*acos(tan(a0)*tan(b));
6 ^7 ^) }9 z+ _) ]0 x        e=atan(B/A);
        h(k)=sqrt(A^2+B^2)*cos(a0+e-a1);
        %T(k)=k/n;
        %a1(k)=a1(k)*180/pi;
        t=t+h(k);
    end
    f=189.8325-t;
' Y# K: p  E, v  ^/ W. c# y$ X6 K
然后用下面程序调用，怎么老出错呢，希望高手能解答一下

n=182;
; @1 B* _$ E8 }% T+ ba1=fminbnd(@H,0,pi/2,[],n)
: t& a0 B0 `- R4 D' r7 l不知道哪里出错了，谢谢帮助了

吖信

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