第八章 欧氏空间
- 格式:doc
- 大小:58.00 KB
- 文档页数:3
MATLAB 软件应用
第八章 欧氏空间
例1:判断向量()2,1,4α=-和()4,4,1β=--是否正交.
解:建立m 文件i1.m 如下
clc
a=[2 -1 4];
b=[-4 -4 1];
c=dot(a,b) %求向量a ,b 的内积 运行结果如下:
c =
a 与
b 内积为0,则a 和b 正交.
例2:将向量组()1,0,1α=,()1,0,1β=-规范正交化
解:建立m 文件i2.m 如下
clc
A=[1 1;0 0;1 -1];
B=sym(orth(A)) %将A 的列向量组正交规范化,
%并以符号的形式输出
B'*B %可以验证两个向量是否规范正交化 运行结果如下
B =
[ -sqrt(1/2), -sqrt(1/2)]
[ 0, 0]
[ -sqrt(1/2), sqrt(1/2)]
ans =
[ 1, 0]
[ 0, 1]
例3 将矩阵400031013A ⎛⎫ ⎪= ⎪ ⎪⎝⎭
规范正交化。
解:建立m 文件i3.m 如下
clc
A=[4 0 0; 0 3 1; 0 1 3];
B=orth(A)
Q=B'*B
运行结果如下:
B =
0 1 0 -985/1393 0 -985/1393 -985/1393 0 985/1393 Q =
1 0 * 0 1 0 * 0 1
例4:对实对称矩阵222254245A -⎛⎫ ⎪=- ⎪ ⎪--⎝⎭
,求正交矩阵U , 使得T U AU 为对角矩阵
解:建立m 文件i4.m 如下
clc
A=[2 2 -2;2 5 -4;-2 -4 5]; %实对称矩阵A
[V,D]=eig(A) %矩阵A 的对角化 V'*V %验证V 是正交阵 运行结果如下:
V =
-963/3230 2584/2889 1/3 -963/1615 -1292/2889 2/3 -963/1292 0 -2/3
D =
1 0 0 0 1 0 0 0 10 ans =
1 * * * 1 0 * 0 1
例5:求向量()2,1,4α=-的模长
解:建立m 文件i5.m 如下
clc
a=[2 -1 4]
b=norm(a)
运行结果如下:
a =
2 -1 4
b =
3524/769
【练习与思考】
1、判别下列矩阵是否正交矩阵
(1)
11
1
23
11
1
22
11
1
32
A
⎛⎫
-
⎪
⎪
⎪
=-
⎪
⎪
⎪
-
⎪
⎝⎭
(2)
184
999
814
999
447
999
A
⎛⎫
--
⎪
⎪
⎪
=--
⎪
⎪
⎪
--
⎪
⎝⎭
2、将下列向量组规范正交化
(1)、()
123111 124 139
ααα
⎛⎫ ⎪= ⎪ ⎪⎝⎭
(2)、()
123111 011 101 110
ααα
-⎛⎫ ⎪
-
⎪=
⎪
-
⎪⎝⎭。