数学建模社区-数学中国
标题:
fminbnd的优化问题
[打印本页]
作者:
羽逐云
时间:
2009-7-23 09:17
标题:
fminbnd的优化问题
function f=H(n,a1)
0 N/ B2 I; o3 v2 D; m
c=23.43/180*pi;
" K3 o( K4 T1 V
a0=32.05*pi/180;
9 z. W, G! I, k$ e0 `
h=zeros(n,1);
2 O- N) {$ E w* M
%a1=zeros(1,n);
; R- i/ F- u( w/ _. q
t=0;
) U. D9 E' c) i- p
for k = 1: n
* K4 C- m1 }% L6 v0 _" R7 y
b=asin(sin(c)*cos(pi*k/n));
: G `8 Q( H7 _) m, B7 L) l/ u
A=sqrt(1-(sec(a0))^2*(sin(b))^2);
8 ` U: t0 n/ p% a
B=sin(b)*acos(tan(a0)*tan(b));
4 v3 {( i! W- V5 F
e=atan(B/A);
7 ~& R& G( e* _/ Z
h(k)=sqrt(A^2+B^2)*cos(a0+e-a1);
% C7 ?* Y# E* V$ _# F9 E5 t
%T(k)=k/n;
4 b' ]. N- k) G0 S! `
%a1(k)=a1(k)*180/pi;
; V! C! T, O6 O1 I, y. j& n
t=t+h(k);
+ {% A8 O! S5 ?+ T* W1 b: J p
end
& F: V' @% D* @2 I9 [) |
f=189.8325-t;
d- P$ p% k f
! R3 E) x$ m- a" Z) k+ L
然后用下面程序调用，怎么老出错呢，希望高手能解答一下
- s& L8 R9 _! j8 j: k
) W& [( A; E7 m/ h o6 ^/ J
n=182;
% T$ @; E1 R* g# G. `
a1=fminbnd(@H,0,pi/2,[],n)
$ T) n% m- E' B
不知道哪里出错了，谢谢帮助了
作者:
吖信
时间:
2009-7-23 11:23
报错？还是结果错了？我刚试试没报错啊~
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5