# 582hw1 - the path of the marble(use plot3 to plot the path...

This preview shows page 1. Sign up to view the full content.

AMATH 582 Winter Quarter 2012 Homework 1: An ultrasound problem DUE: Friday, Jan. 20, 2012 You dog ﬂuﬀy swallowed a marble. The vet suspects that it has now worked its way into the intestines. Using ultrasound, data is obtained concerning the spatial variations in a small area of the intestines where the marble is suspected to be. Unfortunately, ﬂuﬀy keeps moving and the internal ﬂuid movement through the intestines generates highly noisy data. Do you want your dog to live? In order to save your dog’s life you must located and compute the trajectory of the marble. Go to the class webpage and download: Testdata.mat . This contains 20 rows of data for 20 diﬀerent measurements that were taken in time. 1. Through averaging of the spectrum, determine the frequency signature (center frequency) generated by the marble. 2. Filter the data around the center frequency determined above in order to denoise the data and determine
This is the end of the preview. Sign up to access the rest of the document.

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.

Ask a homework question - tutors are online