fminbnd的优化问题

羽逐云

function f=H(n,a1)  
4 g& m" }$ e9 y( M* t' A" cc=23.43/180*pi;
a0=32.05*pi/180;
h=zeros(n,1);
+ [) d5 f( U6 g* W%a1=zeros(1,n);
    t=0;
    for k = 1: n
        b=asin(sin(c)*cos(pi*k/n));
: C# z4 _/ M$ Y& ?$ A        A=sqrt(1-(sec(a0))^2*(sin(b))^2);
        B=sin(b)*acos(tan(a0)*tan(b));
        e=atan(B/A);
5 \8 r) \4 b4 C* g! J        h(k)=sqrt(A^2+B^2)*cos(a0+e-a1);
$ ^; a/ Z4 e; M" ]5 a1 S% D  W        %T(k)=k/n;
        %a1(k)=a1(k)*180/pi;
, z( }; i8 a1 a+ w* t" z        t=t+h(k);
, d  `! U! z  q3 j3 g% f    end
( v; i! C3 E* y# @    f=189.8325-t;
; s% |# u" l. o
然后用下面程序调用，怎么老出错呢，希望高手能解答一下

n=182;
a1=fminbnd(@H,0,pi/2,[],n)
1 U! d7 J  d% {4 a9 n+ C$ ?$ I不知道哪里出错了，谢谢帮助了

吖信

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