武汉理工大学《通信系统课群综合训练与设计》课程设计说明
%抽样判决 & 逆码变换部分- jt=[];
for i = 1 : length(g) if dtt((2*i-1)*15)>0
jt1( (i-1)*30 + 1 : i*30) = 1; else
jt1( (i-1)*30 + 1 : i*30) = 0; end end
subplot(312); plot (jt1 );
axis( [ 0 , 3 *( length(g) + 1), -2.5,2.5] ); grid on;
title('抽样判决波形');
t=0:dt:endtime;%时间范围设定 %抽样
for i=1:length(g) jt2(i)=jt1(30*i); end
subplot(313) stairs(jt2);
axis([0 , ( length(g) + 1), -2.5,2.5]); grid on;
title('数字波形');
%循环码译码过程
%*******************************************************************% n=7;k=4;
p=cyclpoly(n,k) ;
jt3=decode(jt2,n,k,'cyclic',p);
%HDB3码的编码过程
%*******************************************************************%
38
武汉理工大学《通信系统课群综合训练与设计》课程设计说明
k=1; ihdb3=hdb3;
for i=1:endtime/dt*8 if hdb3(i)~=0 if hdb3(i)*k>0 ihdb3(i-3:i)=0; else
ihdb3(i)=1; end k=hdb3(i); end end figure(6) stairs(ihdb3);
axis([0 15 -1 2]);grid on; title('IHDB3'); %PCM译码部分
%*******************************************************************%
for i=1:length(jt3)/8 if (jt3(8*(i-1)+1)>0)
jp1(i)=jt3(8*(i-1)+1)*2^7+jt3(8*(i-1)+2)*2^6+jt3(8*(i-1)+3)*2^5 +jt3(8*(i-1)+4)*2^4+jt3(8*(i-1)+5)*2^3+jt3(8*(i-1)+6)*2^2+jt3(8*(i-1)+7)*2^1+jt3(8*i)-128; else
jp1(i)=(jt3(8*(i-1)+1)*2^7+jt3(8*(i-1)+2)*2^6+jt3(8*(i-1)+3)*2^5+jt3(8*(i-1)+4)*2^4+jt3(8*(i-1)+5)*2^3+jt3(8*(i-1)+6)*2^2+jt3(8*(i-1)+7)*2^1+jt3(8*i))*(-1); end end
figure(7); subplot(211); stem(xc); grid on;
39
武汉理工大学《通信系统课群综合训练与设计》课程设计说明
title('PCM译码信号');
%低通滤波部分
[N,Wn] = buttord( 2*pi*50, 2*pi*150,3,25,'s'); %临界频率采用角频率表示 [b,a]=butter(N,Wn,'s');
[bz,az]=impinvar(b,a,300); %映射为数字的 dtt = filter(bz,az,xc); subplot(212); plot(dtt); grid on;
title('解调后的模拟信号');
%*******************************************************************%%*******************************************************************%
%************************计算失真度**********************% q=abs(dtt)./abs(x); x=x.*q+2.5;
e=var((x-dtt)./x);fprintf('信号失真度为%9.5f\\n',e)
40
武汉理工大学《通信系统课群综合训练与设计》课程设计说明
本科生课程设计成绩评定表
姓 名 专业、班级 课程设计题目： 课程设计答辩或质疑记录： 性 别 成绩评定依据： 最终评定成绩（以优、良、中、及格、不及格评定） 指导教师签字：
年 月 日
41