+++ /dev/null
-function Bin_Sbox =Convert_Sbox_Bin(Sbox)\r
-Bin_Sbox=zeros(256,8);\r
-for it=1:256\r
- Bin_Sbox(it,:)=byte2bit(Sbox(it));\r
- %double(dec2bin(Sbox(it),8))-48;\r
- %byte2bit\r
-end\r
-\r
-\r
+++ /dev/null
-function [RK,Sbox,Pbox]=Dynamickeygenerationnew(Secretkey,counter,strh,h)\r
-x=java.security.MessageDigest.getInstance(strh);\r
-x.update(bitxor(Secretkey,counter));\r
-DK=double(typecast(x.digest,'uint8'));\r
-M=sqrt(h*8);\r
-%%%%%%%%%%%% Creation of P-box (blocks permutation)%%%%%%%%%%%%%%%%%%%%%%%%\r
-% construction of P-box\r
-\r
-%Pbox=PermutationCat(Keyp,M);\r
-\r
-Qp=8*h;\r
-Ap=1:Qp;\r
-Keyp=reshape((DK(1:1:16)),1,[]);\r
-Keyp=mod(Keyp,Qp)+1;\r
-\r
-Pbox=supboxvector2(Ap,Qp,Keyp)\r
-\r
-\r
-\r
-% construction of S-box\r
-Q=2^8;\r
-A=1:Q;\r
-KS=reshape(DK(17:32),1,[]);\r
-Sbox=supboxvector2(A,Q,KS)-1;\r
-\r
-%%%%%%%%%%%%%%%% Construction of round keys%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \r
-r=4;\r
-X=reshape(DK(33:48),1,[]);\r
-X=double(typecast(uint8(X),'uint32'))\r
-P=reshape(DK(49:end),1,[]);\r
-P=double(typecast(uint8(P),'uint32'))\r
-N=2^32;\r
-m=ceil(h/4)\r
-RK= roundKeyiterations(X,P,N,r,Sbox,m);\r
-\r
-\r
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\r
-% %%%%%%%%%%%%% creation of IV (block of h bits) can be employ for CBC, CFB, OFB, CTR mode %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \r
-% s_rk=rc4key(DK(17:1:32));\r
-% IV=reshape(prga( s_rk,h^2),h,h);\r
-% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\r
-% %%%%%%%%%%%%%%%%%%%%creation of the diffusion matrix%%%%%%%%%%%%%%%%%%%%%%%\r
-% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\r
-% sc2=rc4key(DK(end-16:end));\r
-% outd=prga( sc2,(h^2)/4);\r
-% A=reshape(outd,h/2,h/2);\r
-% In=eye(h/2,h/2);\r
-% if (kindd==0)\r
-% G=[A A+In; A-In A ];\r
-% else\r
-% A=mod(A,2);\r
-% G=[A bitxor(A,In) ; bitxor(A,In) A]; \r
-% end
\ No newline at end of file
+++ /dev/null
-function y = bit2byte(x)\r
-% y = bit2byte(x)\r
-% Transform a number of bits into a number of bytes.\r
-\r
-x = double(x);\r
-\r
-x_n = numel(x);\r
-if rem(x_n,8)\r
- x = [x zeros(1,8-rem(x_n,8))];\r
-end\r
-y = zeros(1,ceil(x_n/8));\r
-y_index = 1;\r
-v = 2.^(0:7);\r
-for i=1:8:x_n\r
- y(y_index) = sum(x(i+7:-1:i).*v);\r
- y_index = y_index + 1;\r
-end\r
-y = uint8(y);\r
+++ /dev/null
-function y = byte2bit(x)\r
-% y = byte2bit(x)\r
-% Transform a number of bytes into a number of bits.\r
-\r
-x_n=numel(x);\r
-y=zeros(1,x_n*8);\r
-y_index=1;\r
-for i=1:x_n\r
- for j=8:-1:1\r
-% for j=1:8\r
- y(y_index+8-j) = bitget(x(i),j);\r
- end\r
- y_index = y_index + 8;\r
-end\r
-y = uint8(y);\r
+++ /dev/null
-function dimg =decryptionprocess(img,h,RK,Inv_SubBytes,Inv_Pbox,r)\r
-%%%%%%%%%%%%%% convert image to a vector seq\r
-%seq=double(reshape(img,1,[]));\r
-seq_bin=double(byte2bit(uint8(img))) ;\r
-decr_seq_binary=zeros(1,length(seq_bin));\r
-%%%%%%%%%%%%%%%%%%%% calculate the number of blocks%%%%%%%%%%%%%%%%%%%%%%%%\r
-nb= floor(length(seq_bin)/(8*h));\r
-RK_bin=double(reshape(byte2bit(RK),r,[]));\r
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\r
-%%%%%%%%%%% Convert S-box to binary representation to avoid repeat the\r
-%%%%%%%%%%% byte2bit conversion\r
-Bin_Sbox =double(Convert_Sbox_Bin(Inv_SubBytes));\r
-%%%%%%%%%%%%%%%%%%%% Applied the confusion process (r rounds)\r
-for it=1:nb\r
- temp=seq_bin((it-1)*8*h+1:it*h*8);\r
- for itr=r:-1:1 \r
- temp= temp(Inv_Pbox);\r
- temp=reshape(Bin_Sbox (double(bit2byte(temp))+1,:)',1,[]);\r
- % double(byte2bit(Inv_SubBytes(double(bit2byte(temp))+1)));\r
- %\r
-% size(temp)\r
- %byte2bit(Inv_SubBytes(double(bit2byte(temp))+1));\r
- %temp= Inv_SubBytes(temp+1);\r
- temp=bitxor(temp,RK_bin(itr,1:h*8)); \r
- end\r
- decr_seq_binary((it-1)*8*h+1:it*h*8)=temp;\r
-end\r
-\r
-% Convert to Byte representation\r
-\r
-%valbin=reshape(valbin,Tb/8,8);\r
-dimg=uint8(bit2byte(decr_seq_binary));
\ No newline at end of file
+++ /dev/null
-function dimg =decryptionprocess2_CBC(img,h,RK,Inv_SubBytes,Inv_Pbox,r,IV)\r
-%%%%%%%%%%%%%% convert image to a vector seq\r
-%seq=double(reshape(img,1,[]));\r
-seq_bin=double(byte2bit(uint8(img))) ;\r
-decr_seq_binary=zeros(1,length(seq_bin));\r
-%%%%%%%%%%%%%%%%%%%% calculate the number of blocks%%%%%%%%%%%%%%%%%%%%%%%%\r
-nb= floor(length(seq_bin)/(8*h));\r
-RKB=byte2bit(RK);\r
-RKB=RKB(1:r*h*8);\r
-%RK_bin=double(reshape(byte2bit(RK),r,[]));\r
-RK_bin=double(reshape(RKB,r,[]));\r
-%RK_bin=double(reshape(byte2bit(RK),r,[]));\r
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\r
-%%%%%%%%%%% Convert S-box to binary representation to avoid repeat the\r
-%%%%%%%%%%% byte2bit conversion\r
-Bin_Sbox =double(Convert_Sbox_Bin(Inv_SubBytes));\r
-\r
-iv=[];\r
-for i=1:h\r
- iv=[iv byte2bit(IV(i))];\r
-end\r
-IV=double(iv);\r
-\r
-\r
-%%%%%%%%%%%%%%%%%%%% Applied the confusion process (r rounds)\r
-for it=1:nb\r
- temp=seq_bin((it-1)*8*h+1:it*h*8);\r
- temp_IV=temp;\r
- for itr=r:-1:1 \r
- temp= temp(Inv_Pbox);\r
- temp=reshape(Bin_Sbox (double(bit2byte(temp))+1,:)',1,[]);\r
- temp=bitxor(temp,RK_bin(itr,1:h*8)); \r
- end\r
- temp=bitxor(temp,IV);\r
- decr_seq_binary((it-1)*8*h+1:it*h*8)=temp;\r
- IV= temp_IV;\r
-end\r
-\r
-% Convert to Byte representation\r
-\r
-%valbin=reshape(valbin,Tb/8,8);\r
-dimg=uint8(bit2byte(decr_seq_binary));
\ No newline at end of file
+++ /dev/null
-function eimg =encryptionprocess(img,h,RK,SubBytes,Pbox,r)\r
-%%%%%%%%%%%%%% convert image to a vector seq\r
-%seq=double(reshape(img,1,[]));\r
-seq_bin=double(byte2bit(uint8(img))) ;\r
-encr_seq_binary=zeros(1,length(seq_bin));\r
-%%%%%%%%%%%%%%%%%%%% calculate the number of blocks%%%%%%%%%%%%%%%%%%%%%%%%\r
-nb= floor(length(seq_bin)/(8*h));\r
-RK_bin=double(reshape(byte2bit(RK),r,[]));\r
-Bin_Sbox =double((Convert_Sbox_Bin(SubBytes)));\r
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\r
-%%%%%%%%%%%%%%%%%%%% Applied the confusion process then the bits permutation for r rounds\r
-for it=1:nb\r
- temp=seq_bin((it-1)*8*h+1:it*8*h);\r
- for itr=1:r\r
- temp=bitxor(temp,RK_bin(itr,1:h*8)); \r
- temp= (reshape(Bin_Sbox (double(bit2byte(temp))+1,:)',1,[]));\r
- temp= temp(Pbox);\r
- end\r
- encr_seq_binary((it-1)*8*h+1:it*h*8)=temp;\r
-end\r
-\r
-% Convert to Byte representation\r
-\r
-%valbin=reshape(valbin,Tb/8,8);\r
-eimg=(bit2byte( encr_seq_binary));
\ No newline at end of file
+++ /dev/null
-function eimg =encryptionprocess2_CBC(img,h,RK,SubBytes,Pbox,r,IV)\r
-%%%%%%%%%%%%%% convert image to a vector seq\r
-%seq=double(reshape(img,1,[]));\r
-seq_bin=double(byte2bit(uint8(img))) ;\r
-encr_seq_binary=zeros(1,length(seq_bin));\r
-%%%%%%%%%%%%%%%%%%%% calculate the number of blocks%%%%%%%%%%%%%%%%%%%%%%%%\r
-nb= floor(length(seq_bin)/(8*h));\r
-RKB=byte2bit(RK);\r
-RKB=RKB(1:r*h*8);\r
-%RK_bin=double(reshape(byte2bit(RK),r,[]));\r
-RK_bin=double(reshape(RKB,r,[]));\r
-Bin_Sbox =double((Convert_Sbox_Bin(SubBytes)));\r
-%iv=zeros(1,8*h);\r
-%iv(end-h+1:end)=IV;\r
-%IV=iv;\r
-iv=[];\r
-for i=1:h\r
- iv=[iv byte2bit(IV(i))];\r
-end\r
-IV=double(iv);\r
-\r
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\r
-%%%%%%%%%%%%%%%%%%%% Applied the confusion process then the bits permutation for r rounds\r
-for it=1:nb\r
- temp=seq_bin((it-1)*8*h+1:it*8*h);\r
- temp=bitxor(temp,IV);\r
- for itr=1:r\r
- temp=bitxor(temp,RK_bin(itr,1:h*8)); \r
- temp= (reshape(Bin_Sbox (double(bit2byte(temp))+1,:)',1,[]));\r
- temp= temp(Pbox);\r
- end\r
- encr_seq_binary((it-1)*8*h+1:it*h*8)=temp;\r
- IV=temp;\r
-end\r
-\r
-% Convert to Byte representation\r
-\r
-%valbin=reshape(valbin,Tb/8,8);\r
-eimg=(bit2byte( encr_seq_binary));
\ No newline at end of file
+++ /dev/null
-function inv_perm_tabs= inverse_tables(tab)\r
-le=length(tab);\r
-inv_perm_tabs=zeros(1,le);\r
-for ir = 1 :le\r
- inv_perm_tabs(tab(ir)+1) = ir-1 ;\r
-end\r
-\r
+++ /dev/null
-function inv_perm_tabs= inverse_tables2(tab)\r
-tab2=reshape(tab,1,[]);\r
-le=length(tab2);\r
-inv_perm_tabs=zeros(1,le);\r
-for ir = 1 :le\r
- inv_perm_tabs(tab2(ir)) = ir ;\r
-end\r
-%inv_perm_tabs=reshape(inv_perm_tabs,size(tab,1),size(tab,2));\r
+++ /dev/null
-%clear all\r
-%close all\r
-%clc\r
-\r
-\r
-\r
-%%%%%%%%%%%%%%%% Input data%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\r
-[filename, pathname] = uigetfile( {'*.bmp';}, 'Pick a file');\r
-name=[ pathname filename];\r
-data=imread(name);\r
-%data=imresize(data,[64 64])\r
-\r
-h=2;% number of bytes in a block (Tb=h*8)\r
-r=4;% number of rounds\r
-sk=16 % number of byes in a block cipher\r
-strh='SHA-512';\r
-% mode=0==> ECB, \r
-% mode 1==> CBC\r
-\r
-mode=1\r
-\r
-\r
-\r
-% calculate the length of data\r
-SI=size(data);\r
-%len=SI(1)*SI(2)*SI(3);\r
-%lenBits=len*8;\r
-%data_line=reshape(double(data),1,len); \r
-\r
-figure()\r
-imshow(uint8(data));\r
-\r
-%%%%%%%%%%%%%%%%% Key generation%%%%%%%%%%%%%%%%%%%%%%%%\r
-Secretkey=floor(rand(1,sk)*256);\r
-counter=floor(rand(1,sk)*256);\r
-%%%%%%%%%%%%%%%%%%%%%%%%%Creation of cipher layer's%%%%%%%%%%%%%%%%%%%%\r
-[RK,Sbox,Pbox]=Dynamickeygenerationnew(Secretkey,counter,strh,h);\r
- %%%%%%%%%%%%%%%%%%%%%%%%%%%%% Change in data block%%%%%%%%%%%%%%%%%%%%%%%% \r
-\r
-\r
-\r
-\r
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Encryption Process%%%%%%%%%%%%%%%%%%%%%\r
-tic\r
-if mode==0\r
-eimg =encryptionprocess2(data,h,RK,Sbox,Pbox,r);\r
-elseif mode==1\r
-IV=double(floor(rand(1,h)*256)); \r
-eimg =encryptionprocess2_CBC(data,h,RK,Sbox,Pbox,r,IV); \r
-end\r
-\r
-toc\r
-eimg2=reshape(eimg,size(data));\r
-\r
-figure()\r
-imshow(uint8(eimg2));\r
-\r
-%%%%%%%%%%%%%%%%%%%%% Decryption Side%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\r
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\r
-%%%%%%%%%%%%%%%%%%%% calculate the inverse substitution and diffusion\r
-%%%%%%%%%%%%%%%%%%%% layer%%%%%%%%%%%%%%%%\r
- Inv_Sbox= inverse_tables(Sbox);\r
- Inv_Pbox= inverse_tables2(Pbox);\r
- \r
- \r
- if mode==0\r
- dimg =decryptionprocess2(eimg,h,RK,Inv_Sbox,Inv_Pbox,r);\r
-elseif mode==1\r
-\r
- dimg =decryptionprocess2_CBC(eimg,h,RK,Inv_Sbox,Inv_Pbox,r,IV); \r
- end\r
-\r
-% tic\r
-% dimg =decryptionprocess2(eimg,h,RK,Inv_Sbox,Inv_Pbox,r);\r
-% toc\r
- dimg=reshape(dimg,size(data));\r
- figure()\r
- imshow(uint8(dimg));\r
- \r
+++ /dev/null
-function RK=roundKeyiterations (X,P, N,r,Sbox,m)\r
-% DIfferent size of input block\r
-%m=1;\r
-RK=zeros(r,4*m);\r
-for itr=1:r\r
- for itm=1:m\r
- y=skewtententierperturbed(X(itr),P(itr),N)\r
- %Y(i)=double(typecast(y,'uint8'))\r
- op=double(typecast(uint32(y),'uint8'))\r
- opd(1)=bitxor(bitxor(op(1),op(2)),op(3));\r
- opd(2)=bitxor(bitxor(op(1),op(2)),op(4));\r
- opd(3)=bitxor(bitxor(op(1),op(4)),op(3));\r
- opd(4)=bitxor(bitxor(op(4),op(2)),op(3));\r
- opd=Sbox(opd+1);\r
- RK(itr,(itm-1)*4+1:itm*4)=opd(:);\r
- X(itr)=y;\r
- end\r
-end
\ No newline at end of file
+++ /dev/null
-function y=skewtententierperturbed(x,p,Q)\r
-% x : initial condition\r
-% p: control parameter\r
-% N:precision\r
-if x<=p\r
- y=ceil(Q*(x/p));\r
- else \r
- y=floor(Q*((Q-x)/(Q-p)))+1;\r
-end\r
+++ /dev/null
-function x=supboxvector2(x,Q,p)\r
-% r=1 for each parameters\r
-y=zeros(1,Q);\r
-for i=1:length(p)\r
-%for w=1:r\r
-dinf=find (x<=p(1,i));\r
-dsup=find (x>p(1,i));\r
-y(dinf)=ceil((Q*x(dinf))./p(1,i));\r
-y(dsup)=floor((Q*(Q-x(dsup)))./(Q-p(1,i)))+1; \r
-x=y;\r
-%end\r
-end\r
