【正文】 i)*fc*t)。 //SSB信號上邊帶頻域表達(dá)式//以上是仿真計算部分//以下為繪圖部分//AM信號xset(39。,1)plot(f,abs(S1))title(39。)xlabel(39。)ylabel(39。)mtlb_axis([25,25,0,max(abs(S1))])。window39。AM信號波形39。t39。s(t)39。//DSBSC信號xset(39。,3)plot(f,abs(S2))title(39。)xlabel(39。)ylabel(39。)mtlb_axis([25,25,0,max(abs(S2))])。window39。DSBSC信號波形39。t39。s(t)39。//SSB信號（以上邊帶為例）xset(39。,5)plot(f,abs(S3))title(39。)xlabel(39。)ylabel(39。)mtlb_axis([25,25,0,max(abs(S3))])xset(39。,6)plot(t,s3)title(39。)xlabel(39。)ylabel(39。)mtlb_axis([0,6,3,3]) 產(chǎn)生的波形圖：1） AM信號：2） AM信號頻譜：3） DSBSC信號：4）DSBSC信號的頻譜：5） SSB信號波形：6）SSB信號頻譜： 實驗心得：做這個實驗時的我的理論知識已經(jīng)準(zhǔn)備得比較充分了，所以難點主要在編程方面。這個實驗給了我一種比課本上的理論更深的體會，因為這里的信號頻率有多個，當(dāng)畫出信號波形的時候，我觀察到了直觀的多頻信號形成的調(diào)制信號，這比書本上的公式來的深刻得多。由于之前研究調(diào)幅都是在單頻信號的情況下，而這里是多頻率，所以這對我對調(diào)幅的理解十分有幫助。FM的頻率偏移常數(shù)為5kHz/V。exec 。xdel(winsid())。 //采樣速率T=16。 //采樣點數(shù)dt=1/fs。 //時域采樣點df=1/T。 //頻域采樣點數(shù)fm1=1。fm3=。Kf=5。phi1=(2*%pi)*Kf*cumsum(m1)*dt。phi2=(2*%pi)*Kf*cumsum(m2)*dt。phi3=(2*%pi)*Kf*cumsum(m3)*dt。 //FM信號時域表達(dá)式S=t2f(s,fs)。 //調(diào)制信號表達(dá)式M=t2f(m,fs)。window39。FM信號頻譜39。f39。S(f)39。xset(39。,2)plot(t,s)title(39。)xlabel(39。)ylabel(39。)mtlb_axis([0,3,2,2])。window39。m(t)信號頻譜39。f39。S(f)39。xset(39。,4)plot(t,m)title(39。)xlabel(39。)ylabel(39。)mtlb_axis([8,8,6,8])。記得剛開始畫出信號波形圖的時候我都以為自己弄錯了，因為之前很少觀察符合頻率波的波形，所以一下子難以接受，為什么它是這樣子的呢？在一番觀察驗證之后，我的觀念很快就調(diào)整過來，這對我是一個啟發(fā)，因為我們所研究的真實信號一定是多頻率的，所以當(dāng)我們使用單頻信號研究透了某個調(diào)制解調(diào)原理之后，一定要使用多頻信號再去仿真研究一下，因為這時候所做的工作才是更加符合真實情況的。實驗三通過仿真測量占空比為25%、50%、75%以及100%的單雙極性歸零碼波形及其功率譜。exec 。L=32。 //總采樣點數(shù)M=N/L。 //碼元速率Ts=1/Rb。 //采樣速率T=N/fs。 //系統(tǒng)帶寬t=T/2+[0:N1]/fs。 //頻域采樣點EP1=zeros(1,N)。L0=input(39。)for loop=1:1000 tmp1=zeros(L,M)。 a=(rand(1,M))+0。 //雙極性序列 L1=L*L0。 s1=tmp1(:)39。 //樣本信號的功率譜密度 P1=abs(S1).^2/T。 //雙極性 tmp2([1:L1],:)=ones(L1,1)*b。 S2=t2f(s2,fs)。 //隨機過程的功率譜是各個樣本功率譜的數(shù)學(xué)期望 EP2=EP2*(11/loop)+P2/loop。單極性碼時域圖39。t39。S(t)39。xset(window,2) plot(f,10*log10(EP1+%eps))set(gca(),grid,[1,1]) title(39。)xlabel(39。)ylabel(39。) xset(window,3)plot(t,s2)set(gca(),grid,[1,1]) title(39。)xlabel(39。)ylabel(39。)mtlb_axis([3,3,])。雙極性碼功率譜圖(dB)39。f39。功率39。2） 占空比為25%雙極性歸零碼波形及其功率譜。3） 占空比為50%單極性歸零碼波形及其功率譜。5） 占空比為75%單極性歸零碼波形及其功率譜。當(dāng)信號占空比增大，可以發(fā)現(xiàn)對應(yīng)的信號功率譜變得更加密集。8） 占空比為100%雙極性歸零碼波形及其功率譜。 實驗心得：這次實驗比較簡單，雖然數(shù)據(jù)較多，但是原理相同。實驗四利用編程實現(xiàn)數(shù)字基帶傳輸系統(tǒng)（）1）畫出發(fā)送端輸入碼序列波形和功率譜、發(fā)送濾波器輸出波形和功率譜2）畫出接收端采樣判決后碼序列波形和功率譜、接收濾波器輸出波形和功率譜3）畫出接收濾波器輸出信號眼圖（在升余弦滾將系數(shù)分別為0、1的情況下）4）分別畫出升余弦滾將系數(shù)為0、1，采樣判決點在眼圖最大處的系統(tǒng)的實際誤碼曲線（Pe~s/n0曲線），并在同坐標(biāo)系中畫出理論誤碼曲線5）改變采樣點重復(fù)1）~4）。exec 。 tt=0:1:Na*L。 endset(gca(),auto_clear,on)endfunctionk=input(39。)。N=2^k。M=N/L。 //碼率2MbpsTs=1/Rs。 //時域采樣間隔fs=1/dt。 //頻域采樣間隔T=N*dt。 //系統(tǒng)帶寬t=[T/2+dt/2:dt:T/2]。 //頻域橫坐標(biāo)alpha=input(39。)。 // 升余弦濾波器算法i1=find(abs(f)(1alpha)/(2*Ts)amp。 Hcos(i1)=Ts/2