Simple DES加密，简单的DES加密的Matlab实现

来自信息安全的作业，包含两个文件，读者需要自己创建.只进行了两轮des加密，输入为字符

输出为8位的数字

密文是“ABCD”

输出结果：

 

代码

% 实现3-des 算法
% 2019-5-15
% clear all;
text='ABCD';
% ascii码65 66 67 68
P10=[3 5 2 7 4 10 1 9 8 6];
P8=[6 3 7 4 8 5 10 9];
P4=[2 4 3 1];
K=[1 0 1 0 0 0 0 0 1 0];
IP=[2 6 3 1 4 8 5 7];
IP_len=length(IP);
IPinv=[4 1 3 5 7 2 8 6];
EP=[4 1 2 3 2 3 4 1];
S0=[1 0 3 2;3 2 1 0;0 2 1 3;3 1 0 2];
S1=[0 1 2 3;2 0 1 3;3 2 1 0;2 1 0 3];
% S0ch=['01001110';'11100100';'00100111';'11010010'];
% S1ch=['00011011';'10000111';'11100100';'10010011'];
%% 先获得密钥K1与K2
matrixP10=K(P10);
lenP10=length(matrixP10);
matrixP10_L=matrixP10(1:lenP10/2);
matrixP10_Llen=length(matrixP10_L);
% 左移一位
matrixP10_L=matrixP10_L([2:matrixP10_Llen,1]);
matrixP10_R=matrixP10(lenP10/2+1:lenP10);
matrixP10_Rlen=length(matrixP10_R);
% 右移一位
matrixP10_R=matrixP10_R([2:matrixP10_Rlen,1]);
K1temp=[matrixP10_L,matrixP10_R];
K1=K1temp(P8);
% 左右进行移两位
matrixP10_L=matrixP10_L([3:matrixP10_Llen,1:2]);
matrixP10_R=matrixP10_R([3:matrixP10_Rlen,1:2]);
K2temp=[matrixP10_L,matrixP10_R];
K2=K2temp(P8);
%% 进行F函数的实现   dec2bin() 实现数据到二进制的转换
text_len=length(text);
text_matrix=zeros(text_len,8);      % text_matrix记录文本转换矩阵
standN=8;
for i=1:text_len        % 将文本转化为八位的二进制
    temp=dec2bin(text(i));
    temp_len=length(temp);
    for j=1:standN-temp_len
        temp=['0',temp];
    end
    tempbin=changebin(temp);
    text_matrix(i,:)=tempbin;
end
% 进行两轮循环
count=2;
P4bin=zeros(1,4);   % 存储二进制的P4矩阵
mainK=[K1;K2];
coding=zeros(size(text_matrix));% 初始化加密文本
N=1;
for i=1:text_len
    temprow=text_matrix(i,:);
    temp=temprow(IP);
    for N=1:count       
        IP_L=temp(1:IP_len/2);
        IP_R=temp(IP_len/2+1:IP_len);
        tempIP_R1=IP_R(EP(1:4));
        tempIP_R2=IP_R(EP(5:8));
        EP8=[tempIP_R1,tempIP_R2];  % 完成EP8矩阵的扩展 
        result1=xor(EP8,mainK(N,:));
        S0L=result1(1:4);
        S1R=result1(5:8);
        S0output=S0(2*S0L(1)+S0L(4)+1,2*S0L(2)+S0L(3)+1);
        S1output=S1(2*S1R(1)+S1R(4)+1,2*S1R(2)+S1R(3)+1);
        P4input=[S0output,S1output];    % 完成S矩阵数字到二进制的转换
        P4bin(1)=floor(P4input(1)/2);
        P4bin(2)=mod(P4input(1),2);
        P4bin(3)=floor(P4input(2)/2);
        P4bin(4)=mod(P4input(2),2);
        P4output=P4bin(P4);
        SWinput=xor(IP_L,P4output);
        if N==1         % 完成SW模块
            temp=[IP_R,SWinput];
        end
    end
    tempcoding=[SWinput,IP_R];
    coding(i,:)=tempcoding(IPinv);  
end
disp(coding);

 

% 第二个函数

% 该函数主要是将字符串的 ‘0101‘编程数字0 1 0 1
% 2019-5-15
function temp=changebin(str)
strlen=length(str);
temp=zeros(1,strlen);
for i=1:strlen
    if str(i)=='0'
        temp(i)=0;
    elseif str(i)=='1'
        temp(i)=1;
    end
end