.../中國地質大學〔XX通信原理實驗報告姓名：學號：班級：指導老師：陳振興TOC\o"1-3"\h\z\u29735實驗一338331.實驗目的397622.實驗內容3314163.實驗結果4317314.分析與小結421440實驗二5288161.實驗目的5187972.實驗內容574533.實驗結果56774.分析與小結1032247實驗三10243731.實驗目的10320122實驗內容〔一10279732.1實驗結果12215202.2分析與小結13315443實驗內容〔二13116053.1實驗結果15127183.2分析與小結16實驗一1.實驗目的練習使用matlab中的基本函數,為以后的實驗打下基礎2.實驗內容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*df-df;sf=fft<st>;sf=T/N*fftshift<sf>;function[t,st]=F2T<f,sf>%%傅里葉反變換df=f<2>-f<1>;Fmx=<f<end>-f<1>+df>;dt=1/Fmx;N=length<sf>;T=dt*N;t=0:dt:T-dt;sff=fftshift<sf>;st=Fmx*ifft<sff>;clearall;closeall;T=1;N_sample=128;%%去采樣點個數128dt=T/N_sample;t=0:dt:T-dt;st=[ones<1,N_sample/2>,-ones<1,N_sample/2>];%%產生一個周期的方波信號subplot<211>;plot<t,st>;axis<[01-22]>;xlabel<'t'>;ylabel<'s<t>'>;subplot<212>;[f,sf]=T2F<t,st>;plot<f,abs<sf>>;holdon;axis<[-101001]>;xlabel<'f'>;ylabel<'s<f>'>;3.實驗結果4.分析與小結圖一為一個周期的方波信號,圖二為此方波信號經過傅里葉變換變換到頻域后的圖形。由圖可以看出,該方波信號經過傅里葉變換后能量主要集中在-2到2之間,且該頻率圖形對稱。實驗二1.實驗目的用matlab實現2PSK和QPSK,并做它們的仿真。2.實驗內容t=1;forsnr=0:15%%信噪比從0取到15map=[-11];tr_2=1+fix<rand<1,10^5>>;%%產生隨機的0,1信號tx=map<tr_2+1>;rx=awgn<tx,snr,'measured'>;%%通過高斯信道得到rxrx_2=[rx>0];%%抽樣判決pe=length<find<tr_2-rx_2>>;%%比較接收信號與發射信號的不同ps<t>=pe/length<tr_2>;%%計算誤碼率t=t+1;endsemilogy<0:15,ps,'-bo'>;gridon;3.實驗結果2PSK運行結果：QPSK運行結果：QPSK仿真結果：snr=10;星座圖為：snr=5:星座圖為：2PSK仿真結果：星座圖為：4.分析與小結有實驗結果可以看出,實驗結果和仿真出來的結果在相同信噪比的情況下基本相同；QPSK與2PSK相比,在相同信噪比的情況下,QPSK性能較好于2PSK。實驗三1.實驗目的量化編碼,解碼2實驗內容〔一——將給定音頻文件量化編碼functionwaveL=10000;%%選取采樣點100000個；voicefile=wavread<'voicefile.wav'>;%%打開給定音頻文件wavefile并將其作圖figure<1>plot<voicefile>maxvalue=max<abs<voicefile>>;%%求出音頻的最大振幅M=4;k=log2<M>;quan_set=2*maxvalue/M;%%選取量化間隔tx_quanvalue=zeros<1,L>;%%初始化發送和接收數組rx_quanvalue=zeros<1,L>;fori=1:Lsamplevalue=voicefile<i>;%%讀入音頻文件第i個值forn=1:MIfsamplevalue>=-maxvalue+quan_set*<n-1>&&samplevalue<=-maxvalue+quan_set*ntx_quanvalue<i>=<-2*maxvalue+2*quan_set*n-quan_set>/2;%%取區間中點為抽樣值tx2str=dec2bin<n-1,k>;%%轉化為k位二進制foru=1:ktx_b<u,i>=str2double<char<tx2str<u>>>;%%每一行代表一個數據endendendendfigure<2>plot<tx_quanvalue>%%做出輸出波形量化后的圖形tx=tx_b<:>';rx=reshape<tx,k,L>;fori=1:L%%對量化編碼的數據進行解碼n=0;forj=1:kn=n+rx<j,i>*2^<k-j>;rx_quanvalue<i>=<-2*maxvalue+2*quan_set*<n+1>-quan_set>/2;endendwavwrite<rx_quanvalue<:>,48000,'rx_voicefile.wav'>;%%將解碼后的數據保存入新增文件rx_voicefilefigure<3>plot<rx_quanvalue>2.1實驗結果原噪聲圖形量化編碼后的圖形解碼后的圖形2.2分析與小結由以上的圖形可以看出,原來的音頻信號是一串隨機的信號,經過量化編碼后變成圖二所示的8個量化電平的信號,因為沒有加入噪聲,經過解碼后恢復的信號和編碼是一樣的。3實驗內容〔二——將文件rx_voicefile通過2ASK和QPSK編碼,通過高斯信道,然后解碼%%同上將音頻文件量化編碼functionwave2L=10000;voicefile=wavread<'rx_voicefile.wav'>;figure<1>plot<voicefile>maxvalue=max<abs<voicefile>>;M=4;k=log2<M>;quan_set=2*maxvalue/M;fori=1:Lsamplevalue=voicefile<i>;forn=1:Mifsamplevalue>=-maxvalue+quan_set*<n-1>&&samplevalue<=-maxvalue+quan_set*ntx_quanvalue<i>=<-2*maxvalue+2*quan_set*n-quan_set>/2;tx2str=dec2bin<n-1,k>;%%轉化為k位二進制foru=1:ktx_b<u,i>=str2double<char<tx2str<u>>>;%%每一行代表一個數據endendendEnd%%通過高斯信道然后解碼。snr=100;map=[-1,1];fori=1:Lforu=1:ktx<u,i>=map<tx_b<u,i>+1>;endendrx=awgn<tx<:>,snr,'measured'>;%%經過高斯信道fori=1:k*Lrx<i>=rx<i>>0;%%抽樣判決endfori=1:L%%解碼n=0;forj=1:kn=n+rx<<i-1>*k+j>*2^<k-j>;endrx_quanvalue<i>=<-2*maxvalue+2*quan_set*<n+1>-quan_set>/2;endfigure<2>plot<rx_quanvalue>;wavwrite<rx_quanvalue<:>,48000,'rx_2PSKvoice