分享：用四十种语言分别写一个MD5算法 之20 MATLAB语言MD5算法

分享：用四十种语言分别写一个MD5算法 之20 MATLAB语言MD5算法

function digest = md5(message)

%% Helper Function Definitions

%This function converts a binary representation of a number or vector

%of numbers from a string to a row vector of 1's and 0's

function matrix = binStr2mat(binStr)

matrix = zeros(size(binStr));

for n = (1:numel(binStr))

matrix(n) = str2double(binStr(n));

end

end

%This acts as a "lower <= x <= upper" operator

function trueFalse = inRange(lowerBound,theValue,upperBound)

trueFalse = (lowerBound <= theValue) && (theValue <= upperBound);

end

%This converts a decimal number into its 32-bit binary representation,

%and I'm pretty sure it is little-endian

function binaryRep = to32BitBin(decimal)

binaryRep = binStr2mat(dec2bin(decimal,32));

end

%This converts a decimal number into its 64-bit binary representation,

%and I'm pretty sure it is little-endian

function binaryRep = to64BitBin(decimal)

binaryRep = binStr2mat(dec2bin(decimal,64));

end

%This adds multiple binary numbers together modulo 2^32. The sum of these

%numbers will rap around if the sum of two numbers is >= 2^32

function result = addBinary(varargin)

result = 0;

for l = (1:numel(varargin))

temp = num2str(varargin{l});

temp(temp == ' ') = [];

result = mod(result + bin2dec(temp),2^32);

end

result = to32BitBin(result);

end

%% MD5 Hash Algorithm

%Define constants

%r is the bit-shift amount for each round

r =[7,12,17,22,7,12,17,22,7,12,17,22,7,12,17,22,5,9,14,20,5,9,14,20,5,...

9,14,20,5,9,14,20,4,11,16,23,4,11,16,23,4,11,16,23,4,11,16,23,...

6,10,15,21,6,10,15,21,6,10,15,21,6,10,15,21];

%Use the binary representation of the radian values of sine represented

%as 32-bit integers

k = to32BitBin(floor(abs(sin(1:64)).*(2^32)));

%Initialize the hash variables, convert the hex representation of these

%variables to their 32-bit binary representation

h0 = to32BitBin(hex2dec('67452301'));

h1 = to32BitBin(hex2dec('EFCDAB89'));

h2 = to32BitBin(hex2dec('98BADCFE'));

h3 = to32BitBin(hex2dec('10325476'));

%Convert the ASCII values of the input string to a binary string matrix

message = dec2bin(message);

N = numel(message);

%Convert the message matrix to one long bit stream

message = reshape(transpose(message),1,N);

%Perform the pre-prossesing of appending a 1-bit to the message, then

%enough zeros to make the length of the message congruent to 448 mod

%512. Then append the 64-bit representation of the length of the

%original message string to the end.

message = [binStr2mat([message '1']) zeros(1,abs(mod(N+1,512) - 448)) to64BitBin(N)];

%Process 512-bit chunks of the pre-processed message string

for chunk = (1:(numel(message)/512))

%Stored the index of the first bit in the chunk

blockIndex = ((chunk-1)*512)+1;

%Pull out the 512-bit chunk from the message, then reshape it into

%32, 16-bit strings. Then transpose that string matrix to recover 16,

%32-bit words in the correct order. If we do this directly, by

%reshaping the 512-bit chunk to 16, 32-bit strings the original

%order of the bits will be scrambled.

w = transpose(reshape(message(blockIndex:blockIndex+511),32,16));

%Initialize the hashes for this round

a = h0;

b = h1;

c = h2;

d = h3;

%Process the hahes

%Note: The original MD5 algorithm pseudo-code was written for

%0-based arrays. Therefore, the calculation for "g" has been

%modified for MATLAB's 1-based arrays

for i = (1:64)

if inRange(1,i,16)

f = (b & c) | ((~b) & d);

g = i;

elseif inRange(17,i,32)

f = (d & b) | ((~d) & c);

g = mod( (5*(i-1)) + 1 ,16) + 1;

elseif inRange(33,i,48)

f = xor( b,xor(c,d) );

g = mod( (3*(i-1)) + 5 ,16) + 1;

elseif inRange(48,i,64)

f = xor( c, (b | (~d)) );

g = mod( (7*(i-1)),16 ) + 1;

end

temp = d;

d = c;

c = b;

%The circshift has a -r(i) because normally this function

%rotates right, not left.

b = addBinary(b,circshift( addBinary(a,f,k(i,1:32),w(g,1:32)),[0 -r(i)] ));

a = temp;

end

%Add this chunk's hash to all the previous chunks' hashes

h0 = addBinary(h0,a);

h1 = addBinary(h1,b);

h2 = addBinary(h2,c);

h3 = addBinary(h3,d);

end

%Concatenate the hashes together and convert them to a binary string

digest = num2str([h0 h1 h2 h3]);

%Remove whitespace from the digest, leaving 1's and 0's

digest(digest == ' ') = [];

%Convert the binary representation of the digest to a hexadecimal

digest = lower(bin2hex(digest));

end %md5

Sample:

>> md5('hello csdner')

版权所有：百万商业圈

未经许可不得转载，有任何疑问请与我本人联系 QQ 99923309 Mail：bwsyq@bwsyq.com
更多详情>>

开源：完全自主研发搜索引擎1.0源代码及说明，单机400万网页，任意50词以内的检索不超过 20毫秒

开源：基于百万商业圈.NET开发框架开发的并行带分词的采集器

百万商业圈 .NET 开发框架2.0及开发框架API说明书(BWFW)(含并行计算及中英文分词功能)

分享一点代码（小型C web开发框架），用C语言实现的一个WEB 文件上传（含全部源代码）一

天心天字辈ERP全部PDK源代码到了我手上的后果 - 超越天心之WEB天云

大家看看我的BS的甘特图排程做的如何? 无刷新Ajax甘特图 展示生产排程结果 演示

软件工程概述 - 企业架构 - IT企业做大做强之根本 - 之我见

实践检验得出的真理：asp.net 项目在 linux mono上编译需要修改的只有 3个地方

给大家漏一手本人亲自精心撰写的通用ajax框架,完全兼容 IE FireFox各个版本!(附完整源码及完整范例)

开发了一个中文分词服务器(C语言开发+词库+源代码),最大特色可以让javascript来调用!

用纯C语言写了一个HtmlParse(网页分析器)外带采集功能，大小只有200K(免费+开源+操作示意图)

分享：用九种语言分别写一个web server 之 C语言web server

分享：用九种语言分别写一个web server 之 Delphi语言web server

分享：用九种语言分别写一个web server 之 Fantom语言web server

分享：用九种语言分别写一个web server 之 Perl语言web server

分享：用九种语言分别写一个web server 之 Go语言web server

分享：用九种语言分别写一个web server 之 PicoLisp语言web server

分享：用九种语言分别写一个web server 之 PureBasic语言web server

分享：用九种语言分别写一个web server 之 Python语言web server

分享：用九种语言分别写一个web server 之 Tcl语言web server