shengyuxian.m

cdma 系统级仿真

%*** 设计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);