数学建模社区-数学中国
标题:
fminbnd的优化问题
[打印本页]
作者:
羽逐云
时间:
2009-7-23 09:17
标题:
fminbnd的优化问题
function f=H(n,a1)
9 m- F$ T: P) }3 B+ ^" M* r5 V
c=23.43/180*pi;
% a# c9 n" a& W5 H2 `+ w
a0=32.05*pi/180;
4 O0 h) O; e7 I- w5 f' j) i
h=zeros(n,1);
( F) z1 U) x I) [* @- y; d
%a1=zeros(1,n);
# h. q3 z- P& m1 L% e: N
t=0;
7 S' B( h2 O6 j8 K% v; X' U
for k = 1: n
% {' E+ l$ k- ]4 h% w5 ~" Q
b=asin(sin(c)*cos(pi*k/n));
O9 c' P0 F+ f' H8 G& Y* w
A=sqrt(1-(sec(a0))^2*(sin(b))^2);
- C" K5 i& r0 f# n! C0 q8 a
B=sin(b)*acos(tan(a0)*tan(b));
# A5 J/ |) c7 r/ o0 Y! h
e=atan(B/A);
6 Y0 ]0 X7 h# X* x; M( V/ [" q
h(k)=sqrt(A^2+B^2)*cos(a0+e-a1);
) u- K4 X5 T8 I" ], e
%T(k)=k/n;
7 h. S( T4 u, f: L/ i C3 O/ F. e
%a1(k)=a1(k)*180/pi;
3 I# T& }# L0 I; R2 }/ @
t=t+h(k);
; Z4 `) \8 O2 B, B' Y
end
( s* B! J" c5 r6 J4 Q5 x
f=189.8325-t;
8 J. ?( W% R3 F, T7 [% T+ E+ T" e
C) p4 F% e6 u; [
然后用下面程序调用，怎么老出错呢，希望高手能解答一下
: t: ~' T$ T, w) W3 t
. ]" e5 N4 _: u, @4 l2 y( F
n=182;
g' R9 ^1 F: W9 B6 b6 I% G
a1=fminbnd(@H,0,pi/2,[],n)
0 E2 j, m. _2 b
不知道哪里出错了，谢谢帮助了
作者:
吖信
时间:
2009-7-23 11:23
报错？还是结果错了？我刚试试没报错啊~
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5