Download The draft users manual (pdf

CHAPTER 4. EXAMPLES
36
Terminal egg velocity [mm/s]
9
6
1.5
6
8
5
7
5
3
Density difference [kg/m3]
0.3
4
2
0.5
4
3
2
2.5
1
0
0.1
0
0.5
1
1.5
2
Diameter [mm]
0.4
2.5
1
0.2
3
3.5
4
Figure 4.9: Terminal velocity contours
table 1.1 by the method of least squares. Some array manipulation must first be done to
obtain three column vectors SI, TI and B with corresponding values of salinity, temperature
and viscosity. With MN = 35 = 5 × 7, the normal matrix is simply
A = [ones(MN) TI
TI.^2 SI];
In MATLAB the vector X of coefficients is found by
V = eye(MN);
X = lscov(A,B,V);
% Covariance matrix = identity
in Octave the same task is done by
X = ols(B,A);
4.5
Example 5, Non-constant coefficients
This example illustrates how to use several egg groups on a transient problem with nonconstant coefficients. It also provides a verification of parts of the toolbox by reproducing
results by Westgård (1989). The model is contained in the script runex5. The results are
visualized by the scripts plotex5 and by animation in animex5.
The problem is the first example from Westgård (1989). The bottom depth is 100 m
and a grid size of 2 m is used.