tcm网格编码MATLAB程序文档格式.docx
- 文档编号:22275197
- 上传时间:2023-02-03
- 格式:DOCX
- 页数:15
- 大小:48.78KB
tcm网格编码MATLAB程序文档格式.docx
《tcm网格编码MATLAB程序文档格式.docx》由会员分享,可在线阅读,更多相关《tcm网格编码MATLAB程序文档格式.docx(15页珍藏版)》请在冰豆网上搜索。
tyu=length(z)/2;
s=zeros(1,3*tyu);
fori=1:
s(:
3*i-2)=int(i,2)
end
s(:
3*m+1)=0;
3*(m+1)+1)=0;
forj=1:
tyu
3*j-1)=z(2*j-1)
fork2=1:
3*k2)=z(2*k2)
uu=reshape(s,3,tyu);
kk=uu'
;
(注意:
对添零后卷积编码的处理),将编码处理后的信号进行调制,相位调制实现比较容易,对于通过信道后的解调,有两种实现途径:
①接收信号通过相关器后,将接收到的信号矢量映射到M个可能发送的信号矢量上去,并且选出对应于最大映射的矢量;
②亦可计算接收信号矢量的相位,并从M个可能发送的信号矢量中选出相位最接近的信号。
本实验就是通过第二个途径实现的。
其具体程序如下:
n1=gngauss(sgma);
n2=gngauss(sgma);
f(i)=bin2deci(kk(i,:
));
u=cos(2*pi*fc*t+2*pi*f(i)/8);
if((f(i)>
=3)&
(f(i)<
7))
R(i)=sqrt(Es)*cos(2*pi*f(i)/8)+n1;
H(i)=sqrt(Es)*sin(2*pi*f(i)/8)+n2;
T(i)=pi+atan(H(i)/R(i))
elseiff(i)<
3
T(i)=atan(H(i)/R(i))
else
T(i)=2*pi+atan(H(i)/R(i))
end;
cc=zeros(tyu,8);
dd=zeros(tyu,8);
forj=0:
7
cc(i,j+1)=2*pi*j/8-T(i)
dd=abs(cc)
ifdd(i,1)==min(dd(i,:
))
mm(i)=0;
elseifdd(i,2)==min(dd(i,:
mm(i)=1;
elseifdd(i,3)==min(dd(i,:
mm(i)=2;
elseifdd(i,4)==min(dd(i,:
mm(i)=3;
elseifdd(i,5)==min(dd(i,:
mm(i)=4;
elseifdd(i,6)==min(dd(i,:
mm(i)=5;
elseifdd(i,7)==min(dd(i,:
mm(i)=6;
elsedd(i,8)==min(dd(i,:
mm(i)=7;
out=reshape(mm,tyu,1)
对于解码则是编码的逆过程,在此将解调后的信号经反映射,变换为卷积码的码序列,应用viterbi实现解码:
k1=size(out,1);
输出行数
z=zeros(3,k1)
k1
m1=deci2bin(out(i,1),3)
z(:
i)=m1'
hh=reshape(z,1,3*k1)
ty=zeros(k1-2,2);
k1-2
ty(i,2)=hh(:
3*i-2);
hh(:
2*j-1)=[];
[ww,tt,ee]=VITERBI(g,k,hh)
ty(i,1)=ww(i);
ip=reshape(int,1,2*m);
op=reshape(ty,1,2*m);
当然在此基础上可进行误码率的分析,最简单的方法是将以上ip,op对应元素比较,若不相等,计数器加‘1’,可以使用以下程序如下:
numoferr=0;
forw=1:
2*m
ifip(w)~=op(w)
numoferr=numoferr+1;
对此程序的一点说明:
此程序可实现误码的估计,为了从真正意义上对其优点的认识,可仿真8QPSK的误码性能,可调用以下程序:
function[pb]=smld(snr_in_dB)
N=input(‘number’);
Eb=1;
snr=10^(snr_in_dB/10);
sgma=sqrt(1/(2*snr));
s000=[10];
s001=[sqrt
(2)/2sqrt
(2)/2];
s011=[01];
s010=[-sqrt
(2)/2sqrt
(2)/2];
s110=[-10];
s111=[-sqrt
(2)/2-sqrt
(2)/2];
s101=[0-1];
s100=[sqrt
(2)/2-sqrt
(2)/2];
N,
temp=rand;
if(temp<
0.125),
dsource1(i)=0;
dsource2(i)=0;
dsource3(i)=0;
elseif(temp<
0.25),
dsource3(i)=1
0.375),
dsource2(i)=1;
0.5),
dsource3(i)=1;
0.625),
dsource1(i)=1;
0.75),
0.875),
else
n
(1)=gngauss(sgma);
n
(2)=gngauss(sgma);
if((dsource1(i)==0)&
(dsource2(i)==0)&
dsource3(i)==0),
r=s000+n;
elseif((dsource1(i)==0)&
(dsource3(i)==1)),
r=s001+n;
elseif((dsource1(i)==0)&
(dsource2(i)==1)&
(dsource3(i)==0)),
r=s010+n;
(dsource3(i)==1)),
r=s011+n;
elseif((dsource1(i)==1)&
r=s100+n;
r=s101+n;
r=s110+n;
else((dsource1(i)==1)&
r=s111+n;
c000=dot(r,s000);
c001=dot(r,s001);
c010=dot(r,s010);
c011=dot(r,s011);
c100=dot(r,s100);
c101=dot(r,s101);
c110=dot(r,s110);
c111=dot(r,s111);
c_max=max([c000c001c010c011c100c101c110c111]);
if(c000==c_max),
decis1=0;
decis2=0;
decis3=0;
elseif(c001==c_max),
decis3=1;
elseif(c010==c_max),
decis2=1;
elseif(c011==c_max)
elseif(c100==c_max)
decis1=1;
elseif(c101==c_max)
elseif(c110==c_max)
if(decis1~=dsource1(i)),
if(decis2~=dsource2(i)),
if(decis3~=dsource3(i)),
pb=numoferr/(3*N);
对TCM的仿真调用程序:
function[p]=smld2(snr_in_dB)
N=input('
number'
);
Es=3;
0.5)
0.75)
int=zeros(1,2*N);
foro=1:
int(1,2*o-1)=dsource2(o);
forb=1:
int(1,2*b)=dsource1(b);
rrr=reshape(int,2,N);
eee=rrr'
;
m=size(eee,1);
n)=eee(n,1)
3*i-2)=eee(i,2)
3*N+1)=0;
3*(N+1)+1)=0;
3*j-1)=z(2*j-1);
3*k2)=z(2*k2);
n2=gngauss(sgma)
T(i)=pi+atan(H(i)/R(i));
T(i)=atan(H(i)/R(i));
T(i)=2*pi+atan(H(i)/R(i));
cc(i,j+1)=2*pi*j/8-T(i);
dd=abs(cc);
out=reshape(mm,tyu,1);
z=zeros(3,k1);
m1=deci2bin(out(i,1),3);
hh=reshape(z,1,3*k1);
2*j-1)=[];
[ww,tt,ee]=viterbi(g,k,hh);
ty(i,1)=ww(i)
ip=reshape(eee,1,2*N)
op=reshape(ty,1,2*N)
2*N,
p=numoferr/(2*N)
对误码率图的一点说明,其中‘0’代表不经过信道编码,直接将信号送入8QPSK调制器,解调后对其误码率的仿真,‘*’则代表使用TCM方式后对误码率的仿真,两者对比后与理论之间的差别所产生的原因有以下几点:
①仿真时所取的样点是有局限的,以上两图分别是样点数为200、420点的仿真图,它们的差别就说明了样点数的选择对系统性能的仿真是很重要的;
②解调后的判决途径是不一样的,仅用8QPSK方式是采用以上所提到的途径Ⅰ,而TCM方式是途径Ⅱ,虽说两者是等价的,但对于仿真系统还是稍有区别,不过从以上两图可看出TCM对性能的改善。
对于相同的误码率,其的确有所改善。
(注:
对1000点的仿真)
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- tcm 网格 编码 MATLAB 程序