contents.m

数值方法和MATLAB实现与应用.zip

% Numerical Methods with MATLAB Toolbox
% Version 1.04  02-Feb-2003
%
%  MATLAB Routines and data to accompany the book
%  ``Numerical Methods with MATLAB: Implemenations and Applications''
%  Copyright 1996-2001, Gerald W. Recktenwald.  The contents of the
%  NMM Toolbox may be copied and used without charge.  Repackaging or
%  distribution in any form for a fee is prohibited.
%
% ----------------------------------------------------------------
%
% Contents of "data" directory in the NMM toolbox
%
% airSat.dat      Saturation pressure for air as a function of temperature
%                 Pressure (MPa) in column 1, Temperature (K) in column 2
% airvisc.dat     Viscosity of air verses temperature.  Temperature (degrees C)
%                 in column 1, viscosity (kg/m/s) in column2
% bearing.dat     Wear of journal bearing as function of temperature
% capacitor.dat   Capacitor voltage versus time; Data used to create Figure 9.2
%                 Time (s) in column 1, Voltage (V) in column 2
% chip.dat        Yield data for semiconductors.  First column is serial number
%                 second column is speed (MHz) at which chip qualifies.
%                 NaN entries in second column means chip failed all speed tests.
% CorvRain.dat    Monthly total precipitation data for Corvallis, Oregon, 1890 to 1994
%                 Column 1 is the year, column 2 through 13 is monthly total
%                 precipitation in (1/100)^th of an inch
% cucon1.dat ... cucon3.dat  Thermal conductivity data for Copper versus temperature.
%                 Temperature (C) in column 1, Thermal conductivity (W/m/C) in column 2
% emission.dat    Nitrogen oxide emssion from an internal combustion engine as
%                 function of humidity and atmospheric pressure
% fan9v.dat ... fan13v.dat  Fan curve data at different voltages
% flowsys1.dat    System curve data: flow rate (m^3/s) in column 1, pressure
%                 drop (Pa) in column 2
% flowsys2.dat    System curve data: flow rate (m^3/s) in column 1, pressure
%                 drop (Pa) in column 2
% gc87.dat        Flow data for Glen Canyon Dam in 1987.  First column is hour
%                 second column is flow rate in ft^3/s
% gc87flow.dat    Flow data for Glen Canyon Dam in 1987.  Flow data only in one
%                 column (ft^3/s).  Each row is an hour in the year.
% glycerin.dat    Thermophysical properties of glycerin
% H2Odensity.dat  Density of liquid water as a function of temperature.
% H2Osat.dat      Saturation pressure of water verses temperature.
% H2Ovisc.dat     Viscosity of water verses temperature.  Temperature (degrees C)
%                 in column 1, viscosity (kg/m/s) in column2
% Jtcouple.dat    Calibration data for J type thermocouple wire in the range -50 <= T <= 250 C
% pdxPass.dat     Airline passengers from 1986 to 1996 for the Portland, OR airport
%                 Year in column 1, millions of passengers in column 2
% pdxTemp.dat     Historical average of monthly temperature variation in Portland, OR 
%                 Month in column 1, average high temperature in column 2,
%                 average average temperature in column 3, average low temperature
%                 in column 4.
% pdxThead.dat    Same data as pdxTemp.dat.  First row contains text column headings
% pdxTheadLong.dat  Same data as pdxTemp.dat.  First row contains a text description,
%                 second row contains column headings.
% SiC.dat         Selected properties of Silicon Carbide (SiC)
% sphereCd.dat    Drag coefficient for smooth spheres as a function of Reynolds number
% sprint.dat      100m dash times for C Lewis and B Johnson in 1987 World Championships
% stdatm.dat      Temperature, pressure, and density of the standard atmosphere as
%                 a function of elevation.  Additional documentation in the file.
% Tfield.dat      Temperature field from CFD model of flow over a circuit board
% traffic.dat     Set of times at which vehicles cross monitoring point on a highway
% velocity.dat    Vehicle velocity as function of time for coast-down test
% vprofile.dat    Velocity profile in a round pipe.
% wolfSun.dat     Wolfer sunspot index data with text column headings.  Year in column 1
%                 wolfer index in column two
% xinvpx.dat      Synthetic data used in curve fit using y = c1/x + c2*x 
% xtyt.dat        (x,y) data stored in rows:  first row is x, second row is y
% xy.dat          (x,y) data stored in columms:  first column is x, second column is y
% xy2.dat         (x,y1,y2) data stored in columms:  first column is x, second
%                 and third columns are y1 and y2
% xy5.dat         (x,y1,y2,y3,y4,y5) data stored in columms:  first column is x, second
%                 through sixth columns are y1, y2, ... y5
% xydy.dat        Data set used in least squares fitting; (x,y,deltay) in 3 columns
%
% ----------------------------------------------------------------
%
% Contents of "eigen" directory in the NMM toolbox
%
% eigSort     Eigenvalue/vectors sorted in ascending or descending order
% iterMult    Iterated multiplication of a vector by a matrix: u = A*A*...*A*x0
% poweritInv  Inverse power iterations with shift to find an eigenvalue of A
% powerit     Shifted power method for finding matrix eigenvalues
%
% ----------------------------------------------------------------
%
% Contents of "errors" directory in the NMM toolbox
%
% Archimedes     Perimeter of an n-sided polygon inscribed in a circle.
% bin2flt        Expand a binary representation of floating point mantissa
% demoTaylor     Accuracy of a Taylor Series approximations for f(x) = 1/(1-x)
% epprox         Demonstrate catastrophic cancellation in evaluation of
%                e = exp(1) = lim_{n->infinity} (1 + 1/n)^n
% equalTest      Script to demonstrate equality of two floating point numbers
% expSeriesPlot  Evaluate and plot series representation of exp(x)
% fidiff         Evaluate 1st order finite-difference approximation to d/dx of exp(x)
% halfDiff       Reduce the distance between two numbers until it is set to zero
% linsp1         Generate a vector of equally spaced values, version 1
% linsp2         Generate a vector of equally spaced values, version 2
% linsp3         Generate a vector of equally spaced values, version 3
% lintest        Compare schemes for generating a vector of equally spaced values
% newtsqrtBlank  Incomplete function to use Newton's method to compute the square root of a number
% sinser         Evaluate the series representation of sin(x)
% testSqrt       Test the newtsqrt function for a range of inputs
%
% ----------------------------------------------------------------
%
% Contents of "fit" directory in the NMM toolbox
%
% conductFit   LS fit of conductivity data for Copper at low temperatures
% cuconBasis1  Basis fcns for conductivity model:  1/k = c1/T + c2*T^2
% cuconBasis2  Basis fcns for conductivity model:  1/k = c1/T + c2*T + c3*T^2
% demoFanCurve Multivariate fit of fan data:  dp = f(q,v)
% demoH2OSat   Fit saturation pressure versus temperature for water
% demoPlaneFit Least squares fit of data to a plane: z = f(x,y)
% demoSiCmod   Least squares fit of bulk modulus of SiC versus temperature
% demoTcouple  Linear and quadratic fits to thermocouple data with polyfit
% demoXexp     Demonstrate fit of synthetic data to y = c1*x*exp(c2*x).
% fitnorm      Least squares fitting via solution to the normal equations
% fitqr        Least squares fitting via solution of overdetermined system with QR
% linefit      Least squares fit of data to y = a*x + b
% lineTest     Script that uses linefit to fit four data points to a line
% xexpfit      Least squares fit of data to y = c1*x*exp(c2*x)
% xinvpxBasis  Matrix with columns evaluated with 1/x and x
% xinvpxfit    Least squares fit of synthetic data to y = c1/x + c2*x
%
% ----------------------------------------------------------------
%
% Contents of "integrate" directory in the NMM toolbox
%
% adaptGK            Adaptive numerical integration using Gauss-Kronrod 7-15 rule
% adaptSimpson       Adaptive numerical integration based on Simpson's rule
% adaptSimpsonTrace  Adaptive numerical integration based on Simpson's rule, returns
%                    additional vector of x values where f(x) was evaluated
% compIntRules    Compare trapezoid, simpson, and Gauss-Legendre quadrature rules
% demoAdaptGK     Exercise adaptive Gauss-Kronrod quadrature and compare with quad8
% demoAdaptSimp   Integrate humps(x) with adaptive Simpson's rule
% demoGauss       Use Gauss-Legendre quadrature to integrate x*exp(-x) on [0,5]
% demoQuad        Use built in quad and quad8 to integrate 'humps' function on [0,5]
% demoSimp        Use composite Simpson's rule to integrate x*exp(-x) on [0,5]
% demoTrap        Use composite trapezoidal rule to integrate x*exp(-x) on [0,5]
% expmx2          Evaluate exp(-x^2), where x is a scalar or vector
% gaussKronrod15  Gauss-Kronrod quadrature pair of order 7 and 15
% gaussLagQuad    Gauss-Laguerre quadrature for integrals on [0,infinity)
% gaussQuad       Composite Gauss-Legendre quadrature
% GLagNodeWt      Nodes and weights for Gauss-Laguerre quadrature of arbitrary order
%                 by solving an eigenvalue problem
% GLagTable       Nodes and weights for Gauss-Laguerre quadrature of order n<=25
% GLNodeWt        Nodes and weights for Gauss-Legendre quadrature of arbitrary order
%                 obtained by solving an eigenvalue problem
% GLTable         Nodes and weights for Gauss-Legendre quadrature of order n<=8
% humpInt         Exact value of integral of humps function on [a,b]
% makeGLagTable   Create a table of Gauss-Laguerre nodes and weights
%                 suitable for copy/paste into the GLagTable.m function
% makeGLTable     Create a table of Gauss-Legendre nodes and weights
%                 suitable for copy/paste into the GLTable.m function
% plotSimpInt     Graphical display of composite Simpson rule integration
% plotTrapInt     Graphical display of composite trapezoid rule integration
% quadToInfinity  Integral from 0 to infinity evaluated as sum of integrals
%                 on subintervals of x axis.  Subintervals size increases
%                 geometrically from x = 0 to infinity.  Sum is terminated when
%                 contribution from subintervals is less than tolerance.
% recursiveIndent Demonstration of a recursive function
% simpson         Composite Simpson's rule
% sinxonx         Evaluate (sin(x)/x)^2  with attention to sin(0)/0 = 1
% trapezoid       Composite trapezoid rule
% trapzDat        Composite trapezoid rule for arbitrarily spaced discrete data
% trapzDatTest    Verify trapzDat function for different types of input
% xemx            Evaluate x exp(-x), where x is a scalar or vector
%
% ----------------------------------------------------------------
%
% Contents of "interact" directory in the NMM toolbox
%
% demoContour     Demonstration of contour plot
% demoSubplot     Demonstration of subplot with four sine functions
% demoSurfTypes   Different types of surface plots for z = 2 - x^2 - y^2
% sourceSinkPlot  Script to create source/sink plot
% tempPlot        Script to plot monthly temperature data in pdxTemp.dat
% tempStats       Load data from pdxTemp.dat and compute simple statistics
%                 on temperature data in the file
% Tfield          Load and plot temperature field stored in custom format text file
%
% ----------------------------------------------------------------
%
% Contents of "interp" directory in the NMM toolbox