【正文】 (receive22(index_min:index_max))。subplot(247)。plot(receive23(index_min:index_max))。subplot(248)。plot(receive24(index_min:index_max))。仿真效果圖如下：上圖中，上面4個(gè)小信號(hào)時(shí)發(fā)射機(jī)合并信號(hào)前的4路信號(hào)，下面4個(gè)信號(hào)時(shí)接收機(jī)接收信號(hào)后分解的4個(gè)信號(hào)，從圖中可以看出，這個(gè)信號(hào)“整體趨勢(shì)”是一樣的，不同的是延時(shí)。這個(gè)題目還用到了濾波器等知識(shí)，和題目“信號(hào)的分解（分析）和綜合（合成）”相關(guān)內(nèi)容一樣，在這里不在枚舉。 小結(jié)本題目用到得知識(shí)點(diǎn)如下：濾波器設(shè)計(jì)及應(yīng)用，特別是半帶濾波器，差值和抽取相關(guān)理論及其對(duì)對(duì)信號(hào)頻譜的改變。8 正交頻分復(fù)用（OFDM）仿真 題目?jī)?nèi)涵及說(shuō)明通信系統(tǒng)從信息工程專業(yè)的角度看，在理論方面只有信息處理和信號(hào)處理兩個(gè)主要內(nèi)容，本題目主要是解決信號(hào)處理的。OFDM技術(shù)是一種很有競(jìng)爭(zhēng)力的技術(shù)，其競(jìng)爭(zhēng)力在工程實(shí)踐中也有體現(xiàn)，例如，歐洲數(shù)字電視廣播系統(tǒng)、無(wú)線局域網(wǎng)等。本課題要求學(xué)生在Matlab環(huán)境中，完成OFDM技術(shù)的仿真，加深對(duì)其原理、頻譜等的理解。 授課教師教學(xué)用的Demo OFDM仿真（一）代碼具體如下：%OFDM仿真，但是，這里不調(diào)用FFT函數(shù)和IFFT函數(shù)，而是采用抽樣定理function [errs,ber] = ofdm_con1(carrier_num,snr)msg = randint(1,carrier_num,[1,8])。T = 1。t = 0:1/500:T。%發(fā)送端形成發(fā)射信號(hào)send = zeros(size(t))。for i=1:carrier_num carrier(i,:) = sin(2*pi*i*1/T*t)。 total(i) = sum(carrier(i,:).*carrier(i,:))。 send = send + msg(i)*carrier(i,:)。end%信道作用receive = awgn(send,snr)。%接受端解調(diào)信號(hào)for i = 1:carrier_num temp = sum(receive.*carrier(i,:))。 data(i) = round(temp/total(i))。end%統(tǒng)計(jì)誤碼率[errs,ber] = biterr(msg,data)。 OFDM仿真（二）%OFDM仿真，這里調(diào)用FFT函數(shù)和IFFT函數(shù)function [errs,ber] = ofdm_con2(carrier_num,snr)msg = randint(1,carrier_num,[1,8])。%發(fā)送端形成發(fā)射信號(hào)send = ifft(msg,carrier_num)。%信道作用receive = awgn(send,snr)。%接受端解調(diào)信號(hào)receive = fft(receive)。receive = round(receive)。data = abs(real(receive))。%統(tǒng)計(jì)誤碼率[errs,ber] = biterr(msg,data) 綜合測(cè)試代碼及效果圖下面的代碼完成第一種OFDM的誤碼率曲線的繪制；clear。clc。snr = 0:2:18。num = 64。for i = 1:length(snr) [errs(i),ber(i)] = ofdm_con1(num,snr(i))。endplot(snr,ber)。下面的代碼完成第二種OFDM的誤碼率曲線的繪制；clear。clc。snr = 0:2:30。num = 2048。for i = 1:length(snr) [errs(i),ber(i)] = ofdm_con2(num,snr(i))。endplot(snr,ber)。 小結(jié)OFDM技術(shù)是一種很有競(jìng)爭(zhēng)力的技術(shù)，其競(jìng)爭(zhēng)力在工程實(shí)踐中也有體現(xiàn)，例如，歐洲數(shù)字電視廣播系統(tǒng)、無(wú)線局域網(wǎng)等。本課題要求學(xué)生在Matlab環(huán)境中，完成OFDM技術(shù)的仿真，加深對(duì)其原理、頻譜等的理解。9 小波函數(shù)、尺度函數(shù)以及頻譜 小波和尺度函數(shù)在做課件或者文章的過(guò)程中，常常要得到一個(gè)小波函數(shù)對(duì)于的小波和尺度函數(shù)，具體代碼如下：[a,d,x] = wavefun(39。db239。)。subplot(211)。plot(x,a)。subplot(212)。plot(x,d)?？梢愿淖僿avefun函數(shù)的參數(shù)來(lái)得到不同的小波的母函數(shù)和尺度函數(shù)。 小波和尺度函數(shù)的頻譜function plot_wavelet_psd(wavename,iter)%這個(gè)函數(shù)本質(zhì)就是完成一個(gè)正交小波的尺度函數(shù)和小波函數(shù)的時(shí)域和頻域的顯示。%下面的小波是正交小波，可以調(diào)用這個(gè)函數(shù)%haar%db1......db45%meyr%sym2 ........%coif1.......coif5%上面的db序列小波、sym序列小波、coif序列小波都是正交的、有限緊支小波，而meyr小波則是無(wú)限緊支的小波。wn = wavename。[s,w,x] = wavefun(wn,iter)。subplot(221)。plot(x,s)。hold on。title(39。尺度函數(shù)時(shí)域39。)。subplot(222)。plot(x,w)。hold on。title(39。小波函數(shù)時(shí)域39。)。T = (max(x)min(x))/length(x)。Fs = floor(1/T)/2+1。fft_length = 2^(floor(log2(length(x)))+1)。delt_f = Fs/fft_length。f = 0:delt_f:delt_f*(fft_length/21)。ffts = fft(s,fft_length)。fftw = fft(w,fft_length)。ffts = ffts(1:length(ffts)/2)。fftw = fftw(1:length(fftw)/2)。psd_s = abs(ffts.^2)。psd_w = abs(fftw.^2)。%psd_s = log10(psd_s)。%psd_w = log10(psd_w)。plot_length = floor(1/delt_f)+1。subplot(223)。plot(f(1:plot_length),psd_s(1:plot_length))。hold on。title(39。尺度函數(shù)頻域39。)。subplot(224)。plot(f(1:plot_length),psd_w(1:plot_length))。hold on。title(39。小波函數(shù)頻域39。)。10 小波用于信號(hào)分解、合成 最簡(jiǎn)單的例子clear。load leleccum。a0 = leleccum(1:3920)。[a1,d1] = dwt(a0,39。db139。)。[a2,d2] = dwt(a1,39。db139。)。[a3,d3] = dwt(a2,39。db139。)。[a4,d4] = dwt(a3,39。db139。)。subplot(511)。plot(a0)。subplot(512)。plot(a1)。subplot(513)。plot(a2)。subplot(514)。plot(a3)。subplot(515)。plot(a4)。clear。load leleccum。ls = 3920。a = zeros(4,3920)。a(1,:) = leleccum(1:ls)。subplot(411)。plot(leleccum(1:ls))。for i=2:4 ls = ls/2。 [a(i,1:ls),d(i,1:ls)] = dwt(a(i1,1:2*ls),39。db139。)。 temp = 2^(i+12)。 %temp = 2。 subplot(4,temp,(i1)*temp+1)。plot(a(i,1:ls))。 subplot(4,temp,(i1)*temp+2)。plot(d(i,1:ls))。end 信號(hào)小波分解后合成（保持不變）clc。clear。Fs = 1000。T = 1/Fs。dec_length = 1024。t = 0:T:T*(dec_length1)。s = sin(2*pi*100*t)+sin(2*pi*150*t)+sin(2*pi*200*t)+sin(2*pi*250*t)+sin(2*pi*300*t)。[dec_s,L]= wavedec(s,4,39。db339。)。rec_s = waverec(dec_s,L,39。db339。)。plot_length = 200。subplot(211)。plot(t(1:plot_length),s(1:plot_length))。subplot(212)。plot(t(1:plot_length),rec_s(1:plot_length))。 信號(hào)小波合成后分解（保持不變）clc。clear。Fs = 1000。T = 1/Fs。dec_length = 1041。t = 0:T:T*(dec_length1)。s = sin(2*pi*100*t)+sin(2*pi*150*t)+sin(2*pi*200*t)+sin(2*pi*250*t)+sin(2*pi*300*t)。my_L = [68 68 132 259 514 1024]。rec_s = waverec(s,my_L,39。db339。)。[dec_s,L]= wavedec(rec_s,4,39。db339。)。plot_length = 200。subplot(211)。plot(t(1:plot_length),s(1:plot_length))。title(39。原始信號(hào)39。)。subplot(212)。plot(t(1:plot_length),dec_s(1:plot_length))。title(39。小波合成后又分解的信號(hào)39。)。