基带脉冲成形滤波器.docVIP

基带脉冲成形滤波器 一 基本原理 基带脉冲成形滤波器实际上是一个内插滤波器 上图为一内插器，内插M（码元间差M-1个0），增加数模转换的精度，脉冲成形滤波器h(t),使得码元成形，消除码间干扰，并且能压缩频谱，在基带中h(t)可视为发送滤波器，直接将成形波形发射出去，在频带中经过调制后发射出去。 一般选择发端的脉冲成形滤波器具有根升余弦特性。 升余弦滚降函数为： α=W2/W1为滚降系数。TS为码元间隔。 二 MATLAB仿真 Fd=1; Fs=8; Delay=3; R=0.5; [yf,tf]=rcosine(Fd,Fs,fir,R,Delay); %将原始信号内插后通过升余弦滚降滤波器后的输出 figure(1) plot(yf); grid; xlabel(Time); ylabel(Amplitude); title(升余弦滚降滤波器); x=randint(100,1)*2-1;%原始输入信号为+1，-1码 xt=zeros(1,800); xt(1:8:end)=x; y=filter(yf,tf,xt); yt=y((size(yf)+1)/2:8:end); figure(2); stem(yt(1:40)); title(抽取后输出) grid; figure(3); stem(x(1:40)); title(原始信号输出) figure(4) plot(y(1:100)); title(滤波后输出) grid; 三.硬件实现 方法1.用数字滤波器设计脉冲成形滤波器时关键是抽头系数的确定，可以根据h(t)采样得到。 方法2.利用MATLAB中的FDATOOL工具对升余弦滚降滤波器进行设计。 附：用FPGA实现的波形如下： y = fmmod(x,Fc,Fs,freqdev) - y：调制后的输出的时域信号数据； - x：输入的时域数据； - Fc：调制载波频率。绝对值，Hz。 - Fs：整个系统的采样频率。当然越高越好，能够更清晰的描述整个模拟调制的情况。 - freqdev：每单位的调频带宽（如果输入信号最大值是2，调频带宽就是这个值的2倍）。 这几个参数看起来简单，但是如果没有理解其中含义，初学者把握不住期间关系，就会发现绘制出来的图片和期望相差很大。特别是2007的matlab的帮助文档里面的例子，简直就是让人头晕。 其中的关键点是： X的点数，Fs和Fc间的配置要合理。才能够看到期望的调频图像。具体说，就是： - Fs代表了采样频率，是整个函数的计算基础。 - x里面的每个点，时域上都对应于一个Fs的采样点。 - Fs必须比Fc足够高（10倍），才能够在绘制输出的图像时候，有足够的精度表示出来。 例子如下： x=[0:0.001:1]; x=[x,[1:-0.001:0]]; Fc=800; Fs=100000; fredev=1000; y=fmmod(x,Fc,Fs,fredev); plot(y); % 解调输出 z=fmdemod(y,Fc,Fs,fredev); % 绘制频谱和时域图 lx=length(y); Freqs=Fs; freq=[-Freqs/2:Freqs/lx:Freqs/2-Freqs/lx]; S=fftshift(abs(fft(y))); Sdb=10*log10(S); figure; plot(freq,Sdb); 关于平方根升余弦滚降滤波器 平方根升余弦滚降滤波器在发送端输出的信号是连续信号还是离散信号？ 信道传输要符合NYQUIST不是真定理，要求从发到收信道是升余弦滚降滤波器。实际应用中在发弄个平方根升余弦滚降滤波器，收弄个平方根升余弦滚降滤波器，这个出来是数字的，可以用D/A转换，在调制到射频上发射！ 那像下面发射机的一段函数，delay是自己给定的还是由公式得到的？ 关于升余弦滚降滤波器书上都讲得不全 Fd=1e3 Fc=128e3; % Generate an randm sequence of 0 and 1 % Input binary data of 1KHz n=0:6000 data=100; %发送100个0、1随机序列 rand_data=randn(1,data); for i=1:data if rand_data(i)=0 input(i)=1; else input(i)=0; end end % Series to Parallel Ix=input(1:2:length(input)); Qx=input(2:2:length(input)); % QPSK