电子书_MATLAB程序设计与应用(第二版)

阅阅读读时时：：请请选选择择““视视图图 || 文文档档结结构构图图””，，弹弹出出文文档档中中的的标标题题链链接接。。

数学软件

MATLAB 程序设计与应用

第1章MATLAB系统环境

1.1 MATLAB概貌

1.2 MATLAB环境的准备

1.3 MATLAB操作界面

1.4 MATLAB帮助系统

自上世纪80年代以来，出现了科学计算语言，亦称数学软件。

MATLAB

Mathematica

Mathcad

Maple

LINDO

LINGO

1.1 MATLAB概貌p3

MATLAB 是MATrix LABoratory（矩阵实验室）的缩写。

1984年由Math Works 公司推出，现已成为国际公认的优秀的工程应用开发环境，是影响最大，流行最广的科学计算语言。

1.1.1 MATLAB的发展

2. 例如MATLAB7.6的建造编号是R2008a。说明MATLAB7.6与MATLAB2008a是等同的；

3. 对于建造编号，正规化以后，每年出两个版本。一般来说。a是测试版，b是正式版。a是前半

年出，b是后半年出。

教材采用MATLAB7.0（R14，2004）

实验室采用MATLAB 7.8（R2009a，2009.3，汉化）

1.1.2 MATLAB的主要功能p4

◆数值计算和符号计算功能

◆绘图功能

◆语言体系

◆MATLAB工具箱

(1) 数值计算和符号计算功能

MATLAB以矩阵作为数据操作的基本单位，还提供了十分丰富的数值计算函数。

MATLAB先后和著名的符号计算语言Maple与MuPAD（从MATLAB 2008b开始使用MuPAD）相结合，使得MATLAB具有符号计算功能。

(2) 绘图功能

可以绘制二维和三维图形。

MATLAB提供了两个层次的绘图操作：

●对图形句柄进行的低层绘图操作；

●建立在低层绘图操作之上的高层绘图操作。

(3) 语言体系

MATLAB具有程序结构控制、函数调用、数据结构、输入输出、面向对象等程序语言特征，而且简单易学、编程效率高。

MATLAB是解释性语言，不能脱离MATLAB环境而独立运行。

(4) MATLAB工具箱

MATLAB包含两部分内容：基本部分和各种可选的工具箱。

MATLAB工具箱分为两大类：功能性工具箱和学科性工具箱。

基本部分

构成MATLAB的核心内容，也是使用和构造工具箱的基础。

功能性工具箱

主要用来扩充其符号计算功能、可视建模仿真功能及文字处理功能等。

学科性工具箱

Control System Toolbox 控制系统工具箱

Signal Processing Toolbox 信号处理工具箱

Neural Network Toolbox 神经网络工具箱

Optimization Toolbox 最优化工具箱

Financial Toolbox 金融工具箱

Statistics Toolbox 统计学工具箱

开始→工具箱

MATLAB具备很强的开放性

除内部函数外，所有MATLAB基本文件和各工具箱文件都是可读、可改的源文件，用户可通过对源文件的修改或加入自己编写的文件去构成新的专用工具箱。

1.1.3 MATLAB 功能演示 p5

例1.1 绘曲线图

绘制函数y=2-|x |和y=sin x 的曲线。

例1.2 求方程的根

求方程2x 5

32例1.3 解线性方程组

2328234453923x y z x y z x y z +-=??

++=??++=?

法2：左除运算

法3：符号计算

例1.4 求积分

1

ln(1)x x dx

+?

法2：符号计算

1.2 MATLAB环境的准备p7

1.2.1 MATLAB的安装

安装MATLAB系统，需运行系统自带的安装程序setup.exe，可以按照安装提示依次操作。

1.2.2 MATLAB的启动与退出p8

启动有3种常见方法：

(1) 使用Windows“开始”菜单。

(2) 运行MATLAB（文件夹）系统启动程序matlab.exe。

(3) 利用（桌面）快捷方式。

退出也有3种常见方法：

(1) 在主窗口File菜单中选择Exit MATLAB。

(2) 在命令窗口中键入Exit或Quit命令。

(3) 单击主窗口右上角的关闭按钮。

1.3 MATLAB操作界面p9

1.3.1 主窗口

MATLAB主窗口是MATLAB的主要工作界面。主窗口除了嵌入一些子窗口外，还主要包括菜单栏和工具栏。

?MATLAB的菜单栏

包括6个菜单项:

File 实现有关文件的操作

Edit 用于命令窗口的编辑操作

Debug 用于程序调试

Desktop 用于设置MATLAB集成环境的显示方式（桌面）

Window 用于关闭所有打开的编辑器窗口或选择活动窗口

Help 用于提供帮助信息。

?MATLAB的工具栏

提供了一些命令按钮和一个当前路径列表框。这些命令按钮有对应的菜单命令，但比菜单命令使用起来更快捷、方便。

1.3.2 命令窗口(Command Window)

命令窗口是MATLAB的主要交互窗口，用于输入命令并显示除图形以外的所有执行结果。

>> 为命令提示符

表示MATLAB正在处于准备状态。

在命令提示符后键入命令并回车后，MATLAB就会解释执行所输入的命令，并在命令后面给出计

该窗口既可以内嵌在工作界面内，也可以浮动到工作界面之外。（右上角的按钮）

一般来说，一个命令行输入一条命令，命令行以回车结束。

但一个命令行也可以输入若干条命令，各命令之间以逗号分隔，若前一命令后带有分号，则逗号可以省略。

例

... 为续行符

如果一个命令行很长，一个物理行之内写不下，可以在第一个物理行之后加上续行符并按下回车键，然后接着下一个物理行继续写命令的其他部分。

即把续行符下面的物理行看作该行的逻辑继续。

表1.1命令行编辑中常用的控制及其功能p10

1.3.3 工作空间窗口(Workspace) p10

工作空间是MATLAB用于存储各种变量和结果的内存空间。

在该窗口中显示工作空间中所有的变量，可对变量进行观察、编辑、保存和删除。

1.3.4 当前目录窗口和搜索路径(Cuurent Directory) p11

当前目录

是指MATLAB运行文件时的工作目录，只有在当前目录或搜索路径下的文件、函数可以被运行或

在当前目录窗口中可以显示或改变当前目录，还可以显示当前目录下的文件并提供搜索功能。

cd命令：将用户目录设置成当前目录。

例将用户目录e:\matlab7\work设置为当前目录。

当用户在MATLAB命令窗口输入一条命令后，MATLAB按照一定次序寻找相关的文件。

基本的搜索过程是：

(1) 检查该命令是不是一个变量。

(2) 检查该命令是不是一个内部函数。

(3) 检查该命令是否当前目录下的M文件。

(4) 检查该命令是否MATLAB搜索路径中其他目录下的M文件。

用户可以将自己的工作目录列入MATLAB搜索路径，从而将用户目录纳入MATLAB系统统一管理。

设置搜索路径的方法有：

(1) 用path命令设置搜索路径。

例将用户目录e:\matlab7\work加到搜索路径下。

在命令窗口输入：

(2) 用对话框设置搜索路径

文件→设置路径（或在命令窗口执行pathtool命令）

将出现搜索路径设置对话框。

通过“添加文件夹”或添加包含子文件夹“按钮将指定路径添加到搜索路径列表中。

在修改完搜索路径后，则需要保存搜索路径。

pathdef.m保存搜索路径的信息，通过修改该文件也可以修改搜索路径。

1.3.5 命令历史窗口(Command History)

在默认设置下，历史记录窗口中会自动保留自安装起所有用过的命令的历史记录，并且还标明了使用时间，从而方便用户查询。

通过双击命令可以再次运行历史命令。

要清除这些历史记录，可选择编辑菜单中的Clear Command History命令。

1.3.6 Start菜单p13

在MATLAB主窗口左下角还有一个Start按钮，单击该按钮会弹出一个菜单，选择其中的命令可以快速访问MATLAB的各种工具和查阅MATLAB包含的各种资源。

1.4 MATLAB帮助系统p13

1.4.1 帮助命令p14

MATLAB帮助命令包括help、lookfor以及模糊查询。

1．help命令

在命令窗口中直接输入：

●help显示当前帮助系统中所包含的所有项目，即搜索路径中所有的目录名称。

●help加函数名显示该函数的帮助说明。

2．lookfor命令

help命令只搜索出那些关键字完全匹配的结果。

lookfor命令只对M文件的第一行进行关键字搜索。若在lookfor命令加上-all选项，则可对M文件进行全文搜索。（列出的内容很长，可用Ctrl+C终止）

3．模糊查询

用户只需要输入命令的前几个字母，然后按Tab键，系统就会列出所有以这几个字母开头的命令。

1.4.2 帮助窗口p15

进入帮助窗口可通过3种方法：

(1) 单击MATLAB主窗口工具栏中的○?按钮。

(2) 在命令窗口中输入helpwin、helpdesk或doc。

(3) 选择Help菜单中MATLAB Help命令。

1.4.3 演示系统

●在帮助窗口中选择演示系统(Demos)选项卡，然后在其中选择相应的演示模块。

●在命令窗口输入Demos。

●选择主窗口Help菜单中的Demos子菜单，打开演示系统。

第2章MATLAB数据及其运算

2.1 MATLAB数据的特点

2.2 变量及其操作

2.3 MATLAB矩阵的表示

2.4 MATLAB数据的运算

2.5 字符串

2.6 结构数据和单元数据

2.1 MATLAB数据的特点p17

矩阵

矩阵是MATLAB最基本、最重要的数据对象。

MATLAB的大部分运算或命令都是在矩阵（在复数域上）运算的意义下执行的。

标量对应仅含一个元素的矩阵。

数值数据

双精度数，占64位，转换函数double

单精度数，占32位，转换函数single

带符号整数，转换函数int8,int16,int32

无符号整数，转换函数uint8,uint16,uint32

字符数据，转换函数char

结构体类型（Structure）

单元类型（Cell）

逻辑型，非0为true，0为false

2.2 变量及其操作

2.2.1 变量与赋值

1．变量命名

变量名是以字母开头，后接字母、数字或下划线的字符序列，最多63个字符。

变量名区分字母的大小写。

MATLAB提供的标准函数名以及命令名必须用小写字母。

2．赋值语句

(1) 变量=表达式

(2) 表达式（值赋给预定义变量ans）

其中表达式是用运算符将有关运算量连接起来的式子，其结果是一个矩阵。

例2.1 计算表达式的值 p18

0的计算结果赋给变量x ，然后显示出结果。

其中pi 和i 都是MATLAB 预先定义的变量，分别代表圆周率π和虚数单位。 3. 预定义变量 p18

表2.1 常用的预定义变量及其含义

2.2.2 变量的管理 p19

1．内存变量的显示与删除

显示工作空间中驻留的变量名清单及信息命令： 删除工作空间中的变量命令：clear 变量名 工作空间窗口专门用于内存变量的管理。

在工作空间窗口中显示所有内存变量的属性。

● 当选中某些变量（[Ctrl 或Shift]+单击）后，再单击Delete 按钮，就能删除这些变量。

● 双击变量后，将进入变量编辑器。通过变量编辑器可以直接观察变量中的具体元素，也可修改变量中的具体元素。

● 输入较大矩阵时，可采用变量编辑器。 2．内存变量文件

利用MAT文件可以把当前工作空间中的一些有用变量长久地保留下来，扩展名是.mat。

MAT文件的生成和装入由save和load命令来完成。

常用格式为：

save 文件名[变量名表] [-append][-ascii]

load 文件名[变量名表] [-ascii]

●文件名可以带路径，但不需带扩展名.mat，命令默认对.mat文件进行操作。

●变量名表中的变量个数不限，只要内存或文件中存在即可，变量名之间以空格分隔。当变量名表

省略时，保存或装入全部变量。

●-ascii选项使文件以ASCII格式处理，省略该选项时文件将以二进制格式处理。

●-append选项控制将变量追加到MAT文件中。

2.2.3 数据的输出格式p20

用十进制数表示一个常数，具体可采用日常记数法和科学记数法两种表示方法。

format命令的格式为：

format 格式符

其中格式符决定数据的输出格式

表2.2控制数据输出格式格式符及其含义

例

2.3 矩阵的表示p21

表矩阵操作函数及其含义

2.3.1 矩阵的建立

1. 直接输入法

最简单的建立矩阵的方法是从键盘直接输入矩阵的元素。具体方法如下：

将矩阵的元素用方括号括起来，按矩阵行的顺序输入各元素，同一行的各元素之间用空格或逗号分隔，不同行的元素之间用分号分隔。（可用回车代替分号）

例

对于比较大且比较复杂的矩阵，可以为它专门建立一个M文件。

例2.2利用M文件建立矩阵

(2) 存盘(文件名为mymatrix.m)。

(3) 在命令窗口中输入mymatrix，即运行该M文件，就建立一个名为MYMAT的矩阵。

3. 建立大矩阵（矩阵拼接）

eye(3) 为3阶单位矩阵

ones(3) 为3阶全1矩阵

2.3.2 冒号表达式

用冒号表达式产生行向量，一般格式：

e1:e2:e3

其中e1为初始值，e2为步长，e3为终止值。

用linspace函数产生行向量。调用格式：

linspace(a,b,n)

其中a和b是生成向量的第一个和最后一个元素，n是元素总数。

2.3.3 矩阵的拆分

1．矩阵元素

通过下标引用矩阵的元素，例如

采用矩阵元素的序号来引用矩阵元素。

矩阵元素的序号是相应元素在内存中的排列顺序。在MATLAB中，矩阵元素按列存储，先第一列，再第二列，依次类推。

(1)利用冒号表达式获得子矩阵

A(:,j)表示取A矩阵的第j列全部元素；

A(i,:)表示A矩阵第i行的全部元素；

A(i,j)表示取A矩阵第i行、第j列的元素；

A(i:i+m,:)表示取A矩阵第i～i+m行的全部元素；

A(:,k:k+m)表示取A矩阵第k～k+m列的全部元素；

A(i:i+m,k:k+m)表示取A矩阵第i～i+m行内，并在第k～k+m列中的所有元素；

A(:)将矩阵A每一列元素堆叠起来，成为一个列向量。

end表示某一维的末尾元素下标。

A(end,:)

A(:,end)

(2)利用空矩阵删除矩阵的元素

给变量X赋空矩阵的语句：

X=[ ]

注意，X=[]与clear X不同，clear是将X从工作空间中删除，而空矩阵则存在于工作空间中，只是维数为0。

2.4 MATLAB数据的运算p26

2.4.1 算术运算

1．基本算术运算

MATLAB的基本算术运算有：

＋加

－减

* 乘

/ 右除

\ 左除

^ 乘方

注意，运算是在矩阵意义下进行的。

(1) 矩阵加减运算

假定有两个矩阵A和B，则可以由A+B和A-B实现矩阵的加减运算。

运算规则是：

●若A和B矩阵的维数相同，则可以执行矩阵的加减运算，A和B矩阵的相应元素相加减。

●若A与B的维数不相同，将给出错误信息，提示用户两个矩阵的维数不匹配。

(2) 矩阵乘法

若A为m×n矩阵，B为n×p矩阵，则

C=A*B

为m×p矩阵。

(3) 矩阵除法

有两种矩阵除法运算：

\ 左除

/ 右除

如果A矩阵是非奇异方阵，则A\B和B/A运算可以实现。

A\B等效于inv(A)*B

B/A等效于B*inv(A)

对于含有标量的运算，两种除法运算的结果相同。

如3/4和4\3有相同的值，都等于0.75。

又如设a=[10.5,25]，则a/5=5\a=[2.1000 5.0000]。

对于矩阵运算，一般A\B≠B/A。

(4) 矩阵的乘方

一个矩阵的乘方运算可以表示成

A^x

要求A为方阵，x为标量。

2．点运算

有一种特殊的运算，因为其运算符是在有关算术运算符前面加点，所以叫点运算。点运算符有

.* 点乘

./ 点右除

.\ 点左除

.^ 点乘方

两矩阵进行点运算是指它们的对应元素进行相关运算，要求两矩阵的维数相同。

3. MATLAB常用数学函数p29

●函数的自变量规定为矩阵变量。

●运算法则是将函数逐项作用于矩阵的元素上。

●结果是一个与自变量同维数的矩阵。

表2.3常用数学函数及其含义p29

函数使用说明：

(1) 三角函数以弧度为单位计算。

(3) 用于取整的函数有

fix向0方向取整

floor向-∞方向取整

ceil向+∞方向取整

round四舍五入到最邻近的整数

(4) rem(x,y)和mod(x,y)要求x,y必须为相同大小的实矩阵或为标量。

rem求余数或模运算

mod模运算

●当y≠0时

rem(x,y)=x-y.*fix(x./y)

mod(x,y)=x-y.*floor(x./y)

●当y=0时

rem(x,0)=NaN（非数）

mod(x,0)=x

2.4.2 关系运算p30

6种关系运算符：

< 小于

<= 小于或等于

> 大于

>= 大于或等于

== 等于

～= 不等于

运算法则：

(1) 两个标量的比较

若关系成立，结果为1，否则为0。

(2) 两个同维数矩阵的比较

对相同位置的元素比较。

(3) 标量与矩阵的比较

标量与矩阵的每一个元素比较。

例2.3建立5阶方阵，判断A的元素是否能被3整除。

2.4.3 逻辑运算p31

逻辑运算符：

& 与

| 或

～非

运算法则：

(1)在逻辑运算中，确认非零元素为真，用1表示，零元素为假，用0表示。

(2) a和b是两个标量

●a&b a,b全为非零时，结果为1，否则为0。

●a|b a,b中只要有一个非零，结果为1，否则为0。

●～a 当a是零时，结果为1；当a非零时，结果为0。

(3) 两个运算量是同维矩阵

相同位置上的元素按标量规则逐个进行。

(4) 一个是标量，另一个是矩阵

标量与矩阵中的每个元素之间按标量规则逐个进行。

(5)逻辑非是单目运算符，也服从矩阵运算规则。

(6)在算术、关系、逻辑运算中，算术运算优先级最高，逻辑运算优先级最低。

例2.4在[0, 3π]区间，求y=sin x的值。

要求：

①消去负半波，即(π, 2π)内的函数值置0。

②

2

,

33

ππ

??

?

??

和

78

,

33

ππ

??

?

??

内取均值为sin

3

π

。

表2.4关系逻辑运算函数及其含义p32

例2.5建立矩阵A，然后找出在[10,20]区间的元素的位置。

2.5 字符串p33

字符串是用单撇号括起来的字符序列。

将字符串当作一个行向量，每个元素对应一个字符，其标识方法和数值向量相同。

可以建立多行字符串矩阵。

字符串是以ASCII码形式存储的。

abs和double函数都可以用来获取字符串矩阵所对应的ASCII码数值矩阵。

相反，char函数可以把ASCII码矩阵转换为字符串矩阵。

例2.6建立一个字符串向量，然后对该向量做如下处理：

(1)取第1～5个字符组成的子字符串。

(2)将字符串倒过来重新排列。

(3)将字符串中的小写字母变成相应的大写字母，其余字符不变。

函数eval，其调用格式为：

eval(t)

其中t为字符串。

作用：把字符串的内容作为对应的MATLAB语句来执行。

表2.5字符串处理函数及其含义

(1)若字符串中含有单撇号，则该撇号用两个单撇号表示。

2.6 结构数据和单元数据p34

2.6.1 结构数据

1．结构矩阵的建立与引用

结构矩阵的元素可以是不同的数据类型，它能将一组具有不同属性的数据纳入到一个统一的变量名下进行管理。

建立一个结构矩阵可采用给结构成员赋值的办法。

格式为：

结构矩阵名.成员名=表达式

其中表达式应理解为矩阵表达式。

例建立含有3个元素的结构矩阵

2．结构成员的修改p35

可以根据需要增加或删除结构的成员。

给a中任意一个元素增加成员x4：

a(1).x4='410075';

但其他成员均为空矩阵，可以使用赋值语句给它赋确定的值。

用rmfield函数删除结构的成员。

a=rmfield(a,'x4');

3．关于结构的函数

表2.6结构的函数及其含义p36

2.6.2 单元数据p36

单元矩阵的各个元素是不同类型的数据，用带有大括号下标的形式引用单元矩阵元素。

1．单元矩阵的建立与引用

建立单元矩阵和一般矩阵相似，只是矩阵元素用大括号括起来。

例建立单元矩阵。

2．关于单元的函数

表2.7单元的函数及其含义p37

《MATLAB程序设计与应用(刘卫国)》(第二版) 答案

《MATLAB程序设计与应用（刘卫国）》（第二版）实验一MATLAB运算基础 1.（1） z1=2*sin(pi*85/180)/(1+exp(2)) （2）x=[2,1+2i;-0.45,5]; z2=log(x+sqrt(1+x.^2))/2 （3）a=-3.0:0.1:3.0; z3=0.5*(exp(0.3*a)-exp(-0.3*a)).*sin(a+0.3)+log((0.3+a)/2) （4）t=0:0.5:2.5; z4=(t>=0&t<1).*t.^2+(t>=1&t<2).*(t.^2-1)+(t>=2&t<3).*(t.^2-2*t+1) 2. A=[12,34,-4;34,7,87;3,65,7]; B=[1,3,-1;2,0,3;3,-2,7]; （1）a=A+6*B b=A-B+eye(size(A)) %I=eye(size(A)) （2）c=A*B d=A.*B （3）e=A^3 f=A.^3 （4）g=A/B h=B\A （5）m=[A,B] n=[A([1,3],:);B^2] 3. A=[1,2,3,4,5;6,7,8,9,10;11,12,13,14,15;16,17,18,19,20;21,22,23,24,25]; B=[3,0,16;17,-6,9;0,23,-4;9,7,0;4,13,11]; （1）C=A*B （2）D=C(3:end,2:end) 4.（1） a=100:999; b=rem(a,21)==0; c=find(b); d=length(c) （2）ch='Just as Bianhaiman said,Xiehong is ...'; e=find(ch>='A'&ch<='Z'); ch(e)=[] 实验二 MATLAB矩阵分析与处理 1.E=eye(3);

Matlab编程与应用习题和一些参考答案

Matlab 上机实验一、二 3.求下列联立方程的解???????=+-+-=-+=++-=--+4 1025695842475412743w z y x w z x w z y x w z y x >> a=[3 4 -7 -12;5 -7 4 2;1 0 8 -5;-6 5 -2 10]; >> b=[4;4;9;4]; >> c=a\b 4．设???? ??????------=81272956313841A ，??????????-----=793183262345B ，求C1=A*B’;C2=A’*B;C3=A.*B,并求上述所有方阵的逆阵。 >> A=[1 4 8 13;-3 6 -5 -9;2 -7 -12 -8]; >> B=[5 4 3 -2;6 -2 3 -8;-1 3 -9 7]; >> C1=A*B' >> C2=A'*B >> C3=A.*B >> inv(C1) >> inv(C2) >> inv(C3) 5．设 ?? ????++=)1(sin 35.0cos 2x x x y ，把x=0~2π间分为101点，画出以x 为横坐标，y 为纵坐标的曲线。 >> x=linspace(0,2*pi,101); >> y=cos(x)*(0.5+(1+x.^2)\3*sin(x)); >> plot(x,y,'r') 6．产生8×6阶的正态分布随机数矩阵R1, 求其各列的平均值和均方差。并求该矩阵全体数的平均值和均方差。 （mean var ） a=randn(8,6) mean(a) var(a) k=mean(a) k1=mean(k) i=ones(8,6) i1=i*k1 i2=a-i1 i3=i2.*i2 g=mean(i3) g2=mean(g)

刘卫国版MATLAB程序设计与应用课后实验六八九

实验六 高层绘图操作 %第一题： 程序代码如下： x=linspace(0,2*pi,101); y=(0.5+3*sin(x)./(1+x.^2)).*cos(x); plot(x,y) 01234567 -1 -0.5 0.5 1 1.5 %第二题： %（1） 程序代码如下： x=linspace(-2*pi,2*pi,100); y1=x.^2; y2=cos(2*x); y3=y1.*y2; plot(x,y1,'b-',x,y2,'r:',x,y3,'y--'); text(4,16,'\leftarrow y1=x^2'); text(6*pi/4,-1,'\downarrow y2=cos(2*x)'); text(-1.5*pi,-2.25*pi*pi,'\uparrow y3=y1*y2');

-8 -6 -4 -2 2 4 6 8 -30-20 -10 10 20 30 40 %（2） 程序代码如下： x=linspace(-2*pi,2*pi,100); y1=x.^2; y2=cos(2*x); y3=y1.*y2; subplot(1,3,1);%分区 plot(x,y1); title('y1=x^2');%设置标题 subplot(1,3,2); plot(x,y2); title('y2=cos(2*x)'); subplot(1,3,3); plot(x,y3); title('y3=x^2*cos(2*x)');

-10 10 0510 15202530 35 40y1=x 2 -10 10 -1-0.8 -0.6 -0.4-0.200.20.4 0.6 0.8 1y2=cos(2*x) -10 10 -30-20 -10 10 20 30 40 y3=x 2*cos(2*x) %（3） 程序代码如下： x=linspace(-2*pi,2*pi,20); y1=x.^2; subplot(2,2,1);%分区 bar(x,y1); title('y1=x^2的条形图');%设置标题 subplot(2,2,2); stairs(x,y1); title('y1=x^2的阶梯图'); subplot(2,2,3); stem(x,y1); title('y1=x^2的杆图'); subplot(2,2,4); fill(x,y1,'r');%如果少了'r'则会出错 title('y1=x^2的填充图'); %其他的函数照样做。

matlab程序设计与应用第二版习题答案

matlab程序设计与应用第二版习题答案【篇一：matlab程序设计与应用(第二版)实验答案】 %实验一 matlab运算基础 %第1题 %（1） z1=2*sin(85*pi/180)/(1+exp(2)) %（2） x=[2,1+2i;-0.45,5]; z2=0.5*log(x+sqrt(1+x.^2)) %(3) a=-3.0:0.1:3.0; z3=(exp(0.3*a)-exp(-0.3*a))/2.*sin(a+0.3)+log((0.3+a)/2) %(4) t=0:0.5:2.5; z4=t.^2.*(t=0t1)+(t.^2-1).*(t=1t2)+(t.^2-2*t+1).*(t=2t3) %第2题 a=[12 34 -4;34 7 87;3 65 7]; b=[1 3 -1;2 0 3;3 -2 7]; a+6*b a-b+eye(size(a)) a*b a.*b a^3 a.^3 a/b b\a [a,b] [a([1,3],:);b^2] %第3题 a=[1 2 3 4 5;6 7 8 9 10;11 12 13 14 15;16 17 18 19 20;21 22 23 24 25] b=[3 0 16;17 -6 9;0 23 -4;9 7 0;4 13 11] c=a*b f=size(c) d=c(f(1)-2:f(1),f(2)-1:f(2)) whos %第4题 %(1):

a=100:999; b=rem(a,21); c=length(find(b==0)) %(2): a=lsdhksdlkklsdkl; k=find(a=aa=z); a(k)=[] %实验二 matlab矩阵分析与处理 %第1题 e=eye(3); r=rand(3,2); o=zeros(2,3); s=diag([2,3]); a=[e,r;o,s]; a^2 b=[e,(r+r*s);o,s^2] %第2题 h=hilb(5) p=pascal(5) hh=det(h) hp=det(p) th=cond(h) tp=cond(p) %第3题 a=fix(10*rand(5)) h=det(a) trace=trace(a) rank=rank(a) norm=norm(a) %第4题 a=[-29,6,18;20,5,12;-8,8,5] [v,d]=eig(a) %数学意义略 %第5题方法一 %(1): a=[1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6]; b=[0.95,0.67,0.52]; x=inv(a)*b %(2):

MATLAB程序设计与应用(第二版)

阅阅读读时时：：请请选选择择““视视图图||文文档档结结构构图图””，，弹弹出出文文档档中中的的标标题题链链接接。。数学软件 MATLAB 程序设计与应用

第1章MATLAB系统环境 1.1 MATLAB概貌 1.2 MATLAB环境的准备 1.3 MATLAB操作界面 1.4 MATLAB帮助系统 自上世纪80年代以来，出现了科学计算语言，亦称数学软件。 MATLAB Mathematica Mathcad Maple LINDO LINGO 1.1 MATLAB概貌p3 MATLAB 是MATrix LABoratory（矩阵实验室）的缩写。 1984年由Math Works 公司推出，现已成为国际公认的优秀的工程应用开发环境，是影响最大，流行最广的科学计算语言。 1.1.1 MATLAB的发展

1. 从MATLAB 4.2c开始，每个版本增加了一个建造编号； 2. 例如MATLAB7.6的建造编号是R2008a。说明MATLAB7.6与MATLAB2008a是等同的； 3. 对于建造编号，正规化以后，每年出两个版本。一般来说。a是测试版，b是正式版。a是前半 年出，b是后半年出。 教材采用MATLAB7.0（R14，2004） 实验室采用MATLAB 7.8（R2009a，2009.3，汉化） 1.1.2 MATLAB的主要功能p4 ◆数值计算和符号计算功能 ◆绘图功能 ◆语言体系 ◆MATLAB工具箱 (1) 数值计算和符号计算功能 MATLAB以矩阵作为数据操作的基本单位，还提供了十分丰富的数值计算函数。 MATLAB先后和著名的符号计算语言Maple与MuPAD（从MATLAB 2008b开始使用MuPAD）相结合，使得MATLAB具有符号计算功能。 (2) 绘图功能 可以绘制二维和三维图形。 MATLAB提供了两个层次的绘图操作： ●对图形句柄进行的低层绘图操作； ●建立在低层绘图操作之上的高层绘图操作。 (3) 语言体系 MATLAB具有程序结构控制、函数调用、数据结构、输入输出、面向对象等程序语言特征，而且简单易学、编程效率高。 MATLAB是解释性语言，不能脱离MATLAB环境而独立运行。 (4) MATLAB工具箱 MATLAB包含两部分内容：基本部分和各种可选的工具箱。 MATLAB工具箱分为两大类：功能性工具箱和学科性工具箱。 基本部分 构成MATLAB的核心内容，也是使用和构造工具箱的基础。 功能性工具箱 主要用来扩充其符号计算功能、可视建模仿真功能及文字处理功能等。 学科性工具箱 Control System Toolbox 控制系统工具箱 Signal Processing Toolbox 信号处理工具箱 Neural Network Toolbox 神经网络工具箱 Optimization Toolbox 最优化工具箱 Financial Toolbox 金融工具箱 Statistics Toolbox 统计学工具箱 开始→工具箱 MATLAB具备很强的开放性 除内部函数外，所有MATLAB基本文件和各工具箱文件都是可读、可改的源文件，用户可通过对源文件的修改或加入自己编写的文件去构成新的专用工具箱。

MATLAB程序设计及应用(第二版)课后实验答案

。 Matlab课后实验题答案 《 实验一 MATLAB运算基础 1. 先求下列表达式的值，然后显示MATLAB工作空间的使用情况并保存全部变量。 (1) 0 12 2sin85 1 z e = + (2) 21 ln( 2 z x =，其中 212 0.455 i x + ??=?? -?? (3) 0.30.3 3 0.3 sin(0.3)ln, 3.0, 2.9,,2.9,3.0 22 a a e e a z a a - -+ =++=--

(4) 22 42011 122123t t z t t t t t ?≤

) 3. 设有矩阵A 和B 1234 53 166789101769,11 121314150 23416171819209 7021222324254 1311A B ???? ????-??? ?????==-??? ? ???????????? (1) 求它们的乘积C 。 (2) 将矩阵C 的右下角3×2子矩阵赋给D 。 (3) 查看MATLAB 工作空间的使用情况。 4. 完成下列操作： (1) 求[100,999]之间能被21整除的数的个数。 (2) 建立一个字符串向量，删除其中的大写字母。 解：(1) 结果：

MATLAB程序设计与应用课后习题答案

西安科技大学MATLAB程序设计 专业：信息与计算科学 班级：1001班 学号：1008060129 姓名：刘仲能 2012年6月27日

实验一 2.已知： ??? ? ? ?????-= 765 3 8773443412A ，???? ? ?????--=72 3 302131 B 求下列表达式的值： （1）A+6*B 和A-B+I （其中I 为单位矩阵） （2）A*B 和A.*B （3）A^3和A.^3 （4）A/B 及B\A （5）[A,B]和 [A([1,3],:);B^2]

3.设有矩阵A 和B ????? ?? ? ????????= 2524 23 22 21 2019181716151413121110987654321A ，??????? ? ????????--=1113 4 079423096171603 B （1） 求它们的乘积 C 。 （2） 将矩阵C 的右下角3×2子矩阵赋给D 。 （3） 查看MATLAB 工作空间的使用情况 (1) (2) （3）

4.完成下列操作 （1）求[100,999]之间能被21整除的数的个数。（2）建立一个字符串向量，删除其中的大写字母。 (1)(2)

实验二 3.建立一个5×5矩阵，求它的行列式值、迹、秩和范数。 运行截图： A 矩阵的行列式值、迹、秩分别如下： 范数如下： 4.已知 ???? ? ?????--= 58 8 1252018629A 求A 的特征值及特征向量，并分析其数学意义。 运行截图：

5.下面是一个线性方程组：???? ??????= ???? ? ????????? ? ?????52.067.095.06/15 /14 /15 /14/13 /14/13/12 /1321x x x （1）求方程的解； （2）将方程右边向量元素改为0.53，在求解，并比较的变化和解的相对 变化； （3）计算系数矩阵A 的条件数并分析结论。 （2） 变大，其解中，相对未变化前的的解：x1变大，x2变小，x3变大。 （3） 由于A 矩阵的条件数很大，故当线性方程组中的b 变大时，x 也将发生很大的变

Matlab程序设计与应用第二版刘卫国

实验 T1： %%第一小题 z1=2*sin(85*pi/180)/(1+exp(2)) %%第二小题x=[2,1+2i;-0.45,5]; z2 %%第三小题 a=-3.0:0.1:3.0; 0.3)+log((0.3+a)/2) %%第四题 t=0:0.5:2.5 z4=(t>=0&t<1).*(t.A2)+(t>=1 &t<2).*(t 八2-1)+(t>=2&t<3).*(t.A2-2*t+1) T2： A=[12,34,-4;34,7,87;3,65,7] B=[1,3,-1;2,0,3;3,-2,7] disp ('A+6*B='); disp(A+6*B); disp('A-B+I=');disp(A-B+eye(3)); disp('A*B=');

disp(A*B); disp('A.*B='); disp(A.*B); disp('A A3='); disp(AT); disp('A.A3二')； disp(A.A3); disp('A/B='); disp(A/B); disp('B\A='); disp(B\A); disp('[A,B]='); disp([A,B]); disp('[A([1,3],:);BA2]='); disp([A([1,3],:);BA2]); T3: z=1:25; A=reshape(z,5,5)'; B=[3,0,16;17,-6,9;0,23,-4;9,7,0;4,13,11];C=A*BD=C(3:5,2:3) T4-1: a=100:999; b=find(rem(a,21)==0); c=length(b)

MATLAB程序设计与应用(第二版)实验参考答案

MATLAB程序设计与应用（第二版）实验参考答案 %实验一MATLAB运算基础 %第一题 %（1） z1=2*sin(85*pi/180)/(1+exp(2)) %（2） x=[2,1+2i;-0.45,5]; z2=0.5*log(x+sqrt(1+x.^2)) %(3) a=-3.0:0.1:3.0; z3=(exp(0.3*a)-exp(-0.3*a))/2.*sin(a+0.3)+log((0.3+a)/2) %(4) t=0:0.5:2.5; z4=t.^2.*(t>=0&t<1)+(t.^2-1).*(t>=1&t<2)+(t.^2-2*t+1).*(t>=2&t<3) %第二题 A=[12 34 -4;34 7 87;3 65 7]; B=[1 3 -1;2 0 3;3 -2 7]; A+6*B A-B+eye(size(A)) A*B A.*B A^3 A.^3 A/B B\A [A,B] [A([1,3],:);B^2] %第三题 A=[1 2 3 4 5;6 7 8 9 10;11 12 13 14 15;16 17 18 19 20;21 22 23 24 25] B=[3 0 16;17 -6 9;0 23 -4;9 7 0;4 13 11] C=A*B F=size(C) D=C(F(1)-2:F(1),F(2)-1:F(2)) whos %第四题 %(1): A=100:999; B=rem(A,21); C=length(find(B==0)) %(2): A='lsdhKSDLKklsdkl';

k=find(A>='A'&A<='Z'); A(k)=[] %实验二MATLAB矩阵分析与处理 %第一题 E=eye(3); R=rand(3,2); O=zeros(2,3); S=diag([2,3]); A=[E,R;O,S]; A^2 B=[E,(R+R*S);O,S^2] %第二题 H=hilb(5) P=pascal(5) Hh=det(H) Hp=det(P) Th=cond(H) Tp=cond(P) %第三题: A=fix(10*rand(5)) H=det(A) Trace=trace(A) Rank=rank(A) Norm=norm(A) %第四题: A=[-29,6,18;20,5,12;-8,8,5] [V,D]=eig(A) %数学意义略 %第五题方法一: %(1): A=[1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6]; b=[0.95,0.67,0.52]'; x=inv(A)*b %(2): B=[0.95,0.67,0.53]'; x=inv(A)*B %(3): cond(A) %第五题方法二: A=hilb(4) A(:,1)=[] A(4,:)=[] B=[0.95,0.67,0.52]';

(完整word版)MATLAB程序设计及应用(第二版)课后实验答案

Matlab课后实验题答案 实验一 MATLAB运算基础 1. 先求下列表达式的值，然后显示MATLAB工作空间的使用情况并保存全部变量。 (1) 0 12 2sin85 1 z e = + (2) 21 ln( 2 z x =+，其中 212 0.455 i x + ??=?? -?? (3) 0.30.3 3 0.3 sin(0.3)ln, 3.0, 2.9,,2.9,3.0 22 a a e e a z a a - -+ =++=--L (4) 2 2 4 2 01 112 2123 t t z t t t t t ?≤< ? =-≤< ? ?-+≤< ? ，其中t=0:0.5:2.5 2. 已知：

1234413134787,2033657327A B --???? ????==???? ????-???? 求下列表达式的值： (1) A+6*B 和A-B+I （其中I 为单位矩阵） (2) A*B 和A.*B (3) A^3和A.^3 (4) A/B 及B\A (5) [A,B]和[A([1,3],:);B^2] 解： 3. 设有矩阵A 和B 1234 53 166789101769,11 121314150 23416171819209 7021222324254 1311A B ???? ????-??? ?????==-??? ? ???????????? (1) 求它们的乘积C 。 (2) 将矩阵C 的右下角3×2子矩阵赋给D 。 (3) 查看MATLAB 工作空间的使用情况。

4. 完成下列操作： (1) 求[100,999]之间能被21整除的数的个数。 (2) 建立一个字符串向量，删除其中的大写字母。 解：(1) 结果： (2). 建立一个字符串向量 例如： ch='ABC123d4e56Fg9';则要求结果是： 实验二 MATLAB 矩阵分析与处理 1. 设有分块矩阵33 322322E R A O S ?????? =? ??? ，其中E 、R 、O 、S 分别为单位矩阵、随机矩阵、零矩

电子书_MATLAB程序设计与应用(第二版)

阅阅读读时时：：请请选选择择““视视图图 || 文文档档结结构构图图””，，弹弹出出文文档档中中的的标标题题链链接接。。 数学软件 MATLAB 程序设计与应用

第1章MATLAB系统环境 1.1 MATLAB概貌 1.2 MATLAB环境的准备 1.3 MATLAB操作界面 1.4 MATLAB帮助系统 自上世纪80年代以来，出现了科学计算语言，亦称数学软件。 MATLAB Mathematica Mathcad Maple LINDO LINGO 1.1 MATLAB概貌p3 MATLAB 是MATrix LABoratory（矩阵实验室）的缩写。 1984年由Math Works 公司推出，现已成为国际公认的优秀的工程应用开发环境，是影响最大，流行最广的科学计算语言。 1.1.1 MATLAB的发展 2. 例如MATLAB7.6的建造编号是R2008a。说明MATLAB7.6与MATLAB2008a是等同的；

3. 对于建造编号，正规化以后，每年出两个版本。一般来说。a是测试版，b是正式版。a是前半 年出，b是后半年出。 教材采用MATLAB7.0（R14，2004） 实验室采用MATLAB 7.8（R2009a，2009.3，汉化） 1.1.2 MATLAB的主要功能p4 ◆数值计算和符号计算功能 ◆绘图功能 ◆语言体系 ◆MATLAB工具箱 (1) 数值计算和符号计算功能 MATLAB以矩阵作为数据操作的基本单位，还提供了十分丰富的数值计算函数。 MATLAB先后和著名的符号计算语言Maple与MuPAD（从MATLAB 2008b开始使用MuPAD）相结合，使得MATLAB具有符号计算功能。 (2) 绘图功能 可以绘制二维和三维图形。 MATLAB提供了两个层次的绘图操作： ●对图形句柄进行的低层绘图操作； ●建立在低层绘图操作之上的高层绘图操作。 (3) 语言体系 MATLAB具有程序结构控制、函数调用、数据结构、输入输出、面向对象等程序语言特征，而且简单易学、编程效率高。 MATLAB是解释性语言，不能脱离MATLAB环境而独立运行。 (4) MATLAB工具箱 MATLAB包含两部分内容：基本部分和各种可选的工具箱。 MATLAB工具箱分为两大类：功能性工具箱和学科性工具箱。 基本部分 构成MATLAB的核心内容，也是使用和构造工具箱的基础。 功能性工具箱 主要用来扩充其符号计算功能、可视建模仿真功能及文字处理功能等。 学科性工具箱 Control System Toolbox 控制系统工具箱 Signal Processing Toolbox 信号处理工具箱 Neural Network Toolbox 神经网络工具箱 Optimization Toolbox 最优化工具箱 Financial Toolbox 金融工具箱 Statistics Toolbox 统计学工具箱 开始→工具箱 MATLAB具备很强的开放性 除内部函数外，所有MATLAB基本文件和各工具箱文件都是可读、可改的源文件，用户可通过对源文件的修改或加入自己编写的文件去构成新的专用工具箱。

Matlab程序设计与应用第二版刘卫国

实验一：T1： %%第一小题 z1=2*sin(85*pi/180)/(1+exp(2)) %%第二小题 x=[2,1+2i;-0.45,5]; z2 %%第三小题 a=-3.0:0.1:3.0; 0.3)+log((0.3+a)/2) %%第四题 t=0:0.5:2.5 z4=(t>=0&t<1).*(t.^2)+(t>=1&t<2).*(t .^2-1)+(t>=2&t<3).*(t.^2-2*t+1) T2： A=[12,34,-4;34,7,87;3,65,7] B=[1,3,-1;2,0,3;3,-2,7] disp ('A+6*B='); disp(A+6*B); disp('A-B+I=');disp(A-B+eye(3)); disp('A*B=');

disp(A*B); disp('A.*B='); disp(A.*B); disp('A^3='); disp(A^3); disp('A.^3='); disp(A.^3); disp('A/B='); disp(A/B); disp('B\A='); disp(B\A); disp('[A,B]='); disp([A,B]); disp('[A([1,3],:);B^2]='); disp([A([1,3],:);B^2]); T3: z=1:25; A=reshape(z,5,5)'; B=[3,0,16;17,-6,9;0,23,-4;9,7,0;4,13,11];C=A*BD=C(3:5,2:3) T4-1: a=100:999;

b=find(rem(a,21)==0); c=length(b) T4-2: a=input('请输入一个字符串：','s'); b=find(a>='A'&a<='Z'); a(b)=[]; disp(a); 实验二： T1： E=eye(3),R=rand(3,2),O=zeros(2,3),S=d iag([1,2]); A=[E,R;O,S] disp('A^2='); disp(A^2); disp('[E,R+RS;O,S^2]'); B=[E,R+R*S;O,S^2] T2： H=hilb(5) P=pascal(5)Hh=det(H) Hp=det(P) Th=cond(H)

最新MATLAB程序设计与应用(第二版)课后实验答案

Matlab 课后实验题答案 实验一 MATLAB 运算基础 1. 先求下列表达式的值，然后显示MATLAB 工作空间的使用情况并保存全部变量。 (1) 0 12 2sin851z e =+ (2) 221 ln(1)2z x x = ++，其中2120.45 5i x +??=??-?? (3) 0.30.330.3sin(0.3)ln , 3.0, 2.9,,2.9,3.022a a e e a z a a --+= ++=-- (4) 22 42011 122123t t z t t t t t ?≤=0&t<1).*(t.^2)+(t>=1&t<2).*(t.^2-1)+(t>=2&t<3) .*(t.^2-2*t+1) 2. 已知：

1234413134787,2033657327A B --???? ????==???? ????-???? 求下列表达式的值： (1) A+6*B 和A-B+I （其中I 为单位矩阵） (2) A*B 和A.*B (3) A^3和A.^3 (4) A/B 及B\A (5) [A,B]和[A([1,3],:);B^2] 解： M 文件： A=[12 34 -4;34 7 87;3 65 7];B=[1 3 -1;2 0 3;3 -2 7]; A+6.*B A-B+eye(3) A*B A.*B A^3 A.^3 A/B B\A [A,B] [A([1,3],:);B^2] 3. 设有矩阵A 和B 1234 53 166789101769,11 121314150 23416171819209 7021222324254 1311A B ???? ????-??? ?????==-??? ? ???????????? (1) 求它们的乘积C 。 (2) 将矩阵C 的右下角3×2子矩阵赋给D 。 (3) 查看MATLAB 工作空间的使用情况。 解：. 运算结果： E=(reshape(1:1:25,5,5))';F=[3 0 16;17 -6 9;0 23 -4;9 7 0;4 13 11]; C= E*F H=C(3:5,2:3) C = 93 150 77

MATLAB程序设计与应用

实验一 MATLAB 运算基础 第二题：已知 A=??????????7653877344-3412 B=??? ? ??????72-33021-31 求下列问题： 1) A+6*B 和A-B+I 2) A*B 和A.*B 3) A^3 A.^3 4) A/B B\A 5) [A,B] [A(1,3),:B.^2] 解： >> A=[12 34 -4;34 7 87;3 65 7]; B=[1 3 -1;2 0 3;3 -2 7]; I=[1 0 0;0 1 0;0 0 1]; （1） >> A+6*B ans = 18 52 -10 46 7 105 21 53 49 >> A-B+I ans = 12 31 -3 32 8 84 0 67 1 （2） >> A*B ans = 68 44 62 309 -72 596 154 -5 241 >> A.*B

ans = 12 102 4 68 0 261 9 -130 49 （3） >> A^3 ans = 37226 233824 48604 247370 149188 600766 78688 454142 118820 >> A.^3 ans = 1728 39304 -64 39304 343 658503 27 274625 343 （4） >> A/B ans = 16.4000 -13.6000 7.6000 35.8000 -76.2000 50.2000 67.0000 -134.0000 68.0000 >> B\A ans = 109.4000 -131.2000 322.8000 -53.0000 85.0000 -171.0000 -61.6000 89.8000 -186.2000 （5） >> [A,B] ans =

matlab程序设计与应用(第二版)第三章部分课后答案解析

第三章1. (1A=eye(3 (2A=100+100*rand(5,6 (3A=1+sqrt(0.2*randn(10,50 (4B=ones(size(A (5A+30*eye(size(A (6B=diag(diag(A 2. B=rot90(A C=rot90(A,-1 3. B=inv(A ;A的逆矩阵 C=det(A ;A的行列式的值 D=A*B E=B*A D=E 因此A与A-1是互逆的。 4. A=[4 2 -1;3 -1 2;12 3 0]; b=[2;10;8];

x=inv(A*b x = -6.0000 26.6667 27.3333 5. (1 diag(A ;主对角线元素ans = 1 1 5 9 triu(A ;上三角阵 ans = 1 -1 2 3 0 1 -4 2 0 0 5 2

0 0 0 9 tril(A ;下三角阵ans = 1 0 0 0 5 1 0 0 3 0 5 0 11 15 0 9 rank(A ;秩ans = 4 norm(A ;范数 ans = 21.3005 cond(A ;条件数ans = 11.1739 trace(A ;迹 ans = 16 (2略

6. A=[1 1 0.5;1 1 0.25;0.5 0.25 2] A = 1.0000 1.0000 0.5000 1.0000 1.0000 0.2500 0.5000 0.2500 2.0000 [V,D]=eig(A V = 0.7212 0.4443 0.5315 -0.6863 0.5621 0.4615 -0.0937 -0.6976 0.7103 D = -0.0166 0 0 0 1.4801 0 0 0 2.5365

MATLAB程序设计和应用课后习题答案解析

专业技术资料分享 西安科技大学MATLAB程序设计 专业：信息与计算科学 班级： 1001班 学号：1008060129 姓名：刘仲能 2012年6月27日

实验一 2.已知： ?? ??? ?????-=76538773443412A ，???? ? ?????--=723302131B 求下列表达式的值： （1）A+6*B 和A-B+I （其中I 为单位矩阵） （2）A*B 和A.*B （3）A^3和A.^3 （4）A/B 及B\A （5）[A,B]和 [A([1,3],:);B^2]

3.设有矩阵A 和B ????? ?? ?????????=25242322212019181716151413121110987654321A ，??????? ?????????--= 11134079423096171603B （1） 求它们的乘积C 。 （2） 将矩阵C 的右下角3×2子矩阵赋给D 。 （3） 查看MATLAB 工作空间的使用情况 (1) (2) （3）

4.完成下列操作 （1）求[100,999]之间能被21整除的数的个数。 （2）建立一个字符串向量，删除其中的大写字母。 (1) (2) 实验二 3.建立一个5×5矩阵，求它的行列式值、迹、秩和范数。

运行截图： A 矩阵的行列式值、迹、秩分别如下： 范数如下： 4.已知 ?? ?? ? ?????--=5881252018629A 求A 的特征值及特征向量，并分析其数学意义。 运行截图：

5.下面是一个线性方程组：?? ??? ?????=???? ????????????????52.067.095.06/15/14/15/14/13/14/13/12/1321x x x （1） 求方程的解； （2） 将方程右边向量元素 改为0.53 ，在求解，并比较 的变化和解的相 对变化； （3） 计算系数矩阵A 的条件数并分析结论。 （2） 变大，其解中，相对未变化前的 的解：x1变大，x2变小，x3变大。 （3） 由于A 矩阵的条件数很大，故当线性方程组中的b 变大时，x 也将发生很大的变化，即数值稳定性较差。 实验三 3.硅谷公司员工的工资计算方法如下： （1）工作时数超过120小时者，超过部分加发15%； （2）工作时数低于60小时者，扣发700元； （3）其余按每小时84元计发。

Matlab程序设计及应用实验讲义

Matlab程序设计及应用实验讲义 自编 电子科学与工程系 2013.02

实验一 MATLAB 环境与命令窗口 1、实验目的 1）熟悉MATLAB 的操作环境及基本操作方法； 2）掌握MATLAB 的搜索路径及其设置方法； 3）熟悉MATLAB 帮助信息的查阅方法； 2、实验主要仪器设备和材料 计算机PC 一台 2010a 软件 3、实验内容和原理 1、先建立自己的工作目录，再将自己的工作目录设置到MATLAB 的搜索路径下，再试验用help 命令能否查询到自己的工作目录。 2、在MATLAB 环境下验证下面几个例子，并总结MATLAB 的主要优点。 1）绘制正弦曲线和余弦曲线 2）求方程432379230x x x ++-=的全部根 3）求积分()1 0ln 1x x dx +? 4）求解线性方程组234832245917x y z x y z x y z -+=?? ++=??+-=? 3、利用MATLAB 的帮助功能分别查询inv 、plot 、 max 、round 等函数的功能与用法 4、完成下列操作： 1）在MATLAB 命令窗口输入以下命令： x=0:pi/10:2*pi;

y=sin(x); 2）在工作空间窗口选择变量y，再在工作空间窗口选择绘图菜单命令或在工具栏中单击绘图命令按钮，绘制变量的图形，并分析图形的含义。 5、访问MathsWorks公司的主页，查询有关MATLAB的产品信息。 思考与练习 1、如何启动和退出MATLAB的集成环境？ 2、简述MATLAB的主要功能。 3、如果一个MATLAB命令包含的字符很多，需要分成多行输入，该如何处理？ 4、help命令和look for命令有何区别？ 5、在MATLAB环境下，建立一个变量fac,同时又在当前目录下建立了一个M 文件fac.m，如果需要运行fac.m文件，该如何处理？

MATLAB程序设计与应用(刘卫国编)课后实验答案

实验一 MATLAB 运算基础 1. 先求下列表达式的值，然后显示MATLAB 工作空间的使用情况并保存全部变量。 (1) 012 2sin851z e =+ (2) 21ln(2 z x =+，其中2 120.45 5i x +??=? ?-?? (3) 0.30.330.3sin(0.3)ln , 3.0, 2.9,,2.9,3.022a a e e a z a a --+= ++=-- (4) 2242011 122123t t z t t t t t ?≤

- 2. 已知： 1234413134787,2033657327A B --???? ????==???? ????-???? 求下列表达式的值： (1) A+6*B 和A-B+I （其中I 为单位矩阵） (2) A*B 和A.*B (3) A^3和A.^3 (4) A/B 及B\A

(5) [A,B]和[A([1,3],:);B^2] @ 解： 运算结果：

3. 设有矩阵A和B

1234 53 166789101769,11 121314150 23416171819209 7021222324254 1311A B ???? ????-??? ?????==-??? ? ???????????? (1) 求它们的乘积C 。 (2) 将矩阵C 的右下角3×2子矩阵赋给D 。 ( (3) 查看MATLAB 工作空间的使用情况。 解：. 运算结果： ! 4. 完成下列操作： (1) 求[100,999]之间能被21整除的数的个数。 (2) 建立一个字符串向量，删除其中的大写字母。 解：(1) 结果：