Matlab 与线性代数教案_线性代数及matlab教案

Matlab 与线性代数教案由刀豆文库小编整理，希望给你工作、学习、生活带来方便，猜你可能喜欢“线性代数及matlab教案”。

Matlab 与线性代数

一、Matlab 入门：

1.启动、退出、运行： 2.窗口介绍： 3.基本符号： =：赋值符号

[ ]：数组定义符号 , 区分列 函数参数分隔符;区分行 取消运行显示 % 注释标记

: 具有多种应用功能

4.matlab的变量(区分大小写): 预定义变量: ans

pi 相关命令: format(显示格式 rat long short)

who whos clear

5.M 文件（纯文本文件，扩展名为.m）建立 修改 保存 运行

二、Matlab 与线性代数的基本运算

1.矩阵的输入

数字矩阵：A=[1 2 3;3 2 1]

或 A=[1, 2, 3;3, 2, 1] 或 A=[1 2 32 1]

符号矩阵(显示出来元素之间有逗号): 定义符号变量 sym syms

用法：(1).sym(‘[a,b,c;b,c,a]’)或 sym(‘[a b c;b c a]’)

(2).syms a b c

A=[a b c;b c a]

2.产生特殊矩阵的函数：

zeros(m,n)zeros(n)

ones(m,n)ones(n)eye(n)

magic(n)rand(m,n)randn(n)% 产生（0,1）区间均匀分布的随机矩阵

3.相关命令：

round(A)% 表示对矩阵A中所有元素进行四舍五入 length(A)% 返回A的长度（列数）size(A)% 返回Ａ的尺寸，行数 列数 A(i,j)% 引用矩阵A的第i行第j列元素

4.矩阵的基本运算

(1).+-*.*

(2).转置 A’

(3).方阵的幂：A^3

5.求向量组的极大无关组

A[1,2,3 ]

(1).U=rref(A)% U为A的行最简形

(2).[U,s]=rref(A)% U为A的行最简形, s为首非零元所在列组成的向量

(3).rrefmovie(A)% 返回A的行最简形，且给出每一步化简过程

6.求线性方程组的解

情形1。Ax=b，其中A为n阶可逆阵

法1： x=inv(A)*b 或 x=A^(-1)*b

法2： U=rref([A,b])% 返回值U为矩阵的行最简形，最后一列即为解x。

情形2。Ax=0, 其中A 为m*n 矩阵，R(A)=r

法1：U=rref(A), 选定自由变量，得到一组基础解系

法2：z=null(A)

% z的列向量为Ax=0的一组标准正交基。

情形3。Ax=b, 其中A 为m*n 矩阵, 求通解

U=rref([A,b])从最后一列找特解，前n列找导出组的基础解系，然后按格式写

出Ax=b的通解。（或先写出以U为增广矩阵的同解方程组也可。）

6x13x22x33x44x554x12x2x32x43x54

例子: .4x2x3x2xx0234512xx7x3x2x112345(4).(5).(6).(7).方阵行列式 det(A)方阵的秩 rank(A)方阵的逆 inv(A)或 A^(-1)矩阵的除法 左除 右除/

AB=C

则 A=C/B B=AC 输入：A=[6 3 2 3 4;4 2 1 2 3;4 2 3 2 1;2 1 7 3 2];

b=[5 4 0 1]’;

U=rref([A,b])10得到：U001/2000010000103/417/203/22 603/20取x2,x5为自由变量，令x20,x50得Ax=b的特解*2

60

1/23/410x210分别令和得导出组的基础解系为：10，21

x50107/2013x24x5x112x3x5或：导出组Ax=0的同解方程组：，x2,x5为自由变量，分别令x47x521/23/410x21,x50和x20,x51得导出组的基础解系为：10，21。

07/2017.求矩阵的特征值与特征向量

(1).d=eig(A)% d为矩阵A的特征值构成的向量

(2).[V,D]=eig(A)% D为A 的特征值构成的对角阵，V 的列为A的单位特征向

量，与D中的特征值对应，满足：AVDV8.Schmidt 正交化方法

B=orth(A)% B的列向量为A的列空间的一组标准正交基，换句话说，B的列是

A的列向量的正交标准化，满足B*Beye(rank(A))。

9.用正交变换化二次型为标准形

先写出所给二次型的矩阵A，则A为实对称矩阵，[V,D]=eig(A)% D 为A的特征值构成的对角阵，V的列向量为A的正交单位特征

向量，次序与D的元素对应。满足VAVDVT1'1，即AVVD。

AV。