數(shù)字基帶傳輸系統(tǒng)的matlab仿真

《數(shù)字基帶傳輸系統(tǒng)的matlab仿真》由會員上傳分享，免費在線閱讀，更多相關內容在教育資源-天天文庫。

1、姓名：范浩如學號：107551400879數(shù)字通信作業(yè)：二進制調制最佳接收機的matlab仿真第一章基礎理論知識第一節(jié)知識框架圖第二節(jié)數(shù)字基帶信號傳輸模型第一步:對原始數(shù)字序列進行調制典型的調制方法第二步:調制信號在信道中的傳輸?shù)谌剑赫{制型號的解調（相干解調法）方法一：針對2ASK,2PSK信號進行接收與解調判決門限：2ASK0.5A2PSK0方法二：針對2FSK信號進行接收與解調判決門限:比較兩路信號，誰大選誰?；蛘邔陕沸盘栂鄿p，然后以0為判決門限。大于0時為1（W1代表的），小于0時為0（W0代表的）第二章Matlab程序與仿真：第

2、一節(jié)理想信噪比-誤碼率曲線的產(chǎn)生N=10;k0=2;Th1=0;A=1;%設置基本參量fordb=-5:0.2:15%設置信噪比范圍（取101個點）r=10^(db/10);%對數(shù)化后的信噪比ri0=round((db+6)*5-4);%數(shù)組序號必須從1開始pe0ASKa(i0)=0.5*erfc(sqrt(r/2));%理想條件下Ask的誤碼率(a)非相干pe0ASKb(i0)=0.5*erfc(sqrt(r/4));%理想條件下Ask的誤碼率（b）相干pe0FSK(i0)=pe0ASKa(i0);%理想條件下Fsk的誤碼率pe0PSK(

3、i0)=0.5*erfc(sqrt(r));%理想條件下psk的誤碼率i=1:0.2:21;i1=i-6;%設置i1的范圍從-5到15endsemilogy(i1,pe0PSK,'b',i1，pe0ASKb,'b',i1,pe0FSK,'b');axis([-5151.0e-061.0e0])第二節(jié)2Ask信號的調制解調Matlab程序：Test=[0,1,0,0,0,0,1,1,1,0,0,1,1,1,1,0,1,1,1,1,0,1,1,1,0,0,1,1,1,1,1,0,1,1,1,0,1,1,1,1,1,0,0,0,1,1,1,1,0

4、,1,]fordb=-5:0.2:15sigma=0.5*A*sqrt(N)/10^(db/20);%設置噪聲功率Th=0.5*0.5*A^2*N;%2ASK的判決門限count=0;forj=1:50;%取50個樣本fori=1:Nrand('state',sum(100*clock));%將時鐘作為產(chǎn)生隨機信號的種子NN(i)=sigma*randn(1,1);%產(chǎn)生隨機高斯噪聲s0(i)=0;%輸入端0碼元波形s1(i)=A*sin(2*pi*k0*i/N);%輸入端1碼元波形x0(i)=s0(i)+NN(i);%接收端0碼元波形x1

5、(i)=s1(i)+NN(i);%接收端1碼元波形endsumm0=0;iftest(j)==0fori=1:Nsumm0=summ0+x0(i)*s1(i);endelseiftest(j)==1fori=1:Nsumm0=summ0+x1(i)*s1(i);endendsumma=summ0;Th=0.5*0.5*A^2*N;%2ASK的判決門限ifsumma>Th%和門限值進行比較out(i0,j)=1elseout(i0,j)=0endiftest(j)==out(i0,j)elsecount=count+1;endendpeASK

6、(i0)=count/L;End測試結果如下：原始輸入序列test1-20:0100001110011110111121-30:0111001111101110111141-50:1000111101信噪比db=-5時，out(1,:)錯誤碼元數(shù)：171-20:0101011100100000001021-30:1101001111011110101141-50:1101111111Ps：傳錯的碼加下劃線信噪比db=5時，out(51,:)錯誤碼元數(shù)：11-20:0100101110011110111121-30:0111001111101

7、110111141-50:1000111101信噪比db=15時，out(101,:)錯誤碼元數(shù)：01-20:0100001110011110111121-30:0111001111101110111141-50:1000111101測試集大小為1000時的信噪比-誤碼率曲線第三節(jié)2psk信號的調制解調程序基本上和上面一個相同，只是輸入端碼元改變?yōu)椋簊0(i)=A*sin(2*pi*k0*i/N);%輸入端0碼元S1(i)=-A*sin(2*pi*k0*i/N);%輸入端1碼元門限值改變?yōu)椋篢h=0;測試結果如下：原始輸入序列test1-2

8、0:0100001110011110111121-30:0111001111101110111141-50:1000111101信噪比db=-5時，out(1,:)錯誤碼元數(shù)：8