Matlab在高等数学和线性代数中的应用

Matlab在高等数学和线性代数中的应用

开放实验指导书

曲庆国 编

山东交通学院理学院

2014年9月

⒈目的

本实验旨在向学生介绍一种解决专业问题的快速有效且具有强大功能的科学与工程计算软件。通过本实验，应使学生掌握的内容是：利用Matlab求极限，求导数，求极值，求积分，级数求和；利用Matlab求向量组的最大无关组，求齐次和非齐次线性方程组，求方阵的特征值和特征向量，求正交变换把二次型化成标准型。该实验主要为上机实验，要求学生按要求上机实现相关的程序的设计，自己动手编写程序并验证程序的正确性。

⒉实验任务分解

通过一些实例初步掌握Matlab在高等数学和线性代数中的应用。实验任务可分解为：Matlab在高等数学的应用，Matlab在线性代数中的应用(一)，Matlab在线性代数中的应用(二)。

⒊实验环境介绍

长清校区数学实验室

⒋实验时数

20学时

实验一 Matlab在高等数学的应用（8学时）

实验目的：

1. 2. 3. 4. 5.

利用Matlab求极限； 利用Matlab求导数； 利用Matlab求极值； 利用Matlab求积分； 利用Matlab对级数求和。

实验要求:给出程序和实验结果。 实验相关理论内容:

1 利用Matlab求极限

Matlab求极限的命令为： limit(F,x,a) limit(F,a) limit(F)

limit(F,x,a,'right') limit(F,x,a,'left')

2．利用Matlab求导数

Matlab求导数的命令为： Y = diff(X) Y = diff(X,n) Y = diff(X,n,dim)

3．利用Matlab求极值

先用diff求出导函数，再用solve求出导函数的根。 4．利用Matlab求积分

Matlab求积分的命令为： R = int(S) R = int(S,v) R = int(S,a,b) R = int(S,v,a,b)

5．利用Matlab对级数求和

Matlab级数求和的命令为： r = symsum(s) r = symsum(s,v) r = symsum(s,a,b) r = symsum(s,v,a,b)

实验题目：

1?x2?1?a?一、（1）lim; （2）lim?1??

x?01?cosxx????x?x二、求函数y?lnx?2的三阶导数。 1?x三、求函数f(x)?x3?6x2?8x?1的极值点，并画出函数的图形。 四、（1）求不定积分

?1?211?x2dx;

?（2）求定积分

??cosxcos2xdx.

2?五、求如下级数的和： （1）

2n?1; ?n2n?11. ?2n?1n??（2）

实验二 Matlab在线性代数中的应用(一)（6学时） 实验目的：

1. 利用Matlab求出向量组的一个最大无关组； 2. 利用Matlab求齐次线性方程组； 3. 利用Matlab求非齐次线性方程组。

实验要求:给出程序和实验结果。 实验相关理论内容:

1．向量组的线性相关性

求列向量租A的一个最大无关组，可用命令rref(A)将A化为行最简形，其中单位向量对应的列向量即为最大无关组所含向量。 2. 齐次线性方程组

在Matlab中，函数null用来求解零空间，即满足Ax?0的解空间，实际上是求出解空间的一组基（基础解系）。 3. 非齐次线性方程组

Matlab中解非齐次线性方程组可用使用“\\”。虽然表面上只是一个简单的符号，而它的内部却包含许多自适应算法，如对超定方程组（无解）用最小二乘法，对欠定方程组（多解）它将给出范数最小的一个解。

另外求解欠定方程组（多解）可以使用求矩阵A的行最简形命令rref(A)，求出所有的基础解系。

对超定方程组，可以用pinv求解。Pinv使用范围比“\\”更广泛，其可给出最小二乘解或最小范数解。

实验题目：

一、求下列矩阵列向量组的一个最大无关组。

?1?2?1??242A???2?10??333二、设A??a102?6?6?? 23??34?a2?22?1??14??，B?bb??03?，

a3???2?12?12?????验证a1,a2,a3是

????122????42??R3的一个基，并把b1,b2用这个基线性表示。

三、求方程组的通解

?x1?2x2?2x3?x4?0,??2x1?x2?2x3?2x4?0, ?x?x?4x?3x?0.34?12四、求超定方程组

?2x1?4x2?11,?3x?5x?3,?12 ?x?2x?6,2?1??2x1?x2?7,.五、用最小二乘法解方程组

x1?x2?1,??x?x?2,?13 ?x?x?x?0,?123??x1?2x2?x3??1.六、求解方程组

?x1?x2?x3?x4?0,??x1?x2?x3?3x4?1, ?x?x?2x?x??1/2.34?12

实验三 Matlab在线性代数中的应用(二)（6学时）

实验目的：

1. 利用Matlab求方阵的特征值及相应的特征向量；

2. 利用Matlab求一个正交变换把二次型化成标准型。

实验要求:给出程序和实验结果。

实验相关理论内容:

1、Matlab求方阵的特征值及相应的特征向量

Matlab求特征值及特征向量的命令为： d = eig(A) d = eig(A,B) [V,D] = eig(A)

[V,D] = eig(A,'nobalance') [V,D] = eig(A,B) [V,D] = eig(A,B,flag)

2、Matlab求一个正交变换把二次型化成标准型

应用功能强大的eig命令编程可以实现这一功能。

实验题目：

一、求下列矩阵的特征值及相应特征向量

??35 1 6 26 19 24? 3 32 7 21 23 25?A??22?1????? 1 1 1 1?2?12??31 9 2 22 27 20??122?;B???;?C=? 1 2 3 4? 8 28 33 17 10 15???????? 1 3 6 10?30 5 34 12 14 16?? 1 4 10 ?4 36 29 13 18 11??二、求一个正交变换x?py，把二次型

f?2x1x2?2x1x3?2x1x4?2x2x3?2x2x4?2x3x4

化成标准型。 三、判别二次型

f?2x221?4x22?5x3?4x1x2

的正定性，并求正交变换把二次型化成标准型。

????. 20??

本文来源：https://www.bwwdw.com/article/tm32.html