..中國地質(zhì)大學(xué)〔〕通信原理實驗報告姓名：學(xué)號：班級：指導(dǎo)教師：振興TOC\o"1-3"\h\z\u29735實驗一338331.實驗?zāi)康?97622.實驗容3314163.實驗結(jié)果4317314.分析與小結(jié)421440實驗二5288161.實驗?zāi)康?187972.實驗容574533.實驗結(jié)果56774.分析與小結(jié)1032247實驗三10243731.實驗?zāi)康?0320122實驗容〔一〕10279732.1實驗結(jié)果12215202.2分析與小結(jié)13315443實驗容〔二〕13116053.1實驗結(jié)果15127183.2分析與小結(jié)16實驗一1.實驗?zāi)康木毩?xí)使用matlab中的根本函數(shù)，為以后的實驗打下根底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);Fm*=(f(end)-f(1)+df);dt=1/Fm*;N=length(sf);T=dt*N;t=0:dt:T-dt;sff=fftshift(sf);st=Fm**ifft(sff);clearall;closeall;T=1;N_sample=128;%%去采樣點個數(shù)128dt=T/N_sample;t=0:dt:T-dt;st=[ones(1,N_sample/2),-ones(1,N_sample/2)];%%產(chǎn)生一個周期的方波信號subplot(211);plot(t,st);a*is([01-22]);*label('t');ylabel('s(t)');subplot(212);[f,sf]=T2F(t,st);plot(f,abs(sf));holdon;a*is([-101001]);*label('f');ylabel('s(f)');3.實驗結(jié)果4.分析與小結(jié)圖一為一個周期的方波信號，圖二為此方波信號經(jīng)過傅里葉變換變換到頻域后的圖形。由圖可以看出，該方波信號經(jīng)過傅里葉變換后能量主要集中在-2到2之間，且該頻率圖形對稱。實驗二1.實驗?zāi)康挠胢atlab實現(xiàn)2PSK和QPSK，并做它們的仿真。2.實驗容t=1;forsnr=0:15%%信噪比從0取到15map=[-11];tr_2=1+fi*(rand(1,10^5));%%產(chǎn)生隨機的0,1信號t*=map(tr_2+1);r*=awgn(t*,snr,'measured');%%通過高斯信道得到r*r*_2=[r*>0];%%抽樣判決pe=length(find(tr_2-r*_2));%%比擬接收信號與發(fā)射信號的不同ps(t)=pe/length(tr_2);%%計算誤碼率t=t+1;endsemilogy(0:15,ps,'-bo');gridon;3.實驗結(jié)果2PSK運行結(jié)果：QPSK運行結(jié)果：QPSK仿真結(jié)果：snr=10;星座圖為：snr=5:星座圖為：2PSK仿真結(jié)果：星座圖為：4.分析與小結(jié)有實驗結(jié)果可以看出，實驗結(jié)果和仿真出來的結(jié)果在一樣信噪比的情況下根本一樣；QPSK與2PSK相比，在一樣信噪比的情況下，QPSK性能較好于2PSK。實驗三1.實驗?zāi)康牧炕幋a，解碼2實驗容〔一〕——將給定音頻文件量化編碼functionwaveL=10000;%%選取采樣點100000個；voicefile=wavread('voicefile.wav');%%翻開給定音頻文件wavefile并將其作圖figure(1)plot(voicefile)ma*value=ma*(abs(voicefile));%%求出音頻的最大振幅M=4;k=log2(M);quan_set=2*ma*value/M;%%選取量化間隔t*_quanvalue=zeros(1,L);%%初始化發(fā)送和接收數(shù)組r*_quanvalue=zeros(1,L);fori=1:Lsamplevalue=voicefile(i);%%讀入音頻文件第i個值forn=1:MIfsamplevalue>=-ma*value+quan_set*(n-1)&&samplevalue<=-ma*value+quan_set*nt*_quanvalue(i)=(-2*ma*value+2*quan_set*n-quan_set)/2;%%取區(qū)間中點為抽樣值t*2str=dec2bin(n-1,k);%%轉(zhuǎn)化為k位二進制foru=1:kt*_b(u,i)=str2double(char(t*2str(u)));%%每一行代表一個數(shù)據(jù)endendendendfigure(2)plot(t*_quanvalue)%%做出輸出波形量化后的圖形t*=t*_b(:)';r*=reshape(t*,k,L);fori=1:L%%對量化編碼的數(shù)據(jù)進展解碼n=0;forj=1:kn=n+r*(j,i)*2^(k-j);r*_quanvalue(i)=(-2*ma*value+2*quan_set*(n+1)-quan_set)/2;endendwavwrite(r*_quanvalue(:),48000,'r*_voicefile.wav');%%將解碼后的數(shù)據(jù)保存入新增文件r*_voicefilefigure(3)plot(r*_quanvalue)2.1實驗結(jié)果原噪聲圖形量化編碼后的圖形解碼后的圖形2.2分析與小結(jié)由以上的圖形可以看出，原來的音頻信號是一串隨機的信號，經(jīng)過量化編碼后變成圖二所示的8個量化電平的信號，因為沒有參加噪聲，經(jīng)過解碼后恢復(fù)的信號和編碼是一樣的。3實驗容〔二〕——將文件r*_voicefile通過2ASK和QPSK編碼，通過高斯信道，然后解碼%%同上將音頻文件量化編碼functionwave2L=10000;voicefile=wavread('r*_voicefile.wav');figure(1)plot(voicefile)ma*value=ma*(abs(voicefile));M=4;k=log2(M);quan_set=2*ma*value/M;fori=1:Lsamplevalue=voicefile(i);forn=1:Mifsamplevalue>=-ma*value+quan_set*(n-1)&&samplevalue<=-ma*value+quan_set*nt*_quanvalue(i)=(-2*ma*value+2*quan_set*n-quan_set)/2;t*2str=dec2bin(n-1,k);%%轉(zhuǎn)化為k位二進制foru=1:kt*_b(u,i)=str2double(char(t*2str(u)));%%每一行代表一個數(shù)據(jù)endendendEnd%%通過高斯信道然后解碼。snr=100;map=[-1,1];fori=1:Lforu=1:kt*(u,i)=map(t*_b(u,i)+1);endendr*=awgn(t*(:),snr,'measured');%%經(jīng)過高斯信道fori=1:k*Lr*(i)=r*(i)>0;%%抽樣判決endfori=1:L%%解碼n=0;forj=1:kn=n+r*((i-1)*k+j)*2^(k-j);endr*_quanvalue(i)=(-2*ma*value+2*quan_set*(n+1)-quan_set)/2;endfigure(2)plot(r*_quanvalue);wavwrite(r*_quanvalue(:),48000,'r*_2PSKvoicefile.wav');%%將文件解碼后寫入r*_2PSKvoic