OFDM调制解调仿真.docx
- 文档编号:9826210
- 上传时间:2023-02-06
- 格式:DOCX
- 页数:11
- 大小:226.98KB
OFDM调制解调仿真.docx
《OFDM调制解调仿真.docx》由会员分享,可在线阅读,更多相关《OFDM调制解调仿真.docx(11页珍藏版)》请在冰豆网上搜索。
OFDM调制解调仿真
实验四OFDM调制解调仿真
一、实验目的
(1)了解OFDM调制解调的原理。
(2)学会用星座图分析系统性能。
二、实验容
(1)编写MATLAB程序,实现OFDM系统调制解调。
(2)绘制各步骤图形并分析系统特性。
三、实验代码
1、主代码如下:
clearall;
closeall;
N=input('请输入码元数');
SNR=input('请输入信噪比');
xx=randint(1,4*N);%原序列
figure
(1),stem(xx,'.k');%原序列图形
title('原序列');
B=0;
form=1:
4:
4*N
A=xx(m)*8+xx(m+1)*4+xx(m+2)*2+xx(m+3);
B=B+1;
ee(B)=A;
end
figure
(2),stem(ee,'.b');
title('化为0~15的码元');
yy=star(ee,N);
figure(3),plot(yy,'.r');%映射后的星座图
title('映射后的星座图');
ff=ifft(yy,N);%傅里叶反变换
N1=floor(N*1/4);
N3=floor(N*3/4);
N5=floor(N*5/4);
figure(4),stem(ff,'.m');
title('傅里叶反变换后');
forj=1:
N1%加循环前缀
ss(j)=ff(N3+j);
end
forj=1:
N%变成长度为N*5/4的序列
ss(N1+j)=ff(j);
end
figure(5),stem(ss,'.k');%画出图形
title('加N/4循环前后缀');
%ss=wgn(1,N5,0,10,'dBW','complex');
ss=awgn(ss,SNR);%加入噪声
figure(6),stem(ss,'.m');%加入噪声后的图形
title('加入噪声后');
zz=fft(ss((N1+1):
N5),N);%傅里叶变换
figure(7),plot(zz,'.b');%画图
title('傅里叶变换后');
rr=istar(zz,N);%星座图纠错
figure(8),plot(rr,'.r');%画图
title('纠错后的星座图');
dd=decode(rr,N);%解码
figure(9),stem(dd,'.m');%画图
title('星座图纠错并解码后');
%bb=d2b(dd,N);
bb=d2bb(dd,N);%转化为0/1比特流
figure(10),stem(bb,'.b');
title('转化为0/1比特流后');
2、调用的函数:
(1)星座图映射
%星座图映射
functionyy=star(xx,N)
B=[-3-3*i,-3-i,-1-3*i,-1-i,-3+3*i,-3+i,-1+3*i,-1+i,3-3*i,3-i,1-3*i,1-i,3+3*i,3+i,1+3*i,1+i];
forj=1:
N
yy(j)=B(xx(j)+1);
end
(2)星座图逆映射
%星座图逆映射
functionrr=istar(zz,N)
forj=1:
N
if(mod((floor(real(zz(j)))),2)==0)
zz1(j)=ceil(real(zz(j)));
elsezz1(j)=floor(real(zz(j)));
end
if(mod((floor(imag(zz(j)))),2)==0)
zz1(j)=zz1(j)+ceil(imag(zz(j)))*i;
elsezz1(j)=zz1(j)+floor(imag(zz(j)))*i;
end
rr(j)=zz1(j);
end
(3)十进制转二进制
functionbb=d2b(dd,N)
%十进制转二进制
forj=1:
N*4
bb(j)=;
end
forj=1:
4:
N*4
bb1=dec2bin(dd(floor(j/4)+1),4);
fork=1:
4
bb(4*(j-1)+k)=bb1(k);
end
end
(4)十进制转化为01比特流
%十进制转化为01比特流
functionbb=d2bb(dd,N)
forj=1:
N*4
bb(j)=1;
end
j=1;
while(j<=N*4)
N1=ceil(j/4)
a4=mod(dd(N1),2);
dd(N1)=floor(dd(N1)/2);
a3=mod(dd(N1),2);
dd(N1)=floor(dd(N1)/2);
a2=mod(dd(N1),2);
dd(N1)=floor(dd(N1)/2);
a1=mod(dd(N1),2);
bb(j)=a1;
j=j+1;
bb(j)=a2;
j=j+1;
bb(j)=a3;
j=j+1;
bb(j)=a4;
j=j+1;
end
(5)星座图逆映射
functionyy=decode(rr,N)
forj=1:
N
switch(rr(j))%星座图逆映射
case-3-3*i
yy(j)=0;
case-3-i
yy(j)=1;
case-1-3*i
yy(j)=2;
case-1-i
yy(j)=3;
case-3+3*i
yy(j)=4;
case-3+i
yy(j)=5;
case-1+3*i
yy(j)=6;
case-1+i
yy(j)=7;
case3-3*i
yy(j)=8;
case3-i
yy(j)=9;
case1-3*i
yy(j)=10;
case1-i
yy(j)=11;
case3+3*i
yy(j)=12;
case3+i
yy(j)=13;
case1+3*i
yy(j)=14;
case1+i
yy(j)=15;
otherwisebreak;
end
end
四、
五、
实验结果
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- OFDM 调制 解调 仿真