求一个最优化算法的问题

changyl01

小弟初学matlab
想验证一下袁亚湘提出的最速下降法的新步长
所以要写一个函数： y=(X-X*)A(X-X*)'
要求这个X是1*10000的
A是10000*10000的
因为要用最速下降法
所以要求上面函数的梯度函数（在作步长的时候用）
我不知道怎么写梯度函数
跪求各位大侠帮忙！

[ 本帖最后由 eight 于 2007-1-24 18:48 编辑 ]

happy

最速下下降法程序，其中梯度函数你可以用diff获得

1. function y=F(x)
   
2. y=(x(1)-15)^2+(x(2)-45)^4;
   
3. 
   
4. function y=D_F(x)
   
5. %D_F Summary of this function goes here
   
6. %  Detailed explanation goes here
   
7. y=[0,0];
   
8. y1=2*(x(1)-15),  y2=4*((x(2)-45)^3);
   
9. y=[y1,y2];
   
10. 
    
11. 
    
12. format long
    
13. x=[0,0];
    
14. x=input('请输入初始值[x0,x1]=',[0,0]);
    
15. s=input('请输入初始精度s＝',0.1);
    
16. g=D_F(x);
    
17. g=g.';%转置为列向量；
    
18. while(norm(g)>s)
    
19.     s_k=-g,   r=min_F(x,s_k),   x=x+r.*s_k.';
    
20.     g=D_F(x);
    
21.     g=g.';
    
22. fprintf('The result is %8.5f\n', x);
    
23.     fprintf('The daosu is %8.5f\n', g);   
    
24. end
    
25. axis normal;
    
26. hold on;
    
27. view([1,1,1]);
    
28. for i=-100:100
    
29.    for j=-100:100
    
30. z=[i,j],  z=F(z);
    
31. plot3(i,j,z);
    
32. end
    
33. end
    
34. hold off;

复制代码

changyl01

太感谢了