【正文】 %時域采樣間隔 df=1/(N*dt) %頻域采樣間隔 T=N*dt %時域截短時間 Bs=N*df/2 %系統(tǒng)帶寬 Na=4。 t=[T/2+dt/2:dt:T/2]。 %頻域橫坐標 g1=sin(pi*t/Ts)./(pi*t/Ts)。 g=g1.*g2。 %升余弦脈沖波形的傅氏正變換 figure(1) %set(1,39。,[100,50,300,200]) %設(shè)定窗口位置及大小 %set(2,39。,[350,50,300,200])。t in us39。s(t) in V39。 %加eps是為了躲零 for ii=1:2 % 設(shè)累計平均的循環(huán) a=sign(randn(1,M))。 %產(chǎn)生沖激序列 imp(L/2:L:N)=a/dt。 %升余弦信號的傅氏變換 s=f2t(t2f(imp,df).*G,df)。 P=S.*conj(S)/T。 %累計平均 figure(1) plot(f,30+10*log10(EP),39。)。f (MHz)39。Ps(f) (dBm/MHz)39。 %設(shè)定畫眼圖的坐標的位置 for jj=1:Na*L:NNa*L %Na*L是每次畫眼圖的間隔步長 ，NNa*L是畫眼圖的 %最后一點的起始點。 %畫升余弦信號的眼圖 end end % 對應(yīng)累計平均的循環(huán)實現(xiàn)數(shù)字基帶傳輸系統(tǒng)實驗?zāi)康模豪镁幊虒崿F(xiàn)數(shù)字基帶傳輸系統(tǒng)1）畫出發(fā)送端輸入碼序列波形和功率譜、發(fā)送濾波器輸出波形和功率譜2）畫出接收端采樣判決后碼序列波形和功率譜、接收濾波器輸出波形和功率譜3）畫出接收濾波器輸出信號眼圖（在升余弦滾將系數(shù)分別為0、1的情況下）4）分別畫出升余弦滾將系數(shù)為0、1，采樣判決點在眼圖最大處的系統(tǒng)的實際誤碼曲線（Pe~s/n0曲線），并在同坐標系中畫出理論誤碼曲線5）改變采樣點重復(fù)1）~4）。exec 。k=13。//采樣點數(shù)L=16。//碼元數(shù)Rb=2。 //碼元間隔dt=Ts/L。//采樣速率df=1/(N*dt)。 //截短時間Bs=N*df/2。 //頻域橫坐標t=[T/2+dt/2:dt:T/2]。 //升余弦滾降系數(shù)Hcos=zeros(1,N)。abs(f)=(1+alpha)/(2*Ts))。ii=find(abs(f)=(1alpha)/(2*Ts))。GT=sqrt(Hcos)。 //最佳系統(tǒng)的發(fā)送接收濾波器的傅式變換EPC=zeros(1,N)。EP1=zeros(1,N)。RECT=ones(L,1)。 Eb_N0(l1)=(l11)。 Eb=1。 //信道的噪聲譜密度 sita=n0*Bs。 //誤碼計數(shù) for loop=1:100 a=sign(rand(1,M,normal))。 seq_send=tmp1(:)39。 P1=abs(S1).^2/T。//隨機過程的功率譜是各個樣本的功率譜的數(shù)學(xué)期望 EP11=20+10*log10(EP1+%eps)。//產(chǎn)生沖激序列 sp(1:L:N)=a/dt。 SH=SP.*GT。 //通過發(fā)送濾波器后的信號?。? PB=abs(SH).^2/T。//求平均 EPBB=20+10*log10(EPB+%eps)。 //信道噪聲 SY=SH.*GR。 r=sy+nw。 PC=abs(R).^2/T。//收端 EPCC=20+10*log10(EPC+%eps)。//*取樣?。。? sdp=sign(sd)。 seq_reci=tmp2(:)39。 P2=abs(S2).^2/T。 EP22=20+10*log10(EP2+%eps)。//錯誤累計endPe(l1)=n_err/(M*loop)。g39。 //Pe~Eb/N0曲線畫圖 xlabel(39。)。Pe39。title(Pe~Eb/N0曲線)。 set(gca(),auto_clear,off) k = log10(*erfc(sqrt(eb_n0)))。 xlabel(39。) ylabel(39。) legend(39。,39。)。tt=0:1:Na*L。end set(gca(),auto_clear,on)。 通信原理軟件實驗的整個實驗周期很長，所以我們有充裕的時間去思考自己下一個實驗的內(nèi)容和實現(xiàn)方式。因為可以用自己的筆記本，感覺通原軟件實驗比硬件實驗的機動性更強，而且軟件實驗中避免了很多噪聲和電路上的因素，生成的圖像更加直觀，利于思考和理解。顯然，MatlaB的用途遠遠大于在信號處理上面。31 / 31