数学建模社区-数学中国
标题:
fminbnd的优化问题
[打印本页]
作者:
羽逐云
时间:
2009-7-23 09:17
标题:
fminbnd的优化问题
function f=H(n,a1)
* R# W. E8 e- E! b
c=23.43/180*pi;
8 [% K5 j$ [5 F- y" a2 v. M
a0=32.05*pi/180;
( E) G: b1 d- A" Y9 i& K _& {) W
h=zeros(n,1);
; \8 \1 @+ b4 p
%a1=zeros(1,n);
9 U- ~9 X$ {3 T3 u8 M
t=0;
& J1 |* ~/ m. E2 y
for k = 1: n
7 P/ B5 C$ L( \4 p, e# p1 @' Z
b=asin(sin(c)*cos(pi*k/n));
. \% Q* B/ s% F: z' @0 |$ |: i
A=sqrt(1-(sec(a0))^2*(sin(b))^2);
& Y4 c6 s" D; g E) w: u9 [3 \
B=sin(b)*acos(tan(a0)*tan(b));
* ~9 r7 \& b6 y( T
e=atan(B/A);
7 N S P% r6 r
h(k)=sqrt(A^2+B^2)*cos(a0+e-a1);
; b# y. J# _: n
%T(k)=k/n;
! [" c& j- j0 t
%a1(k)=a1(k)*180/pi;
: O$ _ l+ z$ h9 P
t=t+h(k);
3 R0 \/ k4 R& A
end
9 q5 P2 d1 P! k; Y
f=189.8325-t;
2 d+ I" e% m. L6 ^- A3 p4 L; E, y
. Q6 r( v* T+ J8 E ]! D6 S
然后用下面程序调用，怎么老出错呢，希望高手能解答一下
9 s* f. f- A- _" q+ ?0 r0 m* X+ U
5 n" C/ V1 G3 \, E6 t) M
n=182;
0 z% F4 V/ N* Z; A. o; G8 v
a1=fminbnd(@H,0,pi/2,[],n)
3 W& s7 n. P9 `5 Y& x B) R# A
不知道哪里出错了，谢谢帮助了
作者:
吖信
时间:
2009-7-23 11:23
报错？还是结果错了？我刚试试没报错啊~
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5