DES加密解密matlab實(shí)現(xiàn)

1、% 注意這是一個(gè) function .m 文件， 直接復(fù)制即可function out, bin, hex = DES( Text, key, flag )% key and Text 是字符類型% key ='1011101010100011110100011111001110110100111100111101000110100111'% Text ='1011101010100011110100011111001110110100111100111101000110100111'% DES(Text, key)% when flag = 0 （默認(rèn)）， D

2、ES 加密% when flag = 1， DES 解密 % = initial data = % = import 左移位數(shù)表 = LS = 1 1 2 2 2 2 2 2 1 2 2 2 2 2 2 1 ; % = import PC-1 table= PC1 = 57 49 41 33 25 17 9 1 58 50 42 34 26 18 10 2 59 51 43 35 27 19 11 3 60 52 44 36 63 55 47 39 31 23 15 7 62 54 46 38 30 22 14 6 61 53 45 37 29 21 13 5 28 20 12 4 ; % =

3、 import PC-2 table = PC2 = 14 17 11 24 1 5 3 28 15 6 21 10 23 19 12 4 26 8 16 7 27 20 13 2 41 52 31 37 47 55 30 40 51 45 33 48 44 49 39 56 34 53 46 42 50 36 29 32 ; % = import IP table = IP =58 50 42 34 26 18 10 2 60 52 44 36 28 20 12 4 62 54 46 38 30 22 14 6 64 56 48 40 32 24 16 8 57 49 41 33 25 17

4、 9 1 59 51 43 35 27 19 11 3 61 53 45 37 29 21 13 5 63 55 47 39 31 23 15 7 ; % = import E-table = E = 32 1 2 3 4 5 4 5 6 7 8 9 8 9 10 11 12 13 12 13 14 15 16 17 16 17 18 19 20 21 20 21 22 23 24 25 24 25 26 27 28 29 28 29 30 31 32 1 ; % = import P table = P = 16 7 20 21 29 12 28 17 1 15 23 26 5 18 31

5、10 2 8 24 14 32 27 3 9 19 13 30 6 22 11 4 25 ; % = import IP table = IP_1 = 40 8 48 16 56 24 64 32 39 7 47 15 55 23 63 31 38 6 46 14 54 22 62 30 37 5 45 13 53 21 61 29 36 4 44 12 52 20 60 28 35 3 43 11 51 19 59 27 34 2 42 10 50 18 58 26 33 1 41 9 49 17 57 25 ; % = import S-Box = S1 = 14 4 13 1 2 15

6、11 8 3 10 6 12 5 9 0 7 0 15 7 4 14 2 13 1 10 6 12 11 9 5 3 8 4 1 14 8 13 6 2 11 15 12 9 7 3 10 5 0 15 12 8 2 4 9 1 7 5 11 3 14 10 0 6 13 ; S2 = 15 1 8 14 6 11 3 4 9 7 2 13 12 0 5 10 3 13 4 7 15 2 8 14 12 0 1 10 6 9 11 5 0 14 7 11 10 4 13 1 5 8 12 6 9 3 2 15 13 8 10 1 3 15 4 2 11 6 7 12 0 5 14 9 ; S3

7、 = 10 0 9 14 6 3 15 5 1 13 12 7 11 4 2 8 13 7 0 9 3 4 6 10 2 8 5 14 12 11 15 1 13 6 4 9 8 15 3 0 11 1 2 12 5 10 14 7 1 10 13 0 6 9 8 7 4 15 14 3 11 5 2 12 ; S4 = 7 13 14 3 0 6 9 10 1 2 8 5 11 12 4 15 13 8 11 5 6 15 0 3 4 7 2 12 1 10 14 9 10 6 9 0 12 11 7 13 15 1 3 14 5 2 8 4 3 15 0 6 10 1 13 8 9 4 5

8、 11 12 7 2 14 ; S5 = 2 12 4 1 7 10 11 6 8 5 3 15 13 0 14 9 14 11 2 12 4 7 13 1 5 0 15 10 3 9 8 6 4 2 1 11 10 13 7 8 15 9 12 5 6 3 0 14 11 8 12 7 1 14 2 13 6 15 0 9 10 4 5 3 ; S6 = 12 1 10 15 9 2 6 8 0 13 3 4 14 7 5 11 10 15 4 2 7 12 9 5 6 1 13 14 0 11 3 8 9 14 15 5 2 8 12 3 7 0 4 10 1 13 11 6 4 3 2

9、12 9 5 15 10 11 14 1 7 6 0 8 13 ; S7 = 4 11 2 14 15 0 8 13 3 12 9 7 5 10 6 1 13 0 11 7 4 9 1 10 14 3 5 12 2 15 8 6 1 4 11 13 12 3 7 14 10 15 6 8 0 5 9 2 6 11 13 8 1 4 10 7 9 5 0 15 14 2 3 12 ; S8 = 13 2 8 4 6 15 11 1 10 9 3 14 5 0 12 7 1 15 13 8 10 3 7 4 12 5 6 11 0 14 9 2 7 11 4 1 9 12 14 2 0 6 10

10、13 15 3 5 8 2 1 14 7 4 10 8 13 15 12 9 0 3 5 6 11 ; % = 輸出參數(shù)判斷 = if( 2 = nargin) flag = 0 ; end % = producing the sub cipher code (48 bit) = % = PC-1 = k = zeros(1, 56) ; % 56 bit n = length( PC1 ) ; % n =56 for i = 1: n k( i ) = str2num( key( PC1( i ) ) ) ; % ÐèÒª°Ñ ch

11、ar ת»»Îª number end % = 子密鑰生成 = sub_k = zeros(16, 48) ; % 保存子密鑰 % 評(píng)分密鑰 L = k(1:28) ; R = k(29:56) ; for i = 1:16 n = LS(i) ; % 左移位數(shù) % = 左移 = for j = 1:n L = L(2:end) , L(1) ; R = R(2:end) , R(1) ; end % = PC-2 = k = L, R ; sub_k(i, :) = k( PC2 ) ; end % = 判斷加密還是解密 = if

12、( 1 = flag) sub_k = sub_k(end:-1:1, :) ; end % = DES = % = 初始IP置換 = T = zeros(1, 64) ; % n = length( IP ) ; % n = 64 for i =1: n T( i ) = str2num( Text( IP( i ) ) ) ; end Text = T ; L = zeros(1, 32) ; R = zeros(1, 32) ; for i = 1:32 L(i) = Text(i) ; R(i) = Text(i+32) ; end for i = 1:16 tmp = zeros(1

13、, 48) ; % = E-擴(kuò)展= tmp = R( E ) ; tmp = bitxor( tmp, sub_k(i, :) ) ; % R K 異或運(yùn)算 48bit % = S-盒 = T = zeros(1, 32) ; for j = 1:8 switch j case 1 S = S1 ; case 2 S = S2 ; case 3 S = S3 ; case 4 S = S4 ; case 5 S = S5 ; case 6 S = S6 ; case 7 S = S7 ; case 8 S = S8 ; otherwise disp('´íÎ

14、;ó£º ³¬³ö³ÌÐòËù½ÓÊܵķ¶Î§£¡') ; return ; end ii = tmp( (j-1)*6 +1 )*2+tmp( (j-1)*6 +6 ) ; % h1*2+h6 jj = tmp( (j-1)*6 +2 )*8+tmp( (j-1)*6 +3 )*4+tmp( (j-1)*6 +4 )*2+tmp( (

15、j-1)*6 +5 ) ; % h2*8+h3*4+h4*2+h5 num = S( ii*16 + jj+1 ) ; % S-盒 tmp1 = 0 0 0 0 ; num = dec2bin( num ) ; % dec2bin() 輸出char類型 ii = length(num) ; for jj = 1:ii % ! improtant£¬char to num 類型轉(zhuǎn)換 tmp1( 4-ii+jj ) = str2num( num( jj ) ) ; end T( (j-1)*4+1 : j*4 ) = tmp1 ; end % = P置換 = T = T( P ) ; % 32bit T = bitxor(L, T) ; % L F(R, K) 異或運(yùn)算 % = swap L R = L = R ; R = T ; end % = IP 逆置換= k = R, L ; % 交換L R k = k(IP_1) ; % IP 逆置換 % = output data = % = binary = Text_o = '' ; n = length(