【導讀】交互式程序設計的高科技計算環(huán)境。它將數(shù)值分析、矩陣計算、科學數(shù)據(jù)可視化。C、Fortran）的編輯模式，代表了當今國際科學計算軟件的先進水平[1]。隨著信息技術的發(fā)展，圖像信息被廣泛應用于多媒體通信和計算機系統(tǒng)中，但是圖像數(shù)據(jù)的一個顯著特點就是信息量大。具有龐大的數(shù)據(jù)量，如果不經(jīng)過壓?？s,不僅超出了計算機的存儲和處理能力,而且在現(xiàn)有的通信信道的傳輸速率下，些圖像數(shù)據(jù),必須對其進行壓縮,因此有必要對圖像壓縮編碼進行研究。余、統(tǒng)計冗余、結構冗余、視覺冗余以及知識冗余等。隨著科學技術的發(fā)展,圖像壓縮編碼技術越來越引起人們的關注。謂的圖像壓縮編碼技術就是對要處理的圖像數(shù)據(jù)按一定的規(guī)則進行變換和組合,從而達到以盡可能少的數(shù)據(jù)流來表示盡可能多的數(shù)據(jù)信息。用MATLAB進行了實驗仿真并分析結果得出結論。分量保留低頻分量的做法對于圖像壓縮具有絕對的優(yōu)勢。采用DCT實現(xiàn)對圖像編解碼設計；實現(xiàn)DCT變換和反變換以及JPEG標準量化矩陣的量化、反量化。

　　

【正文】 +1。 else %如果遇到的不是 0，那么 Xpq_coded=（，）的第一個元素存儲之前 0《 信息處理系統(tǒng)綜合訓練 》報告 第 17 頁 的個數(shù)，第二個 元素存儲此時遇到的這個元素 Xpq_coded(press_index,1)=count。 Xpq_coded(press_index,2)=Xpq_estimate(table(index,1),table(index,2))。 count=0。 %過去這個元素之后 0 的個數(shù)就要重新計 press_index=press_index+1。 end if ~count==0 Xpq_coded(press_index,1)=0。 Xpq_coded(press_index,2)=0。 end end %主要思想： %跟壓縮過稱恰好相反 %讀取壓縮文件 %Z 字形游程解碼 %Jpeg 標準量化表反量化 %分塊 DCT 變換 %存儲解壓后的文件 clear all。clc。 global quotiety。 uiload。 h=waitbar(0,39。正在解壓，請稍候 ?? 39。)。 Xpq_decoded=Xpq_coded。 clear Xpq_coded。 [M,N]=size(Xpq_decoded)。 waitbar(1/4,h)。 %游程解碼得到 Xpq_deestimate for index1=1:1:M for index2=1:1:N Xpq_deestimate(8*(index11)+1:8*index1,8*(index21)+1:8*index2)=Run_length_d《 信息處理系統(tǒng)綜合訓練 》報告 第 18 頁 ecoding_function(Xpq_decoded{index1,index2})。 end end waitbar(2/4,h)。 %反量化得到幅頻值 Xpq_resume=blkproc(Xpq_deestimate,[8 8],39。deintegral_function39。)。 waitbar(3/4,h)。 %DCT 逆變換得到 xmn xmn_resume=blkproc(Xpq_resume,[8 8],39。IDCT2_function39。,8)。 xmn_resume=uint8(xmn_resume)。 waitbar(4/4,h)。 close(h)。 figure。imshow(xmn_resume)。 title(39?；謴偷膱D像 39。)。 function xmn_resume=IDCT2_function(Xpq_resume,M) T=dctmtx(M)。 %生成變換矩陣 xmn_resume=T39。*Xpq_resume*T。 % 用變化矩陣的乘法代替效率極低的 for 求和循環(huán) function Xpq_resume=deintegral_function(Xpq_deestimate) integral_table=[16,11,10,16,24,40,51,61。 12,12,14,19,26,58,60,55。 14,13,16,24,40,57,69,56。 14,17,22,29,51,87,80,62。 18,22,37,56,68,109,103,77。 24,35,55,64,81,104,113,92。 49,64,78,87,103,121,120,101。 72,92,95,98,112,100,103,99]。 global quotiety。 Xpq_resume=Xpq_deestimate.*integral_table*quotiety。% 反量化恢復的過程要考慮到量化系數(shù) quotiety 的影響 《 信息處理系統(tǒng)綜合訓練 》報告 第 19 頁 %與 Run_length_coding_function 完全相逆的過程 解碼，此處不再贅述 function Xpq_deestimate=Run_length_decoding_function(Xpq_decoded) table=[1,1。1,2。2,1。3,1。2,2。1,3。1,4。2,3。 3,2。4,1。5,1。4,2。3,3。2,4。1,5。1,6。 2,5。3,4。4,3。5,2。6,1。7,1。6,2。5,3。 4,4。3,5。2,6。1,7。1,8。2,7。3,6。4,5。 5,4。6,3。7,2。8,1。8,2。7,3。6,4。5,5。 4,6。3,7。2,8。3,8。4,7。5,6。6,5。7,4。 8,3。8,4。7,5。6,6。5,7。4,8。5,8。6,7。 7,6。8,5。8,6。7,7。6,8。7,8。8,7。8,8]。 [cols,tmp]=size(Xpq_decoded)。 Xpq_deestimate(1:8,1:8)=0。 depress_index=1。 for index=1:1:cols if Xpq_decoded(index,2)==0 elseif Xpq_decoded(index,1)==0 Xpq_deestimate(table(depress_index,1),table(depress_index,2))=Xpq_decoded(index,2)。 depress_index=depress_index+1。 else depress_index=depress_index+Xpq_decoded(index,1)。 Xpq_deestimate(table(depress_index,1),table(depress_index,2))=Xpq_decoded(index,2)。 depress_index=depress_index+1。 end en