基于Matlab的CDMA通信完整系統(tǒng)分析及仿真

1、推南呼絶爭(zhēng)院課程論文題目：基于MatlabCDMA多址技術(shù)的仿真學(xué)生姓名：蘇未然學(xué)生學(xué)號(hào)：1008030130系另I：電氣信息工程學(xué)院專業(yè)：電子信息工程年級(jí)：10級(jí)指導(dǎo)教師：王麗電氣信息工程學(xué)院制2013年4月i目錄緒論3矚第1章CDMA通信系統(tǒng)分析及仿真41.1 整體仿真框圖4聞1.2 信源5殘1.3 偽隨機(jī)序列生成器61.4擴(kuò)頻7釅1.5編碼和調(diào)制8彈1.6接收端9謀第2章仿真系統(tǒng)10廈礴懇2.1 信源102.2 編碼102.3 擴(kuò)頻112.4 調(diào)制與解調(diào)112.5 誤碼計(jì)算11第3章實(shí)驗(yàn)結(jié)果與總結(jié)12鐃3.2多用戶在相同信道環(huán)境下的仿真12參考文獻(xiàn)13仿真結(jié)果14設(shè)計(jì)的代碼16緒論利用M

2、ATLAB平臺(tái)的SIMULINK可視化仿真功能，結(jié)合CDMA的實(shí)際通信情況，利用MATLAB組建出完整的CDMA通信系統(tǒng)，完成整體設(shè)計(jì)方案，實(shí)現(xiàn)完整的發(fā)送到接收的端到端的CDMA無線通信系統(tǒng)的建模、仿真和分析，介紹了CDMA的主要環(huán)節(jié)（包括擴(kuò)頻技術(shù)、差錯(cuò)控制技術(shù)、調(diào)制技術(shù)、信道等）的參數(shù)設(shè)置。單用戶在不同信道環(huán)境下，信道的噪聲對(duì)結(jié)果影響很大，噪聲功率越大，系統(tǒng)的誤碼率越大。多用戶在相同信道環(huán)境下的仿真時(shí)，由于是采用了噪聲功率為0.01W的信道傳輸環(huán)境，所以在單用戶情況下誤碼率為0，而增加了用戶數(shù)之后，誤碼率也隨之增加?？梢?，信號(hào)在傳輸過程中，除了受到信道噪聲的影響外，還存在多址接入干擾、單頻干

3、擾、窄帶干擾、跟蹤式干擾等，主要是多址接入干擾。20世紀(jì)60年代以來，隨著民用通信事業(yè)的發(fā)展，頻帶擁擠問題日益突出。CDMA（CodeDivisionMultipleAccess，碼分多址）通信，在使用相同頻率資源的情況下，理論上CDMA移動(dòng)網(wǎng)比模擬網(wǎng)容量大20倍，實(shí)際使用中比模擬網(wǎng)大10倍，比GSM要大45倍，所以在通信領(lǐng)域中起著非常重要的作用。CDMA的基本原理是利用互相正交（或盡可能正交）的不同編碼，分配給不同用戶調(diào)制信號(hào)，實(shí)現(xiàn)多用戶同時(shí)使用同一頻率接入系統(tǒng)和網(wǎng)絡(luò)的通信。由于利用互相正交（或盡可能正交）的編碼去調(diào)制信號(hào)，會(huì)將原信號(hào)的頻譜帶寬擴(kuò)展，因此，這種通信方式，又稱為擴(kuò)頻通信。本論文

4、所完成的CDMA通信仿真系統(tǒng)，是結(jié)合CDMA的實(shí)際通信情況，利用MATLAB組建出完整的CDMA通信系統(tǒng)，完成整體設(shè)計(jì)方案，實(shí)現(xiàn)完整的發(fā)送到接收的端到端的CDMA無線通信系統(tǒng)的建模、仿真和分析。買鯛鴯譖曇膚遙閆擷凄第1章CDMA通信系統(tǒng)分析及仿真1.1整體仿真框圖本論文在CDMA通信原理的基礎(chǔ)上，得出CDMA通信系統(tǒng)的仿真框圖（圖2-1）。擴(kuò)頻編碼和調(diào)制JTJUU標(biāo)量化U調(diào)編碼-制矢量化解擴(kuò)解調(diào)和譯碼圖1-1CDMA通信系統(tǒng)的仿真框圖1.2信源二進(jìn)制貝努利序列產(chǎn)生器產(chǎn)生一個(gè)二進(jìn)制序列，并且這個(gè)二進(jìn)制序列中的0和1服從貝努利分布。本文使用4個(gè)二進(jìn)制貝努利信號(hào)發(fā)生器，以子系統(tǒng)形式封裝于輸入信號(hào)中。

5、產(chǎn)生器的產(chǎn)生是由一個(gè)隨機(jī)信號(hào)器與一個(gè)常數(shù)進(jìn)行判決，輸出的二進(jìn)制再進(jìn)行抽樣整形，從而輸出符合參數(shù)設(shè)置的而進(jìn)駐。輸入信號(hào)抽樣的時(shí)間均為1，即碼元寬度為1，選擇產(chǎn)生一維向量。綾鏑鯛駕櫬鶘蹤韋轔糴。1.3偽隨機(jī)序列生成器擴(kuò)頻通信系統(tǒng)中，偽隨機(jī)序列與正交編碼是十分重要的技術(shù)。主要包括m序列，Gold序列，Walsh碼序列等。Walsh碼序列比較復(fù)雜，正交性較好，主要用于CDMAIS-95系統(tǒng)中。而Gold序列可以比m序列產(chǎn)生更多的地址嗎，更適合于大型的通信系統(tǒng)。在本設(shè)計(jì)中，主要是對(duì)簡(jiǎn)單的CDMA系統(tǒng)進(jìn)行仿真，所以選用m序列作為擴(kuò)頻序列，而且有4個(gè)用戶。4個(gè)m序列分別的4級(jí)，5級(jí)，6級(jí)和7級(jí)，周期分別為1

6、5,31，63和127。擴(kuò)頻序列發(fā)生器的主要參數(shù)為生成多項(xiàng)式，試驗(yàn)采用的數(shù)值分別為：11001、110001、1100001、10101011。抽樣時(shí)間設(shè)置為0.1，即碼元寬度為0.1。驅(qū)躓髏彥浹綏譎飴憂錦1.4擴(kuò)頻本文是采用直接序列擴(kuò)頻方式實(shí)現(xiàn)多址接入。在仿真中，將原信號(hào)與偽隨機(jī)序列相乘，從而實(shí)現(xiàn)擴(kuò)頻。但由于輸入信號(hào)和m序列都是單極性的二進(jìn)制數(shù)，所以在進(jìn)入乘法器進(jìn)行擴(kuò)頻之前，還要對(duì)它們進(jìn)行單/雙變換，變成雙極性信號(hào)。圖22分別給出了原信號(hào)波形、擴(kuò)頻序列波形和擴(kuò)頻后的信號(hào)波形。本系統(tǒng)的擴(kuò)頻倍數(shù)為10貓蠆驢繪燈鮒誅髏貺廡a.原信號(hào)波形b.擴(kuò)頻序列波形圖12直接序列擴(kuò)頻方式1.5編碼和調(diào)制1.5.

7、1BCH編碼仿真框圖如23所示。模型采用（7,4）BCH碼，要求送入編碼器的是維數(shù)為4的矢量，編碼器的輸出是維數(shù)為7的矢量，即為每個(gè)信息組添加了3位校驗(yàn)碼元,由圖24得知，只進(jìn)行差錯(cuò)控制編碼，而沒有經(jīng)過擴(kuò)頻的信號(hào)，在給定的高斯信道中傳輸，隨著碼源傳輸?shù)臅r(shí)間增加，誤碼率會(huì)比較高。鍬籟饗逕瑣筆襖鷗婭薔貝努利用列矢肩化悅塊1BCH編碼標(biāo)晴匕模塊1信道圖13BCH碼的仿真框圖0.200.160J20.080.041020304050時(shí)問1/H圖14BCH碼的誤碼率曲線誤碼率計(jì)算公式10Xlog10（power_signal/power_noise）TheAWGNChannel（高斯白噪聲信道）模塊可以

8、在輸入信號(hào)中加入實(shí)信號(hào)噪聲或復(fù)合信號(hào)噪聲。當(dāng)輸入信號(hào)是實(shí)信號(hào)時(shí)，此模塊在輸入信號(hào)中加入實(shí)高斯白噪聲,并輸出實(shí)信號(hào)。當(dāng)輸入信號(hào)是復(fù)合信號(hào)時(shí)，此模塊在輸入信號(hào)中加入復(fù)合的高斯白噪聲，并輸出符合信號(hào)。次模塊從輸入信號(hào)中得到抽樣時(shí)間。這個(gè)模塊用DSPBlockset（數(shù)字信號(hào)處理）模塊中的RandomSource（隨機(jī)信號(hào)源）模塊來產(chǎn)生噪聲。TheInitialseed（初始種子）參數(shù)用來初始化信號(hào)發(fā)生器。TheInitialseed（初始種子）即可以是標(biāo)量也可以是矢量。這個(gè)標(biāo)量或矢量的長(zhǎng)度要與信道匹配。構(gòu)氽頑黌碩飩薺齦話騖1.5 M-PSK仿真圖25給出的是M=16時(shí)M-PSK的仿真框圖，信號(hào)調(diào)制后

9、的頻譜和相位星座圖分別如圖26和圖27所示。本文中4個(gè)調(diào)制器的相數(shù)M分別為16,32,32,40。由星座圖可以得知，將每個(gè)輸入信號(hào)都對(duì)應(yīng)于一個(gè)點(diǎn)，點(diǎn)與點(diǎn)之間的相位差為360°/16=22.5°。1.5信道信道中的噪聲直接影響著信號(hào)的傳輸質(zhì)量。根據(jù)信道中噪聲的特點(diǎn)，可將信道劃分為：加性高斯白噪聲信道（AdditiveWhiteGuassionNoise,AWGN）、二進(jìn)制對(duì)稱信道、多徑瑞利衰落信道和倫琴衰落信道等,我們?cè)诜抡嬷胁捎玫氖茿WGN信道。圖28中的正弦波功率設(shè)置為1W,曲線表明,當(dāng)SNR是-20dB時(shí)，噪聲功率是100W；SNR是20dB時(shí)，噪聲功率是0.01W?？?/p>

10、見，為了得到比較優(yōu)的性能，應(yīng)合理設(shè)置輸入信噪比。輒嶧陽檉籪癤網(wǎng)儂號(hào)澩I町機(jī)序列nrj-ItHiiLlamInteger倍道M-PSK斛調(diào)-1_i_IWXZAM-PSX4"aWGNn十E傾講i卄斤M-FSKlzFirnrRflte起CfllddbiGUli溟碼率il畀J星建分析圖15MPSK仿真模型圖apm_-譽(yù)100-10-20-30101頻率Hz27圖16MPSK信號(hào)頻譜圖-1.5-1-0.50(L?冋相分量1二103102Kio°101圖17MPSK信號(hào)星座圖10'230-20-1001020信噪比/品圖18SNR與噪聲功率的關(guān)系1.6.接收端M-PSK解調(diào)器的

11、參數(shù)設(shè)置與M-PSK調(diào)制器相同，譯碼器的參數(shù)設(shè)置與編碼器相對(duì)應(yīng)。解擴(kuò)過程要求使用的偽隨機(jī)碼與發(fā)送端擴(kuò)頻用的偽隨機(jī)碼不僅碼字相同，而且相位相同。多用戶情況下，由于信號(hào)互相疊加，所以解擴(kuò)后的信號(hào)不再是二進(jìn)制信號(hào)，因?yàn)榇嬖诙嘤脩舾蓴_，所以在進(jìn)入接收端進(jìn)行誤碼統(tǒng)計(jì)之前，必須經(jīng)過濾波和判決。為了減少噪聲影響，在解擴(kuò)之后加入低通濾波器。根據(jù)發(fā)送信號(hào)的頻譜，該濾波器的分子系數(shù)和分母系數(shù)分別設(shè)置為0.00040.00170.00250.00170.0004和1.0000-3.18113.8623-2.11300.4385,初始條件為0。濾波器幅頻、相頻特性曲線如圖29所示。發(fā)送信號(hào)與接收機(jī)恢復(fù)出的信號(hào)同時(shí)送入

12、誤碼儀(Error-ratemeter)模塊進(jìn)行比較，誤碼儀可計(jì)算和顯示誤碼率和誤比特率。堯側(cè)閆繭絳闕絢勵(lì)蜆贅0.20.40.60.81.0500-50-1500-100歸_化頻率a幅頻特性b相頻特性圖1一9濾波器特性第2章仿真系統(tǒng)本系統(tǒng)主要包括信號(hào)源輸入、擴(kuò)頻序列發(fā)生器、擴(kuò)頻解擴(kuò)、編碼譯碼、調(diào)制解調(diào)、信道、接收判決、誤碼計(jì)算等。為了簡(jiǎn)化模型，使用了子系統(tǒng)將其中的一些模塊進(jìn)行了封裝。設(shè)計(jì)中，使用了遞加的方法，可以實(shí)現(xiàn)多個(gè)不同用戶按不同情況接入，圖31為整體仿真框圖。識(shí)饒鎂錕縊灩筧嚌儼淒圖31整體仿真框圖2.1信源本文使用4個(gè)二進(jìn)制貝努利信號(hào)發(fā)生器，以子系統(tǒng)形式封裝于輸入信號(hào)中。產(chǎn)生器的產(chǎn)生是由

13、一個(gè)隨機(jī)信號(hào)器與一個(gè)常數(shù)進(jìn)行判決，輸出的二進(jìn)制再進(jìn)行抽樣整形，從而輸出符合參數(shù)設(shè)置的而進(jìn)駐。輸入信號(hào)抽樣的時(shí)間均為1,即碼元寬度為1,選擇產(chǎn)生一維向量。凍鈹鋨勞臘鍇癇婦脛糴。2.2編碼編碼主要包括m序列，Gold序列，Walsh碼序列等。Walsh碼序列比較復(fù)雜，正交性較好，主要用于CDMAIS-95系統(tǒng)中。而Gold序列可以比m序列產(chǎn)生更多的地址嗎，更適合于大型的通信系統(tǒng)。在本設(shè)計(jì)中，主要是對(duì)簡(jiǎn)單的CDMA系統(tǒng)進(jìn)行仿真，所以選用m序列作試驗(yàn)采用的數(shù)值分別為：11001、110001、1100001、10101011。抽樣時(shí)間設(shè)置為0.1，即碼元寬度為0.1。2.3 擴(kuò)頻本文是采用直接序列擴(kuò)頻

14、方式實(shí)現(xiàn)多址接入。在仿真中，將原信號(hào)與偽隨機(jī)序列相乘，從而實(shí)現(xiàn)擴(kuò)頻。但由于輸入信號(hào)和m序列都是單極性的二進(jìn)制數(shù)，所以在進(jìn)入乘法器進(jìn)行擴(kuò)頻之前，還要對(duì)它們進(jìn)行單/雙變換，變成雙極性信號(hào)。圖22分別給出了原信號(hào)波形、擴(kuò)頻序列波形和擴(kuò)頻后的信號(hào)波形。本系統(tǒng)的擴(kuò)頻倍數(shù)為10。鯊腎鑰詘褳鉀溈懼統(tǒng)庫。2.4 調(diào)制與解調(diào)本文采用MPLK調(diào)制與解調(diào)。2.5 誤碼計(jì)算濾波器幅頻、相頻特性曲線如圖29所示。發(fā)送信號(hào)與接收機(jī)恢復(fù)出的信號(hào)同時(shí)送入誤碼儀(Error-ratemeter)模塊進(jìn)行比較，誤碼儀可計(jì)算和顯示誤碼率和誤比特率。碩癘鄴頏謅攆檸攜驤蘞。第3章實(shí)驗(yàn)結(jié)果與總結(jié)3.1單用戶在不同信道環(huán)境下的仿真由表1可

15、以分析出，在信道傳輸過程中，信道的噪聲對(duì)結(jié)果影響很大，在噪聲功率為100W的情下，誤碼率偏高。降低噪聲功率后，誤碼率也明顯減小?？梢姡趩斡脩羟闆r下，CDMA通信系統(tǒng)的碼率主要取決于信道中的噪聲。閿擻輳嬪諫遷擇楨秘3.2多用戶在相同信道環(huán)境下的仿真仿真條件：用戶數(shù)從1到4，由于本系統(tǒng)采用了遞加的方法，所以可以在同一個(gè)系統(tǒng)中觀察到從1個(gè)用戶一直到4個(gè)用戶的四種不同情況下的仿真。由表2的仿真結(jié)果可以看出，由于是采用了噪聲功率為0.01W的信道傳輸環(huán)境，所以在單用戶情況下誤碼率為0，而增加了用戶數(shù)之后，誤碼率也隨之增加。可見，信號(hào)在傳輸過程中，除了受到信道噪聲的影響外，還存在多址接入干擾、單頻干擾、

16、窄帶干擾、跟蹤式干擾等，主要是多址接入干擾。氬嚕躑竄貿(mào)懇彈瀘頷澩。表31單用戶在不同信道環(huán)境下的仿真用戶數(shù)傳送信息量/bit誤碼數(shù)噪聲功率/W抗干擾方式誤碼率11000444100直擴(kuò)0.4441100034510直擴(kuò)0.345110001611直擴(kuò)0.16111000230.1直擴(kuò)0.0231100000.01直擴(kuò)0表32多用戶在相同信道環(huán)境下的仿真用戶數(shù)傳送信息量/bit誤碼數(shù)噪聲功率/W抗干擾方式誤碼率1100000.01直擴(kuò)0110001830.01直擴(kuò)0.183110001760.01直擴(kuò)0.176110002230.01直擴(kuò)0.223參考文獻(xiàn)1 劉學(xué)勇,編著.通信系統(tǒng)建模與仿真（電

17、子工業(yè)出版社）2 鄧薇，編著.MATALB函數(shù)速查手冊(cè)（人民郵電出版社）3 WilliamC.Y.Lee,著.移動(dòng)通信工程理論和應(yīng)用（第二版）（人民郵電出版社）4 王華奎，李艷萍等編著.移動(dòng)通信原理與技術(shù)（清華大學(xué)出版社）仿真結(jié)果cdKa2000_phlayer/FroMChannel二口Frequency(MHz)HP經(jīng)過信道后的信道頻譜信號(hào)經(jīng)過信道后的星座圖可以看出信號(hào)經(jīng)過信道后，星座點(diǎn)的位置變得雜亂無章，如果不進(jìn)行信道補(bǔ)償，則系統(tǒng)的誤碼率會(huì)大大增加。補(bǔ)償信道后的星座圖可以看出，經(jīng)過信道補(bǔ)償后星座點(diǎn)的位置得到進(jìn)一步集中。設(shè)計(jì)的代碼1. M序列functionmout=mseq(n,taps

18、,inidata,num)%*%n:m序列的階數(shù)n%taps:反饋寄存器的連接位置%inidata:寄存器的初始值序列%num:輸出的m序列的個(gè)數(shù)%mout:輸出的m序列，如果num>1，則每一行為一個(gè)m序列ifnargin<4num=1;endmout=zeros(num,2入n-1);fpos=zeros(n,1);fpos(taps)=1;forii=1:2入n-1mout(1,ii)=inidata(n);%寄存器的輸出值慫temp=mod(inidata*fpos,2);%計(jì)算反饋數(shù)據(jù)諺辭inidata(2:)=inidata(1:n-1);%寄存器移位一次嘰inidat

19、a(1)的值熒紿譏鉦鏌觶鷹緇=temp;%更新第1個(gè)寄存器endifnum>1生成其他m鶼漬螻偉閱ifshiftr>0%如果要輸出多個(gè)m序列，outregi(:,1:shiftr)=inregi(:,v-shiftr+1:v);%循環(huán)移位outregi(:,1+shiftr:v)=inregi(:,1:v-shiftr);elseifshiftr<0outregi(:,1:v+shiftr)=inregi(:,1-shiftr:v);outregi(:,v+shiftr+1:v)=inregi(:,1:-shiftr);end*2、Gold序列：functiongout=go

20、ldseq(m1,m2,num)%*%ml:m序列1%m2:m序列2%num:生成的Gold序列個(gè)數(shù)%gout:生成的Gold序列輸出%*ifnargin<3%如果沒有指定生成的Gold序列個(gè)數(shù)，默認(rèn)為1num=1;endgout=zeros(num,length(m1);forii=1:num%根據(jù)Gold序列生成方法生成Gold序列g(shù)out(ii,:)=xor(m1,m2);m2=shift(m2,1);end%*endoffile其shift函數(shù)代碼同M序列3、仿真部分主程序：%直接序列擴(kuò)頻主程序代碼functionber=dscdma(user,seq)%user:同時(shí)進(jìn)行擴(kuò)頻通

21、信的用戶數(shù)%seq:擴(kuò)頻碼1:M-序列2:Golc序列3:正交Gold序列%ber：該用戶數(shù)下的誤碼率%符號(hào)速率%*初始化部分*sr=25600.0;nSymbol=1000;M=4;br=sr*log2(M);graycode=0132;EbNo=0:2:10;%每種信噪比下發(fā)送的符號(hào)數(shù)%4-QAM調(diào)制裊樣祕(mì)%比特速率倉%Gray編碼規(guī)則綻%Eb/No變化范圍%*脈沖成形濾波器參數(shù)*delay=10;%升余弦濾波器時(shí)延Fs=8;%濾波器過采樣數(shù)櫛rolloff=0.5;%升余弦濾波器滾降因子rrcfilter=rcosine(1,Fs,'fir/sqrt',rolloff,d

22、elay);%設(shè)計(jì)根升余弦濾波器%*擴(kuò)頻碼產(chǎn)生參數(shù)*詩叁撻訥燼憂毀%user=user1;用戶數(shù)%stage=3;%m序列的階數(shù)ptap1=13;%m序列1的寄存器連接方式ptap2=23;%m序列2的寄存器連接方式regi1=111;%m序列1的寄存器初始值陽regi2=111;%m序列2的寄存器初始值%*擴(kuò)頻碼的生成switchseqcase1%M-序列鋇code=mseq(stage,ptap1,regi1,user);case2%Gold序列m1=mseq(stage,ptap1,regi1);m2=mseq(stage,ptap2,regi2);code=goldseq(m1,m2,

23、user);case3%正交Gold序列m1=mseq(stage,ptap1,regi1);m2=mseq(stage,ptap2,regi2);code=goldseq(m1,m2,user),zeros(user,1);頭次頭頭頭頭頭頭夫頭頭頭頭頭夫咲次頭頭咲頭endcode=code*2-1;clen=length(code);%*衰落信道參數(shù)ts=1/Fs/sr/clen;%信道采樣時(shí)間間隔t=(0:nSymbol*Fs*clen-1+2*delay*Fs)*ts;%每種信噪比下的符號(hào)傳輸時(shí)間%fd=160;多普勒頻%移Hz%h=rayleigh(fd,t);仿真開始%*kkic&#

24、39;kkkk'k'kkk':kk'kk'k'kkkk'k'：'kkkk'kforindx=1:length(EbNo)indx%產(chǎn)生各個(gè)用戶的發(fā)射數(shù)據(jù)%Gray編碼%4-QAM調(diào)制%擴(kuò)頻%*發(fā)射端*data=randsrc(user,nSymbol,0:3);data1=graycode(data+1);data1=qammod(data1,M);out=spread(data1,code);out1=rcosflt(out.',sr,Fs*sr,'filter',rrcfilter);%

25、通過脈沖成形濾波器spow=sum(abs(out1).入2)/nSymbol;%計(jì)算每個(gè)用戶信號(hào)功率ifuser>1%用戶數(shù)大于1時(shí)，所有用戶數(shù)據(jù)相加out1=sum(out1.');elseout1=out1.'endsigma=sqrt(0.5*spow*sr/br*10入(-EbNo(indx)/10);緦徑%根據(jù)信噪比計(jì)算高斯白噪聲方差y=;forii=1:usery(ii,:)=out1+sigma(ii).*(randn(1,length(out1)+j*randn(1,length(out1);%加入高斯白噪聲(AWGN)騅憑鈳銘僥張礫陣軫藹。endy=r

26、cosflt(y.',sr,Fs*sr,'Fs/filter',rrcfilter);%通過脈沖成形濾波器進(jìn)行濾波y=downsample(y,Fs);%降采樣forii=1:usery1(:,ii)=y(2*delay+1:end-2*delay,ii);end%4-QAM解調(diào)yd=despread(y1.',code);%數(shù)據(jù)解擴(kuò)鏃demodata=qamdemod(yd,M);demodata=graycode(demodata+1);%Gray編碼逆映射err,ber(indx)=biterr(data,demodata,log2(M);%統(tǒng)計(jì)誤比特率e

27、nd主程序中spread函數(shù)的代碼：%擴(kuò)頻函數(shù)functionout=spread(data,code)switchnargincase0,1%如果輸入?yún)?shù)個(gè)數(shù)不對(duì)，提示錯(cuò)誤error'缺少輸入?yún)?shù)');endhn,vn=size(data);hc,vc=size(code);ifhn>hc%如果擴(kuò)頻碼數(shù)小于輸入的待擴(kuò)頻的數(shù)據(jù)序列，提示錯(cuò)誤error'缺少擴(kuò)頻碼序列');endout=zeros(hn,vn*vc);forii=1:hnout(ii,:)=reshape(code(ii,:).'*data(ii,:),1,vn*vc)臠龍訛驄椏業(yè)變墊,end解擴(kuò)函數(shù)despread的代碼functionout=despread(data,switchnargincode)%如果輸入?yún)?shù)個(gè)數(shù)不對(duì)，提示錯(cuò)誤%信號(hào)解擴(kuò)case0,1error'缺少輸入?yún)?shù)');endhn,vn=size(data);hc,vc=size(code);out=zeros(hc,vn