mirror of
https://github.com/rasbt/python_reference.git
synced 2024-11-30 15:31:12 +00:00
1443 lines
115 KiB
Markdown
1443 lines
115 KiB
Markdown
<link rel="stylesheet" type="text/css" href="../CSS/hover.css" media="screen" />
|
||
|
||
|
||
|
||
|
||
# Numeric matrix manipulation - The cheat sheet for MATLAB, Python NumPy, R, and Julia
|
||
|
||
At its core, this article is about a simple cheat sheet for basic operations on numeric matrices, which can be very useful if you working and experimenting with some of the most popular languages that are used for scientific computing, statistics, and data analysis.
|
||
|
||
<a id='sections'></a>
|
||
|
||
## Sections
|
||
|
||
- [Introduction](#introduction)
|
||
|
||
- [Language overview](#overview)
|
||
|
||
- [MATLAB/Octave](#matlab)
|
||
|
||
- [Python NumPy](#numpy)
|
||
|
||
- [R](#r)
|
||
|
||
- [Julia](#julia)
|
||
|
||
- [Matrix cheat sheet](#cheatsheet)
|
||
|
||
- [Alternative data structures: NumPy matrices vs. NumPy arrays](#numpy_arrays)
|
||
|
||
|
||
|
||
<br>
|
||
|
||
|
||
<a id='introduction'></a>
|
||
|
||
## Introduction
|
||
|
||
[[back to section overview](#sections)]
|
||
<br>
|
||
|
||
Matrices (or multidimensional arrays) are not only presenting the fundamental elements of many algebraic equations that are used in many popular fields, such as pattern classification, machine learning, data mining, and math and engineering in general. But in context of scientific computing, they also come in very handy for managing and storing data in an more organized tabular form.
|
||
Such multidimensional data structures are also very powerful performance-wise thanks to the concept of automatic vectorization: instead of the individual and sequential processing of operations on scalars in loop-structures, the whole computation can be parallelized in order to make optimal use of modern computer architectures.
|
||
|
||
|
||
<br>
|
||
|
||
|
||
|
||
<br>
|
||
<br>
|
||
|
||
<font size=1em>
|
||
<hr>
|
||
**Note:**
|
||
This article originated from an older article with containing a cheat sheet that was just about MATLAB matrices and NumPy arrays. Since then, I added a couple of more rows and doubled the width of the cheat sheet by adding those two other languages R and Julia. Instead of making further modifications, I wanted to keep this old article as is - for future reference and for people who may only be interested in this slimmer version:
|
||
[Moving from MATLAB matrices to NumPy arrays - A Matrix Cheatsheet](http://sebastianraschka.com/Articles/2014_matlab_vs_numpy.html).
|
||
<hr>
|
||
</font>
|
||
|
||
<br>
|
||
<br>
|
||
|
||
<a id='overview'></a>
|
||
|
||
### Language overview
|
||
|
||
[[back to section overview](#sections)]
|
||
<br>
|
||
|
||
Before we **[jump to the actual cheat sheet](#cheatsheet)**, I wanted to give you at least a brief overview of the different languages that we are dealing with.
|
||
|
||
|
||
All four languages, MATLAB/Octave, Python, R, and Julia are dynamically typed, have a command line interface for the interpreter, and come with great number of additional and useful libraries to support scientific and technical computing. Conveniently, these languages also offer great solutions for easy plotting and visualizations.
|
||
|
||
Combined with interactive notebook interfaces or dynamic report generation engines ([MuPAD](http://www.mathworks.com/discovery/mupad.html) for MATLAB, [IPython Notebook](http://ipython.org/notebook.html) for Python, [knitr](http://yihui.name/knitr/) for R, and [IJulia](https://github.com/JuliaLang/IJulia.jl) for Julia based on IPython Notebook) data analysis and documentation has never been easier.
|
||
|
||
<br>
|
||
<br>
|
||
|
||
<a id="matlab"></a>
|
||
|
||
# MATLAB/Octave
|
||
|
||
[[back to section overview](#sections)]
|
||
<br>
|
||
|
||
<br>
|
||
<a href="http://www.mathworks.com/products/matlab/" class="hover_dark_light">
|
||
|
||
</a>
|
||
|
||
<br>
|
||
|
||
[MATLAB](http://www.mathworks.com/products/matlab/) (stands for MATrix LABoratory) is the name of an application and language that was developed by [MathWorks](http://www.mathworks.com/index.html?s_tid=gn_logo) back in 1984. One of its strengths is the variety of different and highly optimized "toolboxes" (including very powerful functions for image and other signal processing task), which makes suitable for tackling basically every possible science and engineering task.
|
||
Like the other languages, which will be covered in this article, it has cross-platform support and is using dynamic types, which allows for a convenient interface, but can also be quite "memory hungry" for computations on large data sets.
|
||
|
||
Even today, MATLAB is probably (still) the most popular language for numeric computation used for engineering tasks in academia as well as in industry.
|
||
|
||
#### GNU Octave
|
||
|
||
<br>
|
||
<a href="http://www.gnu.org/software/octave/" class="hover_dark_light">
|
||
|
||
</a>
|
||
|
||
<br>
|
||
|
||
It is also worth mentioning that MATLAB is the only language in this cheat sheet which is not free and open-sourced. But since it is so immensely popular, I want to mention it nonetheless. And as an alternative there is also the free [GNU Octave re-implementation](http://www.gnu.org/software/octave/) that follows the same syntactic rules so that the code is compatible to MATLAB (except for very specialized libraries).
|
||
|
||
<br>
|
||
|
||
*<font size=1em> This [image](http://commons.wikimedia.org/wiki/File:Matlab_Logo.png) is a freely usable media under public domain and represents the first eigenfunction of the L-shaped membrane, resembling (but not identical to) MATLAB's logo trademarked by MathWorks Inc.</font>
|
||
|
||
<br>
|
||
<br>
|
||
|
||
<a id="numpy"></a>
|
||
|
||
|
||
|
||
|
||
# Python NumPy
|
||
|
||
[[back to section overview](#sections)]
|
||
<br>
|
||
|
||
<br>
|
||
<a href="http://www.numpy.org" class="hover_dark_light"/>
|
||
|
||
</a>
|
||
|
||
<br>
|
||
|
||
Initially, the [NumPy](http://www.numpy.org) project started out under the name "Numeric" in 1995 (renamed to NumPy in 2006) as a Python library for numeric computations based on multi-dimensional data structures, such as arrays and matrices. Since it makes use of pre-compiled C code for operations on its "`ndarray`" objects, it is considerably faster than using equivalent approaches in (C)Python.
|
||
|
||
|
||
Python NumPy is my personal favorite since I am a big fan of the Python programming language. Although similar tools exist for other languages, I found myself to be most productive doing my research and data analyses in [IPython notebooks](http://ipython.org/notebook.html).
|
||
It allows me to easily combine Python code (sometimes optimized by compiling it via the [Cython](http://cython.org) C-Extension or the just-in-time (JIT) [Numba](http://numba.pydata.org) compiler if speed is a concern) with different libraries from the [Scipy stack](http://www.scipy.org/) including [matplotlib](http://matplotlib.org) for inline data visualization (you can find some of my example benchmarks in this [GitHub repository](http://github.com/rasbt/One-Python-benchmark-per-day)).
|
||
|
||
<br>
|
||
<br>
|
||
|
||
<a id="r"></a>
|
||
|
||
# R
|
||
|
||
[[back to section overview](#sections)]
|
||
<br>
|
||
|
||
<br>
|
||
<a href="http://www.r-project.org" class="hover_dark_light"/>
|
||
|
||
</a>
|
||
<br>
|
||
|
||
The [R](http://www.r-project.org) programming language was developed in 1993 and is a modern GNU implementation of an older statistical programming language called [S](http://stat.bell-labs.com/S/), which was developed in the [Bell Laboratories](http://stat.bell-labs.com) in 1976. Since its release, it has a fast-growing user base and is particularly popular among statisticians.
|
||
|
||
R was also the first language which kindled my fascination for statistics and computing. I have used it quite extensively a couple of years ago before I discovered Python as my new favorite language for data analysis.
|
||
Although R has great in-built functions for performing all sorts statistics, as well as a plethora of freely available libraries developed by the large R community, I often hear people complaining about its rather unintuitive syntax.
|
||
|
||
<br>
|
||
<br>
|
||
|
||
<a id="julia"></a>
|
||
|
||
# Julia
|
||
|
||
[[back to section overview](#sections)]
|
||
<br>
|
||
|
||
<br>
|
||
<a href="http://julialang.org" class="hover_dark_light"/>
|
||
|
||
</a>
|
||
|
||
<br>
|
||
|
||
|
||
With its first release in 2012, [Julia](http://julialang.org) is by far the youngest of the programming languages mentioned in this article. a
|
||
While Julia can also be used as an interpreted language with dynamic types from the command line, it aims for high-performance in scientific computing that is superior to the other dynamic programming languages for technical computing thanks to its LLVM-based just-in-time (JIT) compiler.
|
||
|
||
Personally, I haven't used Julia that extensively, yet, but there are some exciting benchmarks that look very promising:
|
||
|
||
<a href="http://julialang.org/benchmarks/">
|
||
|
||
</a>
|
||
|
||
<font size=1em>
|
||
C compiled by gcc 4.8.1, taking best timing from all optimization levels (-O0 through -O3). C, Fortran and Julia use OpenBLAS v0.2.8. The Python implementations of rand_mat_stat and rand_mat_mul use NumPy (v1.6.1) functions; the rest are pure Python implementations.<br>
|
||
|
||
Bezanson, J., Karpinski, S., Shah, V.B. and Edelman, A. (2012), “Julia: A fast dynamic language for technical computing”.
|
||
(Source: [http://julialang.org/benchmarks/](http://julialang.org/benchmarks/), with permission from the copyright holder)
|
||
</font>
|
||
|
||
<br>
|
||
<br>
|
||
<br>
|
||
<a id="cheatsheet"></a>
|
||
<br>
|
||
|
||
|
||
|
||
# Cheat sheet
|
||
|
||
[[back to section overview](#sections)]
|
||
|
||
<a id='cheat_overview'></a>
|
||
<br>
|
||
|
||
|
||
|
||
|
||
### Cheat sheet overview
|
||
|
||
|
||
- [Creating matrices](#creating)
|
||
|
||
- [Accessing matrix elements](#accessing)
|
||
|
||
- [Manipulating shape](#manipulating)
|
||
|
||
- [Basic matrix operations](#basic_ops)
|
||
|
||
- [Advanced matrix operations](#advanced_ops)
|
||
|
||
<br>
|
||
<br>
|
||
<br>
|
||
<br>
|
||
|
||
|
||
**If you are interested in downloading this cheat sheet table for your references, you can find it [here on GitHub](https://github.com/rasbt/python_reference/blob/master/tutorials/matrix_cheatsheet_only.html)**
|
||
|
||
<br>
|
||
<br>
|
||
|
||
|
||
|
||
<table cellpadding="2" cellspacing="0" style="page-break-before: always">
|
||
<col width="383">
|
||
<col width="387">
|
||
<col width="447">
|
||
<col width="469">
|
||
<col width="511">
|
||
<col width="690">
|
||
<tr>
|
||
<td height="59" bgcolor="#0066cc" style="border: none; padding: 0in">
|
||
<p align="center"><font color="#ffffff"><font face="Arial"><font size="5" style="font-size: 18pt"><b>Task</b></font></font></font></p>
|
||
</td>
|
||
<td bgcolor="#0066cc" style="border: none; padding: 0in">
|
||
<p align="center"><font color="#ffffff"><font face="Andale Mono"><font size="5" style="font-size: 18pt"><b>MATLAB/Octave</b></font></font></font></p>
|
||
</td>
|
||
<td bgcolor="#0066cc" style="border: none; padding: 0in">
|
||
<p align="center"><font color="#ffffff"><font face="Andale Mono"><font size="5" style="font-size: 18pt"><b>Python
|
||
NumPy</b></font></font></font></p>
|
||
</td>
|
||
<td bgcolor="#0066cc" style="border: none; padding: 0in">
|
||
<p align="center"><font color="#ffffff"><font face="Andale Mono"><font size="5" style="font-size: 18pt"><b>R</b></font></font></font></p>
|
||
</td>
|
||
<td bgcolor="#0066cc" style="border: none; padding: 0in">
|
||
<p align="center"><font color="#ffffff"><font face="Andale Mono"><font size="5" style="font-size: 18pt"><b>Julia</b></font></font></font></p>
|
||
</td>
|
||
<td bgcolor="#0066cc" style="border: none; padding: 0in">
|
||
<p align="center"><font color="#ffffff"><font face="Arial"><font size="5" style="font-size: 18pt"><b>Task</b></font></font></font></p>
|
||
</td>
|
||
</tr>
|
||
<tr>
|
||
<td colspan="6" height="70" bgcolor="#ffffff" style="border: none; padding: 0in">
|
||
<p align="center"><font face="Arial"><font size="4" style="font-size: 14pt"><b>CREATING
|
||
MATRICES</b></font></font></p>
|
||
</td>
|
||
</tr>
|
||
<tr valign="top">
|
||
<td height="158" bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: none; border-right: 1px solid #000000; padding-top: 0in; padding-bottom: 0in; padding-left: 0in; padding-right: 0.02in">
|
||
<p align="left"><font face="Arial"><font size="2" style="font-size: 10pt"><b>Creating
|
||
Matrices <br>(here: 3x3 matrix)</b></font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: none; border-right: 1px dotted #000000; padding-top: 0in; padding-bottom: 0in; padding-left: 0in; padding-right: 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
A = [1 2 3; 4 5 6; 7 8 9]<br>A =<br> 1 2
|
||
3<br> 4 5 6<br> 7 8
|
||
9</font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: 1px dotted #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
A = np.array([ [1,2,3], [4,5,6], [7,8,9] ])<br><br><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
A<br>array([[1, 2, 3],<br> [4, 5, 6],<br>
|
||
[7, 8, 9]])</font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: 1px dotted #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
A = matrix(c(1,2,3,4,5,6,7,8,9),nrow=3,byrow=T)<br><br>
#
|
||
equivalent to
<br># A = matrix(1:9,nrow=3,byrow=T)
<br><br><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
A<br>[,1] [,2] [,3]<br>[1,] 1 2 3<br>[2,] 4 5 6<br>[3,] 7 8 9</font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: none; padding-top: 0in; padding-bottom: 0in; padding-left: 0.02in; padding-right: 0in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
A=[1 2 3; 4 5 6; 7 8 9]<br>3x3 Array{Int64,2}:<br>1 2 3<br>4 5 6<br>7
|
||
8 9</font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: 1px solid #000000; border-right: 1px solid #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Arial"><font size="2" style="font-size: 10pt"><b>Creating
|
||
Matrices <br>(here: 3x3 matrix)</b></font></font></p>
|
||
</td>
|
||
</tr>
|
||
<tr valign="top">
|
||
<td height="128" bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: none; border-right: 1px solid #000000; padding-top: 0in; padding-bottom: 0in; padding-left: 0in; padding-right: 0.02in">
|
||
<p align="left"><font face="Arial"><font size="2" style="font-size: 10pt"><b>Creating
|
||
an 1D column vector</b></font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: none; border-right: 1px dotted #000000; padding-top: 0in; padding-bottom: 0in; padding-left: 0in; padding-right: 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
a = [1; 2; 3]<br>a =<br> 1<br> 2<br>
|
||
3</font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: 1px dotted #000000; padding: 0in 0.02in">
|
||
<p align="left"><font color="#ff420e"><font face="Andale Mono"><font size="2" style="font-size: 10pt"><b>P></b></font></font></font>
|
||
<font face="Andale Mono"><font size="1" style="font-size: 8pt">a</font></font><font color="#2fff12">
|
||
</font><font color="#000000"><font face="Andale Mono"><font size="1" style="font-size: 8pt">=
|
||
np.array([1,2,3]).reshape(1,3)</font></font></font></p>
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><br></font></font><font color="#ff420e"><font face="Andale Mono"><font size="2" style="font-size: 10pt"><b>P></b></font></font></font>
|
||
<font color="#000000"><font face="Andale Mono"><font size="1" style="font-size: 8pt">b.shape<br>(1,
|
||
3)</font></font></font></p>
|
||
<p align="left"><br>
|
||
</p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: 1px dotted #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
a = matrix(c(1,2,3), nrow=3, byrow=T)<br><br><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
a<br>[,1]<br>[1,] 1<br>[2,] 2<br>[3,] 3</font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: none; padding-top: 0in; padding-bottom: 0in; padding-left: 0.02in; padding-right: 0in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
a=[1; 2; 3]<br>3-element Array{Int64,1}: <br>1<br>2<br>3</font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: 1px solid #000000; border-right: 1px solid #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Arial"><font size="2" style="font-size: 10pt"><b>Creating
|
||
an 1D column vector</b></font></font></p>
|
||
</td>
|
||
</tr>
|
||
<tr valign="top">
|
||
<td height="128" bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: none; border-right: 1px solid #000000; padding-top: 0in; padding-bottom: 0in; padding-left: 0in; padding-right: 0.02in">
|
||
<p align="left"><font face="Arial"><font size="2" style="font-size: 10pt"><b>Creating
|
||
an<br>1D row vector</b></font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: none; border-right: 1px dotted #000000; padding-top: 0in; padding-bottom: 0in; padding-left: 0in; padding-right: 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
b = [1 2 3]<br>b =<br> 1 2 3</font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: 1px dotted #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
b = np.array([1,2,3])<br><br><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
b<br>array([1, 2, 3])</font></font></p>
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt">#
|
||
note that numpy doesn't have <br># explicit “row-vectors”, but
|
||
1-D <br># arrays</font></font></p>
|
||
<p style="margin-bottom: 0in"><font face="Liberation Serif, serif"><font size="3" style="font-size: 12pt"><font color="#ff420e"><font face="Andale Mono"><font size="2" style="font-size: 10pt"><b>P></b></font></font></font><font color="#2fff12">
|
||
</font><font color="#000000"><font face="Andale Mono"><font size="1" style="font-size: 8pt">b.shape</font></font></font></font></font></p>
|
||
<p style="margin-bottom: 0in"><font face="Andale Mono"><font size="1" style="font-size: 8pt">(3,)</font></font></p>
|
||
<p align="left"><br>
|
||
</p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: 1px dotted #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
b = matrix(c(1,2,3), ncol=3)<br><br><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
b<br>[,1] [,2] [,3]<br>[1,] 1 2 3</font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: none; padding-top: 0in; padding-bottom: 0in; padding-left: 0.02in; padding-right: 0in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
b=[1 2 3] <br>1x3 Array{Int64,2}:<br>1 2 3<br><br># note that this
|
||
is a 2D array. <br># vectors in Julia are columns</font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: 1px solid #000000; border-right: 1px solid #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Arial"><font size="2" style="font-size: 10pt"><b>Creating
|
||
an<br>1D row vector</b></font></font></p>
|
||
</td>
|
||
</tr>
|
||
<tr valign="top">
|
||
<td height="128" bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: none; border-right: 1px solid #000000; padding-top: 0in; padding-bottom: 0in; padding-left: 0in; padding-right: 0.02in">
|
||
<p align="left"><font face="Arial"><font size="2" style="font-size: 10pt"><b>Creating
|
||
a <br>random m x n matrix</b></font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: none; border-right: 1px dotted #000000; padding-top: 0in; padding-bottom: 0in; padding-left: 0in; padding-right: 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
rand(3,2)<br>ans =<br> 0.21977 0.10220<br>
|
||
0.38959 0.69911<br> 0.15624 0.65637</font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: 1px dotted #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
np.random.rand(3,2)<br>array([[ 0.29347865, 0.17920462],<br>
|
||
[ 0.51615758, 0.64593471],<br>
|
||
[ 0.01067605, 0.09692771]])</font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: 1px dotted #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
matrix(runif(3*2), ncol=2)<br>[,1] [,2]<br>[1,] 0.5675127
|
||
0.7751204<br>[2,] 0.3439412 0.5261893<br>[3,] 0.2273177 0.223438</font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: none; padding-top: 0in; padding-bottom: 0in; padding-left: 0.02in; padding-right: 0in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
rand(3,2)<br>3x2 Array{Float64,2}:<br>0.36882 0.267725<br>0.571856
|
||
0.601524<br>0.848084 0.858935</font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: 1px solid #000000; border-right: 1px solid #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Arial"><font size="2" style="font-size: 10pt"><b>Creating
|
||
a <br>random m x n matrix</b></font></font></p>
|
||
</td>
|
||
</tr>
|
||
<tr valign="top">
|
||
<td height="128" bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: none; border-right: 1px solid #000000; padding-top: 0in; padding-bottom: 0in; padding-left: 0in; padding-right: 0.02in">
|
||
<p align="left"><font face="Arial"><font size="2" style="font-size: 10pt"><b>Creating
|
||
a<br>zero m x n matrix </b></font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: none; border-right: 1px dotted #000000; padding-top: 0in; padding-bottom: 0in; padding-left: 0in; padding-right: 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
zeros(3,2)<br>ans =<br> 0 0<br> 0
|
||
0<br> 0 0</font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: 1px dotted #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
np.zeros((3,2))<br>array([[ 0., 0.],<br>
|
||
[ 0., 0.],<br> [ 0.,
|
||
0.]])</font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: 1px dotted #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
mat.or.vec(3, 2)<br>[,1] [,2]<br>[1,] 0 0<br>[2,] 0 0<br>[3,] 0 0</font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: none; padding-top: 0in; padding-bottom: 0in; padding-left: 0.02in; padding-right: 0in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
zeros(3,2)<br>3x2 Array{Float64,2}:<br>0.0 0.0<br>0.0 0.0<br>0.0
|
||
0.0</font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: 1px solid #000000; border-right: 1px solid #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Arial"><font size="2" style="font-size: 10pt"><b>Creating
|
||
a<br>zero m x n matrix </b></font></font></p>
|
||
</td>
|
||
</tr>
|
||
<tr valign="top">
|
||
<td height="128" bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: none; border-right: 1px solid #000000; padding-top: 0in; padding-bottom: 0in; padding-left: 0in; padding-right: 0.02in">
|
||
<p align="left"><font face="Arial"><font size="2" style="font-size: 10pt"><b>Creating
|
||
an<br>m x n matrix of ones</b></font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: none; border-right: 1px dotted #000000; padding-top: 0in; padding-bottom: 0in; padding-left: 0in; padding-right: 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
ones(3,2)<br>ans =<br> 1 1<br> 1
|
||
1<br> 1 1</font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: 1px dotted #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
np.ones((3,2))<br>array([[ 1., 1.],<br>
|
||
[ 1., 1.],<br> [ 1., 1.]])</font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: 1px dotted #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
mat.or.vec(3, 2) + 1<br>[,1] [,2]<br>[1,] 1 1<br>[2,] 1 1<br>[3,]
|
||
1 1</font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: none; padding-top: 0in; padding-bottom: 0in; padding-left: 0.02in; padding-right: 0in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
ones(3,2)<br>3x2 Array{Float64,2}:<br>1.0 1.0<br>1.0 1.0<br>1.0
|
||
1.0</font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: 1px solid #000000; border-right: 1px solid #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Arial"><font size="2" style="font-size: 10pt"><b>Creating
|
||
an<br>m x n matrix of ones</b></font></font></p>
|
||
</td>
|
||
</tr>
|
||
<tr valign="top">
|
||
<td height="128" bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: none; border-right: 1px solid #000000; padding-top: 0in; padding-bottom: 0in; padding-left: 0in; padding-right: 0.02in">
|
||
<p align="left"><font face="Arial"><font size="2" style="font-size: 10pt"><b>Creating
|
||
an<br>identity matrix</b></font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: none; border-right: 1px dotted #000000; padding-top: 0in; padding-bottom: 0in; padding-left: 0in; padding-right: 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
eye(3)<br>ans =<br>Diagonal Matrix<br> 1 0
|
||
0<br> 0 1 0<br> 0 0
|
||
1</font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: 1px dotted #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
np.eye(3)<br>array([[ 1., 0., 0.],<br>
|
||
[ 0., 1., 0.],<br> [
|
||
0., 0., 1.]])</font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: 1px dotted #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
diag(3)<br>[,1] [,2] [,3]<br>[1,] 1 0 0<br>[2,] 0 1 0<br>[3,] 0 0
|
||
1</font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: none; padding-top: 0in; padding-bottom: 0in; padding-left: 0.02in; padding-right: 0in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
eye(3)<br>3x3 Array{Float64,2}:<br>1.0 0.0 0.0<br>0.0 1.0 0.0<br>0.0
|
||
0.0 1.0</font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: 1px solid #000000; border-right: 1px solid #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Arial"><font size="2" style="font-size: 10pt"><b>Creating
|
||
an<br>identity matrix</b></font></font></p>
|
||
</td>
|
||
</tr>
|
||
<tr valign="top">
|
||
<td height="217" bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: none; border-right: 1px solid #000000; padding-top: 0in; padding-bottom: 0in; padding-left: 0in; padding-right: 0.02in">
|
||
<p align="left"><font face="Arial"><font size="2" style="font-size: 10pt"><b>Creating
|
||
a<br>diagonal matrix</b></font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: none; border-right: 1px dotted #000000; padding-top: 0in; padding-bottom: 0in; padding-left: 0in; padding-right: 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
a = [1 2 3]<br><br><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
diag(a)<br>ans =<br>Diagonal Matrix<br> 1 0
|
||
0<br> 0 2 0<br> 0 0
|
||
3</font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: 1px dotted #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
a = np.array([1,2,3])<br><br><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
np.diag(a)<br>array([[1, 0, 0],<br> [0,
|
||
2, 0],<br> [0, 0, 3]])</font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: 1px dotted #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
diag(1:3)<br>[,1] [,2] [,3]<br>[1,] 1 0 0<br>[2,] 0 2 0<br>[3,] 0
|
||
0 3</font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: none; padding-top: 0in; padding-bottom: 0in; padding-left: 0.02in; padding-right: 0in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
a=[1, 2, 3] <br><br># added commas because julia<br># vectors are
|
||
columnar<br><br><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
diagm(a)<br>3x3 Array{Int64,2}:<br>1 0 0<br>0 2 0<br>0 0 3</font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: 1px solid #000000; border-right: 1px solid #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Arial"><font size="2" style="font-size: 10pt"><b>Creating
|
||
a<br>diagonal matrix</b></font></font></p>
|
||
</td>
|
||
</tr>
|
||
<tr>
|
||
<td colspan="6" height="128" bgcolor="#ffffff" style="border: none; padding: 0in">
|
||
<p align="center"><font face="Arial"><font size="4" style="font-size: 14pt"><b>ACCESSING
|
||
MATRIX ELEMENTS</b></font></font></p>
|
||
</td>
|
||
</tr>
|
||
<tr valign="top">
|
||
<td height="159" bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: none; border-right: 1px solid #000000; padding-top: 0in; padding-bottom: 0in; padding-left: 0in; padding-right: 0.02in">
|
||
<p align="left"><font face="Arial"><font size="2" style="font-size: 10pt"><b>Getting
|
||
the dimension<br>of a matrix<br>(here: 2D, rows x cols)</b></font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: none; border-right: 1px dotted #000000; padding-top: 0in; padding-bottom: 0in; padding-left: 0in; padding-right: 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
A = [1 2 3; 4 5 6]<br>A =<br> 1 2 3<br>
|
||
4 5 6<br><br><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
size(A)<br>ans =<br> 2 3</font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: 1px dotted #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
A = np.array([ [1,2,3], [4,5,6] ])<br><br><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
A<br>array([[1, 2, 3],<br> [4, 5,
|
||
6]])<br><br><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
A.shape<br>(2, 3)</font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: 1px dotted #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
A = matrix(1:6,nrow=2,byrow=T)<br><br><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
A<br>[,1] [,2] [,3]<br>[1,] 1 2 3<br>[2,] 4 5 6
</font></font></p>
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
dim(A)<br>[1] 2 3</font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: none; padding-top: 0in; padding-bottom: 0in; padding-left: 0.02in; padding-right: 0in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
A=[1 2 3; 4 5 6]<br>2x3 Array{Int64,2}:<br>1 2 3<br>4 5 6<br><br><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
size(A)<br>(2,3)</font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: 1px solid #000000; border-right: 1px solid #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Arial"><font size="2" style="font-size: 10pt"><b>Getting
|
||
the dimension<br>of a matrix<br>(here: 2D, rows x cols)</b></font></font></p>
|
||
</td>
|
||
</tr>
|
||
<tr valign="top">
|
||
<td height="240" bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: none; border-right: 1px solid #000000; padding-top: 0in; padding-bottom: 0in; padding-left: 0in; padding-right: 0.02in">
|
||
<p align="left"><font face="Arial"><font size="2" style="font-size: 10pt"><b>Selecting
|
||
rows </b></font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: none; border-right: 1px dotted #000000; padding-top: 0in; padding-bottom: 0in; padding-left: 0in; padding-right: 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
A = [1 2 3; 4 5 6; 7 8 9]<br><br>% 1st row<br><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
A(1,:)<br>ans =<br> 1 2 3<br><br>% 1st 2
|
||
rows<br><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
A(1:2,:)<br>ans =<br> 1 2 3<br>
|
||
4 5 6</font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: 1px dotted #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
A = np.array([ [1,2,3], [4,5,6], [7,8,9] ])<br><br># 1st row<br><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
A[0,:]<br>array([1, 2, 3])<br><br># 1st 2 rows<br><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
A[0:2,:]<br>array([[1, 2, 3], [4, 5, 6]])</font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: 1px dotted #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
A = matrix(1:9,nrow=3,byrow=T)
<br><br># 1st row<br>
<font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
A[1,]<br>[1] 1 2 3<br><br>
# 1st 2 rows
<br><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
A[1:2,]<br>[,1] [,2] [,3]<br>[1,] 1 2 3<br>[2,] 4 5 6</font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: none; padding-top: 0in; padding-bottom: 0in; padding-left: 0.02in; padding-right: 0in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
A=[1 2 3; 4 5 6; 7 8 9]; <br>#semicolon suppresses output<br><br>#1st
|
||
row<br><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
A[1,:] <br>1x3 Array{Int64,2}:<br>1 2 3<br><br>#1st 2 rows<br><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
A[1:2,:] <br>2x3 Array{Int64,2}:<br>1 2 3<br>4 5 6</font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: 1px solid #000000; border-right: 1px solid #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Arial"><font size="2" style="font-size: 10pt"><b>Selecting
|
||
rows </b></font></font></p>
|
||
</td>
|
||
</tr>
|
||
<tr valign="top">
|
||
<td height="304" bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: none; border-right: 1px solid #000000; padding-top: 0in; padding-bottom: 0in; padding-left: 0in; padding-right: 0.02in">
|
||
<p align="left"><font face="Arial"><font size="2" style="font-size: 10pt"><b>Selecting
|
||
columns</b></font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: none; border-right: 1px dotted #000000; padding-top: 0in; padding-bottom: 0in; padding-left: 0in; padding-right: 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
A = [1 2 3; 4 5 6; 7 8 9]<br><br>% 1st column<br><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
A(:,1)<br>ans =<br> 1<br> 4<br>
|
||
7<br><br>% 1st 2 columns<br><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
A(:,1:2)<br>ans =<br> 1 2<br> 4
|
||
5<br> 7 8</font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: 1px dotted #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
A = np.array([ [1,2,3], [4,5,6], [7,8,9] ])<br><br># 1st column
|
||
(as row vector)<br><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
A[:,0]<br>array([1, 4, 7])<br><br># 1st column (as column
|
||
vector)<br><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
A[:,[0]]<br>array([[1],<br> [4],<br>
|
||
[7]])<br><br># 1st 2 columns<br><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
A[:,0:2]<br>array([[1, 2], <br> [4,
|
||
5], <br> [7, 8]])</font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: 1px dotted #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
A = matrix(1:9,nrow=3,byrow=T)
<br><br># 1st column as row
|
||
vector
<br><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
t(A[,1])<br>[,1] [,2] [,3]<br>[1,] 1 4 7
<br><br># 1st column
|
||
as column vector
<br><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
A[,1]<br>[1] 1 4 7
<br><br># 1st 2 columns
<br><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
A[,1:2]<br>[,1] [,2]<br>[1,] 1 2<br>[2,] 4 5<br>[3,] 7 8</font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: none; padding-top: 0in; padding-bottom: 0in; padding-left: 0.02in; padding-right: 0in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
A=[1 2 3; 4 5 6; 7 8 9];<br><br>#1st column<br><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
A[:,1] <br>3-element Array{Int64,1}:<br>1<br>4<br>7<br><br>#1st 2
|
||
columns<br><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
A[:,1:2] <br>3x2 Array{Int64,2}:<br>1 2<br>4 5<br>7 8</font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: 1px solid #000000; border-right: 1px solid #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Arial"><font size="2" style="font-size: 10pt"><b>Selecting
|
||
columns</b></font></font></p>
|
||
</td>
|
||
</tr>
|
||
<tr valign="top">
|
||
<td height="260" bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: none; border-right: 1px solid #000000; padding-top: 0in; padding-bottom: 0in; padding-left: 0in; padding-right: 0.02in">
|
||
<p align="left"><font face="Arial"><font size="2" style="font-size: 10pt"><b>Extracting
|
||
rows and columns by criteria<br><br>(here: get rows that have
|
||
value 9 in column 3)</b></font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: none; border-right: 1px dotted #000000; padding-top: 0in; padding-bottom: 0in; padding-left: 0in; padding-right: 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
A = [1 2 3; 4 5 9; 7 8 9]<br>A =<br> 1 2
|
||
3<br> 4 5 9<br> 7 8
|
||
9<br><br><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
A(A(:,3) == 9,:)<br>ans =<br> 4 5 9<br>
|
||
7 8 9</font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: 1px dotted #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
A = np.array([ [1,2,3], [4,5,9], [7,8,9]])<br><br><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
A<br>array([[1, 2, 3],<br> [4, 5, 9],<br>
|
||
[7, 8, 9]])<br><br><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
A[A[:,2] == 9]<br>array([[4, 5, 9],<br>
|
||
[7, 8, 9]])</font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: 1px dotted #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
A = matrix(1:9,nrow=3,byrow=T)
<br><br><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
A<br>[,1] [,2] [,3]<br>[1,] 1 2 3<br>[2,] 4 5 9<br>[3,] 7 8
|
||
9
<br><br><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
matrix(A[A[,3]==9], ncol=3)<br>[,1] [,2] [,3]<br>[1,] 4 5 9<br>[2,]
|
||
7 8 9</font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: none; padding-top: 0in; padding-bottom: 0in; padding-left: 0.02in; padding-right: 0in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
A=[1 2 3; 4 5 9; 7 8 9]<br>3x3 Array{Int64,2}:<br>1 2 3<br>4 5 9<br>7
|
||
8 9<br><br># use '.==' for<br># element-wise check<br><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
A[ A[:,3] .==9, :] <br>2x3 Array{Int64,2}:<br>4 5 9<br>7 8 9</font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: 1px solid #000000; border-right: 1px solid #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Arial"><font size="2" style="font-size: 10pt"><b>Extracting
|
||
rows and columns by criteria<br><br>(here: get rows that have
|
||
value 9 in column 3)</b></font></font></p>
|
||
</td>
|
||
</tr>
|
||
<tr valign="top">
|
||
<td height="118" bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: none; border-right: 1px solid #000000; padding-top: 0in; padding-bottom: 0in; padding-left: 0in; padding-right: 0.02in">
|
||
<p align="left"><font face="Arial"><font size="2" style="font-size: 10pt"><b>Accessing
|
||
elements<br>(here: 1st element)</b></font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: none; border-right: 1px dotted #000000; padding-top: 0in; padding-bottom: 0in; padding-left: 0in; padding-right: 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
A = [1 2 3; 4 5 6; 7 8 9]<br><br><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
A(1,1)<br>ans = 1</font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: 1px dotted #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
A = np.array([ [1,2,3], [4,5,6], [7,8,9] ])<br><br><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
A[0,0]<br>1</font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: 1px dotted #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
A = matrix(c(1,2,3,4,5,9,7,8,9),nrow=3,byrow=T)
<br><br>
<font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
A[1,1]<br>[1] 1</font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: none; padding-top: 0in; padding-bottom: 0in; padding-left: 0.02in; padding-right: 0in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
A=[1 2 3; 4 5 6; 7 8 9];<br><br><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
A[1,1]<br>1</font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: 1px solid #000000; border-right: 1px solid #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Arial"><font size="2" style="font-size: 10pt"><b>Accessing
|
||
elements<br>(here: 1st element)</b></font></font></p>
|
||
</td>
|
||
</tr>
|
||
<tr>
|
||
<td colspan="6" height="128" bgcolor="#ffffff" style="border: none; padding: 0in">
|
||
<p align="center"><font face="Arial"><font size="4" style="font-size: 14pt"><b>MANIPULATING
|
||
SHAPE AND DIMENSIONS</b></font></font></p>
|
||
</td>
|
||
</tr>
|
||
<tr valign="top">
|
||
<td height="173" bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: none; border-right: 1px solid #000000; padding-top: 0in; padding-bottom: 0in; padding-left: 0in; padding-right: 0.02in">
|
||
<p align="left"><font face="Arial"><font size="2" style="font-size: 10pt"><b>Converting <br>row
|
||
to column vectors</b></font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: none; border-right: 1px dotted #000000; padding-top: 0in; padding-bottom: 0in; padding-left: 0in; padding-right: 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
b = [1 2 3]
<br><br><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
b = b'<br>b =<br> 1<br> 2<br>
|
||
3</font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: 1px dotted #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
b = np.array([1, 2, 3])<br><br><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
b = b[np.newaxis].T<br># alternatively <br># b =
|
||
b[:,np.newaxis]<br><br><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
b<br>array([[1],<br> [2],<br>
|
||
[3]])</font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: 1px dotted #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
b = matrix(c(1,2,3), ncol=3)<br><br><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
t(b)<br>[,1]<br>[1,] 1<br>[2,] 2<br>[3,] 3</font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: none; padding-top: 0in; padding-bottom: 0in; padding-left: 0.02in; padding-right: 0in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
b=vec([1 2 3])<br>3-element Array{Int64,1}:<br>1<br>2<br>3</font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: 1px solid #000000; border-right: 1px solid #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Arial"><font size="2" style="font-size: 10pt"><b>Converting <br>row
|
||
to column vectors</b></font></font></p>
|
||
</td>
|
||
</tr>
|
||
<tr valign="top">
|
||
<td height="275" bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: none; border-right: 1px solid #000000; padding-top: 0in; padding-bottom: 0in; padding-left: 0in; padding-right: 0.02in">
|
||
<p align="left"><font face="Arial"><font size="2" style="font-size: 10pt"><b>Reshaping
|
||
Matrices<br><br>(here: 3x3 matrix to row vector)</b></font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: none; border-right: 1px dotted #000000; padding-top: 0in; padding-bottom: 0in; padding-left: 0in; padding-right: 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
A = [1 2 3; 4 5 6; 7 8 9]<br>A =<br> 1 2
|
||
3<br> 4 5 6<br> 7 8
|
||
9<br><br><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
total_elements = numel(A)<br><br><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
B = reshape(A,1,total_elements) <br>% or reshape(A,1,9)<br>B
|
||
=<br> 1 4 7 2 5 8
|
||
3 6 9</font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: 1px dotted #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
A = np.array([[1,2,3],[4,5,6],[7,8,9]])<br><br><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
A<br>array([[1, 2, 3],<br> [4, 5, 9],<br>
|
||
[7, 8, 9]])<br><br><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
total_elements = np.prod(A.shape)</font></font></p>
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
B = A.reshape(1, total_elements) <br><br># alternative
|
||
shortcut:<br># A.reshape(1,-1)<br><br><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
B<br>array([[1, 2, 3, 4, 5, 6, 7, 8, 9]])</font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: 1px dotted #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
A = matrix(1:9,nrow=3,byrow=T)
<br><br><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
A<br>[,1] [,2] [,3]<br>[1,] 1 2 3<br>[2,] 4 5 6<br>[3,] 7 8 9
<br><br><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
total_elements = dim(A)[1] * dim(A)[2]<br><br><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
B = matrix(A, ncol=total_elements)<br><br><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
B<br>[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9]<br>[1,] 1 4 7 2
|
||
5 8 3 6 9</font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: none; padding-top: 0in; padding-bottom: 0in; padding-left: 0.02in; padding-right: 0in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
A=[1 2 3; 4 5 6; 7 8 9]<br>3x3 Array{Int64,2}:<br>1 2 3<br>4 5 6<br>7
|
||
8 9<br><br><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
total_elements=length(A)<br>9<br><br><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>B=reshape(A,1,total_elements)<br>1x9
|
||
Array{Int64,2}:<br>1 4 7 2 5 8 3 6 9</font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: 1px solid #000000; border-right: 1px solid #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Arial"><font size="2" style="font-size: 10pt"><b>Reshaping
|
||
Matrices<br><br>(here: 3x3 matrix to row vector)</b></font></font></p>
|
||
</td>
|
||
</tr>
|
||
<tr valign="top">
|
||
<td height="240" bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: none; border-right: 1px solid #000000; padding-top: 0in; padding-bottom: 0in; padding-left: 0in; padding-right: 0.02in">
|
||
<p align="left"><font face="Arial"><font size="2" style="font-size: 10pt"><b>Concatenating
|
||
matrices</b></font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: none; border-right: 1px dotted #000000; padding-top: 0in; padding-bottom: 0in; padding-left: 0in; padding-right: 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
A = [1 2 3; 4 5 6]<br><br><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
B = [7 8 9; 10 11 12]<br><br><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
C = [A; B]<br> 1 2 3<br>
|
||
4 5 6<br> 7
|
||
8 9<br> 10 11 12</font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: 1px dotted #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
A = np.array([[1, 2, 3], [4, 5, 6]])<br><br><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
B = np.array([[7, 8, 9],[10,11,12]])<br><br><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
C = np.concatenate((A, B), axis=0)<br><br><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
C<br>array([[ 1, 2, 3], <br> [ 4,
|
||
5, 6], <br> [ 7, 8, 9], <br>
|
||
[10, 11, 12]])</font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: 1px dotted #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
A = matrix(1:6,nrow=2,byrow=T)<br><br><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
B = matrix(7:12,nrow=2,byrow=T)<br><br><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
C = rbind(A,B)<br><br><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
C<br>[,1] [,2] [,3]<br>[1,] 1 2 3<br>[2,] 4 5 6<br>[3,] 7 8 9<br>[4,]
|
||
10 11 12</font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: none; padding-top: 0in; padding-bottom: 0in; padding-left: 0.02in; padding-right: 0in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
A=[1 2 3; 4 5 6];<br><br><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
B=[7 8 9; 10 11 12];<br><br><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
C=[A; B]<br>4x3 Array{Int64,2}:<br>1 2 3<br>4 5 6<br>7 8 9<br>10
|
||
11 12</font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: 1px solid #000000; border-right: 1px solid #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Arial"><font size="2" style="font-size: 10pt"><b>Concatenating
|
||
matrices</b></font></font></p>
|
||
</td>
|
||
</tr>
|
||
<tr valign="top">
|
||
<td height="256" bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: none; border-right: 1px solid #000000; padding-top: 0in; padding-bottom: 0in; padding-left: 0in; padding-right: 0.02in">
|
||
<p align="left"><font face="Arial"><font size="2" style="font-size: 10pt"><b>Stacking <br>vectors
|
||
and matrices</b></font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: none; border-right: 1px dotted #000000; padding-top: 0in; padding-bottom: 0in; padding-left: 0in; padding-right: 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
a = [1 2 3]<br><br><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
b = [4 5 6]<br><br><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
c = [a' b']<br>c =<br> 1 4<br> 2
|
||
5<br> 3 6<br><br><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
c = [a; b]<br>c =<br> 1 2 3<br>
|
||
4 5 6</font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: 1px dotted #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
a = np.array([1,2,3])<br><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
b = np.array([4,5,6])<br><br><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
np.c_[a,b]<br>array([[1, 4],<br> [2,
|
||
5],<br> [3, 6]])<br><br><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
np.r_[a,b]<br>array([[1, 2, 3],<br> [4,
|
||
5, 6]])</font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: 1px dotted #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
a = matrix(1:3, ncol=3)<br><br><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
b = matrix(4:6, ncol=3)<br><br><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
matrix(rbind(A, B), ncol=2)<br>[,1] [,2]<br>[1,] 1 5<br>[2,] 4
|
||
3
<br><br>
<font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
rbind(A,B)<br>[,1] [,2] [,3]<br>[1,] 1 2 3<br>[2,] 4 5 6</font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: none; padding-top: 0in; padding-bottom: 0in; padding-left: 0.02in; padding-right: 0in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
a=[1 2 3];<br><br><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
b=[4 5 6];<br><br><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
c=[a' b']<br>3x2 Array{Int64,2}:<br>1 4<br>2 5<br>3 6<br><br><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
c=[a; b]<br>2x3 Array{Int64,2}:<br>1 2 3<br>4 5 6</font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: 1px solid #000000; border-right: 1px solid #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Arial"><font size="2" style="font-size: 10pt"><b>Stacking <br>vectors
|
||
and matrices</b></font></font></p>
|
||
</td>
|
||
</tr>
|
||
<tr>
|
||
<td colspan="6" height="114" bgcolor="#ffffff" style="border: none; padding: 0in">
|
||
<p align="center"><font face="Arial"><font size="4" style="font-size: 14pt"><b>BASIC
|
||
MATRIX OPERATIONS</b></font></font></p>
|
||
</td>
|
||
</tr>
|
||
<tr valign="top">
|
||
<td height="253" bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: none; border-right: 1px solid #000000; padding-top: 0in; padding-bottom: 0in; padding-left: 0in; padding-right: 0.02in">
|
||
<p align="left"><font face="Arial"><font size="2" style="font-size: 10pt"><b>Matrix-scalar<br>operations</b></font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: none; border-right: 1px dotted #000000; padding-top: 0in; padding-bottom: 0in; padding-left: 0in; padding-right: 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font> A
|
||
= [1 2 3; 4 5 6; 7 8 9]<br><br><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
A * 2<br>ans =<br> 2 4 6<br>
|
||
8 10 12<br> 14 16
|
||
18<br><br><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
A + 2<br><br><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
A - 2<br><br><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
A / 2</font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: 1px dotted #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
A = np.array([ [1,2,3], [4,5,6], [7,8,9] ])<br><br><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
A * 2<br>array([[ 2, 4, 6],<br>
|
||
[ 8, 10, 12],<br> [14, 16, 18]])<br><br><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
A + 2<br><br><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
A - 2<br><br><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
A / 2</font></font></p>
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt">#
|
||
Note that NumPy was optimized for<br># in-place assignments<br>#
|
||
e.g., A += A instead of <br># A = A + A</font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: 1px dotted #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
A = matrix(1:9, nrow=3, byrow=T)<br><br><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
A * 2<br>[,1] [,2] [,3]<br>[1,] 2 4 6<br>[2,] 8 10 12<br>[3,] 14
|
||
16 18
</font></font></p>
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
A + 2<br><br><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
A - 2<br><br><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
A / 2</font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: none; padding-top: 0in; padding-bottom: 0in; padding-left: 0.02in; padding-right: 0in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
A=[1 2 3; 4 5 6; 7 8 9];<br><br># elementwise operator<br><br><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
A .* 2 <br>3x3 Array{Int64,2}:<br>2 4 6<br>8 10 12<br>14 16 18
|
||
<br><br><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
A .+ 2;<br><br><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
A .- 2;<br><br><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
A ./ 2;</font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: 1px solid #000000; border-right: 1px solid #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Arial"><font size="2" style="font-size: 10pt"><b>Matrix-scalar<br>operations</b></font></font></p>
|
||
</td>
|
||
</tr>
|
||
<tr valign="top">
|
||
<td height="165" bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: none; border-right: 1px solid #000000; padding-top: 0in; padding-bottom: 0in; padding-left: 0in; padding-right: 0.02in">
|
||
<p align="left"><font face="Arial"><font size="2" style="font-size: 10pt"><b>Matrix-matrix<br>multiplication</b></font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: none; border-right: 1px dotted #000000; padding-top: 0in; padding-bottom: 0in; padding-left: 0in; padding-right: 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font> A
|
||
= [1 2 3; 4 5 6; 7 8 9]<br><br><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
A * A<br>ans =<br> 30 36
|
||
42<br> 66 81 96<br>
|
||
102 126 150</font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: 1px dotted #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
A = np.array([ [1,2,3], [4,5,6], [7,8,9] ])<br><br><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
np.dot(A,A) # or A.dot(A)<br>array([[ 30, 36, 42],<br>
|
||
[ 66, 81, 96],<br>
|
||
[102, 126, 150]])</font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: 1px dotted #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
A = matrix(1:9, nrow=3, byrow=T)<br><br><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
A %*% A<br>[,1] [,2] [,3]<br>[1,] 30 36 42<br>[2,] 66 81 96<br>[3,]
|
||
102 126 150</font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: none; padding-top: 0in; padding-bottom: 0in; padding-left: 0.02in; padding-right: 0in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
A=[1 2 3; 4 5 6; 7 8 9];<br><br><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
A * A<br>3x3 Array{Int64,2}:<br>30 36 42<br>66 81 96<br>102 126
|
||
150</font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: 1px solid #000000; border-right: 1px solid #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Arial"><font size="2" style="font-size: 10pt"><b>Matrix-matrix<br>multiplication</b></font></font></p>
|
||
</td>
|
||
</tr>
|
||
<tr valign="top">
|
||
<td height="190" bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: none; border-right: 1px solid #000000; padding-top: 0in; padding-bottom: 0in; padding-left: 0in; padding-right: 0.02in">
|
||
<p align="left"><font face="Arial"><font size="2" style="font-size: 10pt"><b>Matrix-vector<br>multiplication</b></font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: none; border-right: 1px dotted #000000; padding-top: 0in; padding-bottom: 0in; padding-left: 0in; padding-right: 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
A = [1 2 3; 4 5 6; 7 8 9]<br><br><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
b = [ 1; 2; 3 ]<br><br><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
A * b<br>ans =<br> 14<br> 32<br>
|
||
50</font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: 1px dotted #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
A = np.array([ [1,2,3], [4,5,6], [7,8,9] ])<br><br><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
b = np.array([ [1], [2], [3] ])<br><br><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
np.dot(A,b) # or A.dot(b)<br><br>array([[14], [32], [50]])</font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: 1px dotted #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
A = matrix(1:9, ncol=3)<br><br><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
b = matrix(1:3, nrow=3)<br><br>
<font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
t(b %*% A)<br>[,1]<br>[1,] 14<br>[2,] 32<br>[3,] 50</font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: none; padding-top: 0in; padding-bottom: 0in; padding-left: 0.02in; padding-right: 0in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
A=[1 2 3; 4 5 6; 7 8 9];<br><br><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
b=[1; 2; 3];<br><br><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
A*b<br>3-element Array{Int64,1}:<br>14<br>32<br>50</font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: 1px solid #000000; border-right: 1px solid #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Arial"><font size="2" style="font-size: 10pt"><b>Matrix-vector<br>multiplication</b></font></font></p>
|
||
</td>
|
||
</tr>
|
||
<tr valign="top">
|
||
<td height="240" bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: none; border-right: 1px solid #000000; padding-top: 0in; padding-bottom: 0in; padding-left: 0in; padding-right: 0.02in">
|
||
<p align="left"><font face="Arial"><font size="2" style="font-size: 10pt"><b>Element-wise <br>matrix-matrix operations</b></font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: none; border-right: 1px dotted #000000; padding-top: 0in; padding-bottom: 0in; padding-left: 0in; padding-right: 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font> A
|
||
= [1 2 3; 4 5 6; 7 8 9]<br><br><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
A .* A<br>ans =<br> 1 4
|
||
9<br> 16 25 36<br> 49
|
||
64 81<br><br><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
A .+ A<br><br><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
A .- A<br><br><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
A ./ A</font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: 1px dotted #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
A = np.array([ [1,2,3], [4,5,6], [7,8,9] ])<br><br><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
A * A<br>array([[ 1, 4, 9],<br>
|
||
[16, 25, 36],<br> [49, 64, 81]])<br><br><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
A + A<br><br><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
A - A<br><br><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
A / A</font></font></p>
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt">#
|
||
Note that NumPy was optimized for<br># in-place assignments<br>#
|
||
e.g., A += A instead of <br># A = A + A</font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: 1px dotted #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
A = matrix(1:9, nrow=3, byrow=T)
<br><br><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
A * A<br>[,1] [,2] [,3]<br>[1,] 1 4 9<br>[2,] 16 25 36<br>[3,] 49
|
||
64 81
<br><br><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
A + A<br><br><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
A - A<br><br><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
A / A</font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: none; padding-top: 0in; padding-bottom: 0in; padding-left: 0.02in; padding-right: 0in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
A=[1 2 3; 4 5 6; 7 8 9];<br><br><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
A .* A<br>3x3 Array{Int64,2}:<br>1 4 9<br>16 25 36<br>49 64 81
|
||
<br><br><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
A .+ A; <br><br><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
A .- A; <br><br><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
A ./ A;</font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: 1px solid #000000; border-right: 1px solid #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Arial"><font size="2" style="font-size: 10pt"><b>Element-wise <br>matrix-matrix operations</b></font></font></p>
|
||
</td>
|
||
</tr>
|
||
<tr valign="top">
|
||
<td height="164" bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: none; border-right: 1px solid #000000; padding-top: 0in; padding-bottom: 0in; padding-left: 0in; padding-right: 0.02in">
|
||
<p align="left"><font face="Arial"><font size="2" style="font-size: 10pt"><b>Matrix
|
||
elements to power n <br><br>(here: individual elements
|
||
squared)</b></font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: none; border-right: 1px dotted #000000; padding-top: 0in; padding-bottom: 0in; padding-left: 0in; padding-right: 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font> A
|
||
= [1 2 3; 4 5 6; 7 8 9]<br><br><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
A.^2<br>ans =<br> 1 4 9<br>
|
||
16 25 36<br> 49 64 81</font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: 1px dotted #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
A = np.array([ [1,2,3], [4,5,6], [7,8,9] ])<br><br><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
np.power(A,2)<br>array([[ 1, 4, 9],<br>
|
||
[16, 25, 36],<br> [49, 64, 81]])</font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: 1px dotted #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
A = matrix(1:9, nrow=3, byrow=T)<br><br><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
A ^ 2<br>[,1] [,2] [,3]<br>[1,] 1 4 9<br>[2,] 16 25 36<br>[3,] 49
|
||
64 81</font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: none; padding-top: 0in; padding-bottom: 0in; padding-left: 0.02in; padding-right: 0in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
A=[1 2 3; 4 5 6; 7 8 9];<br><br><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
A .^ 2<br>3x3 Array{Int64,2}:<br>1 4 9<br>16 25 36<br>49 64 81</font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: 1px solid #000000; border-right: 1px solid #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Arial"><font size="2" style="font-size: 10pt"><b>Matrix
|
||
elements to power n <br><br>(here: individual elements
|
||
squared)</b></font></font></p>
|
||
</td>
|
||
</tr>
|
||
<tr valign="top">
|
||
<td height="225" bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: none; border-right: 1px solid #000000; padding-top: 0in; padding-bottom: 0in; padding-left: 0in; padding-right: 0.02in">
|
||
<p align="left"><font face="Arial"><font size="2" style="font-size: 10pt"><b>Matrix
|
||
to power n<br><br>(here: matrix-matrix multiplication with
|
||
itself)</b></font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: none; border-right: 1px dotted #000000; padding-top: 0in; padding-bottom: 0in; padding-left: 0in; padding-right: 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font> A
|
||
= [1 2 3; 4 5 6; 7 8 9]<br><br><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
A ^ 2<br>ans =<br> 30 36
|
||
42<br> 66 81 96<br>
|
||
102 126 150</font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: 1px dotted #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
A = np.array([ [1,2,3], [4,5,6], [7,8,9] ])<br><br><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
np.linalg.matrix_power(A,2)<br>array([[ 30, 36,
|
||
42],<br> [ 66, 81, 96],<br>
|
||
[102, 126, 150]])</font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: 1px dotted #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
A = matrix(1:9, ncol=3)<br><br>
# requires the ‘expm’
|
||
package
<br><br><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
install.packages('expm')
<br><br><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
library(expm)
<br><br>
<font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
A %^% 2<br>[,1] [,2] [,3]<br>[1,] 30 66 102<br>[2,] 36 81 126<br>[3,]
|
||
42 96 150</font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: none; padding-top: 0in; padding-bottom: 0in; padding-left: 0.02in; padding-right: 0in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
A=[1 2 3; 4 5 6; 7 8 9];<br><br><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
A ^ 2<br>3x3 Array{Int64,2}:<br>30 36 42<br>66 81 96<br>102 126
|
||
150</font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: 1px solid #000000; border-right: 1px solid #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Arial"><font size="2" style="font-size: 10pt"><b>Matrix
|
||
to power n<br><br>(here: matrix-matrix multiplication with
|
||
itself)</b></font></font></p>
|
||
</td>
|
||
</tr>
|
||
<tr valign="top">
|
||
<td height="172" bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: none; border-right: 1px solid #000000; padding-top: 0in; padding-bottom: 0in; padding-left: 0in; padding-right: 0.02in">
|
||
<p align="left"><font face="Arial"><font size="2" style="font-size: 10pt"><b>Matrix
|
||
transpose</b></font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: none; border-right: 1px dotted #000000; padding-top: 0in; padding-bottom: 0in; padding-left: 0in; padding-right: 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font> A
|
||
= [1 2 3; 4 5 6; 7 8 9]<br><br><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
A'<br>ans =<br> 1 4 7<br> 2
|
||
5 8<br> 3 6 9</font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: 1px dotted #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
A = np.array([ [1,2,3], [4,5,6], [7,8,9] ])<br><br><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
A.T<br>array([[1, 4, 7],<br> [2, 5,
|
||
8],<br> [3, 6, 9]])</font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: 1px dotted #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
A = matrix(1:9, nrow=3, byrow=T)
<br><br><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
t(A)<br>[,1] [,2] [,3]<br>[1,] 1 4 7<br>[2,] 2 5 8<br>[3,] 3 6 9</font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: none; padding-top: 0in; padding-bottom: 0in; padding-left: 0.02in; padding-right: 0in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
A=[1 2 3; 4 5 6; 7 8 9]<br>3x3 Array{Int64,2}:<br>1 2 3<br>4 5 6<br>7
|
||
8 9<br><br><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
A'<br>3x3 Array{Int64,2}:<br>1 4 7<br>2 5 8<br>3 6 9</font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: 1px solid #000000; border-right: 1px solid #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Arial"><font size="2" style="font-size: 10pt"><b>Matrix
|
||
transpose</b></font></font></p>
|
||
</td>
|
||
</tr>
|
||
<tr valign="top">
|
||
<td height="176" bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: none; border-right: 1px solid #000000; padding-top: 0in; padding-bottom: 0in; padding-left: 0in; padding-right: 0.02in">
|
||
<p align="left"><font face="Arial"><font size="2" style="font-size: 10pt"><b>Determinant
|
||
of a matrix:<br> A -> |A|</b></font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: none; border-right: 1px dotted #000000; padding-top: 0in; padding-bottom: 0in; padding-left: 0in; padding-right: 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
A = [6 1 1; 4 -2 5; 2 8 7]<br>A =<br> 6 1
|
||
1<br> 4 -2 5<br> 2 8
|
||
7<br><br><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
det(A)<br>ans = -306</font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: 1px dotted #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font> A
|
||
= np.array([[6,1,1],[4,-2,5],[2,8,7]])<br><br><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
A<br>array([[ 6, 1, 1],<br>
|
||
[ 4, -2, 5],<br> [ 2, 8,
|
||
7]])<br><br><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
np.linalg.det(A)<br>-306.0</font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: 1px dotted #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
A = matrix(c(6,1,1,4,-2,5,2,8,7), nrow=3, byrow=T)<br><br><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
A<br>[,1] [,2] [,3]<br>[1,] 6 1 1<br>[2,] 4 -2 5<br>[3,] 2 8 7<br><br><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
det(A)<br>[1] -306</font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: none; padding-top: 0in; padding-bottom: 0in; padding-left: 0.02in; padding-right: 0in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
A=[6 1 1; 4 -2 5; 2 8 7]<br>3x3 Array{Int64,2}:<br>6 1 1<br>4 -2
|
||
5<br>2 8 7<br><br><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
det(A)<br>-306.0</font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: 1px solid #000000; border-right: 1px solid #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Arial"><font size="2" style="font-size: 10pt"><b>Determinant
|
||
of a matrix:<br> A -> |A|</b></font></font></p>
|
||
</td>
|
||
</tr>
|
||
<tr valign="top">
|
||
<td height="202" bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: none; border-right: 1px solid #000000; padding-top: 0in; padding-bottom: 0in; padding-left: 0in; padding-right: 0.02in">
|
||
<p align="left"><font face="Arial"><font size="2" style="font-size: 10pt"><b>Inverse
|
||
of a matrix</b></font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: none; border-right: 1px dotted #000000; padding-top: 0in; padding-bottom: 0in; padding-left: 0in; padding-right: 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
A = [4 7; 2 6]<br>A =<br> 4 7<br> 2
|
||
6<br><br><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
A_inv = inv(A)<br>A_inv =<br> 0.60000 -0.70000<br>
|
||
-0.20000 0.40000</font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: 1px dotted #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
A = np.array([[4, 7], [2, 6]])<br><br><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
A<br>array([[4, 7], <br> [2,
|
||
6]])<br><br><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
A_inverse = np.linalg.inv(A)<br><br><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
A_inverse<br>array([[ 0.6, -0.7], <br>
|
||
[-0.2, 0.4]])</font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: 1px dotted #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
A = matrix(c(4,7,2,6), nrow=2, byrow=T)<br><br><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
A<br>[,1] [,2]<br>[1,] 4 7<br>[2,] 2 6<br><br><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
solve(A)<br>[,1] [,2]<br>[1,] 0.6 -0.7<br>[2,] -0.2 0.4</font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: none; padding-top: 0in; padding-bottom: 0in; padding-left: 0.02in; padding-right: 0in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
A=[4 7; 2 6]<br>2x2 Array{Int64,2}:<br>4 7<br>2 6<br><br><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
A_inv=inv(A)<br>2x2 Array{Float64,2}:<br>0.6 -0.7<br>-0.2 0.4</font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: 1px solid #000000; border-right: 1px solid #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Arial"><font size="2" style="font-size: 10pt"><b>Inverse
|
||
of a matrix</b></font></font></p>
|
||
</td>
|
||
</tr>
|
||
<tr>
|
||
<td colspan="6" height="202" bgcolor="#ffffff" style="border: none; padding: 0in">
|
||
<p align="center"><font face="Arial"><font size="4" style="font-size: 14pt"><b>ADVANCED
|
||
MATRIX OPERATIONS</b></font></font></p>
|
||
</td>
|
||
</tr>
|
||
<tr valign="top">
|
||
<td height="311" bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: none; border-right: 1px solid #000000; padding-top: 0in; padding-bottom: 0in; padding-left: 0in; padding-right: 0.02in">
|
||
<p align="left"><font face="Arial"><font size="2" style="font-size: 10pt"><b>Calculating
|
||
the covariance matrix <br>of 3 random variables <br><br>(here:
|
||
covariances of the means <br>of x1, x2, and x3)</b></font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: none; border-right: 1px dotted #000000; padding-top: 0in; padding-bottom: 0in; padding-left: 0in; padding-right: 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
x1 = [4.0000 4.2000 3.9000 4.3000 4.1000]’<br><br><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
x2 = [2.0000 2.1000 2.0000 2.1000 2.2000]'<br><br><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
x3 = [0.60000 0.59000 0.58000 0.62000 0.63000]’<br><br><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
cov( [x1,x2,x3] )<br>ans =<br> 2.5000e-02
|
||
7.5000e-03 1.7500e-03<br> 7.5000e-03
|
||
7.0000e-03 1.3500e-03<br> 1.7500e-03
|
||
1.3500e-03 4.3000e-04</font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: 1px dotted #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
x1 = np.array([ 4, 4.2, 3.9, 4.3, 4.1])<br><br><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
x2 = np.array([ 2, 2.1, 2, 2.1, 2.2])<br><br><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
x3 = np.array([ 0.6, 0.59, 0.58, 0.62, 0.63])<br><br><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
np.cov([x1, x2, x3])<br>Array([[ 0.025 , 0.0075 ,
|
||
0.00175],<br> [ 0.0075 , 0.007
|
||
, 0.00135],<br> [ 0.00175,
|
||
0.00135, 0.00043]])</font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: 1px dotted #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
x1 = matrix(c(4, 4.2, 3.9, 4.3, 4.1), ncol=5)<br><br><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
x2 = matrix(c(2, 2.1, 2, 2.1, 2.2), ncol=5)<br><br><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
x3 = matrix(c(0.6, 0.59, 0.58, 0.62, 0.63), ncol=5)
<br><br><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
cov(matrix(c(x1, x2, x3), ncol=3))<br>[,1] [,2] [,3]<br>[1,]
|
||
0.02500 0.00750 0.00175<br>[2,] 0.00750 0.00700 0.00135<br>[3,]
|
||
0.00175 0.00135 0.00043</font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: none; padding-top: 0in; padding-bottom: 0in; padding-left: 0.02in; padding-right: 0in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
x1=[4.0 4.2 3.9 4.3 4.1]';<br><br><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
x2=[2. 2.1 2. 2.1 2.2]';<br><br><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
x3=[0.6 .59 .58 .62 .63]';<br><br><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
cov([x1 x2 x3])<br>3x3 Array{Float64,2}:<br>0.025 0.0075
|
||
0.00175<br>0.0075 0.007 0.00135<br>0.00175 0.00135 0.00043</font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: 1px solid #000000; border-right: 1px solid #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Arial"><font size="2" style="font-size: 10pt"><b>Calculating
|
||
the covariance matrix <br>of 3 random variables <br><br>(here:
|
||
covariances of the means <br>of x1, x2, and x3)</b></font></font></p>
|
||
</td>
|
||
</tr>
|
||
<tr valign="top">
|
||
<td height="309" bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: none; border-right: 1px solid #000000; padding-top: 0in; padding-bottom: 0in; padding-left: 0in; padding-right: 0.02in">
|
||
<p align="left"><font face="Arial"><font size="2" style="font-size: 10pt"><b>Calculating <br>eigenvectors
|
||
and eigenvalues</b></font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: none; border-right: 1px dotted #000000; padding-top: 0in; padding-bottom: 0in; padding-left: 0in; padding-right: 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
A = [3 1; 1 3]<br>A =<br> 3 1<br> 1
|
||
3<br><br><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
[eig_vec,eig_val] = eig(A)<br>eig_vec =<br> -0.70711
|
||
0.70711<br> 0.70711 0.70711<br>eig_val
|
||
=<br>Diagonal Matrix<br> 2 0<br> 0
|
||
4</font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: 1px dotted #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
A = np.array([[3, 1], [1, 3]])<br><br><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
A<br>array([[3, 1],<br> [1, 3]])<br><br><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
eig_val, eig_vec = np.linalg.eig(A)<br><br><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
eig_val<br>array([ 4., 2.])<br><br><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
eig_vec<br>Array([[ 0.70710678, -0.70710678],<br>
|
||
[ 0.70710678, 0.70710678]])</font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: 1px dotted #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
A = matrix(c(3,1,1,3), ncol=2)<br><br><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
A<br>[,1] [,2]<br>[1,] 3 1<br>[2,] 1 3<br><br><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
eigen(A)<br>$values<br>[1] 4 2<br><br>$vectors<br>[,1] [,2]<br>[1,]
|
||
0.7071068 -0.7071068<br>[2,] 0.7071068 0.7071068</font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: none; padding-top: 0in; padding-bottom: 0in; padding-left: 0.02in; padding-right: 0in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
A=[3 1; 1 3]<br>2x2 Array{Int64,2}:<br>3 1<br>1 3<br><br><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
(eig_vec,eig_val)=eig(a)<br>([2.0,4.0],<br>2x2
|
||
Array{Float64,2}:<br>-0.707107 0.707107<br>0.707107 0.707107)</font></font></p>
|
||
</td>
|
||
<td bgcolor="#cfe7f5" style="border-top: none; border-bottom: none; border-left: 1px solid #000000; border-right: 1px solid #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Arial"><font size="2" style="font-size: 10pt"><b>Calculating <br>eigenvectors
|
||
and eigenvalues</b></font></font></p>
|
||
</td>
|
||
</tr>
|
||
<tr valign="top">
|
||
<td height="541" bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: none; border-right: 1px solid #000000; padding-top: 0in; padding-bottom: 0in; padding-left: 0in; padding-right: 0.02in">
|
||
<p align="left"><font face="Arial"><font size="2" style="font-size: 10pt"><b>Generating
|
||
a Gaussian dataset:<br><br>creating random vectors from the
|
||
multivariate normal<br>distribution given mean and covariance
|
||
matrix<br><br>(here: 5 random vectors with <br>mean 0, covariance
|
||
= 0, variance = 2)</b></font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: none; border-right: 1px dotted #000000; padding-top: 0in; padding-bottom: 0in; padding-left: 0in; padding-right: 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt">%
|
||
requires statistics toolbox package<br>% how to install and load
|
||
it in Octave:<br><br>% download the package from: <br>%
|
||
http://octave.sourceforge.net/packages.php<br>% pkg install <br>%
|
||
~/Desktop/io-2.0.2.tar.gz <br>% pkg install <br>%
|
||
~/Desktop/statistics-1.2.3.tar.gz<br><br><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
pkg load statistics<br><br><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
mean = [0 0]<br><br><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
cov = [2 0; 0 2]<br>cov =<br> 2 0<br> 0
|
||
2<br><br><font color="#006699"><font size="2" style="font-size: 10pt"><b>M></b></font></font>
|
||
mvnrnd(mean,cov,5)<br> 2.480150 -0.559906<br>
|
||
-2.933047 0.560212<br> 0.098206
|
||
3.055316<br> -0.985215 -0.990936<br> 1.122528
|
||
0.686977<br> </font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: 1px dotted #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt"><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
mean = np.array([0,0])<br><br><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
cov = np.array([[2,0],[0,2]])<br><br><font color="#ff420e"><font size="2" style="font-size: 10pt"><b>P></b></font></font>
|
||
np.random.multivariate_normal(mean, cov, 5)<br><br>Array([[
|
||
1.55432624, -1.17972629], <br>
|
||
[-2.01185294, 1.96081908], <br>
|
||
[-2.11810813, 1.45784216], <br>
|
||
[-2.93207591, -0.07369322], <br>
|
||
[-1.37031244, -1.18408792]])</font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: 1px dotted #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt">#
|
||
requires the ‘mass’ package<br><br>
<font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
install.packages('MASS')<br><br>
<font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
library(MASS)
<br><br><font color="#0000ff"><font size="2" style="font-size: 10pt"><b>R></b></font></font>
|
||
mvrnorm(n=10, mean, cov)<br>[,1] [,2]<br>[1,] -0.8407830
|
||
-0.1882706<br>[2,] 0.8496822 -0.7889329<br>[3,] -0.1564171
|
||
0.8422177<br>[4,] -0.6288779 1.0618688<br>[5,] -0.5103879
|
||
0.1303697<br>[6,] 0.8413189 -0.1623758<br>[7,] -1.0495466
|
||
-0.4161082<br>[8,] -1.3236339 0.7755572<br>[9,] 0.2771013
|
||
1.4900494<br>[10,] -1.3536268 0.2338913</font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: 1px dotted #000000; border-right: none; padding-top: 0in; padding-bottom: 0in; padding-left: 0.02in; padding-right: 0in">
|
||
<p align="left"><font face="Andale Mono"><font size="1" style="font-size: 8pt">#
|
||
requires the Distributions package from
|
||
https://github.com/JuliaStats/Distributions.jl<br><br><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
using Distributions<br><br><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
mean=[0., 0.]<br>2-element Array{Float64,1}:<br>0.0<br>0.0<br><br><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
cov=[2. 0.; 0. 2.]<br>2x2 Array{Float64,2}:<br>2.0 0.0<br>0.0
|
||
2.0<br><br><font color="#9966cc"><font size="2" style="font-size: 10pt"><b>J></b></font></font>
|
||
rand( MvNormal(mean, cov), 5)<br>2x5 Array{Float64,2}:<br>-0.527634
|
||
0.370725 -0.761928 -3.91747 1.47516<br>-0.448821 2.21904 2.24561
|
||
0.692063 0.390495</font></font></p>
|
||
</td>
|
||
<td bgcolor="#ffffff" style="border-top: none; border-bottom: none; border-left: 1px solid #000000; border-right: 1px solid #000000; padding: 0in 0.02in">
|
||
<p align="left"><font face="Arial"><font size="2" style="font-size: 10pt"><b>Generating
|
||
a Gaussian dataset:<br><br>creating random vectors from the
|
||
multivariate normal<br>distribution given mean and covariance
|
||
matrix<br><br>(here: 5 random vectors with <br>mean 0, covariance
|
||
= 0, variance = 2)</b></font></font></p>
|
||
</td>
|
||
</tr>
|
||
</table>
|
||
|
||
|
||
<br>
|
||
|
||
<br><br><font size=0.6>
|
||
(Thanks to Keith C. Campbell for providing me with the syntax for the Julia language.)
|
||
</font>
|
||
|
||
<br>
|
||
|
||
<a id='numpy_arrays'></a>
|
||
|
||
|
||
<br>
|
||
<br>
|
||
<br>
|
||
<br>
|
||
|
||
### Alternative data structures: NumPy matrices vs. NumPy arrays
|
||
|
||
[[back to section overview](#sections)]
|
||
|
||
Python's NumPy library also has a dedicated "matrix" type with a syntax that is a little bit closer to the MATLAB matrix: For example, the "<code> * </code>" operator would perform a matrix-matrix multiplication of NumPy matrices - same operator performs element-wise multiplication on NumPy arrays.
|
||
Vice versa, the "`.dot()`" method is used for element-wise multiplication of NumPy matrices, wheras the equivalent operation would for NumPy arrays would be achieved via the "<code> * </code>"-operator.
|
||
**Most people recommend the usage of the NumPy array type over NumPy matrices, since arrays are what most of the NumPy functions return.** |