readz.m

来自「基于hht的大地电磁信号资料处理MATLAB方法实现」· M 代码 · 共 21 行

function blockz=readz()
clear all;
clc;
%读取EH-4Z文件数据,存储矩阵格式位第一列为频率，其他列为对应得Z文件值
[filename,pathname]=uigetfile('*.*','Please input a data file');
EHyfile=[pathname,filename];
fid=fopen(EHyfile,'r');
%得到文件号
[f,count]=fscanf(fid,'%f %f');


fclose(fid);

L=length(f);
line=L/15;
blockz=zeros(line,15);
n=1:15:L;
for i=1:line
  blockz(i,:)=f(n(i):(n(i)+14))';
end
end