matlab读取SAC类文件函数

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_readdata\fget_sac.m
     文件     2171841  2006-09-28 03:55  MatSAC_readdata\MYJH.dat
     文件      720632  2006-09-28 03:55  MatSAC_readdata\MYJH.sac
     文件      720632  2006-09-28 03:55  MatSAC_readdata\N.MYJH.Z.sac
     文件         253  2006-09-28 03:55  MatSAC_readdata\newSacHeader.m
     文件         220  2006-09-28 03:55  MatSAC_readdata\rdSac.m
     文件         428  2006-09-28 03:55  MatSAC_readdata\rdSacHead.m
     文件        1530  2006-11-06 06:52  MatSAC_readdata\readme.txt
     文件         861  2006-09-28 03:55  MatSAC_readdata\sac.m
     文件        1626  2006-11-05 10:07  MatSAC_readdata\sacfft.m
     文件        4785  2006-09-28 03:55  MatSAC_readdata\sachdr.m
     文件         243  2006-09-28 03:55  MatSAC_readdata\wtSac.m