%*** 設計4碼元成型濾波器 rafa=0.6的升余弦濾波器 ***%
clc;%清除命令窗口
clear all;%清除工作空間workspace,以釋放內存
nn1=1:48;
x3=[(nn1-0)>=0];%1-48全1
x4=[(nn1-13)>=0];%前12為0，后36為1
shengyuxuan(nn1)=x3-x4;%前12為1，后36為0
nn2=13:48;
rfa=linspace(0,1,36);%生成36點線性向量，在0-1之間且包括0和1
shengyuxuan(nn2)=0.5*(cos(rfa*pi)+1);%36個點的升余弦yuxian
shengyuxuan(83:129)=fliplr(shengyuxuan(2:48));   %fliplr函數使行向量左右反轉
Hshengyuxuan=shengyuxuan;
%for nn3=1:8    %作IDFT變換求得濾波器的時域特性，生成8個點的時域序列作為濾波器系數
%    nn4=1:length(Hshengyuxuan);
%    htest(nn4)=Hshengyuxuan(nn4).*cos(2*pi*(nn4-1)*(nn3-1)/129);  %兩向量對應元素相乘，點乘
%    hsheng1(nn3)=sum(htest(nn4));   
%end;  
%hsheng3=hsheng1/length(Hshengyuxuan);
hsheng1=ifft(Hshengyuxuan);
hsheng3=hsheng1(1:8);
hsheng2=fliplr(hsheng3);
hsheng(1:8)=hsheng2(1:8);
hsheng(9:15)=hsheng3(2:8);