fminbnd的优化问题

羽逐云

function f=H(n,a1)  
5 J! v  l% R2 F3 n( M/ dc=23.43/180*pi;
; }# K; J( ~  x& k$ ]a0=32.05*pi/180;
% A, S% g3 a3 c: [h=zeros(n,1);
8 v* R- A( E4 e0 y( G! V& U%a1=zeros(1,n);
& X( P3 f) f2 F$ F, z$ P    t=0;
- T/ ?1 ~4 C% O3 h, ^; ^    for k = 1: n
        b=asin(sin(c)*cos(pi*k/n));
" f" Z& L  x: U* }. Z) r7 n        A=sqrt(1-(sec(a0))^2*(sin(b))^2);
/ |1 e/ h$ Y" n  M) x2 _        B=sin(b)*acos(tan(a0)*tan(b));
9 z/ q4 ]5 d" J; V5 M        e=atan(B/A);
9 p9 i  z* r- Y+ x/ {  n9 f- i        h(k)=sqrt(A^2+B^2)*cos(a0+e-a1);
* D* Y/ p% R+ a8 o        %T(k)=k/n;
        %a1(k)=a1(k)*180/pi;
        t=t+h(k);
2 L5 u! i$ ~; k$ z    end
    f=189.8325-t;
' A0 Y$ t2 ]: ^' r) |
4 T( p- b( D3 E2 c$ [2 @( u然后用下面程序调用，怎么老出错呢，希望高手能解答一下
/ l! \' b$ ?' `4 }
n=182;
a1=fminbnd(@H,0,pi/2,[],n)
不知道哪里出错了，谢谢帮助了

吖信

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