fminbnd的优化问题

羽逐云

function f=H(n,a1)  
% }1 K+ {# y1 k$ {c=23.43/180*pi;
% c; x. n' g! Z4 A' m6 D( ra0=32.05*pi/180;
9 e/ P/ w2 f' A  N9 _' b$ `h=zeros(n,1);
, u/ K3 F' r5 }, R$ ^- R%a1=zeros(1,n);
* {+ p. d% n$ ]& ]& R7 d# |    t=0;
4 n: X+ \% m2 L    for k = 1: n
        b=asin(sin(c)*cos(pi*k/n));
% h! `7 }$ d4 V/ t8 \        A=sqrt(1-(sec(a0))^2*(sin(b))^2);
        B=sin(b)*acos(tan(a0)*tan(b));
        e=atan(B/A);
  B2 \( V+ ]$ ?! b; c4 t! F        h(k)=sqrt(A^2+B^2)*cos(a0+e-a1);
$ Y( @5 T% v0 ?  r& p+ p; l        %T(k)=k/n;
        %a1(k)=a1(k)*180/pi;
        t=t+h(k);
    end
    f=189.8325-t;

& P% Y  i& F$ c. P/ J然后用下面程序调用，怎么老出错呢，希望高手能解答一下

n=182;
a1=fminbnd(@H,0,pi/2,[],n)
不知道哪里出错了，谢谢帮助了

吖信

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