This preview shows page 1. Sign up to view the full content.
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 doesnt 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 doesnt 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
- Winter '11