3G移動(dòng)通信實(shí)驗(yàn)前3

/３G移動(dòng)通信實(shí)驗(yàn)指導(dǎo)教師：章堅(jiān)武學(xué)生姓名:周云杰學(xué)生學(xué)號(hào)：1２081437實(shí)驗(yàn)日期：2０15/3、４實(shí)驗(yàn)一話務(wù)量呼損率計(jì)算機(jī)仿真【實(shí)驗(yàn)?zāi)康摹考由钤拕?wù)量和呼損率概念的理解;能夠使用C語言（或者M(jìn)aｔｌaｂ)計(jì)算話務(wù)量和呼損率之間的關(guān)系；【實(shí)驗(yàn)內(nèi)容】使用Ｃ語言（或者Ｍatlａb）仿真話務(wù)量、呼損率和信道數(shù)之間的數(shù)值關(guān)系;分析話務(wù)量和呼損率之間的關(guān)系；【實(shí)驗(yàn)設(shè)備】一臺(tái)PC機(jī)【實(shí)驗(yàn)步驟】1.采用話務(wù)量、呼損率和信道數(shù)之間解析表達(dá)式計(jì)算出不同信道數(shù)和不同話務(wù)量條件下的呼損率，畫出三維圖.（可以使用Matｌab三維繪圖語句Ｍesh）分析話務(wù)量、呼損率和信道數(shù)之間的關(guān)系.【實(shí)驗(yàn)結(jié)果】clc;cｌｅarall;A=1：0．5：60；n=1:50;forｉ=1：lenｇtｈ（Ａ）foｒj=1:lenｇth(n）teｍp(i，j）=A(ｉ）^n（j）/facｔoｒial(n(j))；suｍ（ｉ,ｊ)＝1；ｆｏrk=1：ｎ(j)sum(i,j)=ｓum（i,j)+A(i）^k/factｏｒial(k）;endendendＢ=ｔemｐ。/suｍ；mｅsh(n,A，B）；xlａbel（＇ｎ’）；ylabeｌ(’A’);zｌabel(’Ｂ’）; 由圖可得結(jié)論：當(dāng)話務(wù)量不變，呼損率隨著信道數(shù)增加而降低；當(dāng)信道數(shù)不變時(shí)，呼損率隨著話務(wù)量增加而增加。實(shí)驗(yàn)二Ｏｋuｍura－Hata方法計(jì)算機(jī)仿真一、實(shí)驗(yàn)?zāi)康?、加深對(duì)奧村模型的理解;2、能夠使用C語言（或者Ｍatｌａb）利用Okｕｍuｒa—Haｔa方法計(jì)算基本傳輸損耗;二、實(shí)驗(yàn)內(nèi)容1、使用C語言(或者M(jìn)atlab)利用Okｕmuｒa—Hata方法計(jì)算基本傳輸損耗；２、分析Oｋumurａ—Hａｔａ方法的誤差;三、實(shí)驗(yàn)步驟1、采用Ｏｋumura－Hata方法分別計(jì)算大城市市區(qū)地區(qū)準(zhǔn)平滑地形、郊區(qū)和開闊區(qū),基站天線高度是2０0米，手機(jī)天線高度是3米情況下,不同傳播距離和不同載波頻率條件下的傳播損耗中值。畫出相應(yīng)的曲線。2、將計(jì)算結(jié)果和通過奧村模型實(shí)測(cè)測(cè)得的結(jié)果進(jìn)行比較，驗(yàn)證計(jì)算結(jié)果的正確性。３、分析Okumｕra-Hａt(yī)ａ方法在何距離以及何頻率范圍內(nèi)存在較大的誤差。四、實(shí)驗(yàn)結(jié)果ｃlearall;clｏseａll；ｃlc;hb=200；hm＝3;lｂ1=0;lb２＝0；lb3=0;lｂ4＝０；ｆord=［125１030５０６0８０1００]ｆ1=１０0:0．1：３00;f2=300：０。1：３000；lb11=６9.5５+2６。１6＊log10(ｆ1）—１3．82＊log1０（ｈb）—（8.29＊（ｌoｇ１0(1.５4＊hｍ）。^２）—1．１）＋（（44。９-6.55＊log10（hb））＊loｇ10（ｄ））；lb12=６9。5５＋26。16*ｌoｇ１0（ｆ2）－１3。82＊loｇ10（hb)-(3．2*(log1０(11.７5*hm）。^2)－4。9７)＋（（４４．9-6.５5＊loｇ10（ｈb））*ｌog1０(d)）;lb21=lb11—2*（loｇ10(f1/28））。＾2－5．４;lb22=ｌb12-２*(ｌog１0(ｆ2/28)).^2—5.４;lb31=lb11－４。78＊(ｌog10（ｆ1))。^2+18.33*ｌog（f1)－40。98;lb32=ｌb12-4.7８*（lｏg10(f2)）。^2＋1８。３３*log（f2）-4０。9８；ｆ=［f1ｆ2];lb1=［lｂ11lb12］；lb２=［lb２1lb22］；ｌｂ3=[lb31ｌｂ32］；figure（１）;ｈoldoｎ;pｌot（f,lb1,'r'）；tｉtｌe（'大城市’）；xlaｂel（'頻率／MＨz＇）；yｌabel（'損耗中值/ｄB’）；gｒiｄ；ｆiｇｕre（2);ｈoldon；pｌｏｔ(f，lｂ2,’b’)；titｌe（＇郊區(qū)'）；xlabｅｌ（＇頻率/MHz’）;ylabｅl（'損耗中值/dB'）；grｉd；figure（3）；holｄon；plot（f，lｂ3，'g＇）；title（'開闊區(qū)＇）；xｌabｅl（'頻率/MＨz＇）；ｙlaｂｅｌ(’損耗中值/ｄB’)；grｉd；end將計(jì)算結(jié)果和通過奧村模型實(shí)測(cè)測(cè)得的結(jié)果進(jìn)行比較，驗(yàn)證計(jì)算結(jié)果的正確性。對(duì)于奧村模型而言,,其中dB可通過書本曲線簇進(jìn)行數(shù)據(jù)查找。為了驗(yàn)證哈特方法的準(zhǔn)確性，查找數(shù)據(jù)分別代入兩個(gè)公式。取d=10(固定不變），ｆ=50０MHZ,奧村模型=133.9dB，奧村哈特方法=1３5.5dBf=1００0MHZ，奧村模型＝14１．5dB，奧村哈特方法＝１４3.4dＢ;ｆ=1５0０MHZ，奧村模型=147．9ｄB，奧村哈特方法＝１４8dＢ；f＝２000ＭHＺ,奧村模型=１5１。9dB，奧村哈特方法＝１51.２dB;f=２500MHZ,奧村模型=15５．４dＢ，奧村哈特方法=１５3。８dB；f=3０00MＨＺ，奧村模型=1５８。5dＢ，奧村哈特方法=155。8dB。取f＝100０MHZ(固定不變）,d=1KM奧村模型=112.5dB，奧村哈特方法=11３。5dB.d=2KM奧村模型=１21。9dＢ,奧村哈特方法=1２2。5dB。d＝３KＭ奧村模型＝127。0dB,奧村哈特方法＝127。7dB。ｄ＝5KＭ奧村模型=135．４dＢ,奧村哈特方法=１３4.4dB。ｄ=１0KM奧村模型=１42.5dＢ，奧村哈特方法=143．４dB。ｄ=20ＫM奧村模型=151。5dB，奧村哈特方法=152.3dＢ.ｄ=30KM奧村模型=158。8ｄＢ，奧村哈特方法=１57.６dB.ｄ＝40KM奧村模型=1６5.5dＢ，奧村哈特方法=1６１.３dB.d＝5０ＫＭ奧村模型=１７2。4ｄB，奧村哈特方法＝１64。1dB。d＝60KM奧村模型=１7８．０ｄB，奧村哈特方法=166。6dB。d=70KM奧村模型=１82．9dB，奧村哈特方法=168．6dB。對(duì)于其他d的取值和不同的地形，同樣采取取值帶入進(jìn)行比較。通過比較發(fā)現(xiàn)，奧村哈特方法相對(duì)于奧村模型的數(shù)據(jù)誤差并不是很大，基本吻合奧村模型。分析Ｏkumｕｒａ-Haｔａ方法在何距離以及何頻率范圍內(nèi)存在較大的誤差。經(jīng)過粗略的比較可以得出,頻率ｆ在５０0M—2０0０M之間Ｏkumuｒa-Hａｔａ方法存在誤差較小,Ｏkｕmura-Hａt(yī)a方法基本與奧村模型數(shù)據(jù)相吻合.距離d在1KＭ—3０KM之間Ｏkumｕra-Hata方法基本與奧村模型數(shù)據(jù)相吻合,距離超過30KM后，兩者數(shù)據(jù)相差較大。實(shí)驗(yàn)三BPＳＫ調(diào)制解調(diào)仿真【實(shí)驗(yàn)?zāi)康摹考由钜苿?dòng)通信系統(tǒng)中調(diào)制解調(diào)的理解；能夠使用Matlab(或者Ｃ語言）進(jìn)行通信系統(tǒng)調(diào)制解調(diào)基帶仿真，并進(jìn)行性能分析;【實(shí)驗(yàn)內(nèi)容】使用Matlab（或者C語言)仿真BPSＫ調(diào)制解調(diào)過程；比較發(fā)送端星座圖以及不同信噪比下接收端星座圖的不同；仿真AWGN信道下ＢPSK的誤碼率;【實(shí)驗(yàn)設(shè)備】一臺(tái)PC機(jī)【實(shí)驗(yàn)步驟】產(chǎn)生發(fā)送信息比特；對(duì)信息比特進(jìn)行BPSK調(diào)制;將BPSK調(diào)制后的發(fā)送信號(hào)經(jīng)過ＡWＧN信道，從而獲得接收信號(hào)波形;比較發(fā)送端星座圖以及不同信噪比下接收端星座圖的不同；對(duì)接收信號(hào)進(jìn)行ＢＰSK解調(diào)；通過蒙特卡洛方法,仿真系統(tǒng)誤碼率，畫出ＢPSK調(diào)制解調(diào)誤碼率曲線.【實(shí)驗(yàn)報(bào)告】按照要求完成實(shí)驗(yàn)報(bào)告。實(shí)驗(yàn)報(bào)告中要求給出不同信噪比條件下（信噪比分別為—5dB，０dＢ,5dB,１0dB)接收端的星座圖。畫出AWGN信道下ＢPSK調(diào)制解調(diào)誤碼率曲線。【實(shí)驗(yàn)結(jié)果】代碼：cｌc；cleaｒalｌ;n=10００;SNR=［—50５1０］;sｏuｒce=randiｎｔ(１,n）;Iｎ_BＰSK＝source＊２-1;ｆori=1:lｅnｇth（ＳＮR）Y(i，：)=ＡWＧN（Ｉｎ＿ＢPSＫ，SNR(i））；endY_RＥ=Y；Ｙ＿RE（ｆiｎｄ（Ｙ〉０）)＝1；Ｙ＿RE(find（Y<0）)=0;fｏrｉ＝1：ｌｅnｇth(SNR）Y_EＲRO（i，：）=ａｂs（Y_RE(ｉ，:）－source）；Y_ＢＩT（i）=sum(Ｙ_EＲRO(ｉ,：）)／n；endsemiｌogy（SＮＲ,Ｙ_BIT）；gｒidｏｎ；title（'誤碼率’）;figｕｒesｕbplot（４,1，1）；ｐloｔ(Y(１，：），０,’＊＇）;titlｅ（'SNＲ=-5dB＇）；sｕbploｔ（4,1，2)；plot（Y（2,：），0，’＊')；titｌe(’ＳＮR=０ｄB');subｐlot(