matlab数值计算课后答案

第5章 MATLAB数值计算

MATLAB的数学计算＝数值计算＋符号计算

其中符号计算是指使用未定义的符号变量进行运算，而数值计算不允许使用未定义的变量。

注：数值计算有误差。数值计算受到计算机所保留的有效位数的限制，因此每一次运算都会有一定的误差，重复的多次数值计算就可能会造成很大的累积误差。

5.1 矩阵和数组

MATLAB最基本也是最重要的功能就是进行实数或复数矩阵的运算。

矩阵和数组的概念

在MATLAB的运算中，经常要使用标量、向量、矩阵和数组，这几个名称的定义如下： ? 标量：是指1×1的矩阵，即为只含一个数的矩阵。 ? 向量：是指1×n或n×1的矩阵，即只有一行或者一列的矩阵。

? 矩阵：是一个矩形的数组，即二维数组，其中向量和标量都是矩阵的特例，0×0矩阵为空矩阵([])。

? 数组：是指n维的数组，为矩阵的延伸，其中矩阵和向量都是数组的特例。

5.1.1矩阵输入（前已讲“向量的输入”，相同！）

(1) 矩阵元素应用方括号([])括住；

(2) 每行内的元素间用逗号或空格隔开； (3) 行与行之间用分号或回车键隔开； (4) 元素可以是数值或表达式。

1. 通过显式元素列表输入矩阵

c=[1 2;3 4;5 3*2]

精通MATLAB语言，有MATLAB编程问题的朋友，请直接联系我！ 联系方式： QQ：993878382

Email：993878382@

MATLAB博客：/46036333_d.html

如果我QQ不在线，可以将问题发到我的邮箱，或者在我的MATLAB博客留言，第一时间答复你！简单问题请直接留言，复杂问题可以提供解决思路！同探讨，共进步！

还可以为公司和科研单位设计各类算法，欢迎前来咨询！

郑重声明：本资源来源于网络，仅限交流使用，请勿用于商业用途！

Matlab与化工数值计算第一讲简介与基本数学运算隋志军化工学院软件应用教科组联系方式：zhjsui@ 64252169实验16楼605室

化学工程师的任务工厂运行系统特性过程设计数学模型反应特征过程特性设备特性工艺开发数学模型

化学工程专业数学模型类型非线性方程(组)1 f= 4.0 0.4 lg(Re gen f 1 n ' 2 ) 1.2 n' 0.75 n'

常微分方程(组)

f" '+ ff"+β (1 ( f ' ) 2 )= 0偏微分方程(组) u1 2 u1= 0.024 2 F (u1 u 2 ) t x

非线性模型，难以获得解析解，必须采用数值解法模型的数值解法

精通MATLAB语言，有MATLAB编程问题的朋友，请直接联系我！ 联系方式： QQ：993878382

Email：993878382@

MATLAB博客：/46036333_d.html

如果我QQ不在线，可以将问题发到我的邮箱，或者在我的MATLAB博客留言，第一时间答复你！简单问题请直接留言，复杂问题可以提供解决思路！同探讨，共进步！

还可以为公司和科研单位设计各类算法，欢迎前来咨询！

郑重声明：本资源来源于网络，仅限交流使用，请勿用于商业用途！

Matlab与化工数值计算第一讲简介与基本数学运算隋志军化工学院软件应用教科组联系方式：zhjsui@ 64252169实验16楼605室

化学工程师的任务工厂运行系统特性过程设计数学模型反应特征过程特性设备特性工艺开发数学模型

化学工程专业数学模型类型非线性方程(组)1 f= 4.0 0.4 lg(Re gen f 1 n ' 2 ) 1.2 n' 0.75 n'

常微分方程(组)

f" '+ ff"+β (1 ( f ' ) 2 )= 0偏微分方程(组) u1 2 u1= 0.024 2 F (u1 u 2 ) t x

非线性模型，难以获得解析解，必须采用数值解法模型的数值解法

实验6 Matlab数值计算

实验6 Matlab数值计算 disp('随机矩阵为：'); 实验目的：

1、 掌握数据统计与分析的方法； 2、 掌握数据插值和曲线拟合的方法及其应用；

3、 掌握多项式的常用运算。 实验内容：

1.

利用randn函数生成符合正态分布的

10×5随机矩阵A，进行如下操作：

（1） 求A的最大元素和最小元素；

（2） 求A的每行元素的和以及全部元素的和；

（3） 分别对A的每列元素按升序、每行元素按降序排列。

a=randn(10,5); am=max(max(a)); ai=min(min(a)); ah=sum(a,2); az=sum(ah); al=sort(a); ahj=-sort(-a,2);

a

disp('最大值：'); am

disp('最小值：'); ai

disp('每行和：'); ah

disp('总和：'); az

disp('列升序：'); al

disp('行降序：');

ahj

随机矩阵为： a =

-0.4326 -0.1867 6041

-1.6656 0.7258 2573

0.1253 -0.5883 0565

0.2877 2.1832 151

-1.1465

MATLAB语言 simulink

第2章 MATLAB数值计算2.1 变量与数据 2.2 矩阵与数组2.3 矩阵与数组运算 2.4 多项式运算 2.5 字符运算1

MATLAB语言 simulink

第2章 MATLAB数值计算2.1 变量与数据2.1.1 数据数据的表达方式：采用十进制表示 矩阵和数组的概念 ：标量：是指1×1的矩阵，即为只含1个数的矩阵。 向量：是指1×n或n×1的矩阵，即为只含1行或1列的矩阵 矩阵：是1个矩形的数组，即二维数组，其中向量和标量都是矩阵 的特例数组：是指n维数组，为矩阵的延伸，其中矩阵和向量都是 数组的特例。

复数：由实部和虚部组成，用特殊变量“i”和“j”表示虚数的单位2

MATLAB语言 simulink

第2章 MATLAB的数值计算功能2.1 变量与数据2.1.1 变量 变量的命名应遵循如下规则： 变量名必须以字母打头，之后可以是字母、数字 或下划线，如x51483,a_b_c_d_e。

变量名区分字母大小写，如Items,items,itEms 及ITEMS都是不同的变量。3

MATLAB语言 simulink

变量的命名应遵循如下规则：

变量名最多可包含63个字符(对于

数值计算方法配套答案

第一章 绪论

一 本章的学习要求

（1）会求有效数字。

（2）会求函数的误差及误差限。 （3）能根据要求进行误差分析。

二 本章应掌握的重点公式

（1）绝对误差：设x为精确值，x为x的一个近似值，称e??x??x为x的绝对误

差。

??e?（2）相对误差：er??。

x?（3）绝对误差限：??e?x?x。 （4）相对误差限：?r???????x??x??xx?。

??df??（5）一元函数的绝对误差限：设一元函数f?x??0,则??f??? ????x?。?dx??df??（6）一元函数的相对误差限：?r?f???1?????x?。 ??f?dx????f???（7）二元函数的绝对误差限：设一元函数f?x,y??0,则??f??? ??y。?????y????f??1?（8）二元函数的相对误差限：?r?f????????x????f???x??????f???????y??。

???y???

第 - 1 - 页

数值计算方法配套答案

三 本章习题解析

1. 下列各数都是经过四舍五入得到的近似值，（1）试指出它们有几位有效数字，（2）分别

X2?估计A1?XX2X及A2?的相对误差限。

X4??1??3x1??1.

实验二 MATLAB数值计算：二阶电路的时域分析

一、实验目的

在物理学和工程技术上，很多问题都可以用一个或一组常微分方程来描述，因此要解决相应的实际问题往往需要首先求解对应的微分方程(组)。在大多数情况下这些微分方程(组)通常是非线性的或者是超越方程(比如范德堡方程，波导本征值方程等)，很难解析地求解（精确解），因此往往需要使用计算机数值求解（近似解）。MATLAB作为一种强大的科学计算语言，其在数值计算和数据的可视化方面具有无以伦比的优势。在解决常微分方程(组)问题上，MATLAB就提供了多种可适用于不同场合(如刚性和非刚性问题)下的求解器(Solver)，例如ode45，ode15s，ode23，ode23s等等。本次实验将以二阶线性电路-RLC电路和二阶非线性电路-范德堡电路的时域计算为例，了解和学习使用MATLAB作为计算工具来解算复杂的微分方程，以期达到如下几个目的：

1. 熟练使用dsolve函数解析求解常微分方程； 2. 熟练运用ode45求解器数值求解常微分方程；

3. 了解状态方程的概念，能使用MATLAB对二阶电路进行计算和分析；

二、实验预备知识

1．微分方程的概念

未知的函数以及它的某些阶的导数连同自变量都由一已

Matlab程序设计实验报告

班级：08通信

姓名：翁惠柳 学号：222008315222023

实验 一 ： matlab的基本操作

一、实验目的： 1， 2， 3， 4，

熟悉掌握基本矩阵的建立。

掌握一些简单操作、逻辑操作以及关系运算。

熟悉MATLAB的矩阵、数组运算及其函数的使用方法； 掌握定位矩阵元素的方法。

二、实验内容：

1、 利用基本矩阵产生3*3和15*8的单位阵、全1阵、全0阵、均匀分布的随机阵（[-1,1]之间）、正态分布随机阵（均值为1，方差为4）。 2、 利用diag等函数产生下列矩阵： ?0?a?0??2?0?738??2??5 b?0????70??0504??0 ??8?然后利用reshape函数将它们变成行向量。

3、 产生一均匀分布的（-5,5）之间的随机阵（50*2），要求精确到小数点后一位。

4、 编程是实现当??[-?, ?],间隔为1?时求解正弦和余弦的值。

5、 利用rand函数产生（0,1）间均匀分布的10*10的随机矩阵A，然后统计A中大于等于

0.6的元素个数。 6、 利用randn函数产生均值为0，方差为1的10*10正态分布随机矩阵A，然后统

Matlab程序设计实验报告

班级：08通信

姓名：翁惠柳 学号：222008315222023

实验 一 ： matlab的基本操作

一、实验目的： 1， 2， 3， 4，

熟悉掌握基本矩阵的建立。

掌握一些简单操作、逻辑操作以及关系运算。

熟悉MATLAB的矩阵、数组运算及其函数的使用方法； 掌握定位矩阵元素的方法。

二、实验内容：

1、 利用基本矩阵产生3*3和15*8的单位阵、全1阵、全0阵、均匀分布的随机阵（[-1,1]之间）、正态分布随机阵（均值为1，方差为4）。 2、 利用diag等函数产生下列矩阵： ?0?a?0??2?0?738??2??5 b?0????70??0504??0 ??8?然后利用reshape函数将它们变成行向量。

3、 产生一均匀分布的（-5,5）之间的随机阵（50*2），要求精确到小数点后一位。

4、 编程是实现当??[-?, ?],间隔为1?时求解正弦和余弦的值。

5、 利用rand函数产生（0,1）间均匀分布的10*10的随机矩阵A，然后统计A中大于等于

0.6的元素个数。 6、 利用randn函数产生均值为0，方差为1的10*10正态分布随机矩阵A，然后统

MatLab & 数学建模

第五讲 数值计算(二)

一、 线性优化

minCTnx,x?Rs.t.Ax?b vlb?x?vub用命令x=lp(C,A,b,vlb,vub)。 [例] 最小值线性优化 f(x)=-5x1-4x2-6x3 x1-x2+x3≦20 3x1+2x2+4x3≦42 3x1+2x2≦30

(0≦x1, 0≦x2,0≦x3)

First, enter the coefficients: f = [-5; -4; -6] A = [1 -1 1 3 2 4 3 2 0]; b = [20; 42; 30];

lb = [0,0,0]; % x的最小值 [0,0,0] ub = [inf,inf,inf];

Next, call a linear programming routine: x= lp(f,A,b,lb,ub);

Entering x x =

0.0000 15.0000 3.0000

1

实际此命令改为：

x = linprog(f,A,b,Aeq,beq)

x = linprog(f,