丹麦理工大学MATLAB编写的超声声场仿真软件

%  Procedure for calculating the spatial impulse response
%  for an aperture.
%
%  Calling:  [h start_time] = calc_h（Thpoints）; 
%
%  Parameters:  Th     - Pointer to the transducer aperture.  
%               points - Field points. Vector with three columns （xyz） 
%                        and one row for each field point. 
%
%  Return:      h          - Spatial impulse response in m/s.
%               start_time - The time for the first sample in h.
%
%  Version 1.01 October 4 1996 by Joergen Arendt Jensen

function [h start_time] = calc_h （Thpoints）

%  Check the point array

  [mn]=size（points）;
  if （n ~= 3）
    error （‘Points array must have three columns‘）;
    end

%  Call the C-part of the program to show aperture

  [h start_time] = Mat_field （4001Thpoints）;