ppfilter.m

QAM调制函数

function ppfilter
%reference:   PIERRER.CHEVILLAT MEMBER IEEE .
%            'Optinum Fir Transmitter and Receiver Filters for Data
%            Transmission over Band-Limited Channels.'IEEE TRANSACTION ON
%            COMMUNICATION VOL.COM-30,AUGUST 1982

%采用了升余弦滤波器来设计匹配滤波器
%提高了在判决时刻的信噪比，达到了近乎零ISI
clear  all;
N=16;               %窗口长度
alpha=0.4;          %滚降系数
d=N-2*floor(N/2);   %判断n的奇偶性
if d==0
N=N;
else N=N-1;
end
L=N/2;
fs=10000;
fb=fs/2;
Ts=1/fs;
Tb=1/fb;    %抽样比

for n=1:1:N
   nT=(n-1-L)*Ts/Tb;

   if nT==0 
    h(n)=(1-alpha)+4*alpha/pi; 
   elseif (1-(4*alpha*nT)^2)==0 
     h(n)=(pi*(1-alpha)*sin(pi*(1-alpha)*nT)+4*alpha*cos(pi*(1+alpha)*nT)-4*alpha*nT*(1+alpha)*sin(pi*(1+alpha)*nT))/(-32*pi*(alpha*nT)^2);
   else
    h(n)=(sin(pi*(1-alpha)*nT)+4*alpha*nT*cos(pi*(1+alpha)*nT))/(pi*nT*(1-(4*alpha*nT)^2));
    end
end

%matlab库函数求的各系数
b=rcosine(fb,fs,'fir/sqrt',alpha);
figure(1);
plot(h,'r*:'); 
hold on;
plot(b,'g+:');
xlabel('t/T');ylabel('h(t)'); 


%一般情况下在性能比较好的情况下，可以采用升余弦滚降滤波器来代替匹配滤波器.
%匹配滤波器

%接收到的信号为调制在高频上的信号
[yid,yqd]=aa;
yid_f=fft(yid)
figure(2)
subplot(2,1,1);
plot(abs(yid_f));
yid_g=conj(yid_f)
subplot(2,1,2);
plot(abs(yid_g));
t0=2;  %判决时刻
t=t0;
yid_o=yid_g*e^(-j*2*pi*t0)