General class information

Syllabus html pdf

Sample homework report

Tentative calendar

Tentative list of topics

Getting Matlab

Please visit W&M IT software support page and download Matlab from appropriate "Math & Statistics Software" section. They have several available versions. Either one is fine. Since we are learning Matlab, we will not have time to go to fancy toolboxes which Matlab provides/removes with new releases.

Grader

Travis Horrom tshorrom@email.wm.edu

---

Week 01

Homework 01

See assignment in the following file hw01.pdf

Lecture 01 - Computers and programming language introduction

• Class goals and structure
• Computers introduction
• Programming languages overview
• Discretization - the main weakness of computers
• Binary representation

presentation and handouts

Lecture 02 - Matlab as a fancy calculator

• Matlab as a fancy calculator
• Efficient editing
• Help related commands
• Operator precedence

presentation and handouts

Week 02

Homework 02

See assignment in the following file hw02.pdf

do not forget to CC the grader

Lecture 03 - Matrices and plotting.

• Matrices
• Colon (:) operator
• Plotting
• Special array arithmetic operators

presentation and handouts

Lecture 04 - Boolean algebra, conditional statements, loops.

• Boolean algebra
• Comparison operators
• Conditional statements
• Loops
• Series

presentation and handouts

Lecture 05 - Functions and scripts

• Functions
• Scripts
• Saving and loading variables

presentation and handouts

Complimentary code

• hypotenuse.m, myfactorial.m, money_grows.m

Week 03

Homework 03

See assignment in the following file hw03.pdf

do not forget to CC the grader

Lecture 06 - Root finding

• Root finding problem
• Bisection method
• False position method
• Solution convergence

presentation and handouts

Complimentary code

• bisection.m

Lecture 07 - Root finding continued

• Secant method
• Newton-Raphson method
• Numerical derivative of a function
• Ridders method
• Root finding algorithms gotchas
• Root finding algorithms summary

presentation and handouts

Week 04

Homework 04

See assignment in the following file hw04.pdf

Lecture 08 - Sorting

• Bubble sort method
• Quick sort method

presentation and handouts

See what senator (at that time) Obama has to say about the bubble sort method. Skip to the 33rd second.

Complimentary code

• qsort.m

Lecture 09 - Sorting continued

• Heap
• Heap sorting method
• Vector heap representation

presentation and handouts

Week 05

Homework 05

See assignment in the following file hw05.pdf

Lecture 10 - Numerical integration

• Integration problem statement
• The Rectangle method
• Trapezoidal method
• Simpson method

presentation and handouts

Lecture 11 - Numerical integration continued

• Monte Carlo Integration
• Numerical Integration methods error estimates

presentation and handouts

Complimentary code

• picalc.m

Midterm project 01

See the project assignment

Due date: Monday October 4th of 2010 at 1pm.

Please, note that the wording and problem is modified compared to the Wednesday announcement. You still work for NASA, though.

Week 06

Lecture 12 - Random number generators and random processes

• Statistic and probability intro
• Random number generators

presentation and handouts

Complimentary code

lcgrand.m and check_lcgrand.m

Lecture 13 - Monte Carlo simulations

• Coin flipping bets
• The Monty Hall Paradox
• Disease spread

presentation and handouts

Complimentary code

• bet_outcome.m,
• gain_vs_bet_fraction.m,
• doors_game.m,
• doors_game_check.m,
• disease_outcome.m,
• colony_life.m.

Week 07

Homework 06

See assignment in the following file hw06.pdf

Lecture 14 - Optimization problem

• Introduction to optimization
• Matlab built in 1D minimization - fminbnd
• 1D optimization - golden section search

presentation and handouts

Complimentary code

• Calculation of a wire temperature as function of the insulation thickness wire_t.m

Lecture 15 - Multi-D optimization problem

• Multi-D optimization
• Matlab built in Multi-D minimization - fminsearch
• Sample problems

presentation and handouts

Complimentary code

• Simple function which makes a 3D plot of a function of 2 variables plotme3dFofXY.m
• Usage example plot_sinc_3d.m of above function which plots a sample function z=-sinc(r) where r=sqrt(x*x+y*y) fsample_sinc.m.
• Merrit function for the suspended weights problem EconstrainedSuspendedWeights.m

Midterm project 02

See the project assignment

Due date: Monday October 25th of 2010 at 1pm.

Week 08

Lecture 16 - Combinatorial optimization

• Introduction to the combinatorial problem
• Brute force optimization

presentation and handouts

Complimentary code

• Recursive permutation generation permutation.m.
• Combinatoric backpack problem solution backpack.m

Lecture 17 - Simulated annealing/Metropolis optimization

• Backpack problem - brute force approach
• Simulated annealing/Metropolis algorithm

presentation and handouts

Complimentary code

• Backpack problem optimization with binary addition combinatorial search backpack_binary.m
• Backpack problem optimization with Metropolis algorithm backpack_metropolis.m
• Backpack problem solution comparison backpack_check.m

Week 09

Homework 07

See assignment in the following file hw07.pdf

Due date: Monday November 1st of 2010 at 1pm.

Homework 07 data files

• cities_for_combinatorial_search.dat
• cities_for_metropolis_search.dat
• hw7p3.dat
• hw7p4.dat

Lecture 18 - Data reduction and fitting

• Traveling salesman problem
• Data reduction

presentation and handouts

Complimentary code

• Combinatorical solution of travelling salesman traveller_comb.m
• Cities coordinates data set cities.dat
• Fitter function fitter.m
• Fitter function usage example fitter_usage_example.m
• Data file for the above example data_to_fit.dat

Lecture 19 - Data analysis example

Complimentary code

• The main analysis script day_temp_analysis_session.m
  • To execute it in step wise fashion use: echodemo day_temp_analysis_session
• Data file used by the above script yearly_outside_temperatures.dat
• Model of the yearly variations temp_vs_day_simple.m
• Model of the dayly variations temp_vs_day_time.m
• Data averaging with reduction reduce_with_average.m

Week 10

Homework 08

See assignment in the following file hw08.pdf

Due date: Monday November 8th of 2010 at 1pm.

Lecture 20 - Ordinary Differential equations

• ODE - ordinary differential equations
• Euler's method

presentation and handouts

Complimentary code

• Euler's method ODE solver odeeuler.m
• Example of the free fall numerical solution with Euler's method ode_free_fall_example.m
  • free fall forces calculation code free_fall_forces.m
• Example of stability issues in the ODE solvers ode_unstable_example.m
  • forces calculation code ode_unstable_forces.m

Lecture 21 - Ordinary Differential equations continued

• The second-order Runge-Kutta method (RK2)
• The forth-order Runge-Kutta method (RK4)
• Matlab functions for ODEs

presentation and handouts

Complimentary code

• parachute jumper model (time dependent force) ode_parachute_jump_model.m
• predator and prey model shows oscillations ode_predator_prey_model.m
• projectile motion with air drag ode_projectile_with_air_drag_model.m

Midterm project 03

See the project assignment

Due date: Monday November 15th of 2010 at 1pm.

• data file solar_system_data.mat

Week 11

Lecture 22 - System of linear algebraic equations

• Matlab solvers

presentation and handouts

Complimentary code

• the mobile problem mobile.m
• the Wheatstone bridge problem wheatstone_bridge.m

Lecture 23 - Data interpolation

• Nearest neighbor interpolation
• Linear interpolation
• Polynomial fit
• Cubic spline interpolation
• Matlab built in interpolation
• Extrapolation

presentation and handouts

Complimentary code

• nearest neighbor interpolation interp_nearest.m
• linear interpolation interp_linear.m

Week 12

Homework 09

See assignment in the following file hw09.pdf

• data file for the problem hw09pr3.dat.

Lecture 24 - Fourier transform

• Fourier series
• Discrete Fourier transform (DFT)
• Matlab implementation

presentation and handouts

Complimentary code

• Fourier series for |t| fourier_abs.m
• Fourier series for step function fourier_step.m

Lecture 25 - Discrete Fourier Transform and filters

• Discrete Fourier transform
• Nyquist criteria
• DFT filters

presentation and handouts

Complimentary code

• filtersdemo.m
• fourier_frequencies.m
• freq_space_low_pass_filter.m
• voice_record.wav
• niceplotset.m

Week 13

Homework 10

See assignment in the following file hw10.pdf

Due date: Monday November 29th of 2010 at 1pm.

• This homework is optional.

Lecture 26 - High performance computing

• Limits to computing
• Parallel computing
• Matlab parallel tools

presentation and handouts

Complimentary code

• for loop example for_script.m
• parfor loop example parfor_script.m

Week 14

Lecture 27 - Digital filters

• Speed of FFT
• Low-pass filter
• High-pass filter
• Band-pass filter
• Band-stop filter
• Brick wall filters artifacts
• Windowing

presentation and handouts

Complimentary code

• voice record with monochromatic noise voice_record_l27.wav
• voice filtering demo filter_voice_demo_l27.m
• filter application help function filter_application.m

Lecture 28 - Fractals

• Dimension definition
• Fractals
• Affine transformations

presentation and handouts

Complimentary code

• sierpinski_random.m
• sierpinski_regular.m
• fern.m

Final exam project

See the final exam project assignment

Due date: Monday December 6th of 2010 at 2pm.

• sample data file FinalExam_example_of_the_problem_file1.mat with relevant parameters
  • trajectory of the missile for above input file missile_test_trajectory1.mat

Voice messages for problem 4. Each group must decode only the one assigned to them.

• noisy_msg1.wav for Andy and Caitlin
• noisy_msg2.wav for Daniel and Jake
• noisy_msg3.wav for David and Jonathan
• noisy_msg4.wav for Patrick and Kristin
• noisy_msg5.wav for Lauren, Kyle, and Ryan