tmp=floor((abs(x(i))-st)/step);%确定x落在某段的第tmp级内 t=dec2bin(tmp,4)-48;%将tmp转化为四位二进制码 out(5:8)=t(1:4);%段内码为t end
a(i*8-7:i*8)=out(1:8); end
max=512;
/SK调制部分
fs=2000; %采样频率 dt=1/fs; f1=20;
f2=120; %两个信号的频率 g1=a;
g2=~a; %信号反转，和g1反向
g11=(ones(1,200))'*g1; %抽样 --首先生成1行2000列的全1矩阵,然后竖过来.然后与g1进行矩阵乘法
g1a=g11(:)'; %将无数行,无数列的一个矩阵g11,转置成为1行N列 g21=(ones(1,200))'*g2; g2a=g21(:)'; t=0:dt:51.2-dt; t1=length(t); subplot(2,1,2); plot(t,g1a);
title('PCM编码后的信号(局部)') ylabel('幅度') grid on;
axis([49.8 50.3 -0.2 1.2]); fsk1=g1a.*cos(2*pi*f1.*t); fsk2=g2a.*cos(2*pi*f2.*t);
fsk=fsk1+fsk2; %产生的信号 figure(2)
subplot(2,1,1); plot(t,fsk);
title('FSK信号(局部)') ylabel('幅度')
axis([49.8 50.3 -1.2 1.2]);
sn=awgn(fsk,10); %通过高斯信道 figure(2)
subplot(2,1,2);
plot(t,sn); %噪声波形
title('通过高斯以后的FSK信号(局部)') ylabel('幅度')
axis([49.8 50.3 -1 1.2]); %FSK解调
11
b1=fir1(101,[10/800 20/800]);
b2=fir1(101,[90/800 110/800]); %设置带通参数 H1=filter(b1,1,sn);
H2=filter(b2,1,sn); %经过带通滤波器后的信号 figure(3)
subplot(2,1,1); plot(t,H1);
title('经过带通滤波器f1后的波形(局部)') axis([49.8 50.3 -1 1.2]); ylabel('幅度') xlabel('t')
subplot(2,1,2); plot(t,H2);
title('经过带通滤波器f2后的波形(局部)') axis([49.8 50.3 -1 1.2]); ylabel('幅度') xlabel('t') sw1=H1.*H1;
sw2=H2.*H2; %经过相乘器 figure(4) subplot(211); plot(t,sw1);
title('经过相乘器h1后的波形') axis([49.8 50.3 -1 1.2]); ylabel('幅度') subplot(212); plot(t,sw2);
title('经过相乘器h2后的波形') ylabel('幅度')
axis([49.8 50.3 -1 1.2]); xlabel('t')
bn=fir1(101,15/800); %经过低通滤波器 figure(5)
st1=filter(bn,1,sw1); st2=filter(bn,1,sw2); subplot(211); plot(t,st1);
title('经过低通滤波器sw1后的波形') axis([49.8 50.3 -1 1.2]); ylabel('幅度') subplot(212); plot(t,st2);
title('经过低通滤波器sw2后的波形') ylabel('幅度')
12 axis([49.8 50.3 -1 1.2]); xlabel('t') %判决
for i=1:length(t)
if(st1(i)>=st2(i)) st(i)=1; else st(i)=0; end end
figure(6) subplot(211); plot(t,st);
axis([0 10 -1.5 1.5]);
title('经过抽样判决器后的波形') axis([49.8 50.3 -1 1.2]); grid on;
ylabel('幅度')
%对接收到的信号进行反PCM编码 %%抽取其中的PCM编码序列 for i=1:512
pcm(i)=st(i*200); end
%循环将PCM编码变成量化电平数 for i=1:64;
%初始化计算参数 dlm=0; dp=0;
%将2-4位段落码换算成起始电平 if pcm(i*8-6)==1; dlm=dlm+4; end
if pcm(i*8-5)==1; dlm=dlm+2; end
if pcm(i*8-4)==1; dlm=dlm+1; end
%得到由段落码换算来的起始电平 dp=2^(dlm+4);
%先加上最小量化电平的一半 dp=dp+2^(dlm)/2;
%将5-8位段内码结合段落码一起转换成电平数 %第5位
if pcm(i*8-3)==1;
13
dp=dp+2^(dlm+3); end %第6位
if pcm(i*8-2)==1;
dp=dp+2^(dlm+2); end %第7位
if pcm(i*8-1)==1;
dp=dp+2^(dlm+1); end %第8位
if pcm(i*8)==1;
dp=dp+2^(dlm); end
%最后判断第1位码,决定正负 if pcm(i*8-7)==0; dp=-dp; end
jiedp(i)=dp; end
%将量化电平数变成电压 for i=1:64;
dianya(i)=jiedp(i)/2048*6; end N=64;
fs=2800; %抽样频率 t=(0:N-1)/fs;%采样时间s f=100;%模拟信号的频率Hz
figure(6);subplot(2,1,2);plot(t,dianya);grid on
axis([0 1/f*2 -5 5]);title('PCM解码后的模拟信号')
五． 课程设计实验结果
图1为输入模拟信号，PCM编码后的波形
14
图2为曾加高斯噪声的信号
图3为经过带通滤波器的波形
图4为相乘后的波形
15