论文1

武汉理工大学《通信系统课群综合训练与设计》课程设计说明

%抽样判决 & 逆码变换部分- 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