数学建模社区-数学中国
标题:
fminbnd的优化问题
[打印本页]
作者:
羽逐云
时间:
2009-7-23 09:17
标题:
fminbnd的优化问题
function f=H(n,a1)
7 A1 c+ N/ m0 v/ D- D+ E) K1 ] Z( R
c=23.43/180*pi;
7 X7 N9 P5 Y$ \' s
a0=32.05*pi/180;
4 {# ?% o/ t7 ]( r
h=zeros(n,1);
, S t3 a( f! l" s% D
%a1=zeros(1,n);
P3 h0 c4 |4 b
t=0;
2 K G. `/ ~4 B7 o
for k = 1: n
8 A) C- r% l2 d" v- q4 P3 m
b=asin(sin(c)*cos(pi*k/n));
9 @8 G! j. X0 [& `% C1 |) F0 G8 X
A=sqrt(1-(sec(a0))^2*(sin(b))^2);
- ^2 c4 [$ O V
B=sin(b)*acos(tan(a0)*tan(b));
# x- c6 F& b& {" Q: Q" D
e=atan(B/A);
3 p; l9 h+ i: U. _8 j7 z
h(k)=sqrt(A^2+B^2)*cos(a0+e-a1);
1 d3 o) a6 v4 b7 G4 b% Z
%T(k)=k/n;
3 s; d0 O, N4 E8 L# u+ Z
%a1(k)=a1(k)*180/pi;
6 A5 o) q3 l+ P8 h
t=t+h(k);
) D0 `- ?/ ~4 D, p
end
, \' S- \1 ]% t& L" ^- u3 H* _
f=189.8325-t;
* I' ] d: @6 F6 G3 @& P2 W
( ?( @9 A) n( T
然后用下面程序调用，怎么老出错呢，希望高手能解答一下
- ~; Q; U' J2 k4 b6 J1 i- M- N
2 U3 v$ E/ N( R( o9 i
n=182;
' g1 C/ W2 |0 g8 S: P5 K
a1=fminbnd(@H,0,pi/2,[],n)
& f9 C2 C' B5 J6 p! G- ?6 l
不知道哪里出错了，谢谢帮助了
作者:
吖信
时间:
2009-7-23 11:23
报错？还是结果错了？我刚试试没报错啊~
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5