agcvec.m

全数字调制解调器设计的matlab程序 有CIC FIR等滤波器

%本程序生成测试向量。
clc;
clear;

framelength = 5000;
datai=randsrc(1,framelength,[-1 1]);
dataq=randsrc(1,framelength,[-1 1]);

stuff=zeros(1,15);
for i=1:length(datai)
    x((i-1)*16+1)=datai(i);
    x=[x stuff];
end;
data_i=x;

for i=1:length(dataq)
    y((i-1)*16+1)=dataq(i);
    y=[y stuff];
end;
data_q=y;


%平方根升余弦滤波器的系数
R=1.0;
n_T=[-4 4];
rate=16;
T=1;
b = rcosfir(R,n_T,rate,T,'sqrt');

wave_i=filter2(b,data_i,'same');
wave_q=filter2(b,data_q,'same');

I=filter2(b,wave_i,'same');
Q=filter2(b,wave_q,'same');

for i = 1:(framelength-1)
    idecision(i) = I(i*16+1);
    qdecision(i) = Q(i*16+1);
end;

idecision = idecision + sqrt(0.005).*randn(1,(framelength-1));
qdecision = qdecision + sqrt(0.005).*randn(1,(framelength-1));

for i = 1:(framelength-1)
    if idecision(i) > 0
        idecision(i) = round(idecision(i) * (2^9-1));
    else
        idecision(i) = round(idecision(i) * (2^9-1)) + 2^16;
    end;
end;

for i = 1:(framelength-1)
    if qdecision(i) > 0
        qdecision(i) = round(qdecision(i) * (2^9-1));
    else
        qdecision(i) = round(qdecision(i) * (2^9-1) + 2^16);
    end;
end;
     
for i = 1:(framelength - 1)
    idec(i,:) = dec2hex(idecision(i),4);
    qdec(i,:) = dec2hex(qdecision(i),4);
end;

fid=fopen('agcvec.txt','w');
for i=1:(framelength - 1)
    fprintf(fid,'%8.1f> 0 %4s %4s  \n',50.0*(2*i-2),idec(i,:),qdec(i,:));
    fprintf(fid,'%8.1f> 1 %4s %4s  \n',50.0*(2*i-1),idec(i,:),qdec(i,:));
end;
fclose(fid);