close all;
clear all;

Fd= 1; %輸入信號的采樣率
Fs = 8;  %輸出信號的采樣率
Delay = 3;  %濾波器的群時延
R = 0.5;   %滾降因子
[yf,tf] = rcosine(Fd,Fs,'fir',R,Delay);  %yf=b,tf=a;
%plot(yf)
freqz(yf,tf,512);
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');
%hold on;stem(yt(1:40),'r');
x_f=fft(x);
xt_f=fft(xt);
y_f=fft(y);
yt_f=fft(yt);
figure;
%stem(xt_f,'.');
%hold on;stem(y_f,'r');
subplot(211);
plot(abs(xt_f));
subplot(212);
plot(abs(y_f));