【正文】 ?？偸且_好多 M— file 文件進行編寫函數(shù)，最麻煩的是寫完一個子函數(shù)在單獨對其進行驗證就存在一些函數(shù)輸入的問題。當時間過了 2 天之后，才拿起設計要求認真看起來，想起老師在課堂上進行 MATLAB 講解時還是很感興趣的，于是按照老師的要求開始完成本次課程 設計。 對通信系統(tǒng)的理解：因為之前有過課程設計的經(jīng)驗，于是剛開始的時候沒有急于編寫代碼，而是仔細認真研究書本上的知識，進行消化和理解。①定性了分析了整個通信系統(tǒng)；②了解了通信系統(tǒng)研究的意義； 通信 系統(tǒng) 課程設計是培養(yǎng) 我們 綜合運用所學知識 ,發(fā)現(xiàn) 、 提出 、 分析和解決實際問題 ,鍛煉實踐能力的重要環(huán)節(jié) ,是對學生實際工作能力的具體訓練和考察過程 。 課程設計的收獲：本次課程設計結(jié)束了，受益良多。它帶來不僅是知識上的收獲，而且精神上的收獲也是無盡的。①本次課程設計是我在大學寥寥可數(shù)的幾次設計性實驗之一，是和同學們一起奮斗的一場 經(jīng)歷?？吹酵瑢W們的努力，看到老師的辛勤，自己也都不敢懈怠了。②實驗開始的那一刻，告誡自己沒有失敗，老師既然出了設計項目，既然已經(jīng)選擇的這個設計項目，就不許中途退出，這不是游戲，我想那一刻要的是勇氣，是決心吧！③實驗成功的那一刻，收獲的是信心，是對自己本學期以來學習通信專業(yè)課的一個肯定，是對以后的學習的一種督促，也是對自己的未來大學校門外的生活的一個鼓勵。④還有就是在設計肯定會有不甚理解的地方，這就需要我們自己動手查找資料，增加的不僅是查找的資料上面的知識，更多的搜索的能力。 編程體會：①在編寫代碼的時候，不 急于剛開始就對著電腦編寫程序，一定要理解和整理出來你所想要做什么；②從前都沒有畫流程圖的習慣，不過此次試驗之前都有仔細分析，畫流程圖，發(fā)現(xiàn)在代碼編寫的過程中這樣不僅效率高，而且成功率還快；③還有就是因為在程序中過多的嵌套過多 for 循環(huán)語句和 if 選擇語句，所以很容易不小心就會導致代碼錯誤，所以程序中應盡量避免。④還有就是一定要仔細，仔細，再仔細。 武漢理工大學《 學科基礎課群綜合訓練 》 報告 23 5 參考文獻 [1]樊昌信，曹麗娜 . 通信原理（第六版） [M],北京：國防工業(yè)出版社， 2021 [2]程相君，陳生潭 信號與系統(tǒng) [M],西安： 西安電子科技 大學出版社 ,1990 [3]曹志剛，錢壓生 現(xiàn)代通信原理 [M],北京： 清華大學出版社， 1992 [4]劉敏，魏玲 MATLAB 通信仿真與應用 [M]，北京：國防工業(yè)出版社，2021 [5]唐向 .岳恒立 .學峰 ATLAB 及在電子信息類課程中的應用 [M]，北京： 電子工業(yè)出版社， 2021 [6]陳杰 MATLAB 寶典 [M]，北京：電子工業(yè)出版社 [7]張德豐 MATLAB/Simulink 建模與仿真 [M]，北京 電子工業(yè)出版社 武漢理工大學《 學科基礎課群綜合訓練 》 報告 24 6 附錄 ts = 1e3。 t = 0:ts:20*ts 。 x = sin( 2*pi*50*t )+*sin( 2*pi*150*t )。 delta = 。 D( 1+length(t) ) = 0。 for k=1:length(t) e(k) = x(k) D(k)。 e_q(k) = delta * ( 2 * ( e(k) = 0 ) 1 )。 D( k+1 ) = e_q(k) + D(k)。 codeout(k) = ( e_q(k)0 )。 end figure(1) subplot(211)。plot( t, x, 39。o39。 )。axis( [0 20*ts, 2 2] )。title(39。隨機信源 39。)。 subplot(212)。stairs( t, codeout )。axis( [0 20*ts, 2 2] )。title(39。增量編碼 39。)。 %AMI 調(diào)制 xn = codeout*1 。 yn = xn。 num = 0。 for k=1:length(yn) if xn(k) == 1 num = num + 1。 if num/2 == fix(num/2) yn(k) = 1。 else yn(k) = 1。 武漢理工大學《 學科基礎課群綜合訓練 》 報告 25 end end end figure(2) stairs( ( 0:length(yn)1 ), yn )。axis([ 0 length(yn),2 2 ])。title(39。AMI39。)。 %漢明碼編碼 HanMing = encode( xn, 7, 4, 39。hamming39。 )。 figure(3) stairs( ( 0 : 41 ), HanMing )。axis([0 41, 0 ])。title(39。漢明編碼 39。)。 %psk g=(HanMing)。 f=2。 t=0::。 cp=[]。sp=[]。 mod=[]。bit=[]。 for n=1:length(g)。 if g(n)==0。 c1=ones(1,100)。 c2=zeros(1,100)。 else g(n)=1。 c1=ones(1,100)。 c2=ones(1,100)。 end c=sin(2*pi*f*t)。 cp=[cp c1]。 mod=[mod c]。 bit=[bit c2]。 end 武漢理工大學《 學科基礎課群綜合訓練 》 報告 26 psk=cp.*mod。 figure(4) subplot(2,1,1)。stairs(bit,39。LineWidth39。,)。grid on。 title(39。Binary Signal39。)。 axis([0 500 ])。 subplot(2,1,2)。plot(psk,39。LineWidth39。,)。grid on。 title(39。PSK modulation39。)。 axis([0 500 ])。 %對 psk 加噪 pskz=awgn(psk,13)。 figure(5) plot(pskz,39。LineWidth39。,)。grid on。title(39。加噪后信號 39。) %psk 解調(diào) reask=pskz.*mod。 figure(6) plot(reask)。title(39。reask39。)。axis([0 500 ])。title(39。PSK 解調(diào)后的信號 39。) b=fir1(60,2/50,39。low39。)。 figure(7)。 freqz(b,1,512)。title(39。濾波器的幅頻特性，相頻特性 39。) at=filter(b,1,reask)。 figure(8) plot(at)。axis([0 500, 0 ])。title(39。psk 解調(diào) 39。)。 武漢理工大學《 學科基礎課群綜合訓練 》 報告 27 hn=[]。 ren=[]。 decision=length(at)/100。 for k=1:decision i=(2*k1)*50。 hn(i)=at(i)max(at)/2。 ren(k)=(hn(i)0)。 end figure(9) stairs(ren)。axis([0 41, 0 ])。title(39。抽樣判決后的圖 39。) %漢明譯碼 HanMing_decode = decode( HanMing,7 ,4, 39。hamming39。 )。 figure(10) stairs( ( 0 : 23 ), HanMing_decode )。axis([0 23, 0 ])。title(39。漢明譯碼 39。)。 %AMI 譯碼 decodeAMI = abs(yn)。 figure(11) stairs( ( 0 : 20 ), decodeAMI )。axis([0 20, 0 ])。title(39。AMI 譯碼 39。)。 %DM 解調(diào) ts = 1e3。 t = 0:ts:20*ts 。 Dr( 1+length(t) ) = 0。 for k=1:length(t) eq(k) = delta * ( 2*decodeAMI(k) 1 )。 xr(k) = eq(k) + Dr(k)。 武漢理工大學《 學科基礎課群綜合訓練 》 報告 28 end figure(12) subplot(2,1,1)。stairs(t, x)。hold on。 subplot(2,1,1)。plot(t, x)。title(39。對比 39。)。 wp=2*pi*2500。ws=2*pi*3500。Rp=1。As=30。 [N,wc]=buttord(wp,ws,Rp,As,39。s39。)。 [B,A]=butter(N,wc,39。s39。)。 fk=0:800/512:8000。wk=2*pi*fk。 Hk=freqs(B,A,wk)。 ts = 1e3。 t = 0:ts:20*ts 。 Dr( 1+length(t) ) = 0。 for k=1:length(t) eq(k) = delta * ( 2*decodeAMI(k) 1 )。 xr(k) = eq(k) + Dr(k)。 end figure(13) subplot(2,1,1)。stairs(t, xr)。hold on。 subplot(2,1,1)。plot(t, x)。title(39。譯碼源碼對比 39。)。 len=length(codeout)。 errornum=0。 for i=1:len if(codeout(i)~=decodeAMI(i)) errornum=errornum+1。 end 武漢理工大學《 學科基礎課群綜合訓練 》 報告 29 end disp(39。誤碼率為 39。)。 errornum/len