数学建模社区-数学中国

标题: 求数学大神 帮忙解决拟牛顿法求最大值的最优化问题 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求和。
: G! K. C7 n) K% l9 ju(i)(k)=z(i)(k)exp(j*phi(k)),k=0,....,N-1 z(i)(k)是输入矩阵;
: I7 \5 D" T' j. w2 b/ C; @deta(i)=sqrt(sum(abs(f(i)(n)-miu(i))^2)) sum是对n从0到N-1求和;
- E/ d  M% W7 Q, d9 S" m8 Amiu(i)=sum(abs(f(i)(n)))/N,sum是对n从0到N-1求和. |" V% {8 ?% ^: @  ~' K/ c( Q9 ]
C=sum(deta(i)/miu(i))/M,sum是对i从0到M-1求和
- c/ P! ?4 r. l5 o用拟牛顿法求C的最大值时的pih(k)是多少,也就是自变量为phi(k),目标函数为C的最优化问题
! L' Y0 c# ?5 X4 r2 l0 aC关于phi(k)的倒数公式给出 即梯度:8 d$ x8 c4 C/ E" W, ^; d
dC/dphi(k)=sum(lamda(i)*imag(u*(i)(k)*q(i)(k))) sum是对i从0到M-1求和;
& \) N# O* R$ c' P1 V! G8 U其中 lamda(i)=-(1/deta(i)+deta(i)/miu(i)^2)/M/N;i=0,1,....,M-15 i0 ?% y6 u! M4 m/ Y, |5 l
q(i)(k)=sum(f(i)(n)/abs(f(i)(n))exp(-j*2pi*n*k/N)),sum是对n从0到N-1求和
" K" }2 p, }: ~$ |! ou*(i)(k)是u(i)(k)的共轭 ! I& D& [* R) _1 t& P
求大神帮忙解决了这个问题啊
作者: wssl103050    时间: 2012-6-14 13:16

作者: Smile_Again    时间: 2012-6-14 13:27
木有大神吗?




欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) Powered by Discuz! X2.5