Using Matlab commands!

1. Compute the time-harmonic pressure field generated by a rectangular piston by numerically integrating the Raleigh-Sommerfeld integral with the double sum described in class. Use the following parameters: frequency = 2.5MHz, speed of sound = 1500m/s, piston width a = 52, piston height b = 102, where the piston is centered at the origin and located in the xy plane. Sample the aperture with nxprime = 30 and nyprime = 60. Evaluate the pressure in the xy plane at z= (a+b)/(202). Calculate the pressure on a 61 by 61 pressure grid, where xmax = -xmin = 1.5 *a/2 and ymax = -ymin = 1.5 *b/2. Plot the absolute value of the result with the mesh command. 1. Compute the time-harmonic pressure field generated by a rectangular piston by numerically integrating the Raleigh-Sommerfeld integral with the double sum described in class. Use the following parameters: frequency = 2.5MHz, speed of sound = 1500m/s, piston width a = 52, piston height b = 102, where the piston is centered at the origin and located in the xy plane. Sample the aperture with nxprime = 30 and nyprime = 60. Evaluate the pressure in the xy plane at z= (a+b)/(202). Calculate the pressure on a 61 by 61 pressure grid, where xmax = -xmin = 1.5 *a/2 and ymax = -ymin = 1.5 *b/2. Plot the absolute value of the result with the mesh command