数字通信技术实验报告.docVIP

下载本文档

1
0
约8.98千字
约 15页
2019-07-26 发布于江西
举报
版权申诉

数字通信技术实验报告.doc

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

实验一： clear all [x,fs,bits]=wavread(myheart.wav); z=sign(x); Max=max(abs(x)); x1=abs(x/Max); Q=2048*x1; Y=zeros(length(x1),8); for m=1:length(x) if Q(m)128 Q(m)2048 Y(m,2)=1; end if (Q(m)32 Q(m)128) || (Q(m)512 Q(m)2048) Y(m,3)=1; end if (Q(m)16 Q(m)31) || (Q(m)64 Q(m)128) || (Q(m)256 Q(m)512) || (Q(m)1024 Q(m)2048) Y(m,4)=1; end if z(m)0 Y(m,1)=1; elseif z(m)0 Y(m,1)=0; end end N=zeros(1,length(x1)); for m=1:length(x1) N(m)=Y(m,2)*4+Y(m,3)*2+Y(m,4)+1; end a=[0,16,32,64,128,256,512,1024]; b=[1,1,2,4,8,16,32,64]; for m=1:length(x1) q=ceil(Q(m)-a(N(m))/b(N(m))); if q==0 Y(m,(5:8))=[0,0,0,0]; else k=num2str(dec2bin(q-1,4)); Y(m,5)=str2num(k(1)); Y(m,6)=str2num(k(2)); Y(m,7)=str2num(k(3)); Y(m,8)=str2num(k(4)); end End 实验结果： Max =0.3906 a = Columns 1 through 7 0 16 32 64 128 256 512 Column 8 1024 b = 1 1 2 4 8 16 32 64 bits =16 fs =22050 k =111011011 m =168873 q =476 实验二： RGB=imread(8.jpg); I=rgb2gray(RGB); J=dct2(I); imshow(log(abs(J)),[]),colormap(jet(64)),colorbar J(abs(J)10)=0; K=idct2(J); figure,imshow(I) figure,imshow(K,[0 255]) 原图像经过离散余弦变换后原图像的灰度图像经反变换后的图像： 2.输入DCT进行JPEG图像压缩的代码如下： RGB=imread(8.jpg); I=rgb2gray(RGB); I=im2double(I); T=dctmtx(8); B=blkproc(I,[8 8],P1*x*P2,T,T); Mask=[ 1 1 1 1 0 0 0 0 1 1 1 0 0 0 0 0 1 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]; B2=blkproc(B,[8 8], P1*x,Mask); I2= blkproc(B2,[8,8],P1*x*P2,T,T); Subplot(1,2,1); Imshow(I);title(原图像); Subplot(1,2,2); Imshow(I2);title(压缩图像); 3.游程编码： image1=imread(lena.jpg); imshow(image1); image2=image1(:); image2length=length(image2); for i=1:1:image2length if image2(i)=127 image2(i)=255; else image2(i)=0; end end image3=reshape(image2,512,512);