移動(dòng)通信系統(tǒng)課設(shè)_OFDM系統(tǒng)仿真設(shè)計(jì)

1、移動(dòng)通信系統(tǒng)課程設(shè)計(jì)報(bào)告OFDM系統(tǒng)仿真目錄移動(dòng)通信系統(tǒng)課程設(shè)計(jì)報(bào)告1（一）題目要求：2（二）相關(guān)原理：21)OFDM :2)QPSK調(diào)制:3)導(dǎo)頻與均衡:4)循環(huán)前綴:5)分組交織:14（三）基本思路:（四）結(jié)果:1）軟解碼與硬解碼情況下不同信噪比的誤碼率2）不同信噪比下譯碼相位圖：1.5.（五）總結(jié)體會(huì):16（六）分工合作:17（七）程序代碼:17（一）題目要求：1）OFDM128 路傳輸；2）QPSK調(diào)制3）AWGN信道4）3徑或4徑瑞利衰落信道（二）相關(guān)原理：1） OFDM :將信道分成若干正交子信道，將高速數(shù)據(jù)信號(hào)轉(zhuǎn)換成并行的低速子數(shù)據(jù)流，調(diào)制到在每個(gè)子信道上進(jìn)行傳輸。正交信號(hào)可以通

2、過在接收端采用相關(guān)技術(shù)來分開，這樣可以減少子信道之間的相互干擾（ISI）。每個(gè)子信道上的信號(hào)帶寬小于信道的相關(guān)帶寬，因此每個(gè)子信道上可以看成平坦性衰落，從而可以消除碼間串?dāng)_，而且由于每個(gè)子信道的帶寬僅僅是原 信道帶寬的一小部分，信道均衡變得相對(duì)容易。OFDtvt子載;IS頻譜-8000-6000-400020000:2000400060008000'頻率2) QPSK 調(diào)制：將每?jī)蓚€(gè)相連比特組在一起形成雙比特碼元，它的四種狀態(tài)用4個(gè)不同的相位表示；3) 導(dǎo)頻與均衡：在OFDM信息序列中插入已知的導(dǎo)頻序列x n ,通過信道 后將其提取得y n ,做頻域除法得傳輸函數(shù)H z =衛(wèi)習(xí)，再通X

3、z過線性插值后得到每個(gè)信道頻率響應(yīng)，均衡濾波傳輸函數(shù)1Hz4) 循環(huán)前綴：循環(huán)前綴(Cyclic Prefix, CP)是將OFDM符號(hào)尾部的信號(hào)搬移到頭部構(gòu)成的。用來消去碼間干擾，通常取長(zhǎng)度Tg( 為信道沖激響應(yīng)持續(xù)時(shí)間）5）分組交織：為了解決成串的比特差錯(cuò)問題，采用了交織技術(shù)：把一條消息中的相 繼比特分散開的方法，即一條信息中的相繼比特 以非相繼方式發(fā)送，這樣即使在傳輸過程中發(fā)生了成串差錯(cuò) 恢復(fù)成一條相繼比特串的消息時(shí)，差錯(cuò)也就變成單個(gè)（或者 長(zhǎng)度很短）的錯(cuò)誤比特，這時(shí)再用信道糾正隨機(jī)差錯(cuò)的編碼 技術(shù)（FEC消除隨機(jī)差錯(cuò)。糾錯(cuò)數(shù)bmax D t（三）基本思路：信號(hào)產(chǎn)生->卷積編碼-&

4、gt;交織（32）QPSK調(diào)制（96路）>插入導(dǎo)頻信 號(hào)（32路）1才OFDM調(diào)制 （IFFT）插入循環(huán)前 綴（1/4長(zhǎng)度）->并串轉(zhuǎn)換AWGN信道->3徑瑞利信道17串并轉(zhuǎn)換OFD M解調(diào) （FFT）取導(dǎo)頻信道估計(jì)與 均衡并串轉(zhuǎn)換QPSK解調(diào) 解交織卷積譯碼誤碼率統(tǒng)計(jì)說明：1）編碼：使用216卷積碼;相關(guān)代碼：% 卷積編碼 %ori_data_t_ 仁reshape(ori_data_t, numn f,2*cha nn ell); %58*192 ori_data_t_2=zeros( num,2*cha nn ell);%128*192for i = 1:2*cha n

5、n ellseq = ori_data_t_1(:,i)：seq_code = en code216(seq)：ori_data_t_2(:,i) = seq_code;endori_data=resh ape(ori_data_t_2,1,2* nu m*cha nn ell);%1*24576編碼函數(shù)：fun ctio n code=e ncode216(m)%輸入信息序列%g1=1,0,0,0,0,0;g2=1,1,0,0,1,1;trel=poly2trellis(6,40 63);% 定義網(wǎng)格m仁m,0,0,0,0,0,0;code=c onven c(m1,trel);%卷積碼編碼

6、2) 交織：交織深度為32 ；相關(guān)代碼：%編碼交織% ori_data1=reshape(ori_data,deep,(2* nu m*cha nn ell)/deep)：%768*32ori_data2=reshape(ori_data1,1,2* nu m*cha nn ell);%1*245763) QPSK調(diào)制：轉(zhuǎn)換為96路相位信號(hào)(復(fù)數(shù))；代碼：%QPSK調(diào)制 % ori_data_0=reshape(ori_data2,2,channell*num);%拆分成兩行 12288 列ori_data_仁bi2de(ori_data_0','left-msb')&

7、#39;%QPSK 的未調(diào)制數(shù)據(jù)12288列M=4的數(shù)據(jù)de_OFDM_1=modem.pskmod(4);%生成調(diào)制器對(duì)象，設(shè)置 qpsk調(diào)制QPSK_data=modulate(de_OFDM_1,ori_data_1);%1*12288 復(fù)數(shù)形式%cha nn el_data=reshape(QPSK_data,cha nn ell, nu m);%拆分成128行，每行10個(gè)復(fù)數(shù)128*10cha nn el_data=reshape(QPSK_data ,nu m,cha nn ell);%拆分成128行，每行96個(gè)復(fù)數(shù)128*964）插入導(dǎo)頻：32路已知序列；代碼：%插入導(dǎo)頻信號(hào)%p

8、liot_ 仁2*ra nd(128,32);pliot_2=2*rand(128,32);%插入的導(dǎo)頻信號(hào)，實(shí)數(shù)pliot=pliot_1+pliot_2*1i; % 插入的導(dǎo)頻信號(hào)，復(fù)數(shù) chann el_data_ 1= zeros(128,128);for i = 0:31chann el_data_1(:,4*i+1)=pliot(:,i+1);chann el_data_1(:,4*i+2)=cha nn el_data(:,3*i+1);chann el_data_1(:,4*i+3)=cha nn el_data(:,3*i+2);chann el_data_1(:,4*i+4

9、)=cha nn el_data(:,3*i+3);end5)OFDM調(diào)制：使用IFFT;相關(guān)代碼:%OFDM生%IFFT_data=zeros(128,128);for i = 1:128%128*128ifftIFFT_data(:,i)=ifft(cha nn el_data_1(:,i),128);endr_ifft_data=real(IFFT_data);%實(shí)部128*128i_ifft_data=imag(IFFT_data);%虛部128*1286)插入循環(huán)前綴：這里選取長(zhǎng)度為32 ；相關(guān)代碼：綴%r_data=r_ifft_data( num-num/4+1: num,:);

10、r_ifft_data; %保護(hù)間隔取傳送數(shù)據(jù)的四分之一i_data=i_ifft_data( num-n um/4+1: nu m,:);i_ifft_data; %160*1287) 并串轉(zhuǎn)換：轉(zhuǎn)換成1路；代碼：%并串換%r_out=reshape(r_data,1,( num+n um/4)*cha nn el);%1*1600i_out=reshape(i_data,1,( nu m+num /4)*cha nn el);%1*1600data=r_out+i_out.*1i;8) AWGN :加入信噪比-1020dB 的噪聲；代碼：%AWGN 信道 % data_AWGN=awg n

11、(data,SNR,'measured') ;%力口入噪聲9) 串并轉(zhuǎn)換：轉(zhuǎn)換成128路；代碼：%串并換 % %r_rcv_data1=reshape(r_rcv_data, num+nu m/4,cha nn el);i_rcv_data1=reshape(i_rcv_data, num+nu m/4,cha nn el);10) OFDM解調(diào)：FFT快速算法；代碼：%OFDM解調(diào)FFT% de_OFDM_data=zeros(128,128);for i = 1:128de_OFDM_data(:,i)=fft(rcv_data(:,i),128);%128*128 矩陣e

12、nd11) 去導(dǎo)頻：取出過信道后導(dǎo)頻序列yn ；12) 信道估計(jì)：計(jì)算收到每一路導(dǎo)頻序列的Yz，已知導(dǎo)頻序列的Xz，得到信道H z=g4，利用線性插值得到每一Xz路信息序列的Hi z ；13) 均衡：將每一路信號(hào)Xiz Eiz，Eiz1 ，得到均衡后Hiz的信號(hào)頻域Yiz，然后得到y(tǒng)'(n)；代碼：%信道估計(jì)與均衡(新的)depliot=zeros(128,32);for i = 0:31depliot(:,i+1)=de_OFDM_data(:,4*i+1);enddepliot_ 1= zeros(128,32);depliot_2=zeros(32,128);de_OFDM_0=

13、zeros(128,32);for i = 1:32depliot_1(:,i)=pliot(:,i)./depliot(:,i);%depliot_2(i,:)=reshape(depliot_1(:,i),1,128);%x=1:128;%hx=1:128;%de_OFDM_0(:,i)=i nterp1(x,depliot_2(i,:),hx,'pchip');% 線性插值endde_OFDM_1=zeros(128,128);for i= 0:30de_OFDM_1(:,4*i+1)=depliot_1(:,i+1);de_OFDM_1(:,4*i+2)=depliot

14、_1(:,i+1)+(depliot_1(:,i+2)- depliot_1(:,i+1)/4;de_OFDM_1(:,4*i+3)=depliot_1(:,i+1)+(depliot_1(:,i+2)-depliot_1(:,i+1)/2;de_OFDM_1(:,4*i+4)=depliot_1(:,i+1)+(depliot_1(:,i+2)-depliot_1(:,i+1)/4*3;endde_OFDM_1(:,125)=depliot_1(:,32);de_OFDM_1(:,126)=depliot_1(:,32);de_OFDM_1(:,127)=depliot_1(:,32);de

15、_OFDM_1(:,128)=depliot_1(:,32);de_OFDM_2=de_OFDM_data.*de_OFDM_1;% 還原數(shù)據(jù)de_OFDM_3=;for i = 2:4:126de_OFDM_3=de_OFDM_3,de_OFDM_2(:,i:i+2); % 提取出 128*96 數(shù)據(jù)部分end14) 并串轉(zhuǎn)換：變?yōu)?6路信息序列;代碼：%并串轉(zhuǎn)換%de_OFDM_data3=reshape(de_OFDM_3,1,cha nnell*nu m);%128*9615) QPSK解調(diào)：還原調(diào)制前的數(shù)據(jù)；代碼：%QPSK解調(diào) %de_OFDM_data3_real=real(de

16、_OFDM_data3);de_OFDM_data3_imag=imag(de_OFDM_data3);de_OFDM_4=zeros(1,2*le ngth(de_OFDM_data3_real);%求方差%sigma=sum(abs(real(de_OFDM_data3)-real(QPSK_data)+abs(imag(de_OFDM_data3)-imag(QPSK_data)/(2*le ngth(de_OFDM_data3_real);fori=1:le ngth(de_OFDM_data3_real)%求判決到各個(gè)位置的概率%re_f1=n ormpdf(de_OFDM_data

17、3_real(i),-1,sigma);re_O=normpdf(de_OFDM_data3_real(i),0,sigma);re_1=n ormpdf(de_OFDM_data3_real(i),1,sigma);im_f1=n ormpdf(de_OFDM_data3_imag(i),-1,sigma);im_0=n ormpdf(de_OFDM_data3_imag(i),0,sigma);im_1= no rmpdf(de_OFDM_data3_imag(i),1,sigma);de_OFDM_4(2*i-1)=(re_f1+re_0*im_f1/(im_f1+im_1)/(re_f

18、1+re_1+re_0);de_OFDM_4(2*i)=(im1+im_1)/(im1+im_1+im_0);enddata_last_0=de_OFDM_4;%1*2457616) 解交織；代碼：%解交織%data_last_0_0=reshape(data_last_0,(2* nu m*cha nn ell)/deep,deep)'data_last_0_ 仁reshape(data_last_0_0,1,2* nu m*cha nn ell);17) 卷積譯碼：分別使用維特比譯碼，軟解碼，硬解碼方法，得到譯碼結(jié)果；代碼：%卷積碼碼%data_last_ 仁reshape(dat

19、a_last_0_1, num,cha nn ell*2);%128*256data_last_2=zeros( numn f,cha nn ell*2);%58*256for k = 1:2*cha nn ellseqO = data_last_1(:,k)：seq0=1-2.*seq0;seq0_code = decode216(seq0,1)：data_last_2(:,k) = seq0_code;enddata_last = reshape(data_last_2,1,2* numn f*cha nn ell);%1*14848譯碼函數(shù)：function decode=decode2

20、16(m,k)%m為碼字序列，硬判決譯碼輸入 k為1軟判決譯碼輸入k為0trel=poly2trellis(6,40 63);%定義網(wǎng)格tblen=35;%回朔長(zhǎng)度，為約束長(zhǎng)度的5倍if(k=1)%硬判決譯碼hcode=m<0;decode1=vitdec(hcode,trel,tble n,'term','hard');else%軟判決量化譯碼(3比特量化)in dex,hcode=qua ntiz(m,-0.75,-0.5,-0.25,0,0.25,0.5,0.75,7,6,5,4,3,2,1,0);decode1=vitdec(hcode,trel,

21、tble n,'term','soft',3);enddecode=decode1(1:le ngth(decode1)-6);18)誤碼率統(tǒng)計(jì)；代碼：%誤碼率計(jì)% error_ num=sum(abs(data_last-ori_data_t);bit_error_rate(1,(SNR-SNR_mi n)/2+1)= error_ num/N; % 誤碼率sn r(1,(SNR-SNR_mi n)/2+1)=SNR;endstem(s nr,-10*log10(bit_error_rate),'-*','L in eStyle'

22、;,' non e');xlabel('SNR');ylabel('BER');title('Performa nee of OFDM un der the cha nnel AWGN')（四）結(jié)果：1）軟解碼與硬解碼情況下不同信噪比的誤碼率：Performa nee of OFDM un der the cha nnel AWGN-0.5soft decode hard decode-1.5-2681012141618SNR(dB)由上圖可以看出，信號(hào)誤碼性能隨著信噪比增加而增 力口，當(dāng)信噪比達(dá)到10dB的時(shí)候誤碼率已經(jīng)低于0.

23、01，信 噪比達(dá)到16dB時(shí)誤碼率已經(jīng)達(dá)到-5量級(jí)，再之后出現(xiàn)相 關(guān)波動(dòng)；另外，由于此次用于測(cè)試的數(shù)據(jù)只有24576個(gè)，所以誤碼率更高的時(shí)候誤碼率已經(jīng)變?yōu)?。同時(shí)，在低信噪比時(shí)，軟解碼與硬解碼性能差別不大； 在信噪比達(dá)到8dB之后，軟解碼性能明顯比硬解碼好；2) 不同信噪比下譯碼相位圖：Scatter plot2015oSNR=（dB）eluLakaauQo5-20 -10 0 10 20 In-Phase5o5-o2-5-4-3-2-1012345In-PhaseScatter plotprnrardauQIn-Phase由上圖可以看出，在不同信噪比下，相位圖譯碼結(jié)果不一樣；信噪比越大，相位

24、譯碼越清晰；（五）總結(jié)體會(huì)：1）這次課設(shè)時(shí)間周期較長(zhǎng)，一邊學(xué)習(xí)一邊編寫，前后花了好幾周；最開始搭建基本框架比較順利，通過對(duì)課程的理解，成功地理清了思路；2）編寫導(dǎo)頻和均衡部分時(shí)，由于書上講得不是很透徹，網(wǎng)上查閱資料也紛繁復(fù)雜，幾乎陷入了瓶頸；后來慢慢摸索，詢問老師，小組討論，同學(xué)交流；慢慢知道了自己的問題 所在，混淆了頻域和時(shí)域；3）小組討論后將程序分工，成員各自負(fù)責(zé)不同部分，有問題時(shí)一起交流討論，這次團(tuán)隊(duì)協(xié)作非常愉快，并且高效；4）從最開始做的無(wú)從下手到現(xiàn)在基本理解吃透，經(jīng)歷了一段時(shí)間。中間也詢問了很多人，查閱了許多資料；非常感謝 老師的指導(dǎo)，也希望李老師以后講解可以再耐心一些；（六）分工合

25、作：我們小組成員有倪鴻志，王旭，周裔歡；這次我（倪鴻志） 負(fù)責(zé)的是OFDM系統(tǒng)，信道的編寫，王旭負(fù)責(zé)了編碼譯碼部 分，周裔歡寫了 QPSK調(diào)制，剩下的均衡，導(dǎo)頻，軟譯碼是 大家一起商量完成的。（七）程序代碼：1）主程序：clc;clear all;% 參數(shù)設(shè)置 %p=0.5;%設(shè)定p值cha nn el=128; nu m=128;cha nn ell=96;deep=32;%交織深度numnf=num/2-6;%每一組信息位bit,注意復(fù)數(shù)包含兩組,58N=cha nn ell*2* num_inf;SNR_mi n=6;SNR_max=20;echo_ num=8;ts=1e-5;fd=0

26、;tau=ts 3*ts 8*ts;pdb=0 -3 -8 ;%采樣時(shí)間為 1us ,無(wú)多普勒頻移bit_error_rate=zeros(echo_ num,(SNR_max-SNR_mi n)/2+1);sn r=zeros(1,(SNR_max-SNR_mi n) /2+1);% 信噪比矩陣初始化for SNR=SNR_mi n:2:SNR_maxfor echo=1:echo_ num;%初始信號(hào)產(chǎn)生 %ori_data=ra nd(1,N);ori_data_t=ori_data<p;% 生成伯努利信息序列 1*11136%卷積編碼%ori_data_t_ 1=reshape(

27、ori_data_t, numn f,2*cha nn ell); %58*192ori_data_t_2=zeros( nu m,2*cha nn ell);%128*192for i = 1:2*cha nn ellseq = ori_data_t_1(:,i)'seq_code = en code216(seq)'ori_data_t_2(:,i) = seq_code;endori_data=reshape(ori_data_t_2,1,2* nu m*cha nn ell);%1*24576%編碼交織%ori_data1=reshape(ori_data,deep,(

28、2* nu m*cha nn ell)/deep)：%768*32ori_data2=reshape(ori_data1,1,2* nu m*cha nn ell);%1*24576%QPSK調(diào)制%ori_data_0=reshape(ori_data2,2,channell*num);%拆分成兩行 12288 列ori_data_仁bi2de(ori_data_0','left-msb')'%QPSK 的未調(diào)制數(shù)據(jù) 12288 列M=4的數(shù)據(jù)de_OFDM_1=modem.pskmod(4);% 生成調(diào)制器對(duì)象，設(shè)置qpsk調(diào)制QPSK_data=modula

29、te(de_OFDM_1,ori_data_1);%1*12288 復(fù)數(shù)形式%channel_data=reshape(QPSK_data,channell,num);%拆分成 128 行，每行 10個(gè)復(fù)數(shù)128*10channel_data=reshape(QPSK_data,num,channell);%拆分成 128 行，每行 96個(gè)復(fù)數(shù)128*96%插入導(dǎo)頻信號(hào)%pliot_ 仁2*ra nd(128,32);pliot_2=2*rand(128,32);%插入的導(dǎo)頻信號(hào),實(shí)數(shù)pliot=pliot_1+pliot_2*1i; % 插入的導(dǎo)頻信號(hào),復(fù)數(shù) cha nn el_data_

30、仁 zeros(128,128);for i = 0:31cha nn el_data_1(:,4*i+1)=pliot(:,i+1);cha nn el_data_1(:,4*i+2)=cha nn el_data(:,3*i+1);cha nn el_data_1(:,4*i+3)=cha nn el_data(:,3*i+2);cha nn el_data_1(:,4*i+4)=cha nn el_data(:,3*i+3);end%OFDM 信號(hào)產(chǎn)生 %IFFT_data=ifft(cha nn el_data);%128*10ifft 數(shù)據(jù)IFFT_data=zeros(128,12

31、8);for i = 1:128IFFT_data(:,i)=ifft(cha nn el_data_1(:,i),128);%128*128ifft 數(shù)據(jù)endr_ifft_data=real(IFFT_data);% 實(shí)部128*128i_ifft_data=imag(IFFT_data);% 虛部 128*128%加入循環(huán)前綴 %r_data=r_ifft_data( num-nu m/4+1: nu m,:);r_ifft_data; %保護(hù)間隔取傳送數(shù)據(jù)的四分之一i_data=i_ifft_data(num-num/4+1:num,:);i_ifft_data; %160*128%并

32、串轉(zhuǎn)換%r_out=reshape(r_data,1,( num+nu m/4)*cha nn el);%1*1600i_out=reshape(i_data,1,( num+nu m/4)*cha nn el);%1*1600data=r_out+i_out.*1i;%AWGN信道%data_AWGN=awg n(data,SNR,'measured') ;%加入噪聲%瑞利信道%cha n = rayleighcha n(ts,fd,tau,pdb);ru_data = filter(cha n,data_AWGN);%接收號(hào)%r_rcv_data=real(ru_data)

33、;i_rcv_data=imag(ru_data);%串并換%r_rcv_data1=reshape(r_rcv_data, num+nu m/4,cha nn el);i_rcv_data1=reshape(i_rcv_data, num+nu m/4,cha nn el);%去綴% r_rcv_data2=r_rcv_data1( nu m/4+1: num+nu m/4,:); i_rcv_data2=i_rcv_data1( nu m/4+1: num+nu m/4,:); rcv_data=r_rcv_data2+i_rcv_data2*1i;%OFDMFFT%de_OFDM_dat

34、a=zeros(128,128);for i = 1:128de_OFDM_data(:,i)=fft(rcv_data(:,i),128);%128*128 矩陣end%信道估計(jì)與均衡(新的)depliot=zeros(128,32);for i = 0:31depliot(:,i+1)=de_OFDM_data(:,4*i+1);enddepliot_ 1= zeros(128,32);depliot_2=zeros(32,128); de_OFDM_0=zeros(128,32);for i = 1:32depliot_1(:,i)=pliot(:,i)./depliot(:,i);en

35、dde_OFDM_仁zeros(128,128);for i= 0:30de_OFDM_1(:,4*i+1)=depliot_1(:,i+1);de_OFDM_1(:,4*i+2)=depliot_1(:,i+1)+(depliot_1(:,i+2)-depliot_1(:,i+1)/4;de_OFDM_1(:,4*i+3)=depliot_1(:,i+1)+(depliot_1(:,i+2)-depliot_1(:,i+1)/2;de_OFDM_1(:,4*i+4)=depliot_1(:,i+1)+(depliot_1(:,i+2)-depliot_1(:,i+1)/4*3;endde_O

36、FDM_1(:,125)=depliot_1(:,32);de_OFDM_1(:,126)=depliot_1(:,32);de_OFDM 1(:,127)=depliot_1(:,32);de_OFDM_1(:,128)=depliot_1(:,32);de_OFDM_2=de_OFDM_data.*de_OFDM_1; % 還原數(shù)據(jù)de_OFDM_3=;for i = 2:4:126de_OFDM_3=de_OFDM_3,de_OFDM_2(:,i:i+2);% 提取出 128*96 的數(shù)據(jù)部分end%并串轉(zhuǎn)換%de_OFDM_data3=reshape(de_OFDM_3,1,cha n

37、n ell* nu m);%128*96%QPSK調(diào)%h1 = modem.pskdemod(4);%de_QPSK_data=demodulate(h1,de_OFDM_data3);%demodata=de2bi(de_QPSK_data,'left-msb')'de_OFDM_data3_real=real(de_OFDM_data3);de_OFDM_data3_imag=imag(de_OFDM_data3);de_OFDM_4=zeros(1,2*le ngth(de_0FDM_data3_real);%求方差 %sigma=(sum(de_OFDM_4-o

38、ri_data)/(2*le ngth(de_OFDM_data3_real)A2;for i=1:le ngth(de_OFDM_data3_real)%求判決到各個(gè)位置的概率 %re_f1=n ormpdf(de_OFDM_data3_real(i),-1,sigma);re_0=normpdf(de_OFDM_data3_real(i),0,sigma);re_1=n ormpdf(de_OFDM_data3_real(i),1,sigma);im_f1=n ormpdf(de_OFDM_data3_imag(i),-1,sigma);im_0=normpdf(de_OFDM_data3

39、_imag(i),0,sigma);im_1= no rmpdf(de_OFDM_data3_imag(i),1,sigma);de_OFDM_4(2*i-1)=(re_f1+re_O*im_f1/(im_f1+im_1)/(re_f1+re_1+re_O);de_OFDM_4(2*i)=(im_f1+im_1)/(im_f1+im_1+im_O);enddata_last_0=de_OFDM_4;%1*24576% 解交織 %data_last_0_0=reshape(data_last_0,(2* nu m*cha nn ell)/deep,deep)：data_last_0_ 仁 reshape(data_last_0_0,1,2* nu m*cha nn ell);%卷積碼譯碼%data_last_1=