frm_trk.m

来自「这是一个用于语音信号处理的工具箱」· M 代码 · 共 44 行

function [form_trk] = frm_trk(gci_a, ff, Ts, LEN)

nskip=1;
gci_len=length(gci_a);
len=floor(LEN/nskip);
[M,N]=size(ff);

form_trk=zeros(len,N);
kk=1;

for ii=1:gci_len
	if ii == 1
		for jj=1:floor(gci_a(1)/nskip),
			form_trk(kk,:)=zeros(1,N);
			kk=kk+1;
		end
	elseif ii == gci_len
		nn=floor((gci_a(ii)-gci_a(ii-1))/nskip);		
		for jj=1:nn
			form_trk(kk,:)=ff(ii,:);
			kk=kk+1;
		end
		for jj=floor(gci_a(ii)/nskip):len,
			form_trk(kk,:)=zeros(1,N);
			kk=kk+1;
		end
	else
		int=gci_a(ii)-gci_a(ii-1);
		nn=floor(int/nskip);		
		pit=1/(int*Ts);
		if pit > 300 | pit < 50
			for jj=1:nn
				form_trk(kk,:)=zeros(1,N);
				kk=kk+1;
			end
		else
			for jj=1:nn
				form_trk(kk,:)=ff(ii,:);
				kk=kk+1;
			end
		end
	end
end