基于Matlab实现的DES加密

基于Matlab实现的DES加密方法代码。直接运行后输入函数名称DESDECP(a,b,n),其中a，b分别是16进制的明文字符串和密钥，n为1表示加密，为其他数字表示解密。就可以进行DES加密和解密。

function cipher=DESDECP（DKM）

plaintext=D;
key=K;
DB=[];
for i=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=[DBDBi];
end
D=DB;
KB=[];
for i=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=[KBKBi];
end
K=KB;

%64位数据IP置换表
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];
IP0=reshape（IP‘164）;%转换成行向量
D_IP=D（IP0）;%对D初始置换

%E盒为数据扩展
E=[32 1 2 3 4 5;
    4 5 6 7 8 9;
    8 910111213;
   121314151617;
   161718192021;
   202122232425;
   242526272829;
   28293031321];
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];
%P盒置换
P=[1672021291228171152326518311028241432273919133062211425];
PC1=[5749413325179158504234261

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       6639  2012-12-27 17:08  DES\DESDECP.m

     目录          0  2012-12-29 16:12  DES

----------- ---------  ---------- -----  ----

                 6639                    2