【正文】 在這段時間里，我又學到了新的知識。其次，我還要感謝這四年來給予我辛勤教育的老師們，他們不僅教會了我扎實的 專業(yè)知識，還教會了我如何為人處事，在此對他們致以最誠摯的謝意！沒有他們的幫助和教導，就沒有今天的這份論文，再次為他們的付出表示由衷的感謝。 另外，我還要感謝在論文撰寫過程中給我?guī)椭耐瑢W們，特別是和我同宿舍的她們，她們給予了我很大的幫助和支持，是她們無私的付出才成就了我今天的成功，謝謝！ 最后，我要感謝我的家人，是他們始終如一地支持著我，他們是我學習和生活中的摯友，他們是我人生中最寶貴的財富，這份彌足珍貴的情誼將會永遠深藏在我的心中。在此，我對他們表示真摯的感謝，謝謝你們！ 參考文獻 [1] 樊昌信 等，《通信原理》 國防工業(yè)出版社 2020年 [2] 張輝 等，《現(xiàn)代通信原理與技術》 西安電子科技大學出版社 2020年 [3] 曹志剛 等，《現(xiàn)代通信原理》 淸華大學出版社 2020年 [4] 達新宇 等，《通信原理教程》 北京郵電大學出版社 2020年 [5] 李妍 等，《 MATLAB 通信仿真開發(fā)手冊》 國防工業(yè)出版社 2020年 [6] 樓天順 等，《 MATLAB 》 西安電子科技大學出版社 2020年 [7] 趙謙 等，《通信系統(tǒng) 中 MATLAB基礎與仿真應用》 西安電子科技大學出版社 2020年 [8] 徐利民 等，《基于 MATLAB的信號與系統(tǒng)實驗教程》 清華大學出版社 2020年 [9] 別志松，別紅霞 . 系統(tǒng)與通信系統(tǒng)仿真 . 北京：北京郵電大學出版社， 2020. [10] 鄧華 .MATLAB 通信仿真及應用實例詳解 .北京：人民郵電出版社， 2020 [11] MQAM 調(diào)制技術及其在移動通信中的應用 . 移動通信雜志 . 2020 年第一期 附錄 I 將輸入的序列擴展成間隔為 N1個 0的序列 function [out]=sigexpand(d,M) N=length(d)。 out=zeros(M,N)。 out(1,:)=d。 out=reshape(out,1,M*N)。 計算信號的傅里葉變化 function[f,sf]=T2F(t,st)。 dt=t(2)t(1)。 T=t(end)。 df=1/T。 N=length(st)。 f=N/2*df:df:N/2*dfdf。 sf=fft(st)。 sf=T/N*fftshift(sf)。 主函數(shù)代碼 fc=2。 % 載波頻率 N_sample=8。 % 基帶碼元抽樣點數(shù) N=8。 % 碼元數(shù) Ts=1。 % 碼元寬度 dt=Ts/fc/N_sample。 % 抽樣時間間隔 T=N*Ts。 % 信號持續(xù)時間長度 t=0:dt:Tdt。 % 時間向量 Lt=length(t)。 % 時間向量長度 tx1=0。 % 時域波形圖橫坐標起點 tx2=8。 %時域波形圖橫坐標終點 ty1=。 %時域波形圖縱坐標起點 ty2=。 %時域波形圖縱坐標終點 fx1=10。 %功率譜圖橫坐標起點 fx2=10。 %功率譜圖橫坐標終點 fy1=40。 %功率譜圖縱坐標起點 fy2=25。 %功率譜圖縱坐標終點 %產(chǎn)生二進制信源 d=[1,1,1,1,1,1,1,1] dd=sigexpand(d,fc*N_sample)。 % 雙極性 gt=ones(1,fc*N_sample)。 % NRZ波形 d_NRZ=conv(dd,gt)。 % 基帶信號 figure(1)。 subplot(2,3,1)。 plot(t,d_NRZ(1:Lt))。 axis([tx1,tx2,ty1,ty2])。 xlabel(39。時間 (S)39。)。 ylabel(39。幅度 39。)。 title(39。基帶信號時域波形圖 39。)。 grid。 figure(5)。 subplot(2,2,1)。 plot(t,d_NRZ(1:Lt))。 axis([tx1,tx2,ty1,ty2])。 xlabel(39。時間 (S)39。)。 ylabel(39。幅度 39。)。 title(39?；鶐盘枙r域波形圖 39。)。 grid。 [f1,d_NRZf]=T2F(t,d_NRZ(1:Lt))。 figure(5)。 subplot(2,2,2)。 plot(f1,10*log10(abs(d_NRZf).^2/T))。 axis([fx1,fx2,fy1,fy2])。 xlabel(39。頻率 (Hz)39。)。 ylabel(39。功率譜密度 (dB/Hz)39。)。 title(39?；鶐盘柟β首V圖 39。)。 grid。 % 串并轉(zhuǎn)換 d1=[]。 d2=[]。 for i=1:N/2 d1(i)=d(2*(i1)+1)。 d2(i)=d(2*(i1)+2)。 end % 上支路 dd1=sigexpand(d1,2*fc*N_sample)。 gt1=ones(1,2*fc*N_sample)。 d_NRZ1=conv(dd1,gt1)。 figure(1)。 subplot(2,3,2)。 plot(t,d_NRZ1(1:Lt))。 axis([tx1,tx2,ty1,ty2])。 xlabel(39。時間 (S)39。)。 ylabel(39。幅度 39。)。 title(39。串并轉(zhuǎn)換后上支路信號時域波形圖 39。)。 grid。 % 下支路 dd2=sigexpand(d2,2*fc*N_sample)。 d_NRZ2=conv(dd2,gt1)。 figure(1)。 subplot(2,3,5)。 plot(t,d_NRZ2(1:Lt))。 axis([tx1,tx2,ty1,ty2])。 xlabel(39。時間 (S)39。)。 ylabel(39。幅度 39。)。 title(39。串并轉(zhuǎn)換后下支路信號時域波形圖 39。)。 grid。 % 載波 h1t=cos(2*pi*fc*t)。 h2t=sin(2*pi*fc*t)。 figure(1)。 subplot(2,3,4)。 plot(t,h1t)。 axis([tx1,tx2,ty1,ty2])。 xlabel(39。時間 (S)39。)。 ylabel(39。幅度 39。)。 title(39。載波信號時域波形圖 39。)。 grid。 % 上下支路 24電平轉(zhuǎn)換 d11=[]。 % 上支路 d22=[]。 % 下支路 for m=1:N/4。 d11(m)=2*d1(2*m1)+d1(2*m)。 d22(m)=2*d2(2*m1)+d2(2*m)。 end dd11=sigexpand(d11,4*fc*N_sample)。 % 上支路 gt2=ones(1,4*fc*N_sample)。 d_NRZ11=conv(dd11,gt2)。 figure(1)。 subplot(2,3,3)。 plot(t,d_NRZ11(1:Lt))。 axis([tx1,tx2,ty1,ty2])。 xlabel(39。時間 (S)39。)。 ylabel(39。幅度 39。)。 title(39。24電平轉(zhuǎn)換后上支路信號時域波形圖 39。)。 grid。 figure(3)。 subplot(2,2,1)。 plot(t,d_NRZ11(1:Lt))。 axis([tx1,tx2,ty1,ty2])。 xlabel(39。時間 (S)39。)。 ylabel(39。幅度 39。)。 title(39。24轉(zhuǎn)換后上支路信號時域波形圖 39。)。 grid。 dd22=sigexpand(d22,4*fc*N_sample)。 % 下支路 d_NRZ22=conv(dd22,gt2)。 figure(1)。 subplot(2,3,6)。 plot(t,d_NRZ22(1:Lt))。 axis([tx1,tx2,ty1,ty2])。 xlabel(39。時間 (S)39。)。 ylabel(39。幅度 39。)。 title(39。24轉(zhuǎn)換后下支路信號時域波形圖 39。)。 grid。 figure(3)。 subplot(2,2,3)。 plot(t,d_NRZ22(1:Lt))。 axis([tx1,tx2,ty1,ty2])。 xlabel(39。時間 (S)39。)。 ylabel(39。幅度 39。)。 title(39。24轉(zhuǎn)換后下支路信號時域波形圖 39。)。 grid。 % 對上下支路信號進行平滑 pd_NRZ11=d_NRZ11(1:Lt)。 pd_NRZ22=d_NRZ22(1:Lt)。 tao=3/16。 tm=0:dt:8*Ts*taodt。 for w=2:N/4 phxh1=d11(w1)+(d11(w)d11(w1))**(1+sin(pi*tm/8*tao*Ts))。 phxh2=d22(w1)+(d22(w)d22(w1))**(1+sin(pi*tm/8*tao*Ts))。 for k=1:8*tao*fc*N_sample pd_NRZ11((w1)*64+k)=phxh1(k)。 pd_NRZ22((w1)*64+k)=phxh2(k)。 end end % 生成 16QAM信號 s_16qam1=pd_NRZ11 .* h1t。 %生成上支路頻帶信號 s_16qam2=pd_NRZ22 .* h2t。 %生成下支路頻帶信號 figure(2)。 subplot(2,4,1)。 plot(t,s_16qam1)。 axis([tx1,tx2,ty1,ty2])。 xlabel(39。時間 (S)39。)。 ylabel(39。幅度 39。)。 title(39。上支路已調(diào)信號時域波形圖 39。)。 grid。 [f2,s_16qam1f]=T2F(t,s_16qam1)。 figure(2)。 subplot(2,4,5)。 plot(f2,10*log10(abs(s_16qam1f).^2/T))。 axis([fx1,fx2,fy1,fy2])。 xlabel(39。頻率 (Hz)39。)。 ylabel(39。功率譜密度 (dB/Hz)39。)。 title(39。上支路已調(diào)信號功率譜圖 39。)。 grid。 figure(2)。 subplot(2,4,2)。 plot(t,s_16qam2)。 axis([tx1,tx2,ty1,ty2])。 xlabel(39。時間 (S)39。)。 ylabel(39。幅度 39。)。 title(39。下支路已調(diào)信號時域波形圖 39。)。 grid。 [f3,s_16qam2f]=T2F(t,s_16qam2)。