MATLAB超声传感阵列仿真

用于对超声阵列进行仿真，主要参数有声场强度，阵列个数，间距等，自带GUI界面，易于操作

%% TRANSDUCER ARRAY CALCULATION （T_A_C）%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
%   This program was developed to calculate and visualize the radiation pattern of
%   ultrasound transducers especially for 3D Ultrasound Computer Tomography （USCT）.
%   
%   To run the program properly it is necessary to include the following
%   MATLAB files:
%           - T_A_C_GUI.m                                 [= creating GUI and administrating]
%           - T_A_C_GUI.fig                               [= GUI figure]  
%           - calculate_pattern_acoustic_GUI.m            [= calculation of 3D radiation pattern]
%           - calculate_transducer_array_2D_acoustic.m    [= calculation of 2D radiation pattern and -3 dB range]
%           - T_A_C_global_variables_GUI.m                [= init global variables]
%           - visualize_transducer_array_GUI.m            [= visualization of transducer surfaces]
%           - Soundfield_without_approx.m                 [= calculation of Near Field]
%           - calculate_half_power_Beamwidth.m            [= calculation of half power beamwidth]
%   
%   The TAC Package also includes an “Icons“ and “Examples“ folder. 
%
%   Developed by Benedikt Kohout and Luciano.G. Palacios Folla friendly supported by Robin Dapp
%   at the Karlsruhe Institute of Technology Institute for Data Processing and Electronics Karlsruhe Germany.
%   benedikt.kohout@kit.edu
%                                                                                                      03/2012
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%
function [hwb_int hwb_pres] = calculate_half_power_Beamwidth（C_ALL）
%% calculate half power beamwidth bzw. half power angle -intensity

%constants: defined in calculate_pattern

%T_A_C_global_variables_GUI（）;

global lambda;
global c f;
global Phas_0_x;

lambda = c/f;

%  global hwb_int ;
%  global hwb_pres ;

%lambda = c/f;
%res=res_2D/10;

INT （:1） = C_ALL（:1）;
INT （:2） = （C_ALL（:2）.^2）;
Max_Int = sort （max（INT（:2））‘descend‘）;
Max_Int = Max_Int （11）;
INT （:3） = abs（INT（:2）./Max_Int）;

INT （:4） = 10*log10（INT（:3））;

% C_ALL（:7） = 20*log10（C_ALL（:3））;
% 
% figure（）
% plot（INT（:1）.*（180/pi）INT（:4）‘r‘）
% hold on
% plot（C_ALL（:1）.*（180/pi）C_ALL（:7）‘b‘）


if Phas_0_x >= 0 && Phas_0_x <=pi   %search from 0?to 180?
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    counter_hww_max = 2;
    counter_hww = counter_hww_max;
    hww = 4;%0.6; %define start difference value 
    hww_row = counter_hww;
    max_row = counter_hww_max;
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    
    %find #1 max -> set counter_hww
    while counter_hww_max < length（INT（:1）） 

        if abs（INT（counter_hww_max4）- max（INT（:4））） <=eps

            counter_hww = counter_hww_max;  %save counter max
            max_row = counter_hww_max;  %save max ar

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2014-10-20 14:26  TAC_GUI_V1.3\
     文件        8946  2012-03-16 03:00  TAC_GUI_V1.3\calculate_half_power_Beamwidth.m
     文件        3195  2012-03-16 02:59  TAC_GUI_V1.3\calculate_pattern_acoustic_GUI.m
     文件        5855  2012-03-16 02:59  TAC_GUI_V1.3\calculate_transducer_array_2D_acoustic.m
     文件        1533  2012-03-16 03:09  TAC_GUI_V1.3\COPYING.txt
     目录           0  2014-10-20 14:26  TAC_GUI_V1.3\Examples\
     文件         902  2012-02-24 03:33  TAC_GUI_V1.3\Examples\Attenuation_water.att
     文件         703  2009-07-18 21:39  TAC_GUI_V1.3\Examples\chirp_Example.mat
     文件        1788  2012-03-17 02:13  TAC_GUI_V1.3\Examples\Impedance_example.txt
     文件        1345  2012-03-15 21:22  TAC_GUI_V1.3\Examples\TAC_ART_Example.tac
     文件        2271  2012-03-15 03:57  TAC_GUI_V1.3\Examples\TAC_Welcome_Example.tac
     目录           0  2014-11-03 17:26  TAC_GUI_V1.3\Graficos\
     文件      350478  2014-10-10 15:51  TAC_GUI_V1.3\Graficos\141010-15.50.59-D_P_XY-Array.fig
     文件      193598  2014-10-10 16:15  TAC_GUI_V1.3\Graficos\141010-16.15.12-D_P_XY-Array.fig
     文件      192823  2014-10-10 16:15  TAC_GUI_V1.3\Graficos\141010-16.15.25-D_P_entire_arrangement.fig
     文件       61156  2014-10-10 16:15  TAC_GUI_V1.3\Graficos\141010-16.15.33-D_P_C_all_Vertical.fig
     文件       92687  2014-10-10 16:38  TAC_GUI_V1.3\Graficos\141010-16.38.53-Spatial_trans_characteristics.fig
     文件       61166  2014-10-10 17:28  TAC_GUI_V1.3\Graficos\141010-17.28.58-D_P_C_all_Vertical.fig
     文件      101306  2014-10-13 10:12  TAC_GUI_V1.3\Graficos\141013-10.12.16-D_P_C_all_Horizontal.fig
     文件      193598  2014-11-03 17:18  TAC_GUI_V1.3\Graficos\141103-17.18.16-D_P_XY-Array.fig
     文件      204391  2014-11-03 17:26  TAC_GUI_V1.3\Graficos\141103-17.26.43-D_P_XY-Array.fig
     目录           0  2014-10-20 14:26  TAC_GUI_V1.3\Icons\
     文件        1566  2011-09-20 03:01  TAC_GUI_V1.3\Icons\Active_all.bmp
     文件        1566  2011-09-08 20:38  TAC_GUI_V1.3\Icons\En_dis_ampl_trans_change.bmp
     文件        1566  2011-09-08 20:47  TAC_GUI_V1.3\Icons\En_dis_ampl_trans_change_s&w.bmp
     文件        1566  2011-09-06 20:22  TAC_GUI_V1.3\Icons\Fotoaparat_grau.bmp
     文件        2910  2011-09-20 01:52  TAC_GUI_V1.3\Icons\Fotoaparat_grau_all.bmp
     文件        1566  2011-09-20 03:05  TAC_GUI_V1.3\Icons\Inactive_all.bmp
     文件       27654  2012-02-24 02:35  TAC_GUI_V1.3\Icons\KIT.bmp
     文件         774  2012-02-14 04:25  TAC_GUI_V1.3\Icons\Open_attenuation.bmp
     文件         774  2011-12-09 06:00  TAC_GUI_V1.3\Icons\Open_impedance_characteristic.bmp
............此处省略12个文件信息