【文章內(nèi)容簡介】 ，導(dǎo)頻信號可以在時間和頻率兩個方向上插入，在接收端可提取導(dǎo)頻信號估計信道的傳遞函數(shù)。只要導(dǎo)頻信號在時間和頻率方向上間隔對于信道帶寬足夠少，就可以采用二維內(nèi)插濾波的方法來估計傳遞函數(shù)，當(dāng)然也可以采用分離的一維估計。OFDM系統(tǒng)中常用的導(dǎo)頻信號分布方法有導(dǎo)頻信號塊狀分布、梳狀分布和星狀分布三種?？紤]到實現(xiàn)的復(fù)雜度，信道估計準則選用LS估計準則。插值方法有常值內(nèi)插、線性內(nèi)插和DFT插值，常值內(nèi)插一般用在塊狀導(dǎo)頻結(jié)構(gòu)中，是比較簡單的插值方法，本文接下來就來討論LS算法下不同插值方式下對信道的估計；首先線性內(nèi)插是最簡單也是最傳統(tǒng)的內(nèi)插方法之一，它利用兩個導(dǎo)頻信號來進行內(nèi)插估計。時間方向的線性內(nèi)插的公式為： ()其中 , 。同理，可以得到頻率方向的一階線性內(nèi)插的公式為： ()其中 , 。其次是DFT插值，由于信道沖擊響應(yīng)與信道傳輸函數(shù)是傅氏變換對，內(nèi)插可以利用DFT的性質(zhì)。但是DFT插值一般用在基于梳狀導(dǎo)頻的結(jié)構(gòu)中設(shè)信道沖擊向為，…，0，0…0。信道的傳輸函數(shù)為： ()取整數(shù)，且N是M的整數(shù)倍，對信道傳輸函數(shù)在頻．率方向以N/M為間隔進行抽取，得到其中的元素是： ()可以看出，由頻率軸的M個抽樣值可以恢復(fù)信道沖擊響應(yīng)。再進行N點的DFT就可以得到所有子信道的傳輸函數(shù)值。至于常值插入比較簡單就不再贅述。,由圖可以看出同一信噪比下DFT最為理想，線性內(nèi)插效果最差，而常值內(nèi)插介于他倆之間，但在要求同一REB的情況下DFT要求更大信噪比。所以在大信噪比下還是選擇DFT更為理想，如果是在要求小的誤碼率且在小的信噪比下常值內(nèi)插是一種比較簡單而且效果較理想的內(nèi)插恢復(fù)方法。結(jié) 論本文針對目前的研究熱點OFDM技術(shù)進行計算機仿真研究,在OFDM仿真模型的基礎(chǔ)上用MATLAB語言編寫出OFDM發(fā)送、信道及接收整個系統(tǒng)上的仿真圖形,在系統(tǒng)仿真正確的前提下,對在OFDM信道上加上窗函數(shù)前后以及加上循環(huán)前綴后，采用不同的內(nèi)插方法接收信號的改善程度進行了研究，得出預(yù)想的結(jié)果。致 謝參 考 文 獻[1] Erich Cosby. Orthogonal Frequency Division Multiplexing (OFDM) Tutorial and Analysis[M ]. Northern Virginia Center,2001.[2] Mingqi Li ,Qicong Peng, Yubai Li, Performance Evaluation of MC DS CDMA Systems in Multipath Fading Channels[ J ]. 0 7803 7547 5 /02, IEEE , 2002.[ 3 ] 　A Peled, A Ruiz. Frequency domain data transmission using reduced putational plexity algorithms[C ]. In Proc. IEEE Int. Conf. Acoust. , Speech, Signal Processing, 967.[ 4 ] 　R van Nee. OFDM Wireless Multimedia Communications[M ].Rrasad R. Artech House, 1998[ 5 ] 　周正蘭,等. OFDM及其鏈路級平臺的Simulink實現(xiàn)[ J ]. 中國數(shù)據(jù)通信, 2003, (10) : 90 – 92[ 6 ] 　 尹澤明,等. 精通MATLAB6 [M ]. 清華大學(xué)出版社, 2002.[ 7 ] 　蔡濤, 等譯. 無線通信原理與應(yīng)用[M ]. 電子工業(yè)出版社,1999.[ 8 ] 丁玉美,等. 數(shù)字信號處理[M]. 西安電子科技大學(xué)出版社, 2003[9 ] Reiniers : the COFDMbased system for terrestrial television[J]. Electronics amp。 Communication Engineering Journal, 1997,9,(01):2832.[10] 尹長川,羅濤,[M].北京:北京郵電大學(xué)出版社,2004:2045.[11] 王文博,鄭侃. 寬帶無線通信OFDM技術(shù)[M]. 第2版，內(nèi)蒙古:人民郵電出版社,2007:89.附 錄 clear all。close all。IFFT_bin_length = 1024。 % FFT的點數(shù)carrier_count = 200。 % 載波的數(shù)量 bits_per_symbol = 2。 % 每個符號代表的比特數(shù)symbols_per_carrier = 50。 % 每個載波使用的符號數(shù)SNR = 10。 % 信道中的信噪比（dB） baseband_out_length=carrier_count*symbols_per_carrier*bits_per_symbol。%總比特數(shù)carriers = (1:carrier_count) + (floor(IFFT_bin_length/4) floor(carrier_count/2))。conjugate_carriers = IFFT_bin_length carriers + 2。%發(fā)送端%產(chǎn)生隨機二進制數(shù)據(jù):baseband_out = round(rand(1,baseband_out_length))。convert_matrix=reshape(baseband_out,bits_per_symbol,length(baseband_out)/bits_per_symbol)。for k = 1:(length(baseband_out)/bits_per_symbol) modulo_baseband(k) = 0。 for i = 1:bits_per_symbolmodulo_baseband(k)=modulo_baseband(k)+convert_matrix(i,k)*2^(bits_per_symboli)。endend% 串并轉(zhuǎn)換carrier_matrix = reshape(modulo_baseband, carrier_count, symbols_per_carrier)39。% 對每一個載波的符號進行差分編碼carrier_matrix = [zeros(1,carrier_count)。carrier_matrix]。for i = 2:(symbols_per_carrier + 1) carrier_matrix(i,:)=rem(carrier_matrix(i,:)+carrier_matrix(i1,:),2^bits_per_symbol)。end% 把差分符號代碼轉(zhuǎn)換成相位carrier_matrix = carrier_matrix * ((2*pi)/(2^bits_per_symbol))。% 把相位轉(zhuǎn)換成復(fù)數(shù)[X,Y] = pol2cart(carrier_matrix, ones(size(carrier_matrix,1),size(carrier_matrix,2)))。plex_carrier_matrix = plex(X,Y)。% 分配載波到指定的IFFT位置IFFT_modulation = zeros(symbols_per_carrier + 1, IFFT_bin_length)。IFFT_modulation(:,carriers) = plex_carrier_matrix。IFFT_modulation(:,conjugate_carriers) = conj(plex_carrier_matrix)。% 畫出頻域中的OFDM信號代表figure (1)stem(0:IFFT_bin_length1, abs(IFFT_modulation(2,1:IFFT_bin_length)),39。b*39。)grid onaxis ([0 IFFT_bin_length ])ylabel(39。Magnitude39。)xlabel(39。IFFT Bin39。)title(39。OFDM Carrier Frequency Magnitude39。)% figure (2)plot(0:IFFT_bin_length1, (180/pi)*angle(IFFT_modulation(2,1:IFFT_bin_length)), 39。go39。)hold onstem(carriers1, (180/pi)*angle(IFFT_modulation(2,carriers)),39。b*39。)stem(conjugate_carriers1, (180/pi)*angle(IFFT_modulation(2,conjugate_carriers)),39。b*39。)axis ([0 IFFT_bin_length 200 +200])grid onylabel(39。Phase (degrees)39。)xlabel(39。IFFT Bin39。)title(39。OFDM Carrier Phase39。)% 通過IFFT將頻域轉(zhuǎn)化為時域，得到時域信號time_wave_matrix = ifft(IFFT_modulation39。)。time_wave_matrix = time_wave_matrix39。%畫出一個符號周期的時域OFDM信號figure (3)plot(0:IFFT_bin_length1,time_wave_matrix(2,:))grid onylabel(39。Amplitude39。)xlabel(39。Time39。)title(39。OFDM Time Signal, One Symbol Period39。)%畫出每一個載波對應(yīng)的時域信號（分離的OFDM信號）for f = 1:carrier_count temp_bins(1:IFFT_bin_length)=0+0j。 temp_bins(carriers(f))=IFFT_modulation(2,carriers(f))。 temp_bins(conjugate_carriers(f))=IFFT_modulation(2,conjugate_carriers(f))。 temp_time = ifft(temp_bins39。)。 figure(4) plot(0:IFFT_bin_length1, temp_time) hold onendgrid onylabel(39。Amplitude39。)xlabel(39。Time39。)title(39。Separated Time Waveforms Carriers39。)for i = 1:symbols_per_carrier + 1windowed_time_wave_matrix(i,:)=real(time_wave_matrix(i,:)).*hamming(IFFT_bin_length)39。windowed_time_wave_matrix(i,:) = real(time_wave_matrix(i,:))。end%串并轉(zhuǎn)換ofdm_modulation=reshape(windowed_time_wave_matrix39。,1,IFFT_bin_length*(symbols_per_carrier+