test2.m

来自「Galileo BOC 码产生捕获相关仿真结果」· M 代码 · 共 18 行

%a=MySglGen(1,200)
% data=[1.2,1.7,3,5,7,9,56,-123,-2,-1,-9,-4,-3,-2.7,1.5,2.4,0.6,-0.6,-1.4,-2.7];
% partition=linspace(-3,3,6);
 %codebook=[-1:5];
 %[index,data]=quantiz(data,partition,codebook)
 in=-64:1:64;
 N=length(in);
 out=zeros(1,N);
 extract=zeros(1,floor(N/4));
 lpf_b = [1, 0, -4, 0, 16, 26, 16, 0, -4, 0, 1];
 out = filter(lpf_b,1,in);
 for i=1:floor(N/4)
    extract(i)=out(i*4);    
end
subplot(1,2,1);
stem(in,out)
subplot(1,2,2);
stem(extract)