DES算法Matlab代码.doc

DES算法Matlab代码
functionC=DES(D,K,M)
%输入16进制数表示的字符串格式明文D和密钥K,M=1表示加密
%%----------------------------------数据初始化-----------------------------------
%将以字符串形式输入的明文和密钥转换成01数字串
DB=[];
fori=1:16
Di=D(i);
DBi=['0000',dec2bin(hex2dec(Di))];
DBi=DBi(end-3:end);
DBi=[str2num(DBi(1)),str2num(DBi(2)),str2num(DBi(3)),str2num(DBi(4))];
DB=[DB,DBi];
end
D=DB;
KB=[];
fori=1:16
Ki=K(i);
KBi=['0000',dec2bin(hex2dec(Ki))];
KBi=KBi(end-3:end);
KBi=[str2num(KBi(1)),str2num(KBi(2)),str2num(KBi(3)),str2num(KBi(4))];
KB=[KB,KBi];
end
K=KB;

%64位数据IP置换表
IP=[585042342618102;
605244362820124;
625446383022146;
645648403224168;
57494133251791;
595143352719113;
615345372921135;
635547393123157];
IP0=reshape(IP',1,64);%转换成行向量
D_IP=D(IP0);%对D初始置换

%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];
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];