課程設(shè)計(jì)oqpsk通信系統(tǒng)的matlab仿真分析報(bào)告(參考版)
2025-08-05 22:09本頁面
【正文】 拼一個(gè)春夏秋冬！贏一個(gè)無悔人生！早安！—————獻(xiàn)給所有努力的人.學(xué)習(xí)參考。不奮斗就是每天都很容易，可一年一年越來越難。是狼就要練好牙，是羊就要練好腿。Matlab中有很多自帶的函數(shù)，但變量數(shù)據(jù)的格式是有規(guī)定的，如biterr（）規(guī)定變量只能是0或1，如果原來的1，1作為變量的話就會(huì)報(bào)錯(cuò)。在原來所學(xué)的基礎(chǔ)上，用matlab對(duì)oqpsk的調(diào)制和解調(diào)進(jìn)行編程，從調(diào)制解調(diào)以及加入噪聲都是自己一步一步去寫程序或上網(wǎng)查資料，比如如何產(chǎn)生二進(jìn)制序列，如何延時(shí)以及去除延時(shí)，在不斷的摸索當(dāng)中，學(xué)會(huì)了好多函數(shù)的應(yīng)用以及特定的功能實(shí)現(xiàn)，例如reshape（），biterr(),find(),scatterplot(),eyediagriam(),zeros(),ones()等函數(shù)的應(yīng)用，比如產(chǎn)生特定頻率的基帶信號(hào)可以用經(jīng)過sigexpand函數(shù)擴(kuò)展后的信號(hào)與矩形脈沖相卷得到，更重要的是學(xué)會(huì)了數(shù)據(jù)在matlab中的儲(chǔ)存方式，比如本次試驗(yàn)中輸入序列長度為500，采樣率為32，載波頻率為2，則這個(gè)序列在內(nèi)存中是一個(gè)長度為32000的一位橫矩陣。endQe=qe(1:16)。qe(jj)=rt。qe=[]。備注：如果要畫實(shí)際誤碼率曲線可以在程序最后使用循環(huán)，框架如下m=1。 %頻域的相乘[t,st] = F2T(f,yf)。%規(guī)定低通濾波器的寬度hf(bf) = 1。hf = zeros(1,length(f))。s time samplesdf = f(2)f(1)。% 低通濾波器function [t,st]=lpf(f,sf,B)% f: frequency samples % sf: input data spectrum samples% B: lowpass39。sff=fftshift(sf)。%t=T/2:dt:T/2dt。N=length(sf)。Fmx=(f(end)f(1)+df)。sf=T/N*fftshift(sf)。f=N/2*df:df:N/2*dfdf。df=1/T。%計(jì)算信號(hào)的傅里葉變換%Input is the time and the signal vectors,the length of time must greater%than 2%Output is the frequency and the signal spectrumdt=t(2)t(1)。out=reshape(out,1,M*N)。out=zeros(M,N)。)。title(39。)。 ylabel(39。Es/N0(dB)39。)。 %通過一系列計(jì)算可以得到 % %由于 %可以進(jìn)一步得到，semilogy(EsN0dB,Ps,39。figure[num,rt] = biterr(d_NRZJ,qr_qpsk)。s=find(qr_qpsk0)。r=find(d_NRZ0)。q=find(qr_qpsk0)。p=find(d_NRZJ0)。)。title(39。虛部39。)。xlabel(39。axis([4,4,4,4])。x=d_NRZ1(1:fc*N*N_sample)+j*d_NRZ2DLY(1:fc*N*N_sample)。grid。解調(diào)信號(hào)功率譜圖 Pn=10dB39。)。ylabel(39。頻率(Hz)39。axis([fx1,fx2,fy1,fy2])。subplot(2,2,4)。[f10,qr_qpskf]=T2F(t,r_qpsk(1:Lt))。)。title(39。幅度39。)。xlabel(39。plot(t,qr_qpsk(1:Lt))。qr_qpsk=[r_qpsk(2*N_sample*fc+1:N*N_sample*fc),ones(1,2*N_sample*fc)]%去除延時(shí)figure(4)。endrddd=sigexpand(ddd,fc*N_sample)。for s=1:N/2 ddd(2*s1)=dd22DLY(s)。dd22DLY=d_NRZ11DLY(fc*N_sample:2*fc*N_sample:end)。)。title(39。幅度39。)。xlabel(39。subplot(2,4,8)。 %判決dd444=sigexpand(dd44,2*fc*N_sample)。%抽樣判決dd33=r_qpsk21(fc*N_sample:2*fc*N_sample:end)。)。title(39。幅度39。)。xlabel(39。subplot(2,4,7)。[t2,r_qpsk21]=lpf(f9,r_qpsk2f,B1)。grid。下支路通過乘法器信號(hào)時(shí)域波形圖39。)。ylabel(39。時(shí)間(S)39。axis([tx1,tx2,ty1,ty2])。subplot(2,4,6)。% 通過乘法器2r_qpsk2 = y_qpsk .* h2t。)。title(39。幅度39。)。xlabel(39。plot(t,d_NRZ11DLY(1:Lt))。 %上支路信號(hào)延時(shí)figure(3)。 d_NRZ11=conv(dd222,gt1)。 dd22=sign(dd11)。grid。上支路通過低通濾波器信號(hào)時(shí)域波形圖39。)。ylabel(39。時(shí)間(S)39。plot(t1,r_qpsk11)axis([tx1,tx2,ty1,ty2])。figure(3)。B1=1。grid。上支路通過乘法器信號(hào)時(shí)域波形圖39。)。ylabel(39。時(shí)間(S)39。axis([tx1,tx2,ty1,ty2])。subplot(2,4,2)。% 相干解調(diào)% 通過乘法器1r
外語相關(guān)推薦
鄂ICP備17016276號(hào)-1