fminbnd的优化问题

羽逐云

function f=H(n,a1)  
c=23.43/180*pi;
a0=32.05*pi/180;
% u+ }; E# Y/ Z8 Z8 [- m! Fh=zeros(n,1);
' S. b9 }4 ?) i" {2 f2 `- C- z9 L%a1=zeros(1,n);
    t=0;
    for k = 1: n
" }6 q: L5 F" Z# E8 w# S$ S& S6 V        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));
. h7 C" `) j' g2 A  `# J3 k: G) e        e=atan(B/A);
  L" R0 C+ A& y7 v2 ?        h(k)=sqrt(A^2+B^2)*cos(a0+e-a1);
        %T(k)=k/n;
, I3 r0 q- e. S2 G% ~: G        %a1(k)=a1(k)*180/pi;
        t=t+h(k);
    end
4 I/ W+ c$ O) g+ ?% \! B  N    f=189.8325-t;

然后用下面程序调用，怎么老出错呢，希望高手能解答一下

, D( p' F3 i7 k' L) F6 |n=182;
% w2 u4 z0 f. ~8 O( ha1=fminbnd(@H,0,pi/2,[],n)
不知道哪里出错了，谢谢帮助了

吖信

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