hill.m

hill密码加密解密hill密码加密解密 很好用 是信息与编码的实验报告对很多同学都会有帮助的

function hill()
%输入密钥
disp('输入密钥(矩阵)的维数');
n=input('');
disp('输入密钥(矩阵,按行输入)');
key=zeros(n,n);
for j=1:n
    for k=1:n
    key(j,k)=input('');
    end
end
d=det(key);%求矩阵的行列式
if d==0%判断矩阵是否可逆
    error('密钥矩阵不可逆，无法实现Hill密码');
end
i=1;%欧几里德算法
r0=26;
r1=d;
w=zeros(1,26);
while mod(r0,r1)~=0
    w(i)=fix(r0/r1);
    tmp=r0;
    r0=r1;
    r1=mod(tmp,r1);
    i=i+1;
end
i=i-1;r0=1;r1=1;
while i>0
    tmp=r0;
    r0=r1;
    r1=tmp+w(i)*r1;
    i=i-1;
end
%r1为d的逆
if r1==0
    error('密钥矩阵不可逆，无法实现Hill密码');
end
d=r1;
%输入明文
message=input('输入明文 \n','s');
m=size(message);
m=m(2);
if mod(m,n)~=0
    error('输入错误，明文长度应为矩阵维数的倍数');
end
for i=1:m
    if message(i)>='a' && message(i)<='z'
        message(i)=message(i)-97;
    elseif message(i)>='A' && message(i)<='Z'
        message(i)=message(i)-65;
    else
        error('输入错误，应该输入字母');
    end
end
%加密
i=1;
while i<m+1
    A(i:i+n-1)=message(i:i+n-1)*key;
    i=i+n;
end
for i=1:m
    if A(i)>26
        A(i)=mod(A(i),26);
    end
    A(i)=A(i)+97;
end
str=char(A);
fprintf('密文为%s',str);
%解密
k=1;
b=zeros(n-1,n-1);
for i=1:n                     %求逆矩阵
    for j=1:n
        for p=1:n
            for q=1:n
                if p~=i && q~=j
                    if p>i
                        k1=p-1;
                    else
                        k1=p;
                    end
                    if q>j
                        k2=q-1;
                    else
                        k2=q;
                    end
                    b(k1,k2)=key(p,q);
                end
            end
        end
         tmp=(-1)^(i+j)*det(b);
         tmp=mod(tmp,26);
         if tmp<0
             tmp=tmp+26;
         end
         key2(i,j)=tmp;
    end
end
key2=d*key2';
for i=1:n
    for j=1:n
        key(i,j)=mod(key(i,j),26);
        if key(i,j)<0
            key(i,j)=key(i,j)+26;
        end
    end
end
for i=1:m
    A(i)=A(i)-97;
end
i=1;
while i<m+1
    B(i:i+n-1)=A(i:i+n-1)*key2;
    i=i+n;
end
for i=1:m
    if B(i)>26
        B(i)=mod(B(i),26);
    end
    B(i)=B(i)+97;
end
str2=char(B);
fprintf('\n对密文解密后明文为%s\n',str2)