哈夫曼編碼 MATLAB程序

ID:207134 · 發表于 2017-6-2 16:34

大家可以隨便看看。

MATLAB源程序如下:

clc
clear
p=[0.2 0.2 0.1 0.25 0.15 0.1]%p為概率分布，此函數功能是進行哈夫曼編碼
% 此處顯示詳細說明
% h為各個元素的麻子
% e為輸出的平均碼長
if length(find(p<0))~=0
error('概率不應該小于0！')
end
if abs(sum(p)-1)>10e-10
error('概率之和大于1，請檢查輸入！')
end
n=length(p);
q=sort(p);
m=zeros(n-1,n);
for i=1:n-1
[q,e]=sort(q)
m(i,:)=[e(1:n-i+1),zeros(1,i-1)]%由數組l 構建一個矩陣，該矩陣表明概率合并時的順序，用于后面的編碼
q=[q(1)+q(2),q(3:n),1];
end
for i=1:n-1
c(i,:)=blanks(n*n);%c 矩陣用于進行huffman 編碼
end
c(n-1,n)='1';
c(n-1,2*n)='0';
for i=2:n-1
c(n-i,1:n-1)=c(n-i+1,n*(find(m(n-i+1,:)==1))-(n-2):n*(find(m(n-i+1,:)==1))) ;%矩陣c 的第n-i 的第一個元素的n-1 的字符賦值為對應于a 矩陣中第n-i+1 行中值為1 的位置在c 矩陣中的編碼值
c(n-i,1*n)='1';
c(n-i,n+1:2*n-1)=c(n-i,1:n-1); %矩陣c 的第n-i 的第二個元素的n-1 的字符與第n-i 行的第一個元素的前n-1 個符號相同，因為其根節點相同
c(n-i,2*n)='0';
for j=1:i-1
c(n-i,(j+1)*n+1:(j+2)*n)=c(n-i+1,n*(find(m(n-i+1,:)==j+1)-1)+1:n*find(m(n-i+1,:)==j+1)) ;%矩陣c 中第n-i 行第j+1 列的值等于對應于a 矩陣中第n-i+1 行中值為j+1 的前面一個元素的位置在c 矩陣中的編碼值
end
end
c
for i=1:n
hm(i,1:n)=c(1,n*(find(m(1,:)==i)-1)+1:find(m(1,:)==i)*n);%用hm表示最后的huffman 編碼
len(i)=length(find(abs(hm(i,:))~=32));%計算每一個編碼的長度
end
hm
len
H=sum(-p.*log2(p)) ;
q=sort(p)
pingjunmachang=sum(q.*len) %計算平均碼長
xiaolv=H/pingjunmachang*100

復制代碼

clc
clear
fid=fopen( 'C:\Users\yichao\Desktop\新建文本文檔.txt');%打開 txt 文件
[zimu]=fscanf(fid, '%c'); %讀取二進制文件的數據，并將數據存入矩陣
sort_str1=sort(zimu)%按照字符的 ASCII 值對字符串排序
j=1;
str2(1)=sort_str1(1);
for i=1:length(sort_str1)-1 %計算出字符串的種類
if strcmp(sort_str1(i),sort_str1(i+1))~=1
j=j+1;
end
str2(j)=sort_str1(i+1);
end
str2
for i = 1:length(str2) %length 函數獲取字符串長度
str_num =strfind(sort_str1,str2(i)); %strfind(S1,S2): 尋找 S2 是否匹配 S1 ，并返回 S2 的位置
count1(i) = length(str_num);
end
count1
count=count1(1:end);
p=count./sum(count);
sum(-p.*log2(p)) %計算信源熵
if length(find(p<0))~=0
error('概率不應該小于0！')
end
if abs(sum(p)-1)>10e-10
error('概率之和大于1，請檢查輸入！')
end
n=length(p);
q=sort(p);
m=zeros(n-1,n);
for i=1:n-1
[q,e]=sort(q);
m(i,:)=[e(1:n-i+1),zeros(1,i-1)]; %由數組l 構建一個矩陣，該矩陣表明概率合并時的順序，用于后面的編碼
q=[q(1)+q(2),q(3:n),1];
end
for i=1:n-1
c(i,:)=blanks(n*n);%c 矩陣用于進行huffman 編碼
end
c(n-1,n)='1';
c(n-1,2*n)='0';
for i=2:n-1
c(n-i,1:n-1)=c(n-i+1,n*(find(m(n-i+1,:)==1))-(n-2):n*(find(m(n-i+1,:)==1))) ;%矩陣c 的第n-i 的第一個元素的n-1 的字符賦值為對應于a 矩陣中第n-i+1 行中值為1 的位置在c 矩陣中的編碼值
c(n-i,1*n)='1';
c(n-i,n+1:2*n-1)=c(n-i,1:n-1); %矩陣c 的第n-i 的第二個元素的n-1 的字符與第n-i 行的第一個元素的前n-1 個符號相同，因為其根節點相同
c(n-i,2*n)='0';
for j=1:i-1
c(n-i,(j+1)*n+1:(j+2)*n)=c(n-i+1,n*(find(m(n-i+1,:)==j+1)-1)+1:n*find(m(n-i+1,:)==j+1)) ;%矩陣c 中第n-i 行第j+1 列的值等于對應于a 矩陣中第n-i+1 行中值為j+1 的前面一個元素的位置在c 矩陣中的編碼值
end
end
for i=1:n
hm(i,1:n)=c(1,n*(find(m(1,:)==i)-1)+1:find(m(1,:)==i)*n);%用hm表示最后的huffman 編碼
len(i)=length(find(abs(hm(i,:))~=32)); %計算每一個編碼的長度
end
hm
len
H=sum(-p.*log2(p)) ;
q=sort(p)
pingjunmachang=sum(q.*len) %計算平均碼長
xiaolv=H/pingjunmachang*100

復制代碼

所有資料51hei提供下載:

信息論.7z (2.86 KB, 下載次數: 18)

ID:143174 · 發表于 2018-6-24 10:42

謝謝了

帳號		自動登錄	找回密碼
密碼			立即注冊

久久久久久久999_99精品久久精品一区二区爱城_成人欧美一区二区三区在线播放_国产精品日本一区二区不卡视频_国产午夜视频_欧美精品在线观看免费

哈夫曼編碼 MATLAB程序