请写一段matlab代码，用Jiacobi方法计算三对角矩阵A的特征值，A的主对角线元素全为2，低对角线和高对角线元素全为1

要使用Jacobi方法（也称为迭代对角化）来计算具有特定结构的三对角矩阵A的特征值，我们可以编写以下MATLAB代码：

function [values, vectors] = jacobi_eigen(A)
% A = [2 ones(1,N-1); ones(1,N-1) 2]; % 创建三对角矩阵，主对角线为2，低/高对角线为1
% N = size(A,1); % 获取矩阵维度

% 初始化
tol = 1e-6; % 设置收敛误差阈值
max_iter = 1000; % 最大迭代次数
v = ones(size(A)); % 初始猜测向量
D = toeplitz(diag(v)); % 对角矩阵，对应于v的diagonal
old_d = Inf; % 初始估计的特征值

for iter = 1:max_iter
    old_v = v;
    DinvV = inv(D) * v; % 计算D^-1 * v
    v = (A - D + DinvV*DinvV') * v; % Jacobi迭代公式
    v = v / norm(v); % 正则化
    d = dot(DinvV, v); % 当前估计的特征值
    
    % 检查收敛条件
    if abs(d - old_d) < tol
        break;
    end
    old_d = d;
end

% 特征值和对应的特征向量
values = d;
vectors(:,1) = v; % 向量存储在矩阵的第一列
end

这段代码定义了一个名为`jacobi_eigen`的函数，它接受一个三对角矩阵`A`作为输入，然后通过迭代更新特征向量`v`来逼近特征值。注意，你需要手动构造你的三对角矩阵`A`，因为在这里我们没有直接给出具体的矩阵形式。