化工数值计算与matlab课后答案隋志军

精通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

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

第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]

实验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

第二章

%分段函数的写法 syms y x z=int(x*y,x,0,1);

g = evalin(symengine,[' piecewise([y > 1/2,' char(z) '], [y <= 1/2, 1])']) G=int(g,y) g =

piecewise([1/2 < y, y/2], [y <= 1/2, 1]) G =

piecewise([1/2 < y, y^2/4], [y <= 1/2, y]) %习题2-1

class(vpa(sym(3/7+0.1),4)) %习题2-2

findsym(sym('sin(w*t)'),1) findsym(sym('a*exp(-X)'),1) findsym(sym('z*exp(j*th)'),1) %习题2-3 syms x a positive x_1=solve('x^3-44.5') syms x unreal

x_2=solve('x^2-a*x+a^2') syms a unreal

x_2=solve('x^2-a*x+a^2') x_3=solve('x^2-a*x-a^2') syms x real

evalin(symengine,'anames(Properties)') evalin(symengine,'getprop(x)')

%习题2-4

7/3, pi/3 , pi*3^(1/3) a = pi*3^(1/3)

第二章

%分段函数的写法 syms y x z=int(x*y,x,0,1);

g = evalin(symengine,[' piecewise([y > 1/2,' char(z) '], [y <= 1/2, 1])']) G=int(g,y) g =

piecewise([1/2 < y, y/2], [y <= 1/2, 1]) G =

piecewise([1/2 < y, y^2/4], [y <= 1/2, y]) %习题2-1

class(vpa(sym(3/7+0.1),4)) %习题2-2

findsym(sym('sin(w*t)'),1) findsym(sym('a*exp(-X)'),1) findsym(sym('z*exp(j*th)'),1) %习题2-3 syms x a positive x_1=solve('x^3-44.5') syms x unreal

x_2=solve('x^2-a*x+a^2') syms a unreal

x_2=solve('x^2-a*x+a^2') x_3=solve('x^2-a*x-a^2') syms x real

evalin(symengine,'anames(Properties)') evalin(symengine,'getprop(x)')

%习题2-4

7/3, pi/3 , pi*3^(1/3) a = pi*3^(1/3)

课程名称： MATLAB实验

题 目：实验二 MATLAB数值计算与符号运算功能

学生姓名： 专 业： 电子信息工程 班 级： 学 号：

指导教师： 张 静 实验地点： 现代通信实验室 日 期： 2012 年 12月 2 日

实验2 MATLAB数值计算、符号运算功能

一、实验目的

1、掌握建立矩阵、矩阵分析与处理的方法。 2、掌握线性方程组的求解方法。

3、掌握数据统计和分析方法、多项式的常用运算。

4、掌握求数值导数和数值积分、常微分方程数值求解、非线性代数方程数值求解的方法。 5、掌握定义符号对象的方法、符号表达式的运算法则及符号矩阵运算、符号函数极限及导数、符号函数定积分和不定积分的方法。 二、预习要求

（1） 复习4、5、6章所讲内容；

（2） 熟悉MATLAB中的数值计算和符号运算的实现方法和主要函数。 三、实验内容

??29618???1、 已知A?20512，求A的特征值及特征向量，并分析其数学意义。 ?????885??程序：

求A的特征值：A=[-29

matlab 数值积分与微分

第8章 数值积分与数值微分本章目标：近似计算 I 8.1 8.2 8.3 8.4 8.5 8.6 8.7 8.8

b a

f ( x)dx 和 f ( x)

Newton-Cotes公式 复化求积公式 自适应步长求积方法 Gauss求积方法 特殊函数的积分 数值积分的MATLAB函数求解 数值微分 实例解析

matlab 数值积分与微分

8.1 Newton-Cotes公式思 利用插值多项式 P ( x ) f ( x ) 则积分易算。 n 路

插值型积分公式

在[a, b]上取 a x0 < x1 <…< xn b,做 f 的 n 次插值多项式 Ln ( x ) f ( x k )l k ( x ) ,即得到k 0 n

误差 R[ f ]b a

b

a

f ( x )dx f ( x k ) l k ( x )dxk 0

n

Ak

f ( x)dx Ak f ( xk )b a k 0

n

Ak

b a

(xj k

(x x j )k

xj )

dx

由节点 决定， b f ( n 1) ( ) n x f (x) 无关。 a ( x xk

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个字符(对于

吉林工业职业技术学院

顶岗实习报告

系别 ： 化工机械系

实践班级 ： 设备3094 班 学生姓名 ： 隋长军 （学号：34） 实习时间 ： 2011年10月21日～2011年6月30日 指导教师 ： 郭春梅 专业教研室 ： 化工设备维修技术

前 言

实习时间：2011年9月16——2012年6月 实习地点：长春轨道客车股份有限公司 实习人员：吉林工业职业技术学院部分学生 实习目的：学习铁路车辆钳工知识

本次在长春轨道客车股份有限公司，主要从事的是外窗密封岗位工作。 在上岗前公司进行了统一的安全教育培训，包括岗位安全生产注意事项和岗位防火制度，通过培训后了解了钳工安全技术所要求的安全规章制度和发生火灾事故岗位要求的责任规范。

本次顶岗实习期间主要就是对CRH3动车组前风挡玻璃密封，了解开工准备、操作要领和操作注意事项，和客车室密封工作，对司机室前侧