function [ef,result]=des(pf,pw,ind)
%函数形式:[ef,result]=des(pf,pw,ind)
%功能：用密码pw对pf做des加密
%输入：pf 一维向量 64位二进制长 8个英文字母长
%pw 一维向量 56位二进制长 7个英文字母长
%ind 实整数 取值范围1-16 指示返回第i轮加密过程相关中间值及结果
%说明：pf，pw要以字符串的形式输入。
%输出：er 一维向量 加密结果
irt=[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];
%irt 明文初始置换表
endrt=[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];
%ert 明文末置换表
ert=[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];
%ert 明文扩展置换表
iprt=[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];
%ipt 密钥初始置换表
pt=[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];
%p盒置换表
pfb=str2bin(pf,0,0);%明文转二进制
result=['明文',pf,'二进制表示：',pfb, sprintf('\n')];
pwb=str2bin(pw,7,2);%密文转二进制,并加入偶校验位
result=[result,'密码',pw,'二进制表示：',pwb,sprintf('\n')];
pfb=rebit(pfb,irt);%明文初始置换
result=[result,'明文',pf,'初始置换结果',pfb,sprintf('\n')];
pwb=rebit(pwb,iprt);
result=[result,'密钥',pw,'初始置换结果',pwb,sprintf('\n')];
lpfb=pfb(1:32);rpfb=pfb(33:64);
ki=gerki(pwb);%计算16轮子密钥
for i=0:15
if i<9
a1=48; a2=48;
else if i==9
a1=48; a2=55;
else a1=55;a2=55;
end
end
if ind==i+1
result=[result,'L',char(i+a1),':',lpfb,sprintf('\n'),'R',char(i+a1),':',rpfb,sprintf('\n')];
end
templ=lpfb;
lpfb=rpfb;%下一轮新的左数据输入
if ind==i+1
result=[result,'由L',char(i+a2+1),'=R',char(i+a1),'得L',char(i+a2+1),':',lpfb,sprintf('\n')];
end
tempr=rebit(rpfb,ert);%rpfb扩展置换成48位
if ind==i+1
result=[result,'按子密钥的生成过程，求得k',char(i+a1),':',ki(i+1,:),sprintf('\n')];
result=[result,'R',char(i+a1),'扩展成48位：',tempr,sprintf('\n')];
end
rpfb=char(xor(ki(i+1,:)-48,tempr-48)+48);
if ind==i+1
result=[result,'和子密钥k',char(i+a1),'做异或运算得:',rpfb,sprintf('\n')];
end
rpfb=sbox(rpfb);
if ind==i+1
result=[result,'通过s盒处理得:',rpfb,sprintf('\n')];
end
rpfb=rebit(rpfb,pt);
if ind==i+1
result=[result,'从s盒输出后再通过P盒处理得:',rpfb,sprintf('\n')];
end
rpfb=char(xor(templ-48,rpfb-48)+48);
if ind==i+1 result=[result,'从p盒输出后与L',char(i+a1),'异或得R',char(i+a2+1),':',rpfb,sprintf('\n')];
result=[result,'第',char(i+a1+1),'轮','处理得[L',char(i+a2+1),' ','R',char(i+a2+1),']:',[lpfb,rpfb],sprintf('\n\n')];
end
end
pfb=[rpfb lpfb];
pfb=rebit(pfb,endrt);%末置换处理
ef=pfb;
result=[result,pf,'用密钥',pw,'进行des加密后得密文:',ef];
