azure011328

导航

< 2025年3月 >
23 24 25 26 27 28 1
2 3 4 5 6 7 8
9 10 11 12 13 14 15
16 17 18 19 20 21 22
23 24 25 26 27 28 29
30 31 1 2 3 4 5
统计
 

今天学习了共轭梯度法 代码如下

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
function [x_min, f_min] = FR_conjugate_gradient()
 
    % 初始点
    x = [2; 2; 2; 2];
     
    % 初始梯度
    g = gradient(x);
     
    % FR共轭方向
    d = -g;
     
    % 终止准则参数
    epsilon = 1e-6;
     
    % Armijo条件参数
    rho = 0.5;
    c = 0.1;
     
    % 最大迭代次数
    max_iter = 1000;
     
    % FR共轭梯度法迭代
    for iter = 1:max_iter
         
        % 非精确搜索
        alpha = 1;
        while f(x + alpha * d) > f(x) + c * alpha * g' * d
            alpha = rho * alpha;
        end
         
        % 更新x
        x = x + alpha * d;
         
        % 更新梯度
        g_old = g;
        g = gradient(x);
         
        % 终止条件
        if norm(g) < epsilon
            break;
        end
         
        % 更新FR共轭方向
        beta = (g' * g) / (g_old' * g_old);
        d = -g + beta * d;
    end
     
    % 返回最优解和最优值
    x_min = x;
    f_min = f(x);
 
end
 
% 目标函数
function val = f(x)
    val = (x(1) + 10*x(2))^2 + 5*(x(3) - x(4))^2 + (x(2) - 2*x(3))^4 + 10*(x(1) - x(4))^4;
end
 
% 梯度函数
function grad = gradient(x)
    grad = [
        2*(x(1) + 10*x(2)) + 40*(x(1) - x(4))^3;
        20*(x(1) + 10*x(2)) + 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
    ];
end

  

posted on   淮竹i  阅读(6)  评论(0编辑  收藏  举报
相关博文:
阅读排行:
· winform 绘制太阳,地球,月球 运作规律
· AI与.NET技术实操系列(五):向量存储与相似性搜索在 .NET 中的实现
· 超详细:普通电脑也行Windows部署deepseek R1训练数据并当服务器共享给他人
· 【硬核科普】Trae如何「偷看」你的代码?零基础破解AI编程运行原理
· 上周热点回顾(3.3-3.9)
 
点击右上角即可分享
微信分享提示