Geoscience Reference
In-Depth Information
Mises distribution (Section 10.4).
We use two synthetic data sets of directional data to illustrate the
application of this test. We i rst load the data and convert the degrees to
radians.
clear
data_degrees_1 = load('directional_1.txt');
data_degrees_2 = load('directional_2.txt');
data_radians_1 = pi*data_degrees_1/180;
data_radians_2 = pi*data_degrees_2/180;
We then compute the lengths of the resultant vectors.
x_1 = sum(sin(data_radians_1));
y_1 = sum(cos(data_radians_1));
x_2 = sum(sin(data_radians_2));
y_2 = sum(cos(data_radians_2));
mean_radians_1 = atan(x_1/y_1);
mean_degrees_1 = 180*mean_radians_1/pi;
mean_radians_2 = atan(x_2/y_2);
mean_degrees_2 = 180*mean_radians_2/pi;
mean_degrees_1 = mean_degrees_1 + 180
mean_degrees_2 = mean_degrees_2 + 180
R_1 = sqrt(x_1^2 + y_1^2);
R_2 = sqrt(x_2^2 + y_2^2);
mean_degrees_1 =
223.2357
mean_degrees_2 =
200.8121
h e orientations of the resultant vectors are approximately 223° and 201°.
We also need the resultant length for both samples combined, so we combine
both data sets and again compute the resultant length.
data_radians_T = [data_radians_1;data_radians_2];
x_T = sum(sin(data_radians_T));
y_T = sum(cos(data_radians_T));
mean_radians_T = atan(x_T/y_T);
mean_degrees_T = 180*mean_radians_T/pi;
mean_degrees_T = mean_degrees_T + 180;
R_T = sqrt(x_T^2 + y_T^2)
Rm_T = R_T / (length(data_radians_T))
