opdl2.m

这是一个关于hht变换很有用的工具箱

function z=opdl2(nam,n,k)
 
% The function OPDL2 opens and reads 2-D binary file in 'double 64' format.
%
% Example-
% 	x=opdl2('x',1000,9).
%
% Calling sequence-
% z=opdl2(nam,n,k)
%
% Input-
%	nam	- string, that represents the name of the file
%	n	- 1st dimension of the data
%	k	- 2nd dimension of the data
% Output-
%	z	- 2-D data z(n,k) in memory

%----- Open the file
fd=fopen(nam,'rb');

%----- Read the record in memory
z=fread(fd,[n,k],'float64');

%----- Close the file
fclose(fd);
return