Unformatted text preview: the path of the marble. (use plot3 to plot the path once you have it) 3. Where should an intense acoustic wave be focused to breakup the marble at the 20th data measurement. Good luck, and I hope your dog doesn’t die. The following code will help you get started in analyzing the data. It also tells you the spatial and spectral resolution of your ultrasound equipment. (NOTE: the reason for the close all command before isosurface is that isosurface doesn’t seem to clear the previous imagine before plotting a new one) clear all; close all; clc; load Testdata L=15; % spatial domain n=64; % Fourier modes x2=linspace(-L,L,n+1); x=x2(1:n); y=x; z=x; k=(2*pi/(2*L))*[0:(n/2-1) -n/2:-1]; ks=fftshift(k); [X,Y,Z]=meshgrid(x,y,z); [Kx,Ky,Kz]=meshgrid(ks,ks,ks); for j=1:20 Un(:,:,:)=reshape(Undata(j,:),n,n,n); close all, isosurface(X,Y,Z,abs(Un),0.4) axis([-20 20 -20 20 -20 20]), grid on, drawnow pause(1) end...
View Full Document
This note was uploaded on 03/23/2012 for the course AMATH 582 taught by Professor N.k during the Winter '11 term at University of Washington.
- Winter '11