图像通信系统设计与仿真—课程设计论文

1、 科f 课 程 设 计报 告report on curriculum design 题 目： 图像通信系统设计与仿真 学科部、系： 信息学科部电子系 目录目录i摘要ii第一章 图像通信系统的概述11.1 图像通信通信课程设计的任务与要求11.1.1 设计任务11.1.2 设计要求11.2 图像通信通信课程设计流程1第二章 图像的编码与解码22.1 哈弗曼编码的原理2 2.2 哈弗曼解码的原理2第三章 数字信号的调制与解调33.1 数字信号的调制33.1.1 调制的概念及分类33.1.2 调制的原理33.2 数字信号的解调43.2.1 解调的概念43.2.2 解调的原理4第四章 图形用户界面的设

2、计与制作6 4.1 图形用户界面gui的概述6 4.2 图形用户界面（gui）的设计与制作6第五章 总结7参考文献8附录1 图像通信在在matlab中算法实现9摘要图像通信是以图像通信的基本原理、基本技术理论为基础发展的新技术，具体内容包括图像与数字图像通信、图像分析与正交变换、图像压缩编码、静止图像编码、活动图像编码、图像通信中的信道编码与传输技术。本课程设计主要研究图像通信系统的设计和仿真，在matlab中对图像进行读取、压缩和哈弗曼编码，将数字化后的图像进行数字调制，已调后的信号可以传输和交换，到达接收端，经解调和解码还原成图像，从而完成了图像信号传送和接收。利用图形用户界面(gui)，

3、通过窗口、菜单、按键等方式来方便地对本次课程设计的图像通信进行操作，观察。本文通过对图像进行哈弗曼编码，对比ask,fsk,psk三种调制方式进行传输，利用相干解调方式，然后再进行哈弗曼解码，在图形用户界面上实现图像通信操作仿真。关键词：图像通信系统，传输接收，压缩，哈弗曼编码解码，调制，相干解调，图形用户界面 第一章 图像通信系统的概述1.1 图像通信通信课程设计的任务与要求1.1.1 设计任务1) 学习并掌握图像通信的基本原理，学会利用仿真的手段对于实用通信系统的基本理论，基本算法进行实际验证；2) 学习通信系统仿真软件matlab的基本使用方法，解决实际系统出现的问题；3) 通过matl

4、ab进行仿真，并制作图形用户界面(gui)加深对图像通信系统理论的理解；1.1.2 设计要求1) 编码：读入一幅图像，对图像进行哈弗曼编码；2) 调制：采用振幅键控(ask），相移键控（psk），频移键控（fsk）三种方式对编码后产生的数字信号进行调制；3) 解调：采用相干解调方式；4) 解码：哈弗曼解码输出原图像；1.2 图像通信通信课程设计流程 根据系统设计的任务要求，确定该系统的设计流程方框图，如图1.1所示。 图1.1 图像的设计流程方框图 第二章 图像的编码与解码2.1 哈弗曼编码的原理哈弗曼编码过程分为两步骤：第一，消减信源符号数量；第二，对每一个信源符号赋值。基于哈夫曼编码图像压

5、缩的基本原理是在对数据进行编码时可以对那些经常出现的数据指定较少的位数表示,而那些不常出现的数据指定较多的位数表示,每个数据的代码各不相同,这是一种典型的无损编码方式。这些代码都是二进制码,且码字长度是不均匀的、平均码率可以接近信息源熵值的一种编码。 哈夫曼编码是一种常用的压缩编码算法,采用变长码编码,属于无损压缩算法的一种。 由于进行的是无损压缩,所以要扫描图像的所有像素点,压缩过程分为四步:1 扫描统计像素出现的概率并按大小排列；2 建立最优二叉树;3 哈夫曼编码;4 保存编码。 2.2 哈弗曼解码的原理哈夫曼解码过程是文件编码过程的逆过程，由于哈夫曼编码可即时解码，因此只要得到一个码字,

6、则通过查找哈夫曼编码表得到相应的字符，映射过程是编码时映射的逆过程。因此，每从压缩文件中读出一个码字，就从通过查找哈夫曼编码表用字符替换相应的码字，当压缩文件中所有的码字被字符替换掉，也就宣告解压过程完成了。 第三章 数字信号的调制与解调 3.1 数字信号的调制3.1.1 调制的概念及分类用基带信号对载波（通常为余弦或正弦）波形的某些参数（如幅度、相位和频率）进行控制，使这些参数随基带信号的变化而变化。因此，数字信号的调制可分为振幅键控(ask），相移键控（psk），频移键控（fsk）。3.1.2 调制的原理1) m进制振幅键控(mask）基本原理m进制数字振幅调制信号的载波幅度有m种取值，在

7、每个符号时间间隔ts内发送m个幅度中的一种幅度的载波。表达式：a. g(t)为基带信号波形b. ts为符号时间间隔c. an为幅度值2) m制频移键控（mfsk）基本原理表达式： 式中 为载波角频率，共有m种取值。通常可选载波频率fn/2ts，n为正整数,此时m种发送信号相互正交。3) m进制相移键控（mpsk）基本原理表达式 ： 式中a. g(t)-信号包络波形，通常为矩形波，幅度为1b. ts-码元时间宽度c. c-载波角频率d. -第n个码元对应的相位，共有m种取值3.2 数字信号的解调3.2.1 解调的概念解调也叫检波，其作用就是从接收到的已调波中无失真地恢复出调制信号。调幅信号的解调

8、技术可大致分为两类：相关解调和非相关解调。 1) 相干解调（又称同步检波）：需要在接收端知道发射载波的频率和相位（恢复相干载波以用于与接收的已调信号相乘）。其中乘机检波器也叫鉴相器，是将输入的带通信号变为基带信号的下变频电路。2) 非相干解调（又称包络检波）：即不需要知道有关的信息，它的输出信号与输入信号实际包络成正比。3.2.2 解调的原理1. ask信号的解调 2ask相干解调方框图如图所示2. fsk信号的解调2fsk相干解调方框图如图所示3. psk信号的解调2psk相干解调方框图如图所示第四章 图形用户界面的设计与制作 图形用户界面或图形用户接口(graphical user int

9、erface，gui)是指采用图形方式显示的计算机操作环境用户接口。与早期计算机使用的命令行界面相比，图形界面对于用户来说更为简便易用。gui的广泛应用是当今计算机发展的重大成就之一，它极大地方便了非专业用户的使用人们从此不再需要死记硬背大量的命令，取而代之的是可用用通过窗口、菜单、按键等方式来方便地进行操作。而嵌入式gui具有下面几个方面的基本要求：轻型、占用资源少、高性能、高可靠性、便于移植、可配置等特点。 4.1 图形用户界面gui的概述图形用户界面gui是一种包含多种对象的图形窗口，并为gui开发提供一个方便高效的集成开发环境guide。guide主要是一个界面设计工具集，matlab

10、将所有gui支持的控件都集成在这个环境中。guide将设计好的gui保存在一个fig文件中，同时生成m文件框架。fig文件：该文件包括gui图形窗口及其所有fig文件包含序列化的图形窗口对象；m文件：该文件包括gui设计、控制函数以及定义为子函数的用户控件回调函数。guide可以根据gui设计过程直接自动生成m文件框架，优点如下：a. m文件已经包含一些必要的代码；b. 管理图形对象句柄并执行回调函数子程序；c. 提供管理全局数据的途径；d. 支持自动插入回调函数原型。4.2 图形用户界面（gui）的设计与制作1) gui创建包括界面设计和控件编程两部分，主要步骤如下：a. 通过设置guide

11、应用程序的选项来运行guide；b. 使用界面设计编辑器进行界面设计；c. 编写控件行为响应控制（即回调函数）代码。2) 图像通信系统的的图形用户界面gui 第五章 总 结 图像通信系统为目前广泛使用的技术，以图像通信理论为基础，反映了当代图像通信发展的新技术。对基本原理和技术的分析同时涉及模拟图像通信和数字图像通信，并以数字图像的分析为主。具体内容包括图像与数字图像通信、图像分析与正交变换、图像压缩编码、静止图像编码、活动图像编码、图像通信中的信道编码与传输技术、图像编码方法的新进展和数字图像通信系统。本次图像通信系统设计与仿真是对图像通信系统在matlab中的简单应用。通过在matlab中

12、对图像进行读取、压缩和哈弗曼编码，将数字化后的图像进行数字调制，已调后的信号可以传输和交换，到达接收端，经解调和解码还原成图像，从而完成了图像信号传送和接收。利用图形用户界面(gui)，通过窗口、菜单、按键等方式来方便地对本次课程设计的图像通信进行操作，观察。本文通过对图像进行哈弗曼编码，对比ask,fsk,psk三种调制方式进行传输，利用相干解调方式，然后再进行哈弗曼解码，在图形用户界面上实现图像通信操作仿真。让我们深入了解图像通信系统及其技术应用。参考文献1 娄莉图像通信原理与技术m. 北京：清华大学出版社，20102 章毓晋. 图像处理和分析教程m. 北京：人民邮电出版社，2009 3

13、夏良正. 数字图像处理（第2版）.南京：东南大学出版社，2005 4 姚敏. 数字图像处理.北京：机械工业出版社，2005 5 李朝晖. 数字图像处理及应用（含1cd）.北京：机械工业出版社，2005 6 何东健，耿楠，张义宽. 数字图像处理.西安：西安电子科技大学出版社，2004 7 贾永红. 数字图像处理.武汉：武汉大学出版社，2004 附录1 图像通信在在matlab中算法实现1、 哈弗曼编码在matlab中的算法1. function ratio,rowsb,colsb,dictionary,huffcode_i=huffencode(i)2. k = imread(cameraman

14、.tif);3. i = imresize(k,20,20);4. maxi = max(max(i); %求图像灰度级最大值5. h = zeros(1,maxi+1); %直方图矩阵初始化6. rows,cols = size(i);7. rowsb=tentotwo(rows); %调用十进制转换成二进制的程序8. colsb=tentotwo(cols);9. for i = 1:rows10. for j = 1:cols11. h(i(i,j)+1) =h(i(i,j)+1)+1;12. end13. end14. h = h/(rows*cols);15. index, hn =

15、 sortdezero(h); %调用sortdezero函数 16. indexcell = num2cell(index); %huffman树的下标变成元包数组17. lenindex= length(index);18. codeword = cell(1,maxi+1) ; %定义元包数组存放码字19. for i =1:lenindex-120. index1 = indexcell1; 21. index2 = indexcell2; for j =1:length(index1)22. codeword(index1(j) = 0,codewordindex1(j) ;23.

16、end24. for j =1:length(index2)25. codeword(index2(j) = 1,codewordindex2(j) ;26. end 27. indexcell = index1,index2,indexcell3:end;28. hn = sum(hn(1:2),hn(3:end);29. index, hn = sorth(hn);30. indexcell = indexcell(index);31. end32. lavg = 0;33. for i = 1:length(h)34. lavg = lavg+ h(i)* length( codewor

17、di);35. end36. ratio=8/lavg; %计算压缩比37. ii =i(:);38. huffcode_i = ; % 存储图像的huffman编码39. for i = 1:length(ii)40. huffcode_i = huffcode_i,codewordii(i)+1;41. end42. huffcode_i = logical(huffcode_i);43. dictionary=zeros(1,uint16(maxi+1)*16);%字典的传送44. a=;45. for i=1:double(maxi+1)46. num=length(codewordi

18、);47. y=num;48. if num=049. for k=1:450. b=mod(y,2);51. dictionary(16*(i-1)+(5-k)=b;52. y=floor(y/2);53. end54. a=codewordi;55. for j=1:num56. dictionary(16*(i-1)+4+j)=a(j);57. end58. end59. end哈弗曼编码在matlab中的算法如下：1. function di = huffdecode(rowsb,colsb,dictionary, huffcode_i) 2. rows=0; %图像列数3. for

19、i=1:length(rowsb)4. if rowsb(i)=05. rows=rows+2(length(rowsb)-i);6. end7. end8. cols=0; %图像行数9. for i=1:length(colsb)10. if colsb(i)=011. cols=cols+2(length(colsb)-i);12. end13. end 14. codeword = cell(1,256); % 字典解码接收,字典解码存放码字15. for i=1:25616. num=dictionary(16*(i-1)+1)*8+dictionary(16*(i-1)+2)*4+

20、dictionary(16*(i-1)+3)*2+dictionary(16*(i-1)+4);17. if num=018. codeword(i)=dictionary(16*(i-1)+5):(16*(i-1)+4+num);19. end20. if length(dictionary(16*i+1):end)=num33. if huffcode_i(k:k+num-1)=codewordi34. d(q)=i-1; %解码35. k=k+num;36. break;37. end38. end39. end40. end41. end42. di = reshape(d,rows,

21、cols);2、 调制在matlab中的算法如下：1. function eth et bn = digtalmod(fs,fc,an,m,rb,df,modtype)2. tr=fs/rb;3. if m=24. bn = bin2decxie(an,m);5. else 6. bn=an;7. end8. n0=length(bn);9. t=n0/rb;%n个码元的传输时间10. t=0:1/fs:t-1/fs;11. for i=1:length(bn)12. s(i-1)*tr+1:i*tr)=bn(i);13. end14. switch modtype15. case ask11

22、6. et = s.*cos(2*pi*fc*t); %ask17. case fsk218. et = cos(2*pi*(fc+s*df).*t);%fsk19. case psk320. et = cos(2*pi*fc*t+s*2*pi/m); 21. end22. eth=hilbert(et);解调在matlab中的算法如下：1. function bn0 = digitaldemod(eth,fs,fc,rb,m,modtype)2. t=length(eth)/fs;3. t=0:1/fs:t-1/fs;4. tr=fs/rb;5. eth=hilbert(eth);6. switch modtype7. case ask18. des = abs(eth);9. case fsk210. phi=unwrap(angle(eth);11. des=diff(phi);12. case psk313. eth=eth.*exp(-j*2*pi*fc*t