matlab DES加密 直接复制运行

%% 这是个.m文件,直接赋值运行,密钥密文自己改 %% === initial data === %% === 子密钥左移次数=== LS = [1 1 2 2 2 2 2 2 1 2 2 2 2 2 2 1] ; %% === import PC-1置换表=== 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 ] ; %% === import PC-2置换表=== 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 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-扩展表=== 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 10 2 8 24 14 32 27 3 9 19 13 30 6 22 11 4 25 ] ; %% === import IP逆置换表=== 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 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 = [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 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 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 13 15 3 5 8 2 1 14 7 4 10 8 13 15 12 9 0 3 5 6 11 ]; %% === input data === %key_in = input('请输入密钥(字符用单引号括起):') ; % 让用户输入 %Text_in = input('请输入明文(字符用单引号括起):') ; key_in = '加密密钥' % 在此输入你要加密的密钥,四个汉字 Text_in = '加密明文' % 在此输入你要加密的明文,四个汉字 %% === characters convert to binary === key_in = unicode2native(key_in, 'GB2312') % 编码转换,把unicode转为GB2312 Text_in = unicode2native(Text_in, 'GB2312') key_in = dec2bin(key_in) % 转换为二进制 Text_in = dec2bin(Text_in) key = '' ; Text = '' ; len = length(key_in) ; for i = 1:len % 对二进制重新排版,排为一行 key = [key, key_in(i, :)] ; Text = [Text, Text_in(i, :)] ; end key % 输出二进制密钥 Text % 输出二进制明文 %% === 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 ) ) ) ; % 需要把char 转换为number end % == 循环生成子密钥== sub_k = zeros(16, 48) ; % 保存所有的子密钥 % 把有效密钥均分为两份 L = k(1:28) ; R = k(29:56) ; for i = 1:16 n = LS(i) ; % 获得每次循环的左循环移位的次数 % == 左循环移位== 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 %% === 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, 48) ; % == E 扩展== tmp = R( E ) ; tmp = bitxor( tmp, sub_k(i, :) ) ; % R K 的异或运算,得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('错误:超出程序所接受的范围!') ; return ; end