fminbnd的优化问题

[复制链接]

字体大小: 正常放大

羽逐云

2 主题	4 听众	12 积分

升级 7.37%

该用户从未签到

电梯直达

1^#

发表于 2009-7-23 09:17 |只看该作者 |倒序浏览

|招呼Ta 关注Ta

function f=H(n,a1)
c=23.43/180*pi;
a0=32.05*pi/180;
h=zeros(n,1);
%a1=zeros(1,n);
t=0;
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));
      e=atan(B/A);
      h(k)=sqrt(A^2+B^2)*cos(a0+e-a1);
      %T(k)=k/n;
      %a1(k)=a1(k)*180/pi;
      t=t+h(k);
end
f=189.8325-t;

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

n=182;
a1=fminbnd(@H,0,pi/2,[],n)
不知道哪里出错了，谢谢帮助了

zan