通信原理程序.docVIP

三、实验内容及程序分析 I、实验一：发送滤波器长度为N=31，时域抽样频率F0为4 /Ts ，滚降系数分别取0.1、0.5、1，计算并画出此发送滤波器的时域波形和频率特性。以此发送滤波器构成最佳基带系统，计算并画出接收滤波器的输出信号波形和整个基带系统的频率特性。 实验程序： N=32;%抽样点数32 L=4; M=N/L;%码元数 Rs=0.25; Ts=1/Rs; fs=L/Ts; Bs=fs/2; T=N/fs; t=-T/2+[0:N-1]/fs; f=-Bs+[0:N-1]/T; %升余弦滚降alpha=0.5 alpha=0.5; Hcos=zeros(1,N); ii=find(abs(f)(1-alpha)/(2*Ts)abs(f)=(1+alpha)/(2*Ts)); Hcos(ii)=Ts/2*(1+cos(pi*Ts/alpha*(abs(f(ii))-(1-alpha)/(2*Ts)))); ii=find(abs(f)=(1-alpha)/(2*Ts)); Hcos(ii)=Ts; %gen sheng yu xian Hrcos=sqrt(Hcos); ft=zeros(1,N); ft=real(ft2(Hrcos,fs)); %alpha=0.1 alpha=0.1; Hcos1=zeros(1,N); ii=find(abs(f)(1-alpha)/(2*Ts)abs(f)=(1+alpha)/(2*Ts)); Hcos1(ii)=Ts/2*(1+cos(pi*Ts/alpha*(abs(f(ii))-(1-alpha)/(2*Ts)))); ii=find(abs(f)=(1-alpha)/(2*Ts)); Hcos1(ii)=Ts; %gen sheng yu xian Hrcos1=sqrt(Hcos1); ft1=zeros(1,N); ft1=real(ft2(Hrcos1,fs)); %升余弦滚降alpha=1 alpha=1; Hcos2=zeros(1,N); ii=find(abs(f)(1-alpha)/(2*Ts)abs(f)=(1+alpha)/(2*Ts)); Hcos2(ii)=Ts/2*(1+cos(pi*Ts/alpha*(abs(f(ii))-(1-alpha)/(2*Ts)))); ii=find(abs(f)=(1-alpha)/(2*Ts)); Hcos2(ii)=Ts; %gen sheng yu xian Hrcos2=sqrt(Hcos2); ft2=zeros(1,N); ft2=real(ft2(Hrcos2,fs)); %画图 subplot(3,2,2); stem(f,Hrcos1,.); axis([-Bs,Bs,0,max(Hrcos1)]); title(alpha=0.1的根升余弦发送滤波器的频域波形); grid; subplot(3,2,1); stem(t,ft1,.); axis([-T/2,T/2,1.1*min(ft1),1.1*max(ft1)]); grid; title(alpha=0.1的根升余弦的发送滤波器的时域冲激响应); subplot(3,2,4); stem(f,Hrcos,.); axis([-Bs,Bs,0,max(Hrcos)]); title(alpha=0.5的根升余弦发送滤波器的频域波形); grid; subplot(3,2,3); stem(t,ft,.); axis([-T/2,T/2,1.1*min(ft),1.1*max(ft)]); grid; title(alpha=0.5的根升余弦的发送滤波器的时域冲激响应); subplot(3,2,6); stem(f,Hrcos2,.); axis([-Bs,Bs,0,max(Hrcos2)]); title(alpha=1的根升余弦发送滤波器的频域波形); grid; subplot(3,2,5); stem(t,ft2,.); axis([-T/2,T/2,1.1*min(ft2),1.1*max(ft2)]); grid; title(alpha=1的根升余弦的发送滤波器的时域冲激响应); 实验一仿真结果： 2、实验二 根据基带系统模型，编写程序，设计无码间干扰的二进制数字基带传输系统。要求要传输的二进制比特个数、比特速率Rb（可用与Ts的关系表示）、信噪比SNR、滚降系数α是可变的。 1) 生成一个0、1 等概率分布的二进制信源序列（伪随机序列）。可用MATLAB 中的rand 函数生成一组0~1 之间均匀分布的随机序列，如产生的随机数在（0，0.5）区间内，则为0；如果在（0.5，1）区间内，则为1。 2) 基带