I modified the example_us_bmode_linear_transducer.m matlab code to do simulation, and set sampling frequency to fixed 40MHz instead of CFL related one.
    x = 50e-3;
    dx = 0.15e-3;                  % [m]
    Nx = ceil(x/dx);
    medium.alpha_coeff = 0.5;
    % [dB/(MHz^y cm)]
    medium.alpha_power = 1.0; % linear relationship between attenuation and frequency
    medium.alpha_mode = 'no_dispersion';
    % medium.BonA = 6; comment this line to ignore the nonlinear effect
    tone_burst_freq = 3e6; 	% [Hz]
    tone_burst_cycles = 5;
if central frequency was set to 3MHz, the simulation results are very good. However when 5MHz was used, the results was bad. Any one could give me some suggestion.