matlab这是ofdm的时域和频域的子载波图形.docVIP

matｌaｂ这是ＯFDM的时域和频域的子载波图形 ％　ofdm sｐectrum plot cleａr aｌｌ; ｃlose all; clc; Num_Sc = 12;　 % １3 ｓubcａrrｉeｒｓ Ts = 1; 　 　 % 1ｓ Ｆ_ｓｐacｅ = 1/Ts； F　= -Ｆ_sｐace*Num_Ｓｃ／2-4:0．0０1:F_spａce＊Num_Ｓc／2+4； F＿spｅctruｍ　= zｅrｏs(Num_Sc,lengtｈ(F）); for　ｉ = -Nuｍ_Sc/2:1:Num_Sc/2 F_spectｒum(i+Num_Ｓｃ／2＋1,1：end) = sin(2*pi*(F-i*F＿space).＊Ｔs／2)./(２*pi*(Ｆ－i＊Ｆ_sｐace)．*Tｓ/２)； ｅｎd pｌot(F,Ｆ＿spectrｕm） grｉd oｎ %　ofdm sｕｂcarrier ｐｌot clear　aｌl; clｏse alｌ; clc; N = 2５6; M = 6; N_symbol　= １; ％ s_data = ［－3－sqｒt(-１)*3 －3-sｑrt（-１) -3－ｓqｒt(1） －3-sｑｒt（-1）*3]/ｓqrt(10)； ｓ_data ＝ (1+sqrt(-1)）／ｓqrｔ（2）; tx_ｄata = repmａｔ(s_data,1,N）; 　% sｏme initial phａｓe %　Oｒｉginaｌ_bin = rａｎdint（2，N*Ｎ＿symｂol，2,２２2）; % Mｏｄobj = moｄｅm.ｑamｍoｄ(＇M,4，SymｂolＯrｄer,Gray＇，InpuｔTypｅ＇，Bit)； %　tｘ_ｄａta = modulate(Mｏdｏbj,Oriｇｉnａl_ｂｉn）/sqrt（10）； ｆｏr　j = 0：1:N_sｙmbol-1 　 for　k　=0：1:N-1 for n = 0:1:Ｎ－1 　 　　 x＿tｍp(k+1,n+１+j*N) = tx_data（k+1+j*N）＊exp(sｑrt(-1）*2＊pi*k*n/N); 　　 　　 　 end 　end end % ｐloｔ fｏur-sｕbcarriｅｒ figure(１) ｐlｏt(rｅal(x_ｔｍp（1:4，1：ｅnd).）） cｌear；?ｃlc；SNＲ＝1０;　 % 信噪比?fｌ＝1２8; % 设置FFT长度Ns=6;　　 ％设置一个祯结构中OFDM信号的个数?para=1２８;%设置并行传输的子载波个数?sr＝2５０00０； %符号速率bｒ＝sｒ.*2;% 每个子载波的比特率?gl=32 %保护时隙的长度Sigｎal=rand(１，parａ*Ns*2)０.5;%产生0,1 随即序列,符号数为ｐarａ*Ns*2?for ｉ=1:para 　for　j=1:Ns*2 　 SigPaｒa(i,ｊ)＝Signａl(i*j）；%串并变换? 　 ｅndｅｎd%QＰSK调制,将数据分为两个通道for j=1：Ns ich(:,j)=SｉgPａra(:，2＊j－1);? 　　 qch(:，j)＝SｉgPaｒａ(:,２*ｊ);end?ｋmod=１．/sqrt（2);ich1=ｉch.*kｍod;ｑｃh1=qｃh.*kmod;?x=ich1＋qch１．*ｓqrt（－1); %频域数据变时域y=ifft（ｘ）;?ich2=rｅａl(y);?qch２=iｍag（y);?%插入保护间隔?ich３＝[ｉcｈ2(fl-ｇl+1:ｆｌ,:);iｃh2];qch3=[qcｈ2（fl－gl+1:fl，:);qｃh２]；?%并串变换?ｉch４＝ｒeshａpe(ich３,1,(ｆl+gl)＊Ns);?qch4＝ｒeｓhape(qcｈ3,１,(ｆl+gl)*Ns）;%形成复数发射数据?ＴrＤata=ich4+qch4.＊sｑrt(－1);%接收机%加入高斯白噪声ＲeData=ａwgn(TｒDａtａ,ＳＮＲ,ｍeaｓureｄ）;?%接收端%移去保护间隔idａta=ｒeａl(ＲｅData);ｑｄata=imag（ReData）；idａtａ1=reshape(ｉdaｔa，fl+gl,Ｎs）;qdaｔa1=reshａｐe（qdata,fl+gl，Ｎs）；ｉdata２＝idatａ1(ｇl+1:gl+fl,：）； 　ｑｄata２＝qdaｔa1（gl+1:gl+fl,:);　 　　 ?%FＦＴRex=idata2＋qdatａ2＊sqrｔ(-1）;ry=ffｔ(Rex）;?ReIＣhan＝reａl（rｙ);ＲeＱCｈaｎ=imａg