基于DSP的電力線載波OFDM調(diào)制解調(diào)器
1 OFDM原理
2 電力線載波通信系統(tǒng)結(jié)構(gòu)
Homeplug是工業(yè)界第一個電力線家庭網(wǎng)絡(luò)標(biāo)準(zhǔn).系統(tǒng)參考Homeplug采用的頻譜范圍4.5MHz~21MHz,并在Homeplug物理參數(shù)的基礎(chǔ)上確定本系統(tǒng)參數(shù)為:
采樣頻率fs=1/T = 15MHz
數(shù)據(jù)符號時間Td = 256×T=17.07μs
循環(huán)前綴時間Tcp = 172×T=11.47μs
OFDM符號時間Ts = 428×T=28.5μs
數(shù)據(jù)子載波數(shù)為256
子載波間隔Δf=1/Td=0.05858MHz
總子載波占用帶寬 N×Δf=15MHz
由于加入了11.47μs的循環(huán)前綴,系統(tǒng)可以消除11.47μs以內(nèi)的回波干擾.但是同時也付出頻帶利用率僅0.59B/Hz和損失功率2.23dB的代價.考慮到電力線惡劣的通信環(huán)境,付出的代價是值得的.
電力線高速通信系統(tǒng)的系統(tǒng)結(jié)構(gòu)如圖1所示.輸入數(shù)據(jù)在OFDM信號調(diào)制部分依次經(jīng)過串/并變換、IFFT、加入循環(huán)前綴、并/串變換后,輸出調(diào)制后的信號,其頻帶范圍為0~15MHz、數(shù)據(jù)速率為8.97MB.經(jīng)過調(diào)制的信號經(jīng)過數(shù)/模變換和上變頻后,通過系統(tǒng)耦合部分進(jìn)入電力線. 電力線上的信號通過系統(tǒng)耦合部分,輸出的信號通過下變頻、模/數(shù)變換后輸入給OFDM信號解調(diào)部分.在經(jīng)過串/并變換、去除循環(huán)前綴、FFT、并/串變換后,輸出串行數(shù)據(jù)流.
3 OFDM調(diào)制解調(diào)器的硬件實現(xiàn)
S5933是AMCC(Applied Micro Circuits Corporation)公司開發(fā)的32bit PCI控制器,具備強大、靈活的PCI接口功能,適用于高速數(shù)據(jù)傳輸場合.S5933芯片的特點是符合PCI2.1規(guī)范,支持PCI主、從兩種工作方式,支持多種數(shù)據(jù)傳輸方式,適用于不同的數(shù)據(jù)傳輸場合,支持PCI全速傳輸,提供8/16/32bit的Add-On用戶總線,有高低字節(jié)順序調(diào)整功能,支持穿行和并行的BOOT/POST碼功能,160腳PQFP封裝.
DSP部分選用TI公司的TMS320C6201.TMS320C6201有32位的外部存儲接口EMIF,為CPU訪問外圍設(shè)備提供了無縫接口.為了便于多信道數(shù)字信號處理,TMS320C6201配備了多信道帶緩沖能力的串口McBSP.McBSP的功能非常強大,除具有一般DSP串口功能之外,還可以支持T1/E1、ST-BUS、IOM2、SPI、IIS等不同標(biāo)準(zhǔn).TMS32C6201提供的16位主機接口(HPI)使得主機設(shè)備可以直接訪問DSP的存儲空間.通過內(nèi)部或外部存儲空間,主機可以與DSP交換信息,也可以利用HPI直接訪問映射進(jìn)存儲空間的外圍設(shè)備.TMS320C6201的DMA控制器有四個獨立的可編程通道,可以同時進(jìn)行四種不同的DMA操作.
4 OFDM在DSP上的軟件實現(xiàn)
調(diào)制部分的子程序被系統(tǒng)調(diào)用前,發(fā)送的數(shù)據(jù)已裝入數(shù)據(jù)存儲器.子程序被調(diào)用時,數(shù)據(jù)區(qū)的首地址以及長度被作為入口參數(shù)傳遞給子程序.程序執(zhí)行時首先進(jìn)行一系列的配置工作,如配置DSP片內(nèi)外設(shè)以及數(shù)模轉(zhuǎn)換器的各種參數(shù)等.之后,串口中斷產(chǎn)生,中斷服務(wù)程序自動依次讀取發(fā)送存儲器中的內(nèi)容,經(jīng)串口輸出給數(shù)模轉(zhuǎn)換器.然后程序從數(shù)據(jù)存儲區(qū)讀取一幀數(shù)據(jù),并行放入IFFT工作區(qū)的相應(yīng)位置,隨后進(jìn)行IFFT以及加入循環(huán)前綴(即復(fù)制數(shù)據(jù)的后若干位插入到數(shù)據(jù)的前段).所得數(shù)據(jù)存入發(fā)送存儲器以便中斷服務(wù)程序?qū)⑵漭敵?
5 FFT在TMS320C6201上的優(yōu)化算法
表1給出256點Radix2FFT和Radix4FFT在TMS320C6201上所需的指令周期,以及在不同的工作頻率下完成FFT所需的時間. 由表1可以看出,在TMS320C6201上采用Raidx4算法比采用Radix2算法更加高效.并且,為了滿足系統(tǒng)需求,即在17.07μs之內(nèi)完成256個復(fù)數(shù)點的FFT運算,TMS320C6201必須采用200MHz的工作頻率.
為使代碼達(dá)到最大效率,程序?qū)⒈M可能將指令安排為并行執(zhí)行.為使指令并行操作,程序確定指令間的相關(guān)性,即一條指令必須發(fā)生在另一條指令之后.根據(jù)TMS320C6201的數(shù)據(jù)通路和流水線工作方式,在此給出一種高效實現(xiàn)16點Radix4FFT的方法.其基本思想是分解傳統(tǒng)的FFT蝶型算法循環(huán)體,將其分別展開在A、B通路內(nèi)計算兩個FFT蝶型算法.每個蝶型算法分別只分配自己這一側(cè)的寄存器組和功能單元.這樣在循環(huán)體內(nèi)兩個蝶型算法是完全不相關(guān)的,能夠并行執(zhí)行.下面給出基于C.S.Burrus和T.W.Parks的Radix4FFT算法的優(yōu)化算法的代碼實現(xiàn).
void radix4(int n,short x[], short w[])
{
int n1,n2,ie,wa1,wa2,wa3, wb1, wb2,wb3,ia0,ia1,ia2,ia3,ib0,ib1,ib2,ib3,j,k;
short ta,tb,ra1,ra2, rb1,rb2,sa1,sa2,sb1,sb2,coa1,coa2,coa3,cob1,cob2,cob3,sia1,sia2,
sia3,sib1,sib2,sib3;
n2=n;
ie=1;
for(k=n;k>1;k>>=2)
{ //number of stage
n1=n2;
n2>>=2; // distance between input datas
wa1=0;
for(j=0;jwb1=wa1+ie;
wa2=wa1+wa1;
wb2=wb1+wb1; //since heremost of the folow-ering two instructions are parallel
wa3=wa2+wa1;
wb3=wb2+wb1;
coa1=w[wa1*2+1];
cob1=w[wb1*2+1];
sia1=w[wa1*2];
sib1=w[wb1*2];
coa2=w[wa2*2+1];
cob2=w[wb2*2+1];
sia2=w[wa2*2];
sib2=w[wb2*2];
coa3=w[wa3*2+1];
cob3=w[wb3*2+1];
sia3=w[wa3*2];
sib3=w[wb3*2];
wa1=wb1+ie;
for(ia0=j,ib0=j+1;ia0{//loop of two butterflies caculation
ia1=ia0+n2;
ib1=ib0+n2;
ia2=ia1+n2;
ib2=ib1+n2;
ia3=ia2+n2;
ib3=ib2+n2;
ra1=x[2*ia0]+x[2*ia2];
rb1=x[2*ib0]+x[2*ib2];
ra1=x[2*ia0]-x[2*ia2];
rb1=x[2*ib0]-x[2*ib2];
ta=x[2*ia1]+x[2*ia3];
tb=x[2*ib1]+x[2*ib3];
x[2*ia0]=ra1+ta; // x[2*ia0]
x[2*ib0]=rb1+tb; // x[2*ia0]
ra1=ra1-ta;
rb1=rb1-tb;
sa1=x[2*ia0+1]+x[2*ia2+1];
sb1=x[2*ib0+1]+x[2*ib2+1];
sa2=x[2*ia0+1]-x[2*ia2+1];
sb2=x[2*ib0+1]-x[2*ib2+1];
ta=x[2*ia1+1]+x[2*ia3+1];
tb=x[2*ib1+1]+x[2*ib3+1];
x[2*ia0+1]=sa1+ta;
x[2*ib0+1]=sb1+tb;
sa1=sa1-ta;
sb1=sb1-tb;
x[2*ia2]=(ra1*coa2+sa1*sia2)>>15;
x[2*ib2]=(rb1*cob2+sb2*sib2)>>15;
x[2*ia2+1]=(sa1*coa2-ra1*sia2)>>15;
x[2*ib2+1]=(sb1*cob2-rb1*sib2)>>15;
ta=x[2*ia1+1]-x[2*ia3+1];
ra1=ra2+ta;
rb1=rb2+tb;
ra2=ra2-ta;
rb2=rb2-tb;
ta=x[2*ia1]-x[2*ia3];
tb=x[2*ib1]-x[2*ib3];
sa1=sa2-ta;
sb1=sb2-tb;
sa2=sa2+ta;
sb2=sb2+tb;
x[2*ia1]=(ra1*coa1+sa1*sia1) >>15;
x[2*ib1]=(rb1*cob1+sb1*sib1) >>15;
x[2*ia1+1]=(sa1*coa1-ra1*sia1)>>15;
x[2*ib1+1]=(sb1*cob1-rb1*sib1)>>15;
x[2*ia3]=(ra2*coa3+sa2*sia3) >>15;
x[2*ib3]=(rb2*cob3+sb2*sib3) >>15;
x[2*ia3+1]=(sa2*coa3-ra2*sia3)>>15;
x[2*ib3+1]=(sb2*cob3-rb2*sib3)>>15;
}
}
ie =2
}
}
評論