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)
