1 close all;
 2 clear all;
 3 clc;
 4 
 5 m=1.8;              %放大或缩小的高度
 6 n=2.3;              %放大或缩小的宽度
 7 img=imread('lena.jpg');
 8 imshow(img);
 9 [h w]=size(img);
10 imgn=zeros(h*m,w*n);
11 rot=[m 0 0;0 n 0;0 0 1];                                   %变换矩阵for i=1:h*m
12     for j=1:w*n
13         pix=[i j 1]/rot;   
14         
15         float_Y=pix(1)-floor(pix(1)); 
16         float_X=pix(2)-floor(pix(2));
17        
18         if pix(1) < 1　　　　　　　　%边界处理
19             pix(1) = 1;
20         end
21         
22         if pix(1) > h
23             pix(1) = h;
24         end
25         
26         if pix(2) < 1
27             pix(2) =1;
28         end
29         
30         if pix(2) > w
31             pix(2) =w;
32         end
33         
34         pix_up_left=[floor(pix(1)) floor(pix(2))];　　　　%四个相邻的点
35         pix_up_right=[floor(pix(1)) ceil(pix(2))];
36         pix_down_left=[ceil(pix(1)) floor(pix(2))];
37         pix_down_right=[ceil(pix(1)) ceil(pix(2))];     
38     
39         value_up_left=(1-float_X)*(1-float_Y);　　　　　　%计算临近四个点的权重
40         value_up_right=float_X*(1-float_Y);
41         value_down_left=(1-float_X)*float_Y;
42         value_down_right=float_X*float_Y;
43 　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　%按权重进行双线性插值
44         imgn(i,j)=value_up_left*img(pix_up_left(1),pix_up_left(2))+ ...
45                   value_up_right*img(pix_up_right(1),pix_up_right(2))+ ...
46                   value_down_left*img(pix_down_left(1),pix_down_left(2))+ ...
47                   value_down_right*img(pix_down_right(1),pix_down_right(2));        
48     end
49 end
50 
51 figure,imshow(uint8(imgn))