dataread5.m

基于麦克风阵波束形成的谱相减方法在语音增强中的应用

%本程序用于处理 原始采集数据

function [n1,n2,n3,n4,n5]=dataread5(ss,count)
%[c1,c2,c3,c4,c5]=dataread5('d:\matlab6p1\workdat\datn\s5ch6cm8k1.dat',count);
%[n1,n2,n3,n4,n5]=dataread5('d:\matlab6p1\workdat\datn\n5ch6cm8k1.dat',count);
%[n1,n2,n3,n4,n5]=dataread5('d:\matlab6p1\workdat\datn\log1.dat',21760);

%128*168=21504,128*170=21760;

if nargin<=1
  count=21760;
  disp '数据个数未指定，但须是128的倍数！'
end;

%fid=fopen('d:\matlab6p1\workdat\log1.dat','r');
%
fid=fopen(ss,'r');

s=fread(fid,262144*3,'int16');% 这里约定了数据的最大长度，可以改变262144
fclose(fid);
l=length(s)
if (l/5)<count
    disp 'count is too large';
end

%plot(s);
%smean=mean(2:l);
%s(1)=smean;
%s=s-mean(s);figure;plot(s); so=s/max(s);
%c=zeros(1,l/6);

n1=zeros(1,count);n2=zeros(1,count);n3=zeros(1,count);
n4=zeros(1,count);n5=zeros(1,count);n6=zeros(1,count);

n1=s(1:5:l)';
n2=s(2:5:l)';
n3=s(3:5:l)';
n4=s(4:5:l)';
n5=s(5:5:l)';

n1=n1(1:count)-mean(n1);
n2=n2(1:count)-mean(n2);
n3=n3(1:count)-mean(n3);
n4=n4(1:count)-mean(n4);
n5=n5(1:count)-mean(n5);

%figure;plot(c1);figure;plot(c2,'r');
%figure;
%cohere(c1(1:5000),c2(1:5000),1024,8000,[],512,'linear');