求数学大神 帮忙解决拟牛顿法求最大值的最优化问题 matlab

Smile_Again

f(i)(n)=sum(u(i)(k))exp(j*2pi*n*k/N),k=0,....,N-1 sum是对k从0到N-1求和。
u(i)(k)=z(i)(k)exp(j*phi(k)),k=0,....,N-1 z(i)(k)是输入矩阵；
. g  t$ H# c0 w% V# mdeta(i)=sqrt(sum(abs(f(i)(n)-miu(i))^2)) sum是对n从0到N-1求和；
miu(i)=sum(abs(f(i)(n)))/N,sum是对n从0到N-1求和
6 @% ?7 C* @6 q4 HC=sum(deta(i)/miu(i))/M,sum是对i从0到M-1求和
用拟牛顿法求C的最大值时的pih(k)是多少，也就是自变量为phi(k),目标函数为C的最优化问题
+ P7 t, s8 Y4 g) k: x2 B7 L4 EC关于phi(k)的倒数公式给出 即梯度：
dC/dphi(k)=sum(lamda(i)*imag(u*(i)(k)*q(i)(k))) sum是对i从0到M-1求和；
+ T' n: f1 ]. r+ g其中 lamda(i)=-(1/deta(i)+deta(i)/miu(i)^2)/M/N;i=0,1,....,M-1
/ E5 ^: A6 C+ b% y6 _. B# [q(i)(k)=sum(f(i)(n)/abs(f(i)(n))exp(-j*2pi*n*k/N)),sum是对n从0到N-1求和
- ?- e; z8 w, g: Qu*(i)(k)是u(i)(k)的共轭
% @! q) ]- d0 o& [# U& ~: ~4 G求大神帮忙解决了这个问题啊

wssl103050

Smile_Again

木有大神吗？