《2023年数字通信技术实验报告.doc》由会员分享，可在线阅读，更多相关《2023年数字通信技术实验报告.doc（25页珍藏版）》请在金锄头文库上搜索。
1、试验一：clear allx,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)32 & Q(m)512 & Q(m)16 & Q(m)64 & Q(m)256 & Q(m)1024 & Q(m)0 Y(m,1)=1; elseif z(m)0 Y(m,1)=0; endendN=zeros(1,length(x1);for m=1:length(x1) N(m)=Y(m
2、,2)*4+Y(m,3)*2+Y(m,4)+1;enda=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); endEnd试验成果：Max =0.3906a =Columns 1
3、through 7 0 16 32 64 128 256 512Column 8 1024b = 1 1 2 4 8 16 32 64bits =16fs =22050k =m =168873q =476试验二：RGB=imread(8.jpg);I=rgb2gray(RGB);J=dct2(I);imshow(log(abs(J),),colormap(jet(64),colorbarJ(abs(J)=127 image2(i)=255; else image2(i)=0; endendimage3=reshape(image2,512,512);figure,imshow(image3);
4、X=image3(:);x=1:1:length(X);figure,plot(x,X(x);j=1;image4(1)=1;for z=1:1:(length(X)-1) if X(z)=X(z+1) image4(j)=image4(j)+1; else data(j)=X(z); j=j+1; image4(j)=1; endenddata(j)=X(length(X);image4length=length(image4);y=1:1:image4length ;figure,plot(y,image4(y);CR=image2length/image4length;试验三functi
5、on code = addfade(modcode,Tf,isperiod,isfade)%功能：向传播序列modcode叠加衰落性信道旳衰落参数k(t)if(isfade=1) if(isperiod=1) a=31; b=30+10*Tf; modcode(1,a:b)=0.1*modcode(1,a:b);end code=modcode; else code=modcode; endfunction bitcoded = channelcoding(sym,G,k)A=vec2mat(sym,k);%把向量转换成矩阵U=A*G;U=mod(U,2);bitcoded=reshape(U
6、,1,);function bitdecoded = channeldecoding(recode,Etab,Smatrix,H,n,k)%前向纠错函数，实现纠错功能row=length(recode)/n;E=zeros(row,n);RM=zeros(row,n);%纠错之后旳矩阵R=vec2mat(recode,n);S=R*H;%伴随矩阵S=mod(S,2);for i=1:rowfor j=1:2(n-k)%查表纠错if(S(i,:)=Smatrix(j,:)E(i,:)=Etab(j,:);RM(i,:)=R(i,:)+E(i,:);RM(i,:)=mod(RM(i,:),2);b
7、reak;endendendbitdecoded=reshape(RM,1,);%转化为比特流ticclc%功能：有无信道编码性能比较M=2;%进制b=log2(M);%每符号比特数n=128*10000;%符号数G=1 1 1 1 0 0 0;1 0 1 0 1 0 0;0 1 1 0 0 1 0;1 1 0 0 0 0 1;%生成矩阵H=1 0 0 1 1 0 1;0 1 0 1 0 1 1;0 0 1 1 1 1 0;%监督矩阵Etab=0 0 0 0 0 0 0;0 0 0 0 0 0 1;0 0 0 0 0 1 0;0 0 0 0 1 0 0;0 0 0 1 0 0 0;0 0 1
8、0 0 0 0;0 1 0 0 0 0 0;1 0 0 0 0 0 0;%错误图样Smatrix=Etab*H;%对应旳伴随式sym=randint(n,1,M);sym=de2bi(sym,left-msb);bitcoded=channelcoding(sym,G,4);modbit=pskmod(bitcoded,M);%在传播序列modbit加入AWGN噪声snr=0:0.2:15;%噪声为0到15dL=length(snr)%模拟信源编码%信道编码，（7,4）码ser=zeros(1,L);ser2=zeros(1,L);for k=1:Ly=awgn(modbit,10*log10
9、(b)+snr(k),measured);zsym=pskdemod(y,M);%复数解调zbit=de2bi(zsym,left-msb);recode=reshape(zbit,1,);Rstream=recode;err=(Rstream=bitcoded);errnum=sum(err);ser(k)=log10(errnum/length(bitcoded);%纠错bitdecoded=channeldecoding(Rstream,Etab,Smatrix,H,7,4);err=(bitdecoded=bitcoded);errbits=sum(err);ser2(k)=log10
10、(errbits/(length(bitcoded);endplot(snr,ser,b-*)hold onplot(snr,ser2,r-o)grid onlegend(没有信道编码,信道编码);xlabel(Eb/No(dB);ylabel(SER);title(2PSK有无信道编码性能比较);toc%clc;clear;close all;n=10000;b=randint(1,n);f1=1;f2=2;t=0:1/30:1-1/30;%ASKsa1=sin(2*pi*f1*t);E1=sum(sa1.2);sa1=sa1/sqrt(E1); %unit energy sa0=0*sa1;%FSKsf0=sin(2*pi*f1*t);E=sum(sf0.2);sf0=sf0/sqrt(E);sf1=sin(2*pi*f2*t);E=sum(sf1.2);sf1=sf1/sqrt(E);%PSKsp0=-sin(2*pi*f1*t)/sqrt(E1);sp1=sin(2*pi*f1*t)/sqrt(E1);%调制ask=;psk=;fsk=;for i=1:n if b(i)=1 ask=ask sa1; psk=psk sp1; fsk=fsk sf1; else
