ofdm技術(shù)仿真matlab代碼-資料下載頁(yè)
2025-07-13 20:39本頁(yè)面
【正文】 )。pwelch(uoft,[],[],[],Rs)。%Upconverters_tilde=(uoft.39。).*exp(1i*2*pi*fc*t)。s=real(s_tilde)。 %通頻帶信號(hào) (E)figure(9)。plot(t(80:480),s(80:480))。figure(10)。subplot(211)。%plot(ff,abs(fft(((real(uoft).39。).*cos(2*pi*fc*t)),q*FS))/FS)。%plot(ff,abs(fft(((imag(uoft).39。).*sin(2*pi*fc*t)),q*FS))/FS)。plot(ff,abs(fft(s,q*FS))/FS)。subplot(212)。%pwelch(((real(uoft).39。).*cos(2*pi*fc*t)),[],[],[],Rs)。%pwelch(((imag(uoft).39。).*sin(2*pi*fc*t)),[],[],[],Rs)。pwelch(s,[],[],[],Rs)。附錄2：OFDM接收部分代碼%DVB 2K模式接收部分clear all。close all。Tu=224e6。 %有用OFDM符號(hào)持續(xù)時(shí)間T=Tu/2048。 %原始基帶周期G=0。 %允許保護(hù)時(shí)間間隔選擇1/4, 1/8, 1/16或1/32delta=G*Tu。 %保護(hù)間隔持續(xù)時(shí)間Ts=delta+Tu。 %整個(gè)OFDM符號(hào)持續(xù)時(shí)間Kmax=1705。 %子載波數(shù)Kmin=0。FS=4096。 %IFFT/FFT 長(zhǎng)度q=10。 %載波周期與原始基帶周期比f(wàn)c=q*1/T。 %載波頻率Rs=4*fc。 %模擬周期t=0:1/Rs:Tu。tt=0:T/2:Tu。%數(shù)據(jù)產(chǎn)生程序sM = 2。[x,y] = meshgrid((sM+1):2:(sM1),(sM+1):2:(sM1))。alphabet = x(:) + 1i*y(:)。N=Kmax+1。rand(39。state39。,0)。a=1+2*round(rand(N,1)).39。+i*(1+2*round(rand(N,1))).39。A=length(a)。info=zeros(FS,1)。info(1:(A/2)) = [ a(1:(A/2)).39。]。info((FS((A/2)1)):FS) = [ a(((A/2)+1):A).39。]。carriers=FS.*ifft(info,FS)。%UpconverterL = length(carriers)。chips = [ carriers.39。zeros((2*q)1,L)]。p=1/Rs:1/Rs:T/2。g=ones(length(p),1)。dummy=conv(g,chips(:))。u=[dummy。 zeros(46,1)]。[b,aa] = butter(13,1/20)。uoft = filter(b,aa,u)。delay=64。 %接收端重構(gòu)濾波器延遲s_tilde=(uoft(delay+(1:length(t))).39。).*exp(1i*2*pi*fc*t)。s=real(s_tilde)。%OFDM 接收%Downconversionr_tilde=exp(1i*2*pi*fc*t).*s。 %(F)figure(1)。subplot(211)。plot(t,real(r_tilde))。axis([0e7 12e7 60 60])。grid on。figure(1)。subplot(212)。plot(t,imag(r_tilde))。axis([0e7 12e7 100 150])。grid on。figure(2)。ff=(Rs)*(1:(q*FS))/(q*FS)。subplot(211)。plot(ff,abs(fft(r_tilde,q*FS))/FS)。grid on。figure(2)。subplot(212)。pwelch(r_tilde,[],[],[],Rs)。%載波抑制[B,AA] = butter(3,1/2)。r_info=2*filter(B,AA,r_tilde)。 %基帶信號(hào)持續(xù)時(shí)間(G)figure(3)。subplot(211)。plot(t,real(r_info))。axis([0 12e7 60 60])。grid on。figure(3)。subplot(212)。plot(t,imag(r_info))。axis([0 12e7 100 150])。grid on。figure(4)。f=(2/T)*(1:(FS))/(FS)。subplot(211)。plot(ff,abs(fft(r_info,q*FS))/FS)。grid on。subplot(212)。pwelch(r_info,[],[],[],Rs)。%抽樣r_data=real(r_info(1:(2*q):length(t)))... %離散時(shí)間基帶信號(hào)+1i*imag(r_info(1:(2*q):length(t)))。 % (H)figure(5) figure(5)。subplot(212)。stem(tt(1:20),(imag(r_data(1:20))))。axis([0 12e7 100 150])。grid on。figure(6)。f=(2/T)*(1:(FS))/(FS)。subplot(211)。plot(f,abs(fft(r_data,FS))/FS)。grid on。subplot(212)。pwelch(r_data,[],[],[],2/T)。%FFTinfo_2N=(1/FS).*fft(r_data,FS)。 % (I)info_h=[info_2N(1:A/2) info_2N((FS((A/2)1)):FS)]。%限幅for k=1:N,a_hat(k)=alphabet((info_h(k)alphabet)==min(info_h(k)alphabet))。 %jend。figure(7)plot(info_h((1:A)),39。.k39。)。title(39。infoh Received Constellation39。)axis square。axis equal。figure(8)plot(a_hat((1:A)),39。or39。)。title(39。a_hat 4QAM39。)axis square。axis equal。grid on。axis([ ])。subplot(211)。stem(tt(1:20),(real(r_data(1:20))))。axis([0 12e7 60 60])。3
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1