关于升余弦滤波器应用的一个小例子

Apologize

1. clc
   
2. clear all
   
3. close all hidden
   
4. M = 4;                   % The order for modulate of DPSK
   
5. nPacket = 5000;           % The signal length
   
6. x = randint(nPacket,1,M); % Signal for modulate
   
7. 
   
8. h = modem.dpskmod('M',M);      % Creat an object of DPSK modulation
   
9. y = modulate(h,x);        % modulate x get y
   
10. scatterplot(y);
    
11. %% Process raised cosine filter pulse shaping
    
12. Nsamp = 4;               % Oversampling rate
    
13. num = rcosine(1,4);
    
14. ypulse = rcosflt(y,1,4);
    
15. ynoisy = awgn(ypulse,15,'measured');
    
16. ynoisy = ynoisy(3/1*4+1:end-3/1*4);
    
17. ydownsamp = downsample(ynoisy,Nsamp);
    
18. scatterplot(ydownsamp);
    
19. reset(h);
    
20. h = modem.dpskdemod('M',M);
    
21. z = demodulate(h,ydownsamp);
    
22. [num,rt]= symerr(x,z)

复制代码

这里需要注意的是rcosin的参数设置问题，第一个参数是码速率fd，第二个参数采样频率fs，默认的滚降系数是0.5，那么滤波器的时延可以delay=3/fd×fs，这样在降采样之前需要对数据进行处理，去除时延ynoisy = ynoisy(3/1*4+1:end-3/1*4);然进行降采样处理。

希望能对需要的朋友有点参考价值！
滤波器的时延可以用
转自：http://blog.sina.com.cn/s/blog_5def5a660100mhru.html