首页 > 其他分享 >工程数学实验四

工程数学实验四

时间:2023-05-25 10:36:32浏览次数:37  
标签:opt end 工程 iter 实验 alpha 数学 prev grad

function [x_opt, f_opt, iter] = conjugate_gradient()

    % 定义目标函数
    f = @(x) 100*(x(1)^2 - x(2))^2 + (x(1) - 1)^2;

    % 计算目标函数的梯度
    grad_f = @(x) [400*x(1)*(x(1)^2 - x(2)) + 2*(x(1) - 1); -200*(x(1)^2 - x(2))];

    % 定义终止准则
    epsilon = 1e-5;

    % 定义 Armijo 条件参数
    rho = 0.5; % 步长衰减因子
    c = 0.1;  % 条件参数

    % 设置初始点
    x0_list = [0, 0; -1, 1; 2, -2]; % 与上面的最速下降法和牛顿法使用相同的初始点

    for k = 1:size(x0_list, 1)
        % 初始化变量
        iter = 0;
        x_opt = x0_list(k, :)';
        f_opt = f(x_opt);
        grad_norm = norm(grad_f(x_opt));
        p = -grad_f(x_opt);

        % 共轭梯度法迭代过程
        while grad_norm >= epsilon
            % 计算步长
            alpha = armijo_line_search(x_opt, p, f, grad_f, rho, c);

            % 更新变量
            x_opt_prev = x_opt;
            x_opt = x_opt + alpha * p;
            f_opt_prev = f_opt;
            f_opt = f(x_opt);
            grad_prev = grad_f(x_opt_prev);
            grad_opt = grad_f(x_opt);

            % 计算 beta 参数
            beta = (grad_opt' * grad_opt) / (grad_prev' * grad_prev);

            % 计算共轭梯度方向
            p = -grad_opt + beta * p;

            % 更新迭代次数
            iter = iter + 1;
            grad_norm = norm(grad_opt);
        end

        % 输出结果
        fprintf('初始点:%s\n', mat2str(x0_list(k, :)));
        fprintf('最优解:%s\n', mat2str(x_opt));
        fprintf('最优值:%f\n', f_opt);
        fprintf('迭代次数:%d\n', iter);
        disp('----------------------');
    end
end

function alpha = armijo_line_search(x, p, f, grad_f, rho, c)
    % Armijo 非精确线搜索方法
    alpha = 1;
    while f(x + alpha * p) > f(x) + c * alpha * grad_f(x)' * p
        alpha = rho * alpha;
    end
end
% 调用 FR 共轭梯度法函数
[x_opt, f_opt, iter] = conjugate_gradient();

 

标签:opt,end,工程,iter,实验,alpha,数学,prev,grad
From: https://www.cnblogs.com/psh888/p/17430388.html

相关文章

  • 创龙教仪TL6748-PlusTEB教学实验箱
    目 录1.实验箱简介 2.软硬件参数 3.可选摄像头模块 4.开发资料 5.电气特性 6.实验箱机械尺寸 7.产品认证 8.实验箱套件清单 9.技术支持 10.增值服务 更多帮助 附录A教学实验 1. 实验箱简介Ø 基于TITMS320C6748定点/浮点DSPC674x处理器,主频456MHz,高达3648MIPS和274......
  • 工程数学
    function[x_opt,f_opt,iter]=steepest_descent()%定义目标函数f=@(x)100*(x(1)^2-x(2))^2+(x(1)-1)^2;%计算目标函数的梯度grad_f=@(x)[400*x(1)*(x(1)^2-x(2))+2*(x(1)-1);-200*(x(1)^2-x(2))];%定义终止准则eps......
  • 系统工程(二)系统工程方法
    系统工程方法是对为了以最好实现系统为目的,而对系统的结构、元素、信息、反馈进行分析研究所形成的一系列规则。分类:1、霍尔三维结构,适合组织管理大型工程建设项目   三维:逻辑、时间、知识。(whyhow what,为什么做:解决问题的过程,怎样做:工作进程,需要什么:专业科学知识)逻......
  • 我是如何学习软件工程的
    软件工程是一个庞大且不断演化的领域,涉及许多创新。虽然如此,大多数技术——如果不是全部——在软件工程中往往最终汇聚为几个基本原则。工程师更好地学习软件工程的方法是学习软件工程的基础知识,而不是那些不断变化的框架、语言或平台。即使学习基础知识也需要时间和努力。人们经......
  • 2023.5.24——软件工程站立会议(阶段二)
    站立会议内容:1.整个项目预期的任务量:目前已经花的时间:剩余的时间:2.任务看板照片: 3.团队照片: 4.产品状态:最新做好的功能:正在完成中5.燃尽图:......
  • (大二下)课程总结——软件工程
    1)回顾你的课程计划(第一周的计划),你完成的程度如何?请列出具体数据和实际例子。  1.你在这门课的计划是什么?参考一些学校的教学,你对这个课程有什么期待?你打算怎样度过这个课程?    计划就是尽力跟上建民老师的节奏同时,还能主动学习一些其他知识。我的期待是希望能跟上建民......
  • 2023.5.24——软件工程日报
    所花时间(包括上课):6h代码量(行):0行博客量(篇):1篇今天,上午学习,下午学习。我了解到的知识点:1.了解了一些数据库的知识;2.了解了一些python的知识;3.了解了一些英语知识;5.了解了一些Javaweb的知识;4.了解了一些数学建模的知识;6.了解了一些计算机网络的知识;......
  • SDN实验指导
    SDNSDN-Mininet源码安装与验证实验目的1、了解Mininet的基本概念与作用。2、掌握Mininet的源码安装方法。3、Mininet部署后的自我验证。实验环境支持OpenFlow1.0及OpenFlow1.3协议的Mininet源码安装与验证实验拓扑如下图所示。设备名称软件环境硬件环境主机Ub......
  • 电子科技大学实验报告4 文件状态测试
    ##Monitorforfilestatus#cleartest_count=0change_count=0echo-n"Inputfilename:"readFILENAMEif[-f"$FILENAME"];thenFILESIZE=`ls-l${FILENAME}|cut-d""-f5`echo"filesizeis$FILESIZE"......
  • 实验4 函数与异常处理编程
    print(sum)sum=42print(sum)definc(n):sum=n+1print(sum)returnsumsum=inc(7)+inc(7)print(sum)四处sum作用域不同。第一个是内置函数名称,第二个是赋值名称,第三个是局部变量,第四个是全局变量list1=[1,9,8,4]print(sorted(list1))print......