【文章內(nèi)容簡介】 ode)/34,34)。%RS解碼code = gf(code,6)。code = rsdec(code,34,24)。code = gf2mat(code,6)。msg = gf2mat(msg,6)。%統(tǒng)計誤碼率[num,rat] = biterr(code(tblen+1:end),msg(1:endtblen))。temp = msg code。 綜合測試代碼及效果圖下面代碼的運行會繪制一個圖形，上面的三條曲線，分別是（1）單純的64QAM調(diào)制解調(diào)（2）64QAM+卷積碼1（3）64QAM+卷積碼2的曲線。其中卷積碼1和卷積碼2都來源于Matlab現(xiàn)成的Demo，其網(wǎng)格想見代碼或者結(jié)構圖。clear。clc。trel1 = poly2trellis(7,[171 133])。 tblen1 = 48。trel2 = poly2trellis([5 4],[23 35 0。0 5 13])。tblen2 = 34。SNR = 1::8。for i = 1:length(SNR) [num0(i),rat0(i)] = QAM64(SNR(i))。 [num1(i),rat1(i)] = Conv_AND_QAM64(trel1,tblen1,SNR(i))。 %rat = 1/2 [num2(i),rat2(i)] = Conv_AND_QAM64(trel2,tblen2,SNR(i))。 %rat = 2/3endfigure。plot(SNR,rat0,39。r39。,SNR,rat1,39。go39。,SNR,rat2,39。b*39。)。hold on。grid on。legend(39。QAM64 (No Conv)39。,39。QAM64+Conv(2,1,7)39。,39。QAM64+Conv(3,2,5)39。,0)。仿真結(jié)果如下，可以看的出來，64QAM的誤碼率最高，64QAM+卷積碼（3,2,5）次之；64QAM+卷積碼（2,1,7）再次之；下面代碼的運行會繪制一個圖形，上面的三條曲線，分別是（1）單純的64QAM調(diào)制解調(diào)（2）64QAM+RS碼1（3）64QAM+ RS碼2的曲線。%okclc。clear。SNR = 1::8。for i = 1:length(SNR) [Num1(i),Rat1(i)] = QAM64(SNR(i))。 [Num2(i),Rat2(i)] = RS_AND_QAM64(36,24,SNR(i))。 [Num3(i),Rat3(i)] = RS_AND_QAM64(36,20,SNR(i))。endfigure。plot(SNR,Rat1,39。r39。,SNR,Rat2,39。go39。,SNR,Rat3,39。b*39。)。hold on。grid on。legend(39。QAM64 (No RS)39。,39。QAM64+RS(36,24)39。,39。QAM64+RS(36,20)39。,0)。下面代碼的運行會繪制一個圖形，上面的四條曲線，分別是（1）單純的64QAM調(diào)制解調(diào)（2）64QAM+RS碼（3）64QAM+ 卷積碼的曲線。（4）64QAM+ 卷積碼+交織碼+RS碼clc。clear。%有關卷積碼部分trel = poly2trellis(7,[171 133])。 tblen = 48。SNR = 1::8。for i = 1:length(SNR) [Num1(i),Rat1(i)] = QAM64(SNR(i))。 [Num2(i),Rat2(i)] = RS_AND_QAM64(36,24,SNR(i))。 [Num3(i),Rat3(i)] = Conv_AND_QAM64(trel,tblen,SNR(i))。 [Num4(i),Rat4(i)] = RS_Conv_QAM64(SNR(i))。endfigure。plot(SNR,Rat1,39。r39。,SNR,Rat2,39。go39。,SNR,Rat3,39。b*39。,SNR,Rat4,39。rs39。)。hold on。grid on。legend(39。QAM6439。,39。QAM64+RS(36,24)39。,39。QAM64+Conv(2,1,7)39。,39。QAM64+RS(36,24)+Conv(2,1,7)39。,0)。 小結(jié)+，我建議在授課的時候，在給學生講述一邊大概的原理之后，“照本宣科”，講解的內(nèi)容重點講述不同的卷積碼、RS碼對誤碼率的影響，例如，編碼效率和誤碼率之間的關系。5 頻分復用（FDMA）仿真 題目內(nèi)涵及仿真函數(shù)說明FDMA技術是信息工程、通信工程等專業(yè)比較掌握的一類類似科普性質(zhì)的知識點，如果從信號處理的角度理解FDMA技術，超級的簡單。本題目利用簡單的調(diào)制，把幾路基帶信號調(diào)制到不同的頻率點，然后合成，在接收端，利用不同截止頻率的帶通濾波實現(xiàn)信號的復用。本題目為了使一目了然，用到了一個畫頻譜函數(shù)。 授課教師教學用的Demo 畫頻譜函數(shù)function plotspec(x,Ts)N=length(x)。 % length of the signal xt=Ts*(1:N)。 % define a time vector ssf=(N/2:N/21)/(Ts*N)。 % frequency vectorfx=fft(x(1:N))。 % do DFT/FFTfxs=fftshift(fx)。 % shift it for plottingsubplot(2,1,1), plot(t,x) % plot the waveformxlabel(39。seconds39。)。 ylabel(39。amplitude39。) % label the axessubplot(2,1,2), plot(ssf,abs(fxs)) % plot magnitude spectrumxlabel(39。frequency39。)。 ylabel(39。magnitude39。) % label the axes FDMA仿真函數(shù)clear。clc。T = 1/20000。t = 0:T:4999*T。p = hamming(100)。carrier1 = cos(2*pi*500*t)。carrier2 = cos(2*pi*2500*t)。carrier3 = cos(2*pi*4500*t)。carrier4 = cos(2*pi*6500*t)。b = fir1(200,39。low39。)。b1 = fir1(200,39。low39。)。b2 = fir1(200,[ ],39。band39。)。b3 = fir1(200,[ ],39。band39。)。b4 = fir1(200,[ ],39。band39。)。msg1 = randint(1,50,[0,1])。index = find(msg1==0)。msg1(index)=1。msg2 = randint(1,50,[0,1])。index = find(msg2==0)。msg2(index)=1。msg3 = randint(1,50,[0,1])。index = find(msg3==0)。msg3(index)=1。msg4 = randint(1,50,[0,1])。index = find(msg4==0)。msg4(index)=1。base1 = []。base2 = []。base3 = []。base4 = []。for i = 1:50 base1 = [base1 msg1(i)*p39。]。 base2 = [base2 msg2(i)*p39。]。 base3 = [base3 msg3(i)*p39。]。 base4 = [base4 msg4(i)*p39。]。endfigure。plotspec(base1,T)。figure。plotspec(base2,T)。figure。plotspec(base3,T)。figure。plotspec(base4,T)。 base1 = base1.*carrier1。base2 = base2.*carrier2。base3 = base3.*carrier3。base4 = base4.*carrier4。send = base1+base2+base3+base4。figure。plotspec(base1,T)。figure。plotspec(base2,T)。figure。plotspec(base3,T)。figure。plotspec(base4,T)。figure。plotspec(send,T)。receive = send。user1 = filter(b1,1,receive)。user2 = filter(b2,1,receive)。user3 = filter(b3,1,receive)。user4 = filter(b4,1,receive)。figure。plotspec(user1,T)。figure。plotspec(user2,T)。figure。plotspec(user3,T)。figure。plotspec(user4,T)。xt1 = user1.*carrier1。xt2 = user2.*carrier2。xt3 = user3.*carrier3。xt4 = user4.*carrier4。new1 = filter(b,1,xt1)。new2 = filter(b,1,xt2)。new3 = filter(b,1,xt3)。new4 = filter(b,1,xt4)。figure。plotspec(new1,T)。figure。plotspec(new1,T)。figure。plotspec(new1,T)。figure。plotspec(new1,T)。new1 = n