数学建模社区-数学中国
标题:
求数学大神 帮忙解决拟牛顿法求最大值的最优化问题 matlab
[打印本页]
作者:
Smile_Again
时间:
2012-6-13 23:25
标题:
求数学大神 帮忙解决拟牛顿法求最大值的最优化问题 matlab
f(i)(n)=sum(u(i)(k))exp(j*2pi*n*k/N),k=0,....,N-1 sum是对k从0到N-1求和。
# h/ ?3 j3 A4 a# ~3 B* f4 e
u(i)(k)=z(i)(k)exp(j*phi(k)),k=0,....,N-1 z(i)(k)是输入矩阵;
, O! `3 r9 w1 r# i1 j( w L! ~
deta(i)=sqrt(sum(abs(f(i)(n)-miu(i))^2)) sum是对n从0到N-1求和;
: j' W n# J% q" H# |4 O* R2 `
miu(i)=sum(abs(f(i)(n)))/N,sum是对n从0到N-1求和
8 R# ^$ L0 z+ M$ p9 h3 }% ?$ @
C=sum(deta(i)/miu(i))/M,sum是对i从0到M-1求和
1 `0 k$ x* m0 V
用拟牛顿法求C的最大值时的pih(k)是多少,也就是自变量为phi(k),目标函数为C的最优化问题
/ b1 a. `0 i& O3 A% K
C关于phi(k)的倒数公式给出 即梯度:
' _" G6 {0 k2 ]9 {: _* _
dC/dphi(k)=sum(lamda(i)*imag(u*(i)(k)*q(i)(k))) sum是对i从0到M-1求和;
6 T j% K0 ~/ H$ g: h+ Z+ h1 M& U
其中 lamda(i)=-(1/deta(i)+deta(i)/miu(i)^2)/M/N;i=0,1,....,M-1
) g+ l7 W3 ~1 }+ a j0 {& S5 \
q(i)(k)=sum(f(i)(n)/abs(f(i)(n))exp(-j*2pi*n*k/N)),sum是对n从0到N-1求和
' @8 \8 |1 d) L2 s& h
u*(i)(k)是u(i)(k)的共轭
" Z* M0 r' K1 H1 H& N' y. o t
求大神帮忙解决了这个问题啊
作者:
wssl103050
时间:
2012-6-14 13:16
作者:
Smile_Again
时间:
2012-6-14 13:27
木有大神吗?
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5