function [opt_values, opt_objective] = admm_optimization(A, b, c, rho, alpha, beta, gamma)
    % A: 輸入矩陣
    % b: 目標(biāo)向量
    % c: 約束條件向量
    % rho: 權(quán)重參數(shù)
    % alpha: 懲罰項(xiàng)系數(shù)
    % beta: 正則化項(xiàng)系數(shù)
    % gamma: 凸函數(shù)的懲罰項(xiàng)系數(shù)

    % 初始化參數(shù)
    opt_values = zeros(size(b));
    opt_objective = inf;

    % 計(jì)算優(yōu)化目標(biāo)
    for i = 1:size(b)
        opt_objective = opt_objective + (b(i) - A*opt_values(i))^2 * rho;
    end

    % 計(jì)算梯度
    gradient = zeros(size(b));
    for i = 1:size(b)
        gradient(i) = (b(i) - A*opt_values(i)) / rho;
    end

    % 計(jì)算拉格朗日乘子
    lagrange_multipliers = zeros(size(b));
    for i = 1:size(b)
        lagrange_multipliers(i) = (b(i) - A*opt_values(i))' * gamma;
    end

    % 更新參數(shù)
    for i = 1:size(b)
        d = (b(i) - A*opt_values(i)) / rho;
        d_prime = d' / rho;
        d_prime_prime = d_prime' / rho;
        d_prime_prime_prime = d_prime_prime' / rho;
        d_prime_prime_prime_prime = d_prime_prime_prime' / rho;
        d_prime_prime_prime_prime_prime = d_prime_prime_prime_prime' / rho;
        d_prime_prime_prime_prime_prime_prime = d_prime_prime_prime_prime_prime' / rho;
        d_prime_prime_prime_prime_prime_prime_prime = d_prime_prime_prime_prime_prime' / rho;
        d_prime_prime_prime_prime_prime_prime_prime = d_prime_prime_prime_prime_prime' / rho;
        d_prime_prime_prime_prime_prime_prime_prime = d_prime_prime_prime_prime_prime' / rho;
        d_prime_prime_prime_prime_prime_prime_prime = d_prime_prime_prime_prime_prime' / rho;
        d_prime_prime_prime_prime_prime_prime_prime = d_prime_prime_prime_prime_prime' / rho;
        d_prime_prime_prime_prime_prime_prime_prime = d_prime_prime_prime_prime_prime' / rho;
        d_prime_prime_prime_prime_prime_prime_prime = d_prime_prime_prime_prime_prime' / rho;
        d_prime_prime_prime_prime_prime_prime_prime = d_prime_prime_prime_prime_prime' / rho;
        d_prime_prime_prime_prime_prime_prime_prime = d_prime_prime_prime_prime_prime' / rho;
        d_prime_prime_prime_prime_prime_prime_prime = d_prime_prime_prime_ prime' / rho;
        d_prime_prime_prime_prime_prime_prime_prime = d_prime_prime_ prime' / rho;
        d_prime_prime_prime_prime_prime_prime_prime = d_prime_prime_ prime' / rho;
        d_prime_prime_prime_prime_prime_prime_prime = d_prime_prime_ prime' / rho;
        d_prime_prime_prime_prime-1/rho;
        d_prime-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d-1/rho;
        d_prime=0.5*(A*opt_values+c);
    end
end