Huffman编码程序.docVIP

Huffman编码程序 Huffman编码的Matlab实现--用于单导联ECG数据的压缩和解压缩 收藏 dataProcess.m ----主程序 norm2huff.m ----编码 huff2norm.m ----解码 （注意上面两个函数文件的末尾附有子函数） Lead1.mat ----单导联数据 Original_6000_123.mat ----原始数据（12导联） PS：除了打断点以外，Matlab的调试手段同VC6.0，即 F5快速调试至断点，F10单步运行，F11:Step in %%%%%%%%%%%%%%%%%%%%%%%%%dataProcess.m view plaincopy to clipboardprint? clear all, close all, clc; load Lead1.mat; data = Lead1; data = data + 150; % 预处理数据，将负转正 data = data(:); % 预处理，将列向量变行向量 % compress data fprintf(Compresing data ... ) [zipped,info] = norm2huff(data); fprintf(Done!\n) % decompress data fprintf(Decompressing data ... ) unzipped = huff2norm(zipped,info); fprintf(Done!\n) % test it isOK = isequal(data(:),unzipped(:)) whos data zipped unzipped clear all, close all, clc; load Lead1.mat; data = Lead1; data = data + 150; % 预处理数据，将负转正 data = data(:); % 预处理，将列向量变行向量 % compress data fprintf(Compresing data ... ) [zipped,info] = norm2huff(data); fprintf(Done!\n) % decompress data fprintf(Decompressing data ... ) unzipped = huff2norm(zipped,info); fprintf(Done!\n) % test it isOK = isequal(data(:),unzipped(:)) whos data zipped unzipped %%%%%%%%%%%%%%%%%%%%%%%%%norm2huff.m view plaincopy to clipboardprint? function [zipped,info] = norm2huff(vector) % NORM2HUFF Huffman codification (encoder)编码器 % Detailed explanation goes here % vector as a row 行向量 vector = vector(:); % frequency 【①求各采样值出现的频率】 f = frequency(vector); % simbols presents in the vector are simbols = find(f~=0); % first value is 1 not 0!!!【②这些采样值中共有多少个不重复的非零值】 f = f(simbols); %【③各采样值的频率】 % sort using the frequency [f,sortindex] = sort(f); % f=f(sortindex) simbols = simbols(sortindex); % 【④频率值由小到大排序，为Huffman编码作准备】 len = length(simbols); % 求向量的长度 simbols_index = num2cell(1:len); % Convert numeric array to cell array codeword_tmp = cell(len,1); % Construct cell array: len * 1;数组初始值为空 【⑤codewor