matlab读取SAC类文件函数

function [tdataSAChdr] = fget_sac（filename）
%[tdataSAChdr] = fget_sac（filename）

% read sac into matlab 
% written by Zhigang Peng
% program called
% [head1 head2 head3 data]=sac（filename）;
% [SAChdr]=sachdr（head1 head2 head3）;

% Updated Mon Jul 30 11:21:24 PDT 2001

if nargin <1 error（‘ERROR!! No input file name‘）; end

[head1 head2 head3 data]=sac（filename）;
[SAChdr]=sachdr（head1 head2 head3）;
t = [SAChdr.times.b:SAChdr.times.delta:（SAChdr.data.trcLen-1）*SAChdr.times.delta+SAChdr.times.b]‘;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2016-01-07 19:42  MatSAC_readdata
     文件         514  2006-09-28 03:55  MatSAC_readdatafget_sac.m
     文件     2171841  2006-09-28 03:55  MatSAC_readdataMYJH.dat
     文件      720632  2006-09-28 03:55  MatSAC_readdataMYJH.sac
     文件      720632  2006-09-28 03:55  MatSAC_readdataN.MYJH.Z.sac
     文件         253  2006-09-28 03:55  MatSAC_readdata
ewSacHeader.m
     文件         220  2006-09-28 03:55  MatSAC_readdata
dSac.m
     文件         428  2006-09-28 03:55  MatSAC_readdata
dSacHead.m
     文件        1530  2006-11-06 06:52  MatSAC_readdata
eadme.txt
     文件         861  2006-09-28 03:55  MatSAC_readdatasac.m
     文件        1626  2006-11-05 10:07  MatSAC_readdatasacfft.m
     文件        4785  2006-09-28 03:55  MatSAC_readdatasachdr.m
     文件         243  2006-09-28 03:55  MatSAC_readdatawtSac.m