nbhkdz.com冰点文库

mathlab作业二


MATLAB 作业二
1、 请将下面给出的矩阵 A 和 B 输入到 MATLAB 环境中,并将它们转换成符号矩阵。若某 一矩阵为数值矩阵,另以矩阵为符号矩阵,两矩阵相乘是符号矩阵还是数值矩阵。

?5 ?2 ? ?6 ? A??3 ?10 ? ?7 ?4 ?

7 6 5 1 6 5? 5 5 0 ?3 ?3 3 1 0 0 1

4? 2 5 4 ? ? ?1 4 2 0 6 4 4? 2 1 1 ? ? 9 6 3 6 6 2? , B ? ? 3 5 1 5 ?4 7 6 0 0 7 7? 1 0 1 ? ? 2 4 4 0 7 7? ? ?3 ?4 ?7 3 ? 1 ?10 7 ?6 8 6 7 2 1 7? ? ?

1 2 6 3 2 2 2 4 1 0

7 8 8 1

3? 5? ? 6? ? 2? 1? ? 12 ? 5? ?

一矩阵为数值矩阵,另以矩阵为符号矩阵,两矩阵相乘是符号矩阵。 2、 利用 MATLAB 语言提供的现成函数对习题 1 中给出的两个矩阵进行分析,判定它们是

否为奇异矩阵,得出矩阵的秩、行列式、迹和逆矩阵,检验得出的逆矩阵是否正确。 对 A: >> det(A) ans = -3.5432e+04 >> rank(A) ans =7 >> trace(A) ans = 27

对 B: >> det(B) >> rank(B)

ans = -2.6326e-26 ans =5

>> trace(B) ans = 26 由 A 和 B 的行列式的值可以看出,矩阵 A 为非奇异矩阵,而矩阵 B 为奇异矩阵,而由 norm 的结果可以说明对 A 直接求逆的结果比较准确, 而对矩阵 B 直接求逆的偏差较大, 所以对奇异矩阵不适合用 inv 函数直接求逆。

3、 试求出习题 1 中给出的 A 和 B 矩阵的特征多项式、 特征值与特征向量, 并对它们进行 LU 分解。 对 A:

对 B:

4、试求下面齐次方程的基础解系。

6 x1 ? x2 ? 4 x3 ? 7 x4 ? 3x5 ? 0 ? ? ?2 x1 ? 7 x2 ? 8 x3 ? 6 x4 ? 0 ? ? ? ?4 x1 ? 5 x2 ? x3 ? 6 x4 ? 8 x5 ? 0 ? ?34 x ? 36 x ? 9 x ? 21x ? 49 x ? 0 1 2 3 4 5 ? ? 26 x ? 12 x ? 27 x ? 27 x ? 17 x ? 1 2 3 4 5 ?0 ?

基础解系 x:

5、试求下面线性代数方程的解析解与数值解,并检验解的正确性。

?2 ?10 ? ?8 ? ? ?5

?9 3 ?2 ?1? ? ?1 ?4 0 ? ? ? ?3 ?8 ?4 ? ?1 10 5 0 ? ? X ?? ?0 3 3? ?2 ?4 ?6 3 ? ? ? ? ?6 ?6 ?8 ?4 ? ? 9 ?5 3 ?

解析解求解过程:

通解:

检验:

数值解求解过程:

通解:

检验:

由以上过程可以对原方程组求解,从 norm 的结果可以看出用解析的方法求得的结果的精 确度更高。

6、试判定下面的线性代数方程是否有解。

?16 2 3 13? ?1 ? ? 5 11 10 8 ? ? ? ? ? X ? ? 3? ? 9 7 6 12? ?4? ? ? ? ? ? 4 14 15 1 ? ?7 ?

通过求方程组的系数矩阵和增广矩阵的秩,因为两个秩不等故原方程组无解。 7、求解能转换成多项式方程的联立方程,并检验得出的高精度数值解(准解析解)的精度。

? x ? x2 ? 1 ? 0 1) ? 2 2 ?( x1 ? 2) ? ( x2 ? 0.5) ? 1 ? 0
2 1

? x 2 y 2 ? zxy ? 4 x 2 yz 2 ? xz 2 ? 2) ? xy 3 ? 2 yz 2 ? 3x3 z 2 ? 4 xzy 2 ? y 2 x ? 7 xy 2 ? 3xz 2 ? x 4 zy ?

(1)

(2)

? 10x1 ? x2 ? 2 x3 ? 72 ? 8 、 用 Jacobi 、 Gauss-Seidel 迭 代 法 求 解 方 程 组 ? ? x1 ? 10 x2 ? 2 x3 ? 83 , 给 定 初 值 为 ? ? x ? x ? 5 x ? 42 3 ? 1 2

x(0) ? (0,0,0)T 。
jocobi 迭代法:

Gauss-Seidel 迭代法:

两种方法的计算结果一致,但是可以看出后一种方法的迭代次数更少。 9、 取? ? 1.4, x
(0)

?1 ? 2 x1 ? x2 ? ? (1,1,1) , 用超松弛法解方程组 ? ? x1 ? 2 x2 ? x3 ? 0 ? ? x ? 2 x ? 1.8 2 3 ?
T


《应用计算方法教程》matlab作业二

《应用计算方法教程》matlab作业二_理学_高等教育_教育专区。作业六 6-1 试验目的计算特征值,实现算法 试验内容:随机产生一个 10 阶整数矩阵,各数均在-5 和 5...

Matlab第二次作业

《数学实验》报告 实验名称 学院 Matlab 第二次作业 冶金与生态工程学院 专业班级 姓学名号 2013 年 11 月 一、 【实验目的】 编辑比较复杂的程序,学会用 ...

MATLAB作业2 参考答案

MATLAB作业2 参考答案_院校资料_高等教育_教育专区。中科院 matlab作业MATLAB 作业二参考答案 1、 试求出如下极限。 (1) lim x 2 y ? xy 3 1 ? cos( x...

MATLAB第二次作业

《数学实验》报告 实验名称 MATLAB 绘图第一次实验 学院 高等工程师学院 专业班级 姓学名号 4 1 2 0 4 0 4 4 2014 年 5 月 数学实验 MATLAB 孙金胜 ...

mathlab作业二

MATLAB 作业二 1、 请将下面给出的矩阵 A 和 B 输入到 MATLAB 环境中,并将它们转换成符号矩阵。若某 一矩阵为数值矩阵,另以矩阵为符号矩阵,两矩阵相乘是符号...

matlab作业二

8页 2下载券 MATLAB作业2参考答案 5页 免费 matlab第二、三章作业 暂无评价 3页 2下载券m​a​t​l​a​b​作​业​二 ...

matlab1-8章课后作业

MATLAB 基础教程 1~8 章作业 Matlab 第一章 1.阐述 Matlab 的功能 Matlab ...>> vecR=vec' vecR = 4 5 2 8 4 7 2 64 2 57 2 45 7 43 2 5 ...

mathlab作业一

MATLAB 作业一 1、 试求出如下极限。 (1) lim ( x ? 2) x ? 2 ( x ? 3) x ?3 x 2 y ? xy 3 1 ? cos( x 2 ? y 2 ) lim lim ,(...

matlab第二次作业答案

6页 1下载券 matlab上机作业答案参考 25页 2下载券 Matlab第六章作业答案 暂无评价 2页 1下载券 matlab程序设计(第二版)... 9页 1下载券喜欢...

matlab第二次作业

matlab第二次作业_理学_高等教育_教育专区。《数学实验》报告 实验名称 学院 MATLAB 绘制曲线 机械工程学院 物流 1401 孙宇 41445011 专业班级 姓学名号 2015 年 ...