Scripting 5.7. Vectors and matrices

### 1. What is a vector?

A numeric vector is an array of numbers, regarded as a single object. For instance, the squares of the first five integers can be collected in the vector { 1, 4, 9, 16, 25 }. In a Praat script, you can put a vector into a variable whose name ends in a number sign ("#"):

``` squares# = { 1, 4, 9, 16, 25 }```

After this, the variable squares# contains the value { 1, 4, 9, 16, 25 }. We say that the vector squares# has five dimensions, i.e. it contains five numbers.

Whereas in Scripting 3.2. Numeric variables we talked about a numeric variable as being analogous to a house where somebody (the numeric value) could live, a numeric vector with five dimensions can be seen as a street that contains five houses, which are numbered with the indexes 1, 2, 3, 4 and 5, each house containing a numeric value. Thus, the street squares# contains the following five houses: squares# [1], squares# [2], squares# [3], squares# [4] and squares# [5]. Their values (the numbers that currently live in these houses) are 1, 4, 9, 16 and 25, respectively.

To list the five values with a loop, you could do:

``` writeInfoLine: "Some squares:"``` ``` for i from 1 to 5``` ```    appendInfoLine: "The square of ", i, " is ", squares# [i]``` ``` endfor```

Instead of the above procedure to get the vector squares#, with a pre-computed list of five squares, you could compute the five values with a formula, as in the example of Scripting 5.6. Arrays and dictionaries. However, in order to put a value into an element of the vector, you have to create the vector first (i.e., you have to build the whole street before you can put something in a house), so we start by creating a vector with five zeroes in it:

``` squares# = zero# (5)```

After this, squares# is the vector { 0, 0, 0, 0, 0 }, i.e., the value of each element is zero. Now that the vector (street) exists, we can put values into (populate) the five elements (houses):

``` for i from 1 to 5``` ```    squares# [i] = i * i``` ``` endfor```

After this, the variable `squares#` has the value { 1, 4, 9, 16, 25 }, as before, but now we had the computer compute the squares.

### 2. Creating a vector

You can create a vector in many ways. The first way we saw was with a vector literal, i.e. a series of numbers (or numeric formulas) between braces:

``` lengths# = { 1.83, 1.795, 1.76 }```

The second way we saw was to create a series of zeroes. To create a vector consisting of 10,000 zeroes, you do

``` zero# (10000)```

Another important type of vector is a series of random numbers. To create a vector consisting of 10,000 values drawn from a Gaussian distribution with true mean 0.0 and true standard deviation 1.0, you could do

``` noise# = randomGauss# (10000, 0.0, 1.0)```

To create a vector consisting of 10,000 values drawn from a uniform distribution of real numbers with true minimum 0.0 and true maximum 1.0, you use

``` randomUniform# (10000, 0.0, 1.0)```

To create a vector consisting of 10,000 values drawn from a uniform distribution of integer numbers with true minimum 1 and true maximum 10, you use

``` randomInteger# (10000, 1, 10)```

Vectors can also be created by some menu commands. For instance, to get vectors representing the times and pitch frequencies of the frames in a Pitch object, you can do

``` selectObject: myPitch``` ``` times# = Get times of frames``` ``` pitches# = Get values in frames```

### 3. Turning a vector into a number

For the vector defined above, you can compute the sum of the five values as

``` sum (squares#)```

which gives 55. You compute the average of the five values as

``` mean (squares#)```

which gives 11. You compute the standard deviation of the values as

``` stdev (squares#)```

which gives 9.669539802906858 (the standard deviation is undefined for vectors with fewer than 2 elements). The center of gravity of the distribution defined by regarding the five values as relative frequencies as a function of the index from 1 to 5 is computed by

``` center (squares#)```

which gives 4.090909090909091 (for a vector with five elements, the result will always be a number between 1.0 and 5.0). You compute the inner product of two equally long vectors as follows:

``` other# = { 2, 1.5, 1, 0.5, 0 }``` ``` result = inner (squares#, other#)```

which gives 1*2 + 4*1.5 + 9*1 + 16*0.5 + 25*0 = 25. The formula for this is ∑i=15 squares[i] * other[i], so that an alternative piece of code could be

``` result = sumOver (i to 5, squares# [i] * other# [i])```

### 4. Converting vectors to vectors

``` a# = squares# + 5 ; adding a number to each element of a vector```

causes a# to become the vector { 6, 9, 14, 21, 30 }.

``` b# = a# + { 3.14, 2.72, 3.16, -1, 7.5 } ; adding two vectors of the same length```

causes b# to become the vector { 9.14, 11.72, 17.16, 20, 37.5 }.

``` c# = b# / 2 ; dividing each element of a vector```

causes c# to become the vector { 4.57, 5.86, 8.58, 10, 18.75 }.

``` d# = b# * c# ; elementwise multiplication```

causes d# to become the vector { 41.7698, 68.6792, 147.2328, 200, 703.125 }.

A vector can also be given to a menu command that returns another vector. For instance, to get a vector representing the pitch frequencies at 0.01-second intervals in a Pitch object, you can do

``` selectObject: myPitch``` ``` tmin = Get start time``` ``` tmax = Get end time``` ``` times# = sequence_by_centre# (tmin, tmax, 0.01)``` ``` pitches# = Get values at times: times#, "hertz", "linear"```