Hi,
I want to put one wire in scattering medium, is this the right way?
sound_speed_map(Nx/2, Ny/2,1:end) = 5900;
density_map(Nx/2, Ny/2,1:end) = 7850/1.5;
% assign to the medium inputs
medium.sound_speed = sound_speed_map;
medium.density = density_map;
Hi shock,
You first need to define the maps to be the correct dimension, and assign the background value. Then you can assign the wire using your code, e.g.,
Nx = 50;
Ny = 50;
Nz = 50;
c0 = 1500;
cwire = 5900;
sound_speed_map = c0 * ones(Nx, Ny, Nz);
sound_speed_map(Nx/2, Ny/2, 1:end) = cwire;
voxelPlot(sound_speed_map);Keep in mind if you are defining a scatterer right to the edge of the grid to place the perfectly matched layer (PML) outside your simulation by setting the optional input 'PMLInside', false.
Brad.
Hi Bradley
Thank you so much!
You must log in to post.