閫氫俊鍘熺悊matlab璇剧▼璁捐�℃姤鍛�-matlab璇剧▼璁捐�℃姤鍛� - 鐧惧害鏂囧簱

三、源程序

clear clc

a=20; %a表示元素个数，可以修改码元为任意个数如1024

r=unidrnd(2,1,a-1); %随机产生离散均匀二值分布（1，2），共a-1个

rn=r-1; %每个元素减1变为（0，1）二值随机数列 rn=[1,rn]; %跟书上一致，从1开始随机产生离散均匀二值分布（1，2），共a个

ori=rn; %将原序列保存起来，便于后面解码后比较 ori

l=1; for k=1:a-3

if abs(rn(k))==1

m=rn(k); %m记录可能存在的连续4个0前面非零码元的符号

end if rn(k)==0

if rn(k+1)==0 if rn(k+2)==0

if rn(k+3)==0 %检测是否连续4个码元都是0 rn(k+3)=2*m; %用2表示书上的V end end end end end

rn ; %插入V(2)后 checked l=1; %l=1表示两个V之间有偶数个非零码，0是偶数，所以初始为1，l=-1表示两个V之间有奇数个非零码 for s=1:a

if abs(rn(s))==2 for d=s+1:a if abs(rn(d))==1 %V之间的非零码只有+1，-1遇到一个一使l变一次号，表示偶数或奇数 l=-l; end

if abs(rn(d))==2 if l~=1

l=1; %不需要插入B时若检测到V也应把计数器清零

5

%最初版本的升级处1 else %检测到下一个V时，若为偶数个，则插入B

rn(d-3)=3*(-(rn(d-4))/(abs(rn(d-4)))); %3代表B,插入不带符号的B

rn(d:a)=-rn(d:a);

end %V后的符号再交替 end end end end

rn %到此处完成了插入不带符号的B %for s=1:a

% if abs(rn(s))==3 %找到B

%rn(s)=rn(s)*(-(rn(s-1))/(abs(rn(s-1)))); %B（3）与前一个位置的带符号的归一值相乘再取反，实现符号B(3)的极性与前一非符号的相反

%B后面第三个就是接下来的V,从它开始非零码正负号交替变化 %end %end

HDB=rn; %给B(3)添加了符号，并且实现了V后的符号再交替

HDB

%以上便实现了HDB3的编码

%下面进行解码 for k=1:a

if abs(rn(k))==2

rn(k-3)=0;rn(k)=0; %每个V（2）都变成0，V前面第三个有可能是B(3)有可能是0也恢复为0 end

if abs(rn(k))==1 rn(k)=1; end end

rn %解码后的恢复序列

rn-ori %解码与原码比较全为0则解码正确

6

s=menu('通信原理','2ASK','2PSK','2FSK','2DPSK') switch s

case 1,scolor='2ASK'; n=8;N=100;K=4; a=randint(1,n); bita=[];sl=[];

bitRate=1e3;fc=1e3;%载频1KHZ t=linspace(0,1/bitRate,N); for i=1:length(a) if a(i)==0

bit1=zeros(1,N); else

bit1=ones(1,N); end

bita=[bita,bit1]; c=sin(2*pi*t*fc); sl=[sl c]; end

figure(1);

subplot(K,1,1);plot(bita,'LineWidth',1.5),title('基带信号'),grid on;axis([0,N*length(a),-2.5,2.5]);

tz=bita*6.*sl;

subplot(K,1,2);plot(tz,'LineWidth',1.5);title('ASK调制后信号');grid on;

signal=awgn(tz,80,'measured');

subplot(K,1,3);plot(signal,'LineWidth',1.5),grid on;title('信号+噪声')

Fs=3e3;

[b,a]=ellip(4,0.1,40,[999.9,1000.1]*2/Fs);%设计IIR带通滤波器，阶数为4，通带纹波0.1，阻带衰减40DB

sf=filter(b,a,signal);%信号通过该滤波器 figure(2); K1=4;

subplot(K1,1,1);plot(sf,'LineWidth',1.5),grid on;title('BPF')

signal2=abs(sf); %乘同频同相sin

subplot(K1,1,2);plot(signal2,'LineWidth',1.5),grid on;title('全波整流器');

Fs=3e3;%抽样频率400HZ

7

[b,a]=ellip(4,0.1,40,[50]*2/Fs);%设计IIR低通滤波器 sf1=filter(b,a,signal2);%信号通过该滤波器，输出信号sf

subplot(K1,1,3);plot(sf1,'LineWidth',1.5),grid on;title('LPF');

sf2=[];

LL=fc/bitRate*N; i=LL/2; bitb=[];

while (i<=length(sf1)) %判决 sf2=[sf2,sf1(i)>=0.001]; i=i+LL; end

for i=1:length(sf2) if sf2(i)==0

bit1=zeros(1,N); else

bit1=ones(1,N); end

bitb=[bitb,bit1]; end

figure(1);

subplot(K,1,4); plot(bitb,'LineWidth',1.5),grid on;title('解调后信号'); axis([0,N*length(sf2),-2.5,2.5]);

case 2,scolor='2PSK';

l=linspace(0,pi,50);% 数据初始化 t=linspace(0,9*pi,450); b=1:1:9; out=1:1:450; f=1:1:450; g=1:1:450;

w1=2 %正弦波f1的频率,可以根据自己想要的频率在此改写 %正弦波f2的频率,可以根据自己想要的频率在此改写 f1=sin(w1*l); figure(1);

f2=sin(w1*l+pi); figure(1);

subplot(2,1,1),plot(l,f1),axis([0 pi -1.2 1.2]),xlabel('t'),ylabel('f1');%画出f1信号波形 subplot(2,1,2),plot(l,f2),axis([0 pi -1.2 1.2]),xlabel('t'),ylabel('f2');%画出f2信号波形 a=[0 1 0 0 0 1 1 0 1]

8