matlab 核岭回归代码

以下是 Matlab 中实现核岭回归的示例代码：

% 加载数据
load data.mat

% 定义核函数
kernel = @(x, y) exp(-norm(x-y).^2/(2*sigma^2));

% 构建核矩阵
K = zeros(n, n);
for i=1:n
    for j=1:n
        K(i,j) = kernel(X(i,:), X(j,:));
    end
end

% 计算岭回归系数
alpha = (K + lambda*eye(n)) \ y;

% 预测
m = size(X_test, 1);
y_pred = zeros(m, 1);
for i=1:m
    k = zeros(n, 1);
    for j=1:n
        k(j) = kernel(X_test(i,:), X(j,:));
    end
    y_pred(i) = k'*alpha;
end

其中，`data.mat` 文件中包含训练集数据 `X`（大小为 $n \times d$，表示 $n$ 个样本，每个样本有 $d$ 个特征）、对应的标签 `y`（大小为 $n \times 1$）、测试集数据 `X_test`（大小为 $m \times d$，表示 $m$ 个样本，每个样本有 $d$ 个特征）、核函数的参数 `sigma` 和岭回归的参数 `lambda`。其中，核函数的定义可以根据具体场景进行修改。