Fd= 1; %輸入信號的采樣率
Fs = 8;  %輸出信號的采樣率
Delay = 3;  %濾波器的群時延
R = 0.5;   %滾降因子
[yf,tf] = rcosine(Fd,Fs,'fir',R,Delay);
plot(yf)
grid
xlabel('Time')
ylabel('Amplitude');
% impz(rrcfilter,1); % 另外一種畫沖激響應的方法
x=rand(100,1);  %產生隨機信號
xt=zeros(1,800);  %對輸入信號后面補零
xt(1:8:end)=x;
y=filter(yf,tf,xt);  %經過升余弦濾波器
yt=y((length(yf)+1)/2:8:end);    %下采樣，去除濾波器的抽頭時延
figure;stem(x(1:40),'.');
hold on;stem(yt(1:40),'r');