Matlab in Chemical Engineering at CMU

Contents

Interpolation of data

when we have data at two points but we need data in between them we use interpolation. Suppose we have the points (4,3) and (6,2) and we want to know the value of y at x=4.65, assuming y varies linearly between these points. we use the interp1 command to achieve this.

x = [4 6]
y = [3 2]

interp1(x,y,4.65)

x =

     4     6


y =

     3     2


ans =

    2.6750

multiple interpolation values

you can interpolate several values

interp1(x,y,[4.65 5.01 4.2 9])
% note that you cannot interpolate outside the region where data is
% defined. That is why the value at x=9 is NaN.

ans =

    2.6750    2.4950    2.9000       NaN

more sophisticated interpolation.

the default interpolation method is simple linear interpolation between points. Other methods exist too, such as fitting a spline to the data and using the spline representation to interpolate from.

x = [1 2 3 4];
y = [1 4 9 16]; % y = x^2

xi = [ 1.5 2.5 3.5]; % we want to interpolate on these values
y1 = interp1(x,y,xi);
y2 = interp1(x,y,xi,'spline')

y2 =

    2.2500    6.2500   12.2500

Compare these methods by plotting each set of interpolated values

the data we interpolated was constructed from y = x^2

figure
hold on
ezplot('x^2',[1 4])
plot(xi,y1,'rs',xi,y2,'bd')
legend('data','linear interpolation','spline interpolation')
hold off

summary: in this case the spline interpolation is a little more accurate than the linear interpolation. That is because the underlying data was polynomial in nature, and a spline is like a polynomial. That may not always be the case, and you need some engineering judgement to know which method is best.

% categories: Basic Matlab
% tags: math

Published with MATLAB® 7.10

Contents

Integrating functions in Matlab

it is a good idea to add these lines to the top of your m-files. They keep things looking neat!

close all  % close all figure windows that are open
clear all  % clear all the variables currently stored in memory
clc        % clear the commands in the command window

Problem statement

find the integral of a function f(x) from a to b i.e.

In matlab we use numerical quadrature to achieve this with the `quad` command. f = @(x) some_function_of_x quad(f,a,b)

as a specific example, lets integrate

from x=0 to x=1. You should be able to work out that the answer is 1/3.

f = @(x) x.^2;  %note we must use .^ to ensure element-wise squaring
a = 0;
b = 1;

integrand = quad(f,a,b)

integrand =

    0.3333

double integrals

we use the `dblquad` command

Integrate over

i.e.

f = @(x,y) y*sin(x)+x*cos(y);
x1=pi; x2=2*pi;
y1=0; y2=pi;
integrand = dblquad(f, x1,x2,y1,y2)

integrand =

   -9.8696

triple integrals

we use the `triplequad` command Integrate over the region

f = @(x,y,z) y*sin(x)+z*cos(x);
x1=0; x2=pi;
y1=0; y2=1;
z1=-1; z2=1;
integrand = triplequad(f, x1,x2,y2,y2,z1,z2)

% categories: Basic matlab
% tags: math

integrand =

     0

Published with MATLAB® 7.10