GardnerForMatlab.docVIP

GardnerForMatlab ff1=fopen(get(handles.filename,string),r); % 以只读方式读取文件 %%%%% rr=fread(ff1,uint16); %%%% 读取二进制文件 %%%% 以无符号(有的数据是有符号的)16位数据读取文件(.dat文件可以设置低位或高位在前) % 设定全局变量。hn为FIR滤波器系数 global hn N FT Hw f0; %% 加上采集频率为100k FT=100000; Input_Ts=1/FT; rr=2*bitand(rr,4095)/4095-1; %%%% 去除控制位 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%% 输入序列 Sample_N=10*1024; %% 样值个数 200k*1024 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%% 输入信号，I路和Q路采样数据 begin00=10001; % 初始位置，去掉前面的错误采集数据16401 Input_In=rr(begin00:2:Sample_N+begin00-1); Input_Qn=rr(begin00+1:2:Sample_N+begin00); %%%%%% 计算经过滤波器后的输出信号序列DATA_I和DATA_Q %%% 是否需要对总的信号进行滤波而不是对I路和Q路分别滤波，其实两种滤波方法计算的结果时一样的。 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%% 对整个信号进行滤波 % Initial=zeros(1,N); % x=[Initial Input_signal]; % %%%%% 利用矩阵（数组）的方法计算 % for mm=1:Sample_N % DATA_Signal(mm)=sum(hn.*x(N+mm:-1:mm)); % end % DATA_I=real(DATA_Signal); % DATA_Q=imag(DATA_Signal); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%% 对I路和Q路分别滤波 % Initial=zeros(1,N); % x1=[Initial Input_In]; % x2=[Initial Input_Qn]; % %%%%% 利用矩阵（数组）的方法计算 % for mm=1:Sample_N/2 % DATA_I(mm)=sum(hn.*x1(N+mm:-1:mm)); % DATA_Q(mm)=sum(hn.*x2(N+mm:-1:mm)); % end %% 如果不考虑滤波 % DATA_I=Input_In; % DATA_Q=(Input_Qn+0.0383)*1.0656; DATA_I=Input_In; DATA_Q=Input_Qn; %DATA_I(10078)=-DATA_I(10078); %DATA_Q(10078)=-DATA_Q(10078); %DATA_Q(10080)=-DATA_Q(10080); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%% 画出y的频谱 % %%%% 计算输入和输出的频谱 % %%%% % % 由于序列x的傅立叶变换是周期imiga为2*pi的周期函数，因此这里的频率应为归一化频率 % n=[0:Sa