fminbnd的优化问题

羽逐云

function f=H(n,a1)  
c=23.43/180*pi;
a0=32.05*pi/180;
h=zeros(n,1);
: p: M3 K% W4 @* p6 P%a1=zeros(1,n);
    t=0;
0 K5 U4 h" `, z; K" j4 y    for k = 1: n
        b=asin(sin(c)*cos(pi*k/n));
        A=sqrt(1-(sec(a0))^2*(sin(b))^2);
        B=sin(b)*acos(tan(a0)*tan(b));
* [3 i: p4 X- O' l2 \        e=atan(B/A);
        h(k)=sqrt(A^2+B^2)*cos(a0+e-a1);
$ b/ H+ L# f3 C( N1 M1 q        %T(k)=k/n;
& R9 ?# j' V+ U$ i7 W        %a1(k)=a1(k)*180/pi;
# J/ n" S; t  ^8 t% j        t=t+h(k);
# @( r  U$ s; X2 b0 f    end
! a2 r/ p" v! p0 \    f=189.8325-t;

然后用下面程序调用，怎么老出错呢，希望高手能解答一下
0 w) Z' k/ I( K3 [4 ~
& W' Z: ]/ N5 f( e% L% Q- H( }n=182;
a1=fminbnd(@H,0,pi/2,[],n)
不知道哪里出错了，谢谢帮助了

吖信

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