fminbnd的优化问题

羽逐云

function f=H(n,a1)  
c=23.43/180*pi;
( S5 _& l: T' f/ Ra0=32.05*pi/180;
4 d( S, p  u, d; c/ J+ uh=zeros(n,1);
%a1=zeros(1,n);
    t=0;
- _) o* \4 J4 `' c5 w    for k = 1: n
' t/ B& f) G& V) H+ R: k( Z4 F        b=asin(sin(c)*cos(pi*k/n));
. y3 {. x& l# r1 l        A=sqrt(1-(sec(a0))^2*(sin(b))^2);
: ~+ l0 |, K, a! E        B=sin(b)*acos(tan(a0)*tan(b));
        e=atan(B/A);
- l" B) o/ L0 J/ K# u        h(k)=sqrt(A^2+B^2)*cos(a0+e-a1);
        %T(k)=k/n;
2 M: r9 F  L* v6 ^        %a1(k)=a1(k)*180/pi;
        t=t+h(k);
    end
1 q* S0 n7 |8 a6 q- c    f=189.8325-t;

! k* j3 Y7 x" E4 q然后用下面程序调用，怎么老出错呢，希望高手能解答一下
  l' j) W. h. c5 m
n=182;
a1=fminbnd(@H,0,pi/2,[],n)
! t( x" G. b( |  u9 Z$ N不知道哪里出错了，谢谢帮助了

吖信

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