5月21日学习进度

合集 - 学习进度（第四学期）(57)

1.3月19日学习进度2024-03-192.第一次结对作业3月27日开发进度2024-03-273.第一次结对作业4月11日开发进度2024-04-114.4月18日冲刺记录2024-04-185.4月19日学习进度2024-04-196.4月19日冲刺记录2024-04-197.4月20日学习进度2024-04-208.4月20日冲刺记录2024-04-209.4月21日冲刺记录2024-04-2110.第一次结对作业4月10日开发进度2024-04-1011.4月23日学习进度2024-04-2312.4月23日冲刺记录2024-04-2313.4月22日冲刺记录2024-04-2214.4月25日冲刺记录2024-04-2515.4月26日冲刺记录2024-04-2616.4月27日学习记录2024-04-2717.4月28日学习记录2024-04-2818.4月24日冲刺记录2024-04-2419.4月22日学习进度2024-04-2220.4月30日学习进度2024-04-3021.5月6日学习进度2024-05-0622.5月8日学习进度2024-05-0823.5月10日学习进度2024-05-1024.6月4日学习进度2024-06-0425.5月11日学习进度2024-05-1126.5月12日学习进度2024-05-1227.5月13日学习进度2024-05-1328.5月15日学习进度2024-05-1529.5月16日学习进度2024-05-1630.5月17日学习进度2024-05-1731.5月18日学习进度2024-05-1832.5月19日学习进度2024-05-1933.5月20日学习进度2024-05-20

34.5月21日学习进度2024-05-21

35.5月22日学习进度2024-05-2236.5月23日学习进度2024-05-2337.5月24日学习进度2024-05-2438.5月25日学习进度2024-05-2539.5月26日学习进度2024-05-2640.5月27日学习进度2024-05-2741.5月28日学习进度2024-05-2842.5月29日学习进度2024-05-2943.5月30日学习进度2024-05-3044.5月31日学习进度2024-05-3145.6月1日学习进度2024-06-0146.6月2日学习进度2024-06-0247.6月3日学习进度2024-06-0348.6月5日学习进度2024-06-0549.6月6日学习进度2024-06-0650.6月7日学习进度2024-06-0751.6月11日学习进度2024-06-1152.6月12日学习进度（个人总结）2024-06-1253.6月14日学习进度2024-06-1454.6月17日学习进度2024-06-1755.人月神话读后感12024-05-1056.人月神话读后感22024-05-2357.人月神话读后感32024-06-01

收起

一.所花时间

0.5h

二.代码量

30行

三.博客量

1篇

四.了解到的知识点

实验二：最速下降法程序设计
一、实验目的
通过最速下降法的程序设计，为今后的约束优化方法的学习和编程奠定基础；掌握负梯度方向的定义和最速下降法的迭代公式 ；通过此次实验，进一步巩固最速下降法的基本原理和思想。
二、实验内容
（1）求解无约束优化问题： ；
（2）终止准则取 ；
（3）完成最速下降法（负梯度法）的MATLAB编程、调试；
（4）要求选取多个不同的初始点，并给出迭代次数，最优函数值等相关信息，有能力的同学尝试画出最优值随迭代次数变化的曲线图；
（5）按照模板撰写实验报告，要求规范整洁。
三、算法步骤、代码、及结果
1. 算法步骤
（1）给出初始点x 0，置iter=0，最大迭代次数max_iter=10000，精度e=1e-6，目标 函数f及其梯度grad
（2）若终止条件满足(| | f | | < e)，则迭代停止
（3）计算d k = −grad
（4）在d k方向上利用一维精确线性搜索Amrijo搜索求步长ak
（5）x k + 1 = x k + α k d k ， k = k + 1，转步骤2

2.代码

function f=fun(x)
f= (x(1)+10*x(2))^2+5*(x(3)-x(4))^2+(x(2)-2*x(3))^4+10*(x(1)-x(4))^4;


function gf=gfun(x)
gf = [2*(x(1)+10*x(2))+40*(x(1)-x(4))^3;  20*(10*x(2)+x(1))+4*(x(2)-2*x(3))^3;  10*(x(3)-x(4))-8*(x(2)-2*x(3))^3;  -10*(x(3)-x(4))-40*(x(1)-x(4))^3];


function [k,x,val] = grad(fun,gfun,x0_list,epsilon)

%功能：用最速下降法求解无约束优化问题

%输入：funf，gfun分别是目标函数和梯度，x0是初始点，epsilon是容许误差

%输出：x,val分别是近似最优解和最优值，k是迭代次数
maxk=10000; %最大迭代次数
value = zeros(maxk, 1);
for i = 1:length(x0_list)-1
    x0 = x0_list(:,i);

    beta=0.5;  sigma=0.4;
    k=0;
    fprintf('初始点 (%g; %g; %g; %g)\n', x0(1), x0(2), x0(3), x0(4));
    while(k<maxk)

        gk=feval(gfun,x0); %计算梯度
    
        dk=-gk;
    
        if(norm(gk)<epsilon), break; end    %检验终止准则
        m=0; mk=0;
    
        %用Amrijo搜索技术确定步长
        while(m<20) %最大迭代次数
    
           if(feval(fun,x0+beta^m*dk)<=feval(fun,x0)+sigma*beta^m*gk'*dk)
               mk=m;  break;
           end
           m=m+1;
        end
        x0=x0+beta^mk*dk;
        k=k+1;
    value(k,1)=feval(fun,x0);
    end
    x=x0;
    val=feval(fun,x0);
    

    fprintf('迭代次数: %d\n', k);
    fprintf('最优点: (%g; %g; %g; %g)\n', x(1), x(2), x(3), x(4));
    fprintf('最优函数值: %g\n', val);
    fprintf('\n');

    % 迭代曲线图
    figure;
    plot(value);
    title('迭代曲线');
    xlabel('迭代次数');
    ylabel('最优值');

end


% 命令行输入
% x0=[-1,-2,3;  1,7,10;  3,-4,5;  3,6,4];
% grad('fun','gfun',x0,1e-6);

3.结果