weisuijijiaozhi.m

这是一个通信中交织器的程序 自己用过了还可以使用 通过将图像中的集中错误分散化来实现的

clear all;
clc;
K=input('K=');
svar=0.6;
x=round(svar*rand(1,K));
y=x;
disp('原始输出序列为：');
disp(x);
figure;
subplot(2,2,1);
stem(x);title('原始输出序列');xlabel('个数');ylabel('值');

%信息序列没经过交织器直接进入信道
nvar=0.165;
noise=nvar*randn(1,K);
z12=x+noise;
uvar=0.3;
for t=round(0.4*K):round(0.6*K)
    wp(t)=u(t);
end
ganrao=uvar*wp;
ganrao(1:(round(0.4*K)-1))=0;
ganrao((round(0.6*K)+1):K)=0;
z=z12+ganrao;

z2=round(z);
disp('经过信道后信息序列为：');
z3=z2(1:K);
disp(z3);
subplot(2,2,2);
stem(z3,'r');title('经过信道后信息序列为');xlabel('个数');ylabel('值');

%计算没经过交织器的误码率
j=0;
for i=1:K
    if(z3(i)~=x(i))
    j=j+1; 
    end
end
Pe1=j/K;
st=sprintf('没经过交织器的误码率为：%.20f',Pe1);
disp(st);

%序列经过伪随机交织器
%(1)产生伪随机序列
N=1;hk=1;fg=1;
while(N<=K)
ga=1+(K-1)*rand(1,1);
xuanze=round(ga);
if(x(xuanze)~=pi)
FG(fg)=xuanze;
jiaozhixulie(hk)=x(xuanze);
x(xuanze)=pi;
hk=hk+1;
fg=fg+1;
N=N+1;
end
end

%(2)进入信道
pattern22=jiaozhixulie+noise;
jiaozhixulie=pattern22+ganrao;

%(3)恢复原始信息序列
huifu=zeros(1,K);
for kh=1:hk-1
huifu(FG(kh))=jiaozhixulie(kh);
end
huifu=round(huifu);
disp('恢复后的信息序列为：');
disp(huifu);
subplot(2,2,3);
stem(huifu,'g');title('经过信道后交织信息序列');xlabel('个数');ylabel('值');

%计算经过交织器的误码率
ji=0;
for ij=1:K
    if(huifu(ij)~=y(ij))
    ji=ji+1; 
    end
end
Pe2=ji/K;
st2=sprintf('经过交织器的误码率为：%.20f',Pe2);
disp(st2);