? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 數組運算
http://wenku.baidu.com/link?url=4WxXul51iVOqyNlyYapmMZmebo9fg-4J-GK_DfdSjw7JW1ELv17FvnmW9n5f7jib6n11AeJ9tQjYuoAZzYzeow2xGe0SZscY2Ua39FmHSv3
例程：>>A=magic(3)? A =
8?????1?????6 ? ? ?
3?????5?????7??????
4?????9?????2
>>A(2:3,3:-1:1)
ans?=
7?????5?????3??????
2?????9?????4
>>A(:,end)
ans?=??????6??????7??????2
>>A(1,end-1)
ans?=??????1
>>A([2?1?3?3],[1?1?2?2?1])
ans?=
3?????3?????5?????5?????3?????
?8?????8?????1?????1?????8??
?4?????4?????9?????9?????4?
?4?????4?????9?????9?????4
A(ii,jj)：其中ii和jj可以是一維向量、標量、“:”號或者“end”
大家對下標估計比較熟悉，由于在C語言中接觸過，但是我這里需要強調的是，Matlab的下標是可以多行多列同時引用的，而像C語言等一次只能引用一個，比如
A(2:3,3:-1:1)表示引用數組中的2~3行，3~1列對應的元素
A(:,end)表示引用最后一列元素，“:”表示所有列或行，“end”表示最后一列或列，“end-n”表示倒數第n行或列
A(1,end-1)表示引用第1行倒數第2個元素
A([2 1 3 3],[1 1 2 2 1])表示引用按兩個向量引用指定的元素，即A中的第2,1,3,3行和第1,1,2,2,1列對應的元
1、向量的創建
1）直接輸入：
行向量：a=[1,2,3,4,5]
列向量：a=[1;2;3;4;5]
2）用“:”生成向量
a=J:K生成的行向量是a=[J,J+1,…,K]
a=J:D:K生成行向量a=[J,J+D,…,J+m*D],m=fix((K-J)/D)
3）函數linspace用來生成數據按等差形式排列的行向量
x=linspace(X1,X2):在X1和X2間生成100個線性分布的數據，相鄰的兩個數據的差保持不變。構成等差數列。
x=linspace(X1,X2,n):在X1和X2間生成n個線性分布的數據，相鄰的兩個數據的差保持不變。構成等差數列。
4）函數logspace用來生成等比形式排列的行向量
X=logspace(x1,x2)在x1和x2之間生成50個對數等分數據的行向量。構成等比數列，數列的第一項x(1)=10x1,x(50)=10x2
X=logspace(x1,x2,n)在x1和x2之間生成n個對數等分數據的行向量。構成等比數列，數列的第一項x(1)=10x1,x(n)=10x2
注：向量的的轉置：x=(0,5)’
2）函數eye，生成單位矩陣
eye(n)?:生成n*n階單位E
eye(m,n):生成m*n的矩陣E，對角線元素為1，其他為0
eye(size(A))：生成一個矩陣Ａ大小相同的單位矩陣
eye(m,n,classname):對角線上生成的元素是1，數據類型用classname指定。其數據類型可以是：duoble、single、int8、uint8、int16、uint16、int32、uint32。
3）函數ones用ones生成全1的矩陣
ones(n)?:生成n*n的全1矩陣
ones(m,n)?:生成m*n的全1矩陣
ones(size(A))?:生成與矩陣A大小相同的全1矩陣
ones(m,n,p,…)生成m*n*p*….的全1的多維矩陣
ones(m,n,…,classname)制定數據類型為classname
4）函數zeros函數zeros生成全0矩陣
zeros(n):生成n*n的全0矩陣
zeros(m,n:)生成m*n的全0矩陣
zeros(size(A)):生成與矩陣A大小相同的全0矩陣
zeros?(m,n,p,…)生成m*n*p*….的全0的多維矩陣
zeros?(m,n,…,classname)指定數據類型為classname
5）函數rand函數rand用來生成[0,1]之間均勻分布的隨機函數，其調用格式是：
Y=rand:生成一個隨機數
Y=rand(n):生成n*n的隨機矩陣
Y=rand(m,n):生成m*n的隨機矩陣
Y=rand(size(A)):生成與矩陣A大小相同的隨機矩陣
Y=rand(m,n,p,…):生成m*n*p*…的隨機數多維數組
6）函數randn函數rand用來生成服從正態分布的隨機函數，其調用格式是：
Y=randn:生成一個服從標準正態分布的隨機數
Y=randn(n):生成n*n的服從標準正態分布的隨機矩陣
Y=randn(m,n):生成m*n的服從標準正態分布的隨機矩陣
Y=randn(size(A)):生成與矩陣A大小相同的服從標準正態分布的隨機矩陣
Y=randn(m,n,p,…):生成m*n*p*…的服從標準正態分布的隨機數多維數組
3、矩陣元素的提取與替換
1）單個元素的提取
如：a=[1,2,3;3,4,5],運行后：
a?=
1?????2?????3
3?????4?????5
輸入b=a(1,2)
b?=
2
2）提取矩陣中某一行的元素，
如：a=[1,2,3;3,4,5],運行后：
a?=
1?????2?????3
3?????4?????5
輸入b=a(1,:)
b?=
1?????2?????3
3）提取矩陣中某一列：
如：a=[1,2,3;3,4,5],運行后：
a?=
1?????2?????3
3?????4?????5
輸入b=a(:,1)
b?=
1
3
4）提取矩陣中的多行元素
如：a=[1,2,3;3,4,5],運行后：
a?=
1?????2?????3
3?????4?????5
輸入b=a([1,2],:)
b?=
1?????2?????3
3?????4?????5
5）提取矩陣中的多列元素
如：a=[1,2,3;3,4,5],運行后：
a?=
1?????2?????3
3?????4?????5
輸入b=a(:,[1,3])
b?=
1?????3
3?????5
6）提取矩陣中多行多列交叉點上的元素
如：a=[1,2,3;3,4,5],運行后：
a?=
1?????2?????3
3?????4?????5
輸入b=a([1,2],[1,3])
b?=
1?????3
3?????5
7）單個元素的替換：
如：a=[1,2,3;3,4,5],運行后：
a?=
1?????2?????3
3?????4?????5
輸入：a(2,3)=-1
a?=
1?????2?????3
3?????4????-1
4、矩陣元素的重排和復制排列
1）矩陣元素的重排
B=reshape(A,m,n):返回的是一個m*n矩陣B，矩陣B的元素就是矩陣A的元素，若矩陣A的元素不是m*n個則提示錯誤。
B=reshape(A,m,n,p):返回的是一個多維的數組B，數組B中的元素個數和矩陣A中的元素個數相等
B=reshape(A,…,[],…):可以默認其中的一個維數
B=reshape(A,siz)?:由向量siz指定數組B的維數，要求siz的各元素之積等于矩陣A的元素個數
2）矩陣的復制排列函數是repmat
B=repmat(A,n):返回B是一個n*n塊大小的矩陣，每一塊矩陣都是A
B=repmat(A,m,n):返回值是由m*n個塊組成的大矩陣，每一個塊都是矩陣A。
B=repmat(A,[m,n,p,…]):返回值B是一個多維數組形式的塊，每一個塊都是矩陣A
