中国地质大学通信原理结课报告.docVIP

中国地质大学（武汉） 通信原理实验报告 姓名： 学号： 班级： 指导老师：陈中兴 实验一 实验目的 练习使用matlab中的基本函数，为今后的实验打下基础 实验内容 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;%%去采样点个数128 dt=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)); 实验结果 解析与小结 图一为一个周期的方波信号，图二为此方波信号经过傅里叶变换变换到频域后的图形。由 图可以看出，该方波信号经过傅里叶变换后能量主要集中在-2到2之间，且该频率图形对称。 实验二 实验目的 用matlab实现2PSK和QPSK，并做它们的仿真。 实验内容 t=1; forsnr=0:15%%信噪比从0取到15 map=[-11]; tr_2=1+fix(rand(1,10^5));%%产生随机的0,1信号 tx=map(tr_2+1); rx=awgn(tx,snr,measured);%%经过高斯信道获取rx rx_2=[rx0];%%抽样判决 pe=length(find(tr_2-rx_2));%%比较接收信号与发射信号的不相同 ps(t)=pe/length(tr_2);%%计算误码率 t=t+1; end semilogy(0:15,ps,-bo); gridon; 实验结果 2PSK运行结果： QPSK运行结果： QPSK仿真结果： snr=10; 星座图为： snr=5: 星座图为： 2PSK仿真结果： 星座图为： 解析与小结 有实验结果可以看出，实验结果和仿真出来的结果在相同信噪比的情况下基真相同；QPSK与2PSK对照，在相同信噪比的情况下，QPSK性能较好于2PSK。 实验三 实验目的 量化编码，解码 实验内容（一） ——将给定音频文件量化编码 functionwave L=10000;  %%采用采样点  100000个； voicefile=wavread();%%打开给定音频文件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:L samplevalue=voicefile(i);  %%读入音频文件第  i个值 forn=1:M Ifsamplevalue=-maxvalue+quan_set*(n-1)samplevalue= -maxvalue+quan_set*n tx_quanvalue(i)=(-2*maxvalue+2*quan_set*n-quan_set)/2;%%取区间中点为 抽样值 tx2str=dec2bin(n-1,k);%%转变成k位二进制 foru=1:k tx_b(u,i)=str2double(char(tx2str(u)));%%每一行代表一个数据 end end end end figure(2) plot(tx_quanvalue)  %%做出输出波形量化后的图形 tx=tx_b(:); rx=reshape(tx,k,L); fori=1:L  %%对量化编码的数据进行解码 n=0; forj=1:k n=n+rx(j,i)*2^(k-j); rx_quanvalue(