请用matlab语言写出图像双三次插值代码

好的，以下是matlab语言实现的双三次插值代码：

function [outimg] = bicubic_interp(inimg, scale)

% 获取输入图像的尺寸
[h, w, c] = size(inimg);

% 计算输出图像的尺寸
hout = round(h * scale);
wout = round(w * scale);

% 初始化输出图像
outimg = zeros(hout, wout, c);

% 计算插值系数矩阵
A = [-0.5, 1.5, -1.5, 0.5;
      1, -2.5, 2, -0.5;
     -0.5, 0, 0.5, 0;
      0, 1, 0, 0];
  
% 对于每个通道（RGB）
for channel = 1 : c
    % 对于输出图像的每个像素
    for i = 1 : hout
        for j = 1 : wout
            % 计算输入图像的坐标
            x = i / scale;
            y = j / scale;
            
            % 计算四个相邻像素的坐标
            x1 = floor(x);
            y1 = floor(y);
            x2 = x1 + 1;
            y2 = y1 + 1;
            
            % 处理边界像素
            if x1 < 1
                x1 = 1;
            end
            if x2 > h
                x2 = h;
            end
            if y1 < 1
                y1 = 1;
            end
            if y2 > w
                y2 = w;
            end
            
            % 计算插值系数
            dx = x - x1;
            dy = y - y1;
            px = [dx^3, dx^2, dx, 1];
            py = [dy^3; dy^2; dy; 1];
            Q = [inimg(x1, y1, channel), inimg(x1, y2, channel), inimg(x2, y1, channel), inimg(x2, y2, channel)];
            coef = px * A * Q * A' * py;
            
            % 将插值结果保存到输出图像中
            outimg(i, j, channel) = coef;
        end
    end
end

% 将输出图像的像素值限制在 [0, 1] 范围内
outimg(outimg < 0) = 0;
outimg(outimg > 1) = 1;
end

这段代码实现了对输入图像进行双三次插值的功能，其中 `inimg` 是输入图像，`scale` 是放大倍数，`outimg` 是输出图像。