用尺度变换法求f(x) = 60 - 10x1 - 4x2 + x1^2 + x2^2 - x1*x2的MATLAB编程

在MATLAB中，尺度变换法通常用于特征向量的计算或者数据预处理，比如主成分分析（PCA）。对于函数`f(x)`，这是一个二次多项式，我们并不直接用尺度变换来求解，而更多是用来表示一个二维空间中的曲面。如果你想要对其进行梯度下降或其他优化操作，你可以编写一个函数来计算函数值及其梯度。

function [f, grad] = compute_function_and_gradient(x)
% 定义函数f(x)
f = 60 - 10*x(1) - 4*x(2) + x(1)^2 + x(2)^2 - x(1)*x(2);

% 计算偏导数，即梯度grad
grad = [-10 - 4 + 2*x(1) - x(2); -4 + 2*x(2) - x(1)];
end

% 示例使用该函数
x_start = [0; 0]; % 初始点
learning_rate = 0.01; % 学习率
num_iterations = 100; % 迭代次数

for i = 1:num_iterations
    [f_value, gradient] = compute_function_and_gradient(x_start);
    x_start = x_start - learning_rate * gradient;
    disp(['Iteration ', num2str(i), ': f-value = ', num2str(f_value)]);
end