毕业设计（论文 ）基于二维Arnold变换非等长数字图像置乱算法研究.doc

《毕业设计（论文 ）基于二维Arnold变换非等长数字图像置乱算法研究.doc》由会员分享，可在线阅读，更多相关《毕业设计（论文 ）基于二维Arnold变换非等长数字图像置乱算法研究.doc（24页珍藏版）》请在三一办公上搜索。

1、摘 要随着网络技术的发展，越来越多的图像信息在网络中传输，使得图像信息的安全性尤为重要。保障图像信息安全性最常用的方法就是图像置乱，图像置乱有许多置乱算法，Arnold变换因其简单而广为使用。现在对于Arnold变换通常应用于正方形图像置乱的研究。而对于高宽不等的矩形图像，现有的方法是将Arnold变换推广后对高宽满足一定条件的图像进行置乱，但它不适用于所有高宽不等的图像；另一种方法是将矩形图像按长边扩充为正方形图像后用Arnold变换进行置乱，但该算法改变了图像的属性。针对现有算法存在的问题提出了一种新颖的方法，该方法采用划分的方式使得非等长图像可应用Arnold变换进行图像置乱。对非等长图

2、像只需根据高与宽计算出划分时的重合区域，按短边为边长对非等长图像划分为几个正方形图像，然后对每个正方形图像依次运用Arnold变换进行置乱。置乱后的图像从反方向依次对划分的正方形图像运用Arnold反变换将图像恢复。最后，在MATLAB下对算法进行了仿真实验。实验结果表明，该算法具有很好的置乱效果，同时该算法可广泛应用于非等长图像的置乱变换过程。关键词：数字图像；图像置乱；Arnold变换；非等长图像ABSTRACTWith the development of network technology, more and more image information is transmitted

3、 in the Internet, which makes the security of image information particularly important. The most common method to ensure the security of image information is image scrambling. There are many hashing algorithms, of which Arnold transform is widely used for its simplicity. Now Arnold transform is appl

4、ied in the study of image scrambling of square. For the rectangular images with different height and width, the current approach is to scramble images which meet the required height and width. But it does not apply to all images especially with different height and width. Another method is to expand

5、 rectangular images to square by lengthen the edge. Then using Arnold transform to scramble the square image. But this hashing algorithm has changed the attribute of the image. Due to the problems existing in the current algorithm, a novel method has been put forward, which applies a way of dividing

6、 tubes to make images use Arnold transform. For those images with different length, we only need to calculate superposition region, and then use Arnold transform to scramble each square image. The scrambled image should be restored by using reverse Arnold transform from reverse direction. In the end

7、, the author carried out simulation experiment under MATLAB. The results show that this algorithm has very good effect on scrambling. At the same time, it can be widely used in the process of scrambling transformation of images with different lengths.Key words: Digital Image; Image Scrambling; Arnol

8、d Transformation; The Tubes image目 录第0章 引 言10.1 图像置乱技术概述10.1.1 图像置乱以及图像置乱技术10.1.2 图像置乱技术的应用10.1.3 图像置乱的意义20.2 本文的研究意义及结构30.2.1 本文的研究意义30.2.2 本论文的结构安排3第1章 基于二维ARNOLD变换的非等长数字图像置乱算法41.1 二维Arnold变换41.2 现有算法41.2.1 算法一：基于扩充图像的方法41.2.2 算法二：基于二维Arnold变换推广式的方法51.2.2.1 二维Arnold变换的推广51.2.2.2 推广式的周期存在性问题51.3 现有

9、算法中的不足61.4 基于区域划分方法的非等长数字图像置乱61.4.1 划分算法61.4.2 置乱恢复8第2章 仿真实验112.1 算法的有效性实验112.2 置乱恢复实验112.3 算法安全性实验122.4 与现有算法一比较实验12第3章 结束语14参考文献15谢 辞16附 录17第0章 引 言0.1 图像置乱技术概述0.1.1 图像置乱以及图像置乱技术 图像是人类描述客观世界的最有效的手段之一，凡是记录在纸上的,拍摄在照片上的,显示在屏幕上的所有具有视觉效果的画面都可以称之为图像。数字图像可以用一个矩阵来表示，矩阵的元素所在的行和列,就是图像显示在计算机屏幕上各像素点的坐标,元素的数值就是

10、该像素的灰度值。图像置乱,顾名思义,就是把图像打乱，从而隐藏原始图像的真实内容。对于数字化的图像，置乱过程不仅可以在数字图像的空间域（色彩空间、位置空间）上进行，而且可以在数字图像的频域上进行。数字图像位置空间上的置乱就是几何变换，即只改变像素的几何位置，而不改变像素的灰度值，因此图像置乱一次的效果取决于几何变换作用的范围大小。对于数字化的图像置乱过程不仅可以在数字图像的空间域（色彩空间、位置空间）上进行,还可以在数字图像的频域上进行。数字图像置乱即是对数字图像的一种加密方法，它使得合法使用者可以自由控制算法的选择、参数的选择以及使用随机数技术，这就加大了攻击者非法破译的难度。空间域的图像置乱

11、是利用某种算法将一幅图像各像素的次序打乱，但像素的总个数不变，直方图不变，使一幅图像变得“面目全非”。数字图像置乱技术,可以看作数字图像加密的一种途径,也就是改变图像中像素的位置或者像素的颜色,将原始图像变换成一个无内容、无纹理、无形状的杂乱无章的新图像的过程。0.1.2 图像置乱技术的应用随着网络技术和多媒体的飞速发展，数字图像信息作为一种重要的资料，被广泛地在网络上传播。人们通过这种新的媒体互相传递信息，简单快捷，但随之而来的是这些信息的安全隐患，越来越多的图像信息被泄露、纂改和假冒，给个人隐私、公司利益，甚至国家安全都带来了重大的威胁。对数字图像的安全性这个领域的主要研究方向有数字图像置

12、乱技术、数字图像分存技术、数字图像隐藏技术和数字水印技术。因此，网络信息加密成为首要研究的问题，而保障图像信息安全最常用、最直接的方法就是加密当前对图形的加密所采用的方法主要是对图像像素点位置进行置乱,或对图像像数值进行置混，或将置乱和置混结合起来。0.1.3 图像置乱的意义数字图像信息安全是数学、密码学、信息论、计算机视觉以及其它计算机应用技术的多学科交叉的研究课题。数字图像置乱作为图像信息加密的一种方法或图像信息隐藏的预处理手段，越来越受到关注，它起源于早期的经典加密学理论和电视图像应用技术。经典密码学对一维数据流提了较好的加密算法，但对数字图像的处理却不能忽视其固有的一些特殊性质，如二维

13、的自相似性、相关性等。图像是含有大量信息的信息载体，且包含大量的冗余信息，作为直观的信息表达方式，具有很大的迷惑性。如果我们对数字化图像做一些“扰乱”，得到一幅完全杂乱无章、面目全非的图像，那么即使非法截获者注意到它，如果不知道如何恢复，对它也无能为力，这在一定程度上增强了图像的安全性。置乱技术就是这样一项研究课题，它从一维的单表密码扩展而来，应用到二维图像平面、甚至三维图像色彩空间中，扰乱图像的组成部分，破坏图像的自相关性，使其所要表达的真实信息无法直观地得到，即使计算机用“穷举法”计算各种组合，也要耗费大量的时间，在一定程度上保护了图像信息。通过置乱数字图像可以达到两个目的：第一是加密处理

14、，图像置乱变换是一种基于内容的图像加密方法，就像不知道加密密钥对加密过的信息进行解密一样，如果不知道置乱所采用的算法，同样难以恢复原始图像；第二是作为信息隐藏的预处理手段可以增强图像伪装的鲁棒性,首先，将图像置乱后，将得到一幅杂乱无章的图像,这样的图像无内容、无纹理、无形状，从中无法读取有意义的信息,那么,将这样一幅图像嵌入到另一幅普通图像时就不容易引起那幅图像在内容、纹理、形状上的太大改变，甚至不会发生改变，这样人眼就不易识别，从而逃出了第三方的视线；其次，由于秘密图像是置乱后的图像，根据上述图像的“三无”特征，第三方难以对其进行内容、纹理、形状等的统计分析，即便他们截取到了秘密图像，也是无

15、能为力的。如果第三者企图对秘密图像进行反置乱，这也是非常困难的，由于图像置乱有很多种方法，每种方法又可以使用不同置乱模板算法，设置不同的参数，使用者有很大的自由度，他可以根据自己的想法得到不同的结果。相反，这给企图截获秘密信息的第三方带来了很大的困难，使他们需要耗费巨大的计算量来穷举测试各种可能性；最后，如果第三方反置乱不成，在隐蔽载体上作了恶意修改，合法接收者在对接受到的秘密图像进行反置乱的过程中，就使得第三方在图像上所涂、画的信息被分散到画面的各个地方，从而形成了点状的随机噪声,不会对视觉效果造成较大的影响。0.2 本文的研究意义及结构0.2.1 本文的研究意义其中Arnold变换算法简单

16、且置乱效果好而在图像中应用最广泛。目前Arnold 变换的研究主要集中在正方形图像置乱的方面,针对非正方形的置乱尚很少。目前已有的方法要么将非正方形图像扩展为正方形图像后应用Arnold 变换，要么将Arnold 变换推广针对非高宽比例满足一定条件的图像进行置乱。本文提出了一种Arnold 变换应用于非等长图像置乱的方法，该方法在不改变图像形状的条件下对任意高宽比例的长方形图像应用Arnold 变换进行置乱过程,扩展了2维图像置乱变换在非等长图像上的应用。0.2.2 本论文的结构安排本文的主要结构包括三个方面,第一部分：介绍了Arnold变换公式,及其置乱过程。其次,阐述了Arnold变换公式

17、不能直接用于非等长图像的原因。第二部分，介绍了现有算法并分析了现有算法存在的问题。第三部分，提出了给予划分思想的非等长图像的Arnold 变换置乱方法,并对其置乱恢复所需的反变换算法进行了阐述。第四部分：给出了实验仿真结果。对本文所提出的划分思想及划分公式进行了验证。第1章 基于二维Arnold变换的非等长数字图像置乱算法1.1 二维Arnold变换二维Arnold变换也就是通常所说的Arnold变换,俗称“猫脸变换”,它是V.I Arnold提出的。定义1 设图像矩阵大小为,像素坐标为且.若映射为满足式（1）,则成为二维等长Arnold变换，简称Arnold变换。 （1）称为Arnold变换

18、矩阵，为数字图像矩阵的阶。将Arnold变换应用在数字图像上,可以通过像素坐标的改变而改变图像灰度值的布局，把数字图像看作一个矩阵，则经Arnold变换变换后的图像会变得“混乱不堪”,这就是图像置乱的概念。1.2 现有算法对于二维Arnold变换非正方形图像现在只有两种算法。算法一将非等长图像扩展为正方形图像,算法二是将Arnold变换推广后,找到一种高宽满足一定条件的非等长图像可运用Arnold变换进行图像置乱处理。1.2.1 算法一：基于扩充图像的方法一般地，Arnold变换都是对正方形图像进行处理。对于非等长图像，由于长度和宽度不相等，分别对长度和宽度的模运算就不一样了，所以在此处Arn

19、old变换就不再是一一映射。因此,孔涛在文献8中提出把矩形图像转化为正方形图像，其基本思想是：把原来的矩形图像扩充成正方形图像,然后再用Arnold变换进行图像置乱。如果长度大于宽度,则按长度的大小扩充成正方形图像；如果宽度大于长度,则按宽度大小扩充成正方形图像，见图1.1。 图1.11.2.2 算法二：基于二维Arnold变换推广式的方法1.2.2.1 二维Arnold变换的推广李永逵,冯乔生等人在文献10中将二维等长Arnold变换推广以适用非等长图像。定义2 设图像矩阵大小为，像素坐标为，且，。若映射为满足式（2），则称为二维非等长图像置乱变换。 （2）1.2.2.2 推广式的周期存在性

20、问题通过二维等长Arnold变换和二维非等长图像置乱变换的定义发现,对于非等长图像若想直接运用Arnold变换将图像置乱，就必须使得周期存在,而对于非等长Arnold变换周期存在必须满足一定的条件。定理1 二维非等长置乱变换周期存在，当且仅当像素矩阵所有元素，按式（2）一次迭代后，映射为新像素矩阵的不同元素。定理2 设像素矩阵为维，对于所有不同为0的整数，若式（3）和式（4）不同时成立,则二维非等长图像置乱变换周期存在。 （3） （4）定理3 二维非等长Arnold变换周期存在的充要条件是M/N为整数或N/M为奇数。邵利平,覃征等人在文献9-10中对上述三个定理进行了详细的证明,这里就不在进行

21、证明了。1.3 现有算法中的不足算法一把非等长图像扩充为等长图像则改变了图像的属性。算法二当所取M,N值满足充要条件时周期都存在。而当所取M,N值不满足充要条件时，都求不出周期,当然更无法进行图像置乱变换了。当M和N值不满足充要条件时,便可能存在两点同时映射到一点，这样变换周期就不存在了。1.4 基于区域划分方法的非等长数字图像置乱由于Arnold变换只能用于正方形图像置乱，对于不满足一一映射关系的非正方形图像并不能直接运用Arnold变换实现图像置乱,为了使的非正方形图像能用Arnold变换对图像实现置乱,就必须将的非等长图像划分成多个正方形图像，再运用Arnold变换对图像进行置乱,如何对

22、非正方形图像进行是下面将要讨论的问题。1.4.1 划分算法1、对于的非等长图像，若，当时，把的非等长图像划分为一个和的两部分图像（见图1.2）,运用Arnold变换将第一部分的图像置乱后，发现第二部分的图像将无法运用Arnold变换进行置乱（见图1.3），为了将这一部分图像运用Arnold变换置乱,只有将的图像补成的正方形图像。如果在图像后面将其补成的正方形图像，则就改变了图像的属性,所以从前面的正方形图像中截取一部分图像将的图像补成的正方形图像。这样就可以运用Arnold变换对第二部分图像进行置乱,从而达到全局置乱的效果。（见图1.4）。2、当时，将图像划分成和的两部分图像（见图1.5）。同

23、样第一部分可以直接运用Arnold变换,第二部分则为上述1的情形，那么可以用上述1的方法对图像进行置乱（见图1.6）。我们发现图像置乱后有明显的区域性。这是由于第一部分的点在置乱后还在第一部分，不会出现在第二部分。同理，第二部分的点在置乱后也还在第二部分中而不可能出现在第一部分中。为了使得的图像在置乱后没有区域性，现把第一部分置乱后的一些点放在第二部分进行置乱。这样的非等长图像以第二部分和前一部分有一定的重合区域的方式进行划分，直到把的非等长图像置乱完全为止（见图1.7）。NNM-NNNM-NNNM-NN 图1.2 图1.3 图1.4NN NNN N NM-N NM-NN 图1.5 图1.6

24、图1.73、当时，为了将的非等长图像进行全局置乱且没有区域性，就按照上述2的方法，先从短边的始端划分一个的正方形图像运用Arnold变换将图像置乱，接着从第二部分开始每部分将和前一部分有一定的重合区域的方式进行划分，依次运用Arnold变换讲话份额正方形图像置乱，直到把的非等长图像置乱完全为止（见图1.8）。MNNN 图1.8用上述方法对图像进行划分,重合区域有很多种情况，即原图像可以划分出正方形图像的个数不同。划分出的正方形个数越少，图像置乱所需时间也就越少。为使得图像划分成尽可能少的正方形,从而大量减少图像置乱所需要的时间，当把非等长图像划分为个正方形图像时，置乱所需的时间是最少。以短边N

25、为边长的划分成个正方形图像。设每个正方形图像之间重合长度为X，则可列方程为： 其中表示的是不超过M/N的最大整数，表示进行正方形所需的长,表示一共所需的重合的长度。求得 。1.4.2 置乱恢复二维Arnold变换的周期性是一个很好的性质,当反复应用Arnold变换时,在某一时刻就能恢复原图。但是对于本文所提的方法有无周期性还不知道,从而我们选择了另一种方法恢复非等长图像，即Arnold反变换。非等长图形对于本文所提出的划分方法利用反变换进行恢复。根据反变换对置乱后图像从较长边的末端开始划分正方形区域，将从反方向也就是若把的图像看作是的矩阵，则将从第列开始对图像进行划分，重合长度还是。逐一对每个

26、正方形进行恢复。刘芳在文献12中对反变换进行了阐述,其证明过程如下：因，由矩阵的性质可知，A可逆,。由公式(1)，根据矩阵运算的性质，我们来推测Arnold反变换的公式为,其中mod代表取余数,floor表示向负无穷方向取整,a、b、c为整数。证明：根据公式(1),Arnold变换的等价形式可以写成:这意味着,Z为整数则即 由上可知， 无论和是正整数、负整数或是0,都有即所以 综上所述，Arnold反变换的公式为： (5)第2章 仿真实验2.1 算法的有效性实验分别选取图像一为150 200（jpg图）、图像二为150 400（jpg图）、图像三为400 250（jpg图）进行实验,实验结果如

27、图2.1：图2.1置乱次数原 始 图 像置 乱 1 次置 乱 2 次置 乱 3 次置 乱 4 次图像一图像二图像三通过实验可以看到图像的置乱刚开始具有区域性，但在多次置乱后区域性将逐渐消失，置乱效果逐步达到最佳。通过选取三组不同高宽的矩形进行置乱，发现图像一开始分为两部分置乱，图像二分为三部分置乱，图像三分为两部分置乱,且到最后都达到了不错的置乱效果。从而说明对于任意高宽不等的矩形图像,都可通过本文提出的划分方法将图像置乱。2.2 置乱恢复实验选取图像一和图像二置乱后的图像进行Arnold反变换实验,实验结果如图2.2：图2.2置乱次数反 变 换 1 次反 变 换 2 次反 变 换 3 次反

28、变 换 4 次反 变 换 5 次图像一图像二实验表明通过Arnold反变换从图像的较长边的末端开始恢复,可以回到原图像,说明该算法具有可恢复性。2.3 算法安全性实验选取图像一置乱后的图像，使用错的重合度对图像恢复，实验结果如图2.3：图2.3置乱次数反 变 换 1 次反 变 换 2 次反 变 换 3 次反 变 换 4 次反 变 换 5 次图像一通过实验表明,在对图像进行恢复时,若使用错误的重合度,那么图像将无法恢复到原图像,从而说明该算法是安全的。2.4 与现有算法一比较实验选取图像一，利用扩展为正方形图像的算法和本文提出的划分算法分别对图像一进行实验,实验结果如图2.4：图2.4置乱次数原

29、 始 图 像置 乱 1 次置 乱 2 次置 乱 3 次反变换到原图像扩充算法存储空间6.15KB10.0KB15.7KB19.4KB6.73KB置乱次数原 始 图 像置 乱 1 次置 乱 2 次置 乱 3 次反变换到原图像划分算法存储空间6.15KB10.6KB14.5KB15.5KB6.15KB实验表明：用Arnold变换进行处理后的结果,我们发现将非等长图像进行划分后，将图像置乱变换后比将非等长图像扩展为等长图像进行置乱后所占用的存储空间要少得多,且置乱效果好，并且划分图像的算法并没有改变图像的属性。第3章 结束语传统的二维Arnold变换图像置乱一般用于正方形图像置乱,而对于高宽不等的矩

30、形图像,需将其扩展为正方形图像进行置乱。对此问题,本文提出了一种新颖的2维非等长图像置乱变换方法,将其按照短边N为边长,以为重合度进行划分正方形，然后利用Arnold变换进行图像置乱，并运用Arnold反变换将其还原。同时利用了数学思想中的方程思想进行了严谨的推理。实验表明了所提出的新的二维非等长Arnold变换算法可以将非等长图像进行置乱并且恢复。该算法克服了将非等长图像扩展为等长图像后需要增加存储空间的缺点,并且该算法的置乱效果比扩展为正方形图像算法的置乱效果要好。下一步将讨论划分算法的周期性。参考文献1 丁玮，闫伟齐,等.基于Arnold变换的数字图像置乱技术J.计算机辅助设计与图形学学

31、报.2001.13(4):338-3412 王冬梅.奇数阶幻方变换数字图像的准周期J.浙江工业大学学报.2005. 33（3）:292-2943 丁玮,齐东旭.数字图像变换及信息隐蔽与伪装技术J.计算机学报.1998. 21(9):838-8434 林雪辉，蔡利栋.基于Hilbert曲线的数字图像置乱方法研究J.中国体视学与图像分析.2004.9（4）:224-2275 齐东旭.矩阵变换及其在图像信息隐蔽中的应用研究J.北方工业大学学报.1999.11(1):24-286 丁玮,闫伟齐等.基于生命游戏的数字图像置乱与数字水印技术J.北方工业大学学报.2000.12(1):1-57 邹建成,李国

32、富等.广义Gray码及其在数字图像置乱中的应用J.高校应用数学学报（A辑）.2002.17(3):363-3708 孔涛. Arnold反变换的一种新算法J.软件学报.2004.15（10）:1558-15649 邵利平,覃征.二维非等长图像置乱变换J.电子学报.2007.35(7):1209- 129410 李永逵,冯乔生等.二维Arnold变换及非等长图像置乱变换J.计算机工程与设计. 2009.30(13):3133-313511 齐东旭,邹建成等.一类新的置乱变换及其在图像信息隐藏中的应用.中国科学(E辑).2000.30(5):440-44712 刘芳.一种新的Arnold反变换在数

33、字水印中的应用J.全国图象图形学学术会议.2006.7:1-6谢 辞经过两个月的努力，终于完成了我的毕业论文。在没有做毕业设计以前觉得毕业设计只是对这几年来所学知识的单纯总结，但是通过这次做毕业设计发现自己的看法有点太片面。毕业设计不仅是对前面所学知识的一种检验，而且也是对自己能力的一种提高。通过这次毕业设计使我明白了自己原来知识还比较欠缺。自己要学习的东西还太多，以前老是觉得自己什么东西都会，什么东西都懂，有点眼高手低。通过这次毕业设计，我才明白学习是一个长期积累的过程，在以后的工作、生活中都应该不断的学习，努力提高自己知识和综合素质。我的毕业论文的完成是在我的指导老师郭琳琴老师的帮助下完成

34、的。她从选题的确定、论文资料的收集、论文框架的确定、开题报告准备及论文初稿与定稿中对字句的斟酌倾注的大量心血，郭琳琴老师治学严谨，在论文写作上对我认真指导、严格要求。在此对指导老师- 郭琳琴老师表示感谢！ 在这里，还要特别感谢大学四年学习期间给我诸多教诲和帮助的数学系的各位老师，你们给予我的指导和教诲我将永远记在心里！师恩永铭，友情难忘。感谢和我一起生活四年的室友，是你们让我们的寝室充满快乐与温馨，愿我们以后的人生都可以充实、多彩与快乐！感谢我的兄弟姐妹们-我的同学和朋友们，他们给我提供了不少的参考资料和相关建议。感谢答辩组老师们的辛勤劳动。更要感谢我的父母，他们虽然不在我的身边，但四年来，他

35、们无时无刻不在关注着我一点一滴的成长。附 录Arnold变换的程序：a=imread(24.jpg);ax=double(a);bx=zeros(200,200);for i=1:200 for j=1:200 ik=i+j; jk=i+2*j; ik=mod(ik,200); if ik=0 ik=200; end jk=mod(jk,200); if jk=0 jk=200; end bx(ik,jk)=ax(i,j); endendimshow(bx,);b=uint8(bx);imwrite(b,25.jpg);图2.1程序：a=imread(2_3.jpg);tic;ax=doubl

36、e(a);cx=zeros(150,200);cx=ax(1:150,1:200);for h=1:50:51dx=zeros(150,150);dx=cx(:,h:h+149);bx=zeros(150,150); for i=1:150 for j=1:150 ik=i+j; jk=i+2*j; ik=mod(ik,150); if ik=0 ik=150; end jk=mod(jk,150); if jk=0 jk=150; end bx(ik,jk)=dx(i,j); end endfor i=1:150 for j=1:150 cx(i,j+h-1)=bx(i,j); endend

37、endimshow(cx,);c=uint8(cx);toc;imwrite(c,2_4.jpg);图2.2反变换程序：a=imread(2_8.jpg);ax=double(a);cx=zeros(150,200);cx=ax(1:150,1:200);for h=200:-50:150dx=zeros(150,150);dx=cx(:,h-149:h);bx=zeros(150,150); for i=1:150 for j=1:150 ik=2*i-j; jk=-i+j; ik=mod(ik,150); if ik=0 ik=150; end jk=mod(jk,150); if jk=

38、0 jk=150; end bx(ik,jk)=dx(i,j); end endfor i=1:150 for j=1:150 cx(i,j+h-150)=bx(i,j); endendendimshow(cx,);c=uint8(cx);imwrite(c,3_7.jpg);将非等长图像扩展为等长图像的程序：a=imread(2.jpg);ax=double(a);bx=zeros(200,200);for i=1:150 for j=1:200 bx(i,j)=ax(i,j); endendcx=zeros(200,200);for i=1:200 for j=1:200 ik=i+j; jk=i+2*j; ik=mod(ik,200); if ik=0 ik=200; end jk=mod(jk,200); if jk=0 jk=200; end cx(ik,jk)=bx(i,j); endendimshow(cx,);b=uint8(cx);imwrite(b,22.jpg);