升余弦滚降滤波器仿真测试.doc

word通信原理仿真作业 班 级 1401014 学 号140101400 姓 名 任课教师 升余弦滤波器仿真测试一、 实验要求利用Matlab做出一组升余弦滚降滤波器的冲激响应，滚降系数为0，0.5，0.75和1，并通过FFT求出其幅频特性二、 实验原理1. 无码间串扰的时域条件假如想要消除码间串扰，应有：由于an是随机的，要想通过在接收滤波器输出的信号抽样信号中的各项相互抵消使码间串扰为0是不行的，这就需要对基带传输系统的总传输特性h(t)的波形提出要求如果相邻码元的前一个码元的波形到达后一个码元抽样判决时刻已经衰减到0，就能满足要求但是，这样的波形不易实现，因为现实中的h(t)波形有很长的“拖尾〞，也正是由于每个码元的“拖尾〞造成了对相邻码元的串扰这就是消除码间串扰的根本思想只要基带传输系统的冲激响应波形h(t)仅在本码元的抽样时刻上有最大值，并在其他码元的抽样时刻上均为0，如此可消除码间串扰所以应满足下式：由此我们可以得到基带传输特性应满足的频域条件：此条件称为奈奎斯特第一准如此2. 由此准如此可设计出理想低通滤波器：但理想低通滤波器存在着问题：理想矩形特性的物理实现极为困难。

理想的冲激响应h(t) 的“尾巴〞很长，尾部摆幅较大，衰减缓慢，对位定时的要求严格，要求抽样时刻严格对准零点当定时存在偏差时，偏离零点，可能出现严重的码间串扰3. 解决方法——引入滚降滚降系数：理论传输特性：理论冲击响应：三、 试验流程1. 确定根本参数码元速率为1000Bd采样速率为 10000Hz输入到响应峰值之间的延迟为5码元时隙数滚降系数分别为0, 0.5, 0.75, 1〔循环执行〕Fd=1e3; % 输入数字序列的采样率即码元速率Fs=Fd*10; %采样频率 此式保证了Fs/Fd为正整数delay=5; %输入到响应峰值之间的延迟〔单位是码元时隙数〕2. 运用rcosine函数进展升余弦滤波器设计num = rcosine(Fd,Fs, 'fir/normal',r,delay);其中‘fir/normal’用于FIR滚升余弦滤波器设计3. 制作冲击响应图每次用不同的颜色标识冲击响应曲线k=[rand(),rand(),rand()];%每个循环改变一次RGB颜色 figure(1); plot(t,num,'Color',k); axis([0 0.01 -0.3 1.1]); xlabel('t'); ylabel('h(t)'); title('冲击响应');hold on;4. 使用快速傅里叶变换制作传输特性曲线Hw=abs(fft(num,1000)); %fft快速傅里叶变换 N=1000 abs求得振幅 f=(1:Fs/1000:Fs)-1; %频率分辨率为Fs/N=10 figure(2); plot(f,Hw,'Color',k); axis([0 1500 0 12]); xlabel('f'); ylabel('H(w)'); title('传输特性'); hold on;四、 实验结果仿真传输特性：仿真冲击响应：理论传输特性：理论冲击响应：通过理论与仿真的升余弦滤波器的冲击响应与传输特性之间的比拟，我们可以发现：当时理论值与仿真值非常吻合，但当时仿真的特性曲线不为矩形。

这是因为理想FIR数字滤波器不是因果稳定的，无法设计但可以用窗函数法进展逼近时升余弦滤波器实为矩形窗函数过渡带很短，但会产生吉普斯效应，即通带与阻带的波纹五、 程序附录与注释〔注：浅绿色的为原有注释，深绿色为新加注释即我对程序的理解与分析〕clcclear allclose allFd=1e3; % 输入数字序列的采样率即码元速率Fs=Fd*10; %采样频率 此式保证了Fs/Fd为正整数delay=5; %输入到响应峰值之间的延迟〔单位是码元时隙数〕for r=[0, 0.5, 0.75, 1] %滚降系数 num = rcosine(Fd,Fs, 'fir/normal',r,delay); %‘fir/normal’用于FIR滚升余弦滤波器设计 fir有限脉冲响应〔Finite Impulse Response〕 normal常规的升余弦函数 t=0:1/Fs:1/Fs*(length(num)-1); k=[rand(),rand(),rand()];%每个循环改变一次RGB颜色 figure(1); plot(t,num,'Color',k); axis([0 0.01 -0.3 1.1]); xlabel('t'); ylabel('h(t)'); title('冲击响应'); hold on; Hw=abs(fft(num,1000)); %fft快速傅里叶变换 N是DFT变换区间长度N=1000 abs求得振幅 f=(1:Fs/1000:Fs)-1; %频率分辨率为Fs/N=10 figure(2); plot(f,Hw,'Color',k); axis([0 1500 0 12]); xlabel('f'); ylabel('H(w)'); title('传输特性'); hold on;end文案大全。