matlab代码耗时很长，如何优化？？

在做图像处理的大作业，写的双线性插值放大的函数运行速度很慢（100*100的放大两倍大概就得10秒左右吧），变换矩阵的地方修改了一点还是耗时不少。附件里是 .m文件，哪位大神写过这方面的程序求看一下，万分感谢！！！
附代码：

function fangda_bilinear(im_before,size_N)       %函数，输入原始图像和放大倍数[m,n],分别为高度和宽度放大。
N=size(im_before);                               %原始图像的大小
im_after=zeros(size_N(1)*N(1),size_N(2)*N(2));   %放大后矩阵的大小
trans_matrix=[size_N(1) 0 0;0 size_N(2) 0;0 0 1];    %变换矩阵，插值用
for k=1:1:3
for i=1:1:size_N(1)*N(1)                     
   for j=1:1:size_N(2)*N(2)                           
     zb_cz=[i j 1]/trans_matrix;     %插值位置的坐标（在原始图像上插值，有小数）
     if zb_cz(1)<1
         zb_cz(1)=1;                 %上界
     end
      if zb_cz(2)<1                  %左界
         zb_cz(2)=1;
      end


      zb_cz_upleft=[floor(zb_cz(1)) floor(zb_cz(2))];  %左上
      zb_cz_upright=[floor(zb_cz(1)) ceil(zb_cz(2))];  %右上
      zb_cz_downleft=[ceil(zb_cz(1)) floor(zb_cz(2))]; %左下
      zb_cz_downright=[ceil(zb_cz(1)) ceil(zb_cz(2))]; %右下
      %以上找到离最近的四个点坐标
      x=zb_cz(2)-floor(zb_cz(2));
      y=zb_cz(1)-floor(zb_cz(1));

     %归一化坐标，简化计算方法参考ppt25页
      im_after(i,j,k)=(im_before(zb_cz_upright(1),zb_cz_upright(2),k)-im_before(zb_cz_upleft(1),zb_cz_upleft(2),k))*x...
                     +(im_before(zb_cz_downleft(1),zb_cz_downleft(2),k)-im_before(zb_cz_upleft(1),zb_cz_upleft(2),k))*y...
                     +(im_before(zb_cz_downright(1),zb_cz_downright(2),k)+im_before(zb_cz_upleft(1),zb_cz_upleft(2),k)...
                     -im_before(zb_cz_downleft(1),zb_cz_downleft(2),k)-im_before(zb_cz_upright(1),zb_cz_upright(2),k))*x*y...
                     +im_before(zb_cz_upleft(1),zb_cz_upleft(2),k);

   end
end
end
imshow(uint8(im_after));
end