A a abs 绝对值, 模
acos 反余弦
acosh 反双曲余弦
acot 反余切
acoth 反双曲余切
acsc 反余割
acsch 反双曲余割
all 所有元素均非零则为真alpha 透明控制angle 相角
ans 最新表达式的运算结果any 有非零元则为真area 面域图
asec 反正割
asech 反双曲正割
asin 反正弦
asinh 反双曲正弦
atan 反正切
atan2 四象限反正切
atanh 反双曲正切
autumn 红、黄浓淡色
axis 轴的刻度和表现
B b
bar 直方图
binocdf 二项分布概率
binopdf 二项分布累积概率binornd 产生二项分布随机数组blanks 空格符号
bode 给出系统的对数频率曲线bone 蓝色调浓淡色阵
box 坐标封闭开关
break 终止最内循环
brighten 控制色彩的明暗butter ButterWorth低通滤波器C c
caxis (伪)颜色轴刻度
cd 设置当前工作目录
cdf2rdf 复数对角型转换到实块对角型
ceil 朝正无穷大方向取整cell 创建单元数组char 创建字符串数组或者将其
他类型
变量转化为字符串数组
charfcn Maple函数
Children 图形对象的子对象
clabel 等高线标注
class 判别数据类别
clc 清除指令窗中显示内容
clear 从内存中清除变量和函数
clf 清除当前图形窗图形
close 关闭图形窗
collect 合并同类项
Color 图形对象色彩属性
colorbar 显示色条
colorcube 三浓淡多彩交错色
colordef 定义图形窗色彩
colormap 设置色图
comet 彗星状轨迹图
comet3 三维彗星动态轨迹线图compass 射线图;主用于方向和速度
cond 矩阵条件数
conj 复数共轭
continue 将控制转交给外层的for或while循环
contour 等高线图
contourf 填色等高线图
conv 卷积和多项式相乘
cool 青和品红浓淡色图copper 线性变化纯铜色调图corrcoef 相关系数
cos 余弦
cosh 双曲余弦
cot 余切
coth 双曲余切
cov 协方差矩阵
csc 余割
csch 双曲余割
cumsum 元素累计和cumtrapz 梯形法累计积分
D d
dblquad 二重(闭型)数值积分指令
deconv 解卷和多项式相除del2 计算曲率
demos 演示函数
det 行列式的值
diag 创建对角阵,抽取对角向量
diff 求导数,差分和近似微分
digits 控制符号数值的有效数字位数
dir 列出目录清单
dirac 单位冲激函数
disp 显示矩阵和文字内容disttool 概率分布计算交互界面doc 列出指定工具包中所有函
数名
docsearch 进行多词条检索double 把符号常数转化为16 位相对精度的浮点数值对象drawnow 刷新屏幕
dsolve 求解符号常微分方程
E e
edit 矩阵编辑器,打开M文件
Ei maple 指数积分
eig 矩阵特征值和特征向量
end 数组的最大下标,结束for,while,if 语句
eps 浮点相对误差EraseMode 图形对象属性
error 显示错误信息
exit 关闭MATLAB
exp 指数
expand 对指定项展开
expm 矩阵指数
eye 单位阵
ezcontour 画等位线ezcontourf 画填色等位线ezmesh 画网线图
ezmeshc 画带等位线的网线图ezplot 绘制符号表达式的二维图形
ezplot3 画三维曲线
ezpolar 画极坐标曲线
ezsurf 画曲面图
ezsurfc 画带等位线的曲面图
F f
factor 进行因式或因子分解false 按指定大小创建全0逻辑数组feather 从X轴出发的复数向量图,羽毛图
feval 函数宏指令
fill 多边形填色图
find 寻找非数单下标标识findsym 确认表达式中自由符号变量
fix 朝零方向取整
flag 红-白-蓝-黑交错色图
fliplr 矩阵的左右翻转
flipud 矩阵的上下翻转
floor 朝负无穷大方向取整fminbnd 非线性函数在某区间中极小值fminsearch 单纯形法求多元函数极值点指令
for (end) 按规定次数重复执行语句
format 设置数据输出格式fourier Fourier变换
fsolve 解非线性方程组的最简单格式
function 函数文件头functions 观察函数句柄内涵function handle 函数句柄funfun 数值泛函函数和ODE解算器
funm 计算一般矩阵函数
fzero 单变量函数的零点
G g
gallery 产生测试矩阵
gca 获得当前轴的柄
gcf 获得当前图的柄
general 通用指令
get 获得图柄
getframe 获得影片动画图象的帧
ginput 用鼠标在图上获取数据
global 定义全局变量
gradient 梯度
gray 线性灰度
grid on/off 画坐标网格线
H h
heaviside 单位阶跃函数
help 在线帮助指令helpbrowser 超文本文档帮助信息
helpdesk 超文本文档帮助信息helpwin 打开在线帮助窗hidden 网线图消隐开关
hist 统计频数直方图
histfit 带拟曲线的统计频数直方图
hold on/off 图形的保持horner 转换成嵌套形式
hot 黑-红-黄-白交错色图
hsv 饱和色彩图
I i
i, j 虚数单位
if end 条件执行语句
if-else-end 程序分支控制ifourier Fourier反变换ilaplace Laplace反变换imag 复数虚部
image 图像
impulse 给出系统的冲激响应ind2sub 据单下标换算出全下标
inf 或Inf 无穷大
inline 创建内联函数
input 提示键盘输入
int 计算积分
int2str 整数转换为字符串
inv 矩阵的逆
invhilb Hilbert矩阵法求逆阵isa 判断指定变量类别
ischar 若是字符串则为真isempty 若是空矩阵则为真isfinite 若是有限数则为真isglobal 若是全局变量则为真ishandle 是否图柄
isinf 若是无穷大则为真isletter 串中是字母则为真islogical 若是逻辑数则为真isnan 若为非数则为真isnumeric 若是数值则为真isolate maple的特殊指令isprime 是否质数
isreal 若是实数矩阵则为真isspace 串中是空格则为真iztrans Z反变换
J j
jacobian Jacobian 矩阵
jet 变异HSV色图
jordan Jordan分解
K k
keyboard 键盘获得控制权
L l
laplace Laplace变换
legend 形成图例说明
length 确定数组长度
light 灯光控制
lighting 设置照明模式
limit 求极限
line 创建线对象
LineStyle 图形线对象属性-线型
LineWidth 图形属性-线宽linmod2 从SIMULINK模型得到系统的状态方程
linspace 线性等分向量
load 从磁盘调入数据变量
Location 图形对象属性-位置
log 自然对数
log10 常用对数
log2 以2为底的对数
logical 将数值转化为逻辑值
logspace 对数等分向量
lookfor 关键词检索
M m
magic 魔方阵
maple 进入MAPLE工作空间计算
Marker 图形对象属性-点形状MarkerEdge-
Color 图形对象属性-点边界色彩
MarkerFace-
Color 图形对象属性-点域色彩MaekerSize 图形对象属性-点大小material 对象材质
max 最大值
md 创建目录
mean 平均值
mesh 三维网线图
meshgrid 用于三维曲面的分格线坐标
mfun 对MAPLE中若干经典特殊函数的数值计算
mfunlist MAPLE经典特殊函数列表
mhelp 查阅Maple中的库函数及其调用方法
min 最小值
minreal 状态方程最小实现mkdir 创建目录
mod 模数求余
more 命令窗口分页输出的控制开关
movie 播放影片动画
moviein 影片动画内存初始化mtaylor Taylor级数展开
N n
NaN或nan 非数
nargin 函数输入量的个数nargout 函数输出量的个数ndims 数组的维数
norm 矩阵或向量范数normcdf 正态分布累计概率normpdf 服从N分布的随机变量取值x的概率密度
normrnd 产生服从N分布的随
机数组
notebook 创建或打开M-book
文件
null 零空间
num2str 把数值转换为字符串
numden 提取公因式
O o
ode45 高阶法解微分方程
ones 全1 数组
optimset
orth 值空间
P p
pack 合并工作内存中的碎块
pascal Pascal 矩阵
path 控制MATLAB的搜索路径pathtool 修改搜索路径
pause 暂停
pcolor 用颜色反映数据的伪色图
peaks 产生peaks图形数据
pi 3.14….
pie 饼形统计图
pink 淡粉红色图阵
plot 直角坐标下线性刻度曲线plot3 三维直角坐标曲线图plotyy 双纵坐标图
polar 极坐标曲线图
poly 特征多项式,由根创建多项式
poly2sym 将多项式转换为符号多项式
polyfit 多项式拟合
polyval 求多项式的值polyvalm 求矩阵多项式的值pow2 2的幂
pretty 习惯方式显示
prism 光谱色图阵
prod 元素积
Q q
quad 低阶法数值积分
quadl 高阶法数值积分
quit 退出MATLAB
quiver 二维箭头图;主用于场强、流向
R r
rand 均匀分布随机数组randn 正态分布随机数组random 产生各种分布随机数
组
randsrc 产生均布数组
rank 秩
real 复数实部
realmax 最大浮点数
realmin 最小正浮点数
rem 求余数
repmat 铺放模块数组reshape 矩阵变维
residue 求部分分式表达return 返回
roots 求多项式的根
rose 频数扇形图;主用于统计rot90 矩阵逆时针旋转90度rotate 旋转指令
round 四舍五入取整
rref 转换为行阶梯形
S s
save 把内存变量存入磁盘
sec 正割
sech 双曲正割
set 设置图形对象属性shading 图形渲染模式
shg 显示图形窗
sign 函数符号,符号函数simple 运用各种指令化简符号表达式
simplify 恒等式简化simulink 打开SIMULINK集成环境
sin 正弦
sinh 双曲正弦
size 确定数组大小
slice 切片图
实验五数据处理、多项式计算 数值微积分与方程数值求解 实验要求： 为达到理想的实验效果，同学们务必做到： (1)实验前认真准备，要根据实验目的和实验内容，复习好实验中可能要用 到的命令，想好编程的思路，做到胸有成竹，提高上机效率。 (2)实验过程中积极思考，要深入分析命令、程序的执行结果以及各种屏幕 信息的含义、出现的原因并提出解决办法。 (3)实验后认真总结，要总结本次实验有哪些收获，还存在哪些问题，并写 出实验报告。实验报告应包括实验目的、实验内容、流程图（较大程序）、程序（命令）清单、运行结果以及实验的收获与体会等内容。 同学们在上机过程中会碰到各种各样的问题，分析问题和解决问题的过程就是积累经验的过程。只要同学们按照上面3点要求去做，在学完本课程后就一定会有很大的收获。 一、实验目的 1. 掌握数据统计和分析的方法。 2. 掌握数值插值与曲线拟合的方法及其应用。 3. 掌握多项式的常用运算。 二、实验内容 1. 利用MATLAB提供的rand函数生成30000个符合均匀分布的随机数，然后检验随机数的性质： (1) 均值和标准方差。 (2) 最大元素和最小元素。 (3) 大于的随机数个数占总数的百分比。 2. 将100个学生5门功课的成绩存入矩阵P中，进行如下处理：
(1) 分别求每门课的最高分、最低分及相应学生序号。 (2) 分别求每门课的平均分和标准方差。 (3) 5门课总分的最高分、最低分及相应学生序号。 (4) 将5门课总分按从大到小顺序存入zcj 中，相应学生序号存入xsxh 。 提示：上机调试时，为避免输入学生成绩的麻烦，可用取值范围在[45,95]之间的随机矩阵来表示学生成绩。 3. 某气象观测得某日6:00~18:00之间每隔2h 的室内外温度（0C ）如实验表1所示。 实验表1 室内外温度观测结果（0C ） 时间h 6 室内温度t1 室外温度t2 试用三次样条插值分别求出该日室内外6:30~18:30之间每隔2h 各点的近似温度（0C ）。 4. 已知lgx 在[1,101]区间10个整数采样点的函数值如实验表2所示。 实验表2 lgx 在10个采样点的函数值 x 1 11 21 31 41 51 61 71 81 91 101 lgx 0 试求lgx 的5次拟合多项式p(x)，并绘制出lgx 和p(x)在[1,101]区间的函数曲线。 5. 有3个多项式P1(x)=x4+2x3+4x2+5，P2(x)=x+2，P3(x)=x2+2x+3，试进行下列操作： (1) 求P(x)=P 1(x)+P 2(x)P 3(x)。 (2) 求P(x)的根。 (3) 当x 取矩阵A 的每一元素时，求P(x)的值。其中 ： 1 1. 2 1.40.752 3.505 2.5A --????=?????? (4) 当以矩阵A 为自变量时，求P(x)的值。其中A 的值与第(3)题相同。 1. 利用MATLAB 提供的rand 函数生成30000个符合均匀分布的随机数，然后检验随机
实验四 数值微积分与方程数值求解 以及符号运算 一、实验目的 1. 掌握求数值导数和数值积分的方法。 2. 掌握代数方程数值求解的方法。 3. 掌握常微分方程数值求解的方法。 4. 掌握定义符号对象的方法。 5. 掌握符号表达式的运算法则以及符号矩阵运算。 6. 掌握求符号函数极限及导数的方法。 7. 掌握求符号函数定积分和不定积分的方法。 8. 掌握级数求和的方法。 二、实验内容 1. 求函数在指定点的数值导数。 23 2()123,1,2,302 6x x x f x x x x x == 2. 用数值方法求定积分。 (1) 210I π =?的近似值。 3. 分别用3种不同的数值方法解线性方程组（左除、求逆、LU 分解）。 65254941334221 39211 x y z u x y z u x y z u x y u +-+=-??-+-=??++-=??-+=? 4. 求代数方程的数值解。 (1) 3x +sin x -e x =0在x 0=1.5附近的根。 5. 求函数在指定区间的极值。 (1) 3cos log ()x x x x x f x e ++=在(0,1)内的最小值。 (2) 332 12112122(,)2410f x x x x x x x x =+-+在[0,0]附近的最小值点和最小值。 6. 已知x=6,y=5，利用符号表达式求 z =
提示：定义符号常数x=sym(‘6’)，y=sym(‘5’)。 7. 分解因式。 (1) x 4-y 4 (2) 5135 8. 化简表达式。 21212 483(1)sin cos cos sin (2)21x x x ββββ++-+ 9. 用符号方法求下列积分。 48 2 ln 22400(1)(2)11(3)(4)(1)1 x x dx x x x dx e e dx x +∞+++++???10. 级数符号求和。 (1) 计算101 121n S n ==-∑。 (2) 求级数211n n n x ∞-=∑的和函数，并求215n n n ∞ =∑之和。 11. 求微分方程组的通解。 233453442dx x y z dt dy x y z dt dz x y z dt ?=-+???=-+???=-+?? 实验四 数值微积分与方程数值求解 以及符号运算 一、实验目的 1. 掌握求数值导数和数值积分的方法。 2. 掌握代数方程数值求解的方法。
示波器的使用和数据保存 1.示波器的参数 " Number of axes" 项用于设定示波器的Y 轴数量，即示波器的输入信号端口的个数，其预设值为"1" ，也就是说该示波器可以用来观 察一路信号，将其设为"2" ，则可以同时观察两路信号，并且示波器的图标也自动变为有两个输入端口，依次类推，这样一个示波器可以同时观察多路信号。 "Time range" (时间范围) ，用于设定示波器时间轴的最大值，这一般可以选自动(auto) ，这样X 轴就自动以系统仿真参数设置中的起始和终止时间作为示披器的时间显示范围。 第三项用于选择标签的贴放位置。 第四项用于选择数据取样方式，其中Decimation 方式是当右边栏设为"3" 时，则每3 个数据取一个，设为"5" 时，则是5 中取1 ，设的数字越大显示的波形就越粗糙，但是数据存储的空间可以减少。一般该项保持预置值"1" ，这样输入的数据都显示，画出的波形较光滑漂亮。如果取样方式选Sample time 采样方式，则其右栏里输入的是采样的时间间隔，这时将
按采样间隔提取数据显示。该页中还有一项"Floating scope" 选择，如果在它左方的小框中点击选中，则该示波器成为浮动的示波器，即没有输入接口，但可以接收其他模块发送来的数据。 示波器设置的第二页是数据页，这里有两项选择。第一项是数据点数，预置值是5000 ，即可以显示5000个数据，若超过5000 个数据，则删掉前面的保留后面的。也可以不选该项，这样所有数据都显示，在计算量大时对内存的要求高一些。如果选中了数据页的第二项"Save data to workspace" ，即将数据放到工作间去，则仿真的结果可以保存起来，并可以用MATLAB 的绘图命令来处理，也可以用其他绘图软件画出更漂亮的图形。 在保存数据栏下，还有两项设置，第一项是保存的数据命名(Variable name) ，这时给数据起一个名，以便将来调用时识别。第二项是选择数据的保存格式(Format) ，该处有3 种选择:Arrary格式适用于只有一个输入变量的情况;Structure with time 和Structure 这两种格式适用于以矢量表示的多个变量情况，并且前者同时保存数值和时间，后者仅保存数值。用Arrary 格式保存的变量，为了以后可以用
6 用抛物线法求飞f(x)=0.25x^4-4/3*x^3+5/2*x^2-2*x的极小值点解： >> f=inline('0.25*x^4-4/3*x^3+5/2*x^2-2*x','x'); x0=0;x1=-6;x2=6; c=0.001; f0=f(x0); f1=f(x1); f2=f(x2); b1=(f1-f0)/(x1-x0); b2=(f2-f1)/(x2-x1); a2=(b2-b1)/(x2-x0); a1=b1-(x1+x0)*a2; x3=0.5*(x1+x0-b1/a2); f3=f(x3); while abs(x2-x0)>c f0=f(x0); f1=f(x1); f2=f(x2); b1=(f1-f0)/(x1-x0); b2=(f2-f1)/(x2-x1); a2=(b2-b1)/(x2-x0); a1=b1-(x1+x0)*a2; x3=0.5*(x1+x0-b1/a2); f3=f(x3); if (f0-f3)<0 if (x0-x3)<0 x2=x3;f2=f3; else x1=x3;f1=x3; end else if (x0-x3)<0 x1=x0;f1=f0; x0=x3;f0=f3; else x2=x0;f2=f0; x0=x3;f0=f3; end end end x3=x0;f3=f0; x3 = 2.0000 f3 = -0.6667 所以：极小值点为x=2.0000 极小值为：f3=-0.6667
9 编写用”成功——失败“法求飞f(x)= x^5+2*x^4-4*x^3+x^2+x+2的极小值点的计算程序。并求解。 解： 编写程序并运行如下： >> f=inline('x^5+2*x^4-4*x^3+x^2+x+2','x') h=0.01; c=0.0001 x0=0 f1=f(x0); while 1 f2=f(x0+h); if f2 地贫患者的基因筛查问题 摘要 地中海贫血（简称“地贫”）是全球广为流行、危害极为严重的遗传性溶血性疾病，全世界至少有亿人携带地中海贫血的致病基因。医学上通过大人群的基因筛查来预防地贫患儿的出生。 本文应用统计学原理，对病人以及健康人的110个基因进行分析，采用Fisher判别模型建立判别标准和多元统计模型spss 软件进行筛选。 问题一，利用费希尔模型判别待测者是否患有地贫，以编号1~20地贫患者的样本，编号21~40健康人员的样本，分别作为模版建立模型，用mathlab软件求解得到待测组的患病者编号41~60个是待筛查人员的样本。 问题二，为确定“地贫”样本与“健康”样本在基因链上的区别。以及癌症样本中是否有子类。我们用1~20数据为标准化并确立相关系数矩阵，求出相关矩阵的特征值和特征向量，然后通过前m 个 主成分的累计贡献率满足 % 85 ) 1 /( ) 1 (≥ ∑ = ∑ =k i k k i k λ λ 来确定贡献率矩阵，从而得出各种基因的权 值，又利用初始特征值需大于 1，再运用逐步剔除法得出关键基因关键字：地贫患者的基因 Fisher判别筛查相关系数矩阵 1 问题重述 化验指标能够协助医生诊断。人们到医院就诊时，诊断就诊人员是否患肾炎时通常要化验人体内各种元素含量。表是确诊病例的化验结果，其中1－30号病例是已经确诊为肾炎病人的化验结果；31－60号病例是已经确定为健康人的结果。表是就诊人员的化验结果。 1.根据表中的数据，提出一种或多种简便的判别方法，判别属于患者或健康人的方 法，并检验你提出方法的正确性。 2.按照1提出的方法，判断表中的30名就诊人员的化验结果进行判别，判定他（她） 们是肾炎病人还是健康人。 3.能否根据表的数据特征，确定哪些指标是影响人们患肾炎的关键或主要因素，以 便减少化验的指标。 4.根据3的结果，重复2的工作。 5.对2和4的结果作进一步的分析。 2 问题分析 问题解决的关键是如何正确判断正常人与患者之间的差异，利用所给数据，可以选择用医学统计方法[1]中的判别分析法[2]进行分析。从题目给出的表中可以得出以下信息：1）表中分别给出正常人与患者各30组数据，每组数据各包含7种元素（Zn、Cu、Fe、Ca、Mg、K、Na）在人体中的含量。通过对这些数据进行分析，可以从中找出数据差异，根据判别法确定判别标准。利用所得判别标准，与就诊人员的化验结果比较可以判 论文关键词: 图像分割边缘检测模糊理论遗传算法 Matlab 论文摘要：分割的目的是将图像划分为不同区域。图像分割算法一般是基于亮度值的两个基本特性之一：不连续性和相似性。第一类性质的已用途径是基于亮度的不连续变化分割图像，比如图像的边缘。第二类的主要应用途径是依据事先制订的准则将图像分割为相似的区域。门限处理、区域生长、区域分离和聚合都是这类方法的实例。遗传算法具有简单、鲁棒性好和本质并行的突出优点。其在应用领域取得的巨大成功，引起了广大学者的关注。在图像分割领域，遗传算法常用来帮助确定分割阈值。 本文介绍讨论了几种目前广泛应用的图像边缘检测、图像阈值分割的各种算法，并给出了对比分析;对遗传算法的基本概念和研究进展进行了综述;给出了标准遗传算法的原理、过程、实验结果及分析. 实验结果表明，本文提出的遗传分割算法优于传统分割算法。 第一章绪论 1.1 图像分割综述 图像分割就是指把图像分成各具特性的区域并提取出感兴趣目标的技术和过程。这里所说的特性可以是灰度、颜色、纹理等，而目标可以对应单个区域，也可以对应多个区域。图像分割是数字图像处理中的一项关键技术，它使得其后的图像分析，识别等高级处理阶段所要处理的数据量大大减少，同时又保留有关图像结构特征的信息。而且，在数字图像处理工程中，一方面，图像分割是目标表达的基础，对特征测量有重要的影响；另一方面，图像分割是自动目标识别的关键步骤，图像分割及其基于分割的目标表达、特征提取和参数测量等将原始图像转化为更抽象更紧凑的形式，分割中出现的误差会传播至高层次处理阶段，因此分割的精确程度是至关重要的。只有通过细致精细的图像分割，才能使得更高层的图像分析和理解成为可能。因此，图像分割是由图像处理进到图像分析的关键步骤，在图像工程中占据重要的位置。 1.2 图像分割的研究意义与发展现状 作为计算机视觉和图像处理中的难点和热点之一，图像分割的研究受到了研究工作者的高度重视，对图像分割进行了深入、广泛的研究。作为一种重要的图像技术，图像分割在不同领域中有时也用其它名称：如目标轮廓(object delineation)技术，阈值化(thresholding)技术，图像区分或求差(image discrimination)技术，目标检测(target detection)技术，目标识别(target recognition)技术，目标跟踪(target tracking)技术等，但这些技术本身或其核心实际上也就是图像分割技术。图像分割作为图像处理、分析的一项基本内容，其应用非常广泛，几乎出现在有关图像处理的所有领域，并涉及各种类型的图像。在工业自动化、在线产品检验、生产程控、文件图像处理、遥感图像、保安监视、以及军事、体育、农业等行业和工程中，图像分割都有着广泛的应用。例如：在遥感图像中，合成孔径雷达图像中目标的分割、遥感云图中不同云系和背景分布的分割等；在医学应用中，脑部 MR 图像分割成灰质（GM）、白质（WM）、脑脊髓（CSF）等脑组织和其它脑组织区域（NB）等；在交通图像分析中，把车辆目标从背景中分割出来等；在面向对象的图像压缩和基于内容的图像检索中将图像分割成不同的对象区域等。在各种图像应用中，只要需对图像目标进行提取，测量等都离不开图像分割。 自 20 世纪 70 年代至今，已提出上千种各种类型的分割算法。如：门限法、匹配法、区域生长法、分裂－合并法、水线法、马尔可夫随机场模型法、多尺度法、小波分析法、数学形态 附：计算使用的MATLAB程序代码及其含义 %不同i0对应车速为40、50、60时的发动机转数 clear clc close all r=0.4275; %车轮半径 u0=[40 50 60]; %速度 i0=[4.22 4.68 4.88 5.01 5.28 5.48]; %主减速比 n=(i0'*u0)/(0.377*r); %对应发动机转数 n %加速燃油消耗率 clear clc close all yita=0.89;G=9160*9.8;f=0.02;CDA=4.64;a=0.1704;Iw1=1.798;Iw2=3.598;r=0.4 275; m=9160;If=0.218;ig=1;i0=[4.22,4.68,4.88,5.01,5.28,5.48]; for x=1:1:6 n=1000:100:2800; %等距取数 ua5=(0.377*r.*n)/i0(x); %求转数对应车速 F5=f*G+CDA*(ua5.^2)/21.15; %求阻力 P_fw5=F5.*ua5./(yita*3.6*1000); %阻力功率 ua0=50:1:60; %50到60公里加速 ua1=50:1:60; delta=1+(Iw1+Iw2)/(m*r^2)+(If*ig^2*i0(x)^2*yita)/(m*r^2); %求旋转质量换算系数 P0=(G*f.*ua0./3600+CDA.*ua0.^3/76140+(delta*m.*ua0/3600)*a)/yita; %初速度为ua0时的阻力功率 P=(G*f.*ua1/3600+CDA.*ua1.^3/76140+(delta*m.*ua1/3600)*a)/yita; %末速度为ua1时的阻力功率 b1=0.0133.*P.^2-2.0715.*P+287.47; %燃油消耗率 Qt=P.*b1./(367.1.*8); %燃油消耗 i1=size(Qt); %求Qt向量长度 i=i1(2); Qt1=Qt(2:i-1); dt=1/(3.6*a); %速度每增加1km/h所需要的时间 q=(Qt(1)+Qt(i))*dt./2+sum(Qt1)*dt; %加速段燃油消耗量 Page 1/9 0020002100220023002400250026002700280029002A002B002C002D002E002F 0030003100320033003400350036003700380039003A003B003C003D003E003F 0040004100420043004400450046004700480049004A004B004C004D004E004F 0050005100520053005400550056005700580059005A005B005C005D005E005F 0060006100620063006400650066006700680069006A006B006C006D006E006F 0070007100720073007400750076007700780079007A007B007C007D007E007F 00A000A100A200A300A400A500A600A700A800A900AA00AB00AC00AD00AE00AF 00B000B100B200B300B400B500B600B700B800B900BA00BB00BC00BD00BE00BF 00C000C100C200C300C400C500C600C700C800C900CA00CB00CC00CD00CE00CF 00D000D100D200D300D400D500D600D700D800D900DA00DB00DC00DD00DE00DF 00E000E100E200E300E400E500E600E700E800E900EA00EB00EC00ED00EE00EF 00F000F100F200F300F400F500F600F700F800F900FA00FB00FC00FD00FE00FF 010901310141014201520153016001610178017D017E0192019F01A901C001C1 01C3030203030304030503060307030A030B030C0311032303240328032C032D 032E032F033003310332033303340335033603370338033F0347035C035D035E 线性规划及其LINDO 求解 姓名 任淑丽 学号 102093116 班级 10计算科学 说明： （1）要求在Notebook 环境下完成； （2）完成后请以姓名(班级学号)实验报告名称命名并存盘; （3）在下一次上课之前由学习委员收齐打包以附件形式后发到 lvximing@https://www.doczj.com/doc/fa7358908.html, ； （5）请自行保存备份,以备后用。 【LINDO 程序】 max 50x1+100x2 ST x1+x2<=300 2x1+x2<=400 x2<=250 end 【计算结果】 OBJECTIVE FUNCTION VALUE 1) 27500.00 1212122max 501003002400.250 0(1,2)LINDO j z x x x x x x st x x j =++≤?? +≤?? ≤??≥=? 一、已知线性规划用编程计算其最优解和最优值。 VARIABLE VALUE REDUCED COST X1 50.000000 0.000000 X2 250.000000 0.000000 ROW SLACK OR SURPLUS DUAL PRICES 2) 0.000000 50.000000 3) 50.000000 0.000000 4) 0.000000 50.000000 NO. ITERATIONS= 0 LP OPTIMUM FOUND AT STEP 0 OBJECTIVE FUNCTION VALUE 1) 27500.00 VARIABLE VALUE REDUCED COST X1 50.000000 0.000000 X2 250.000000 0.000000 ROW SLACK OR SURPLUS DUAL PRICES 2) 0.000000 50.000000 实验二 MATLAB 结构程序设计与绘图操作 实验要求： 为达到理想的实验效果，同学们务必做到： (1) 实验前认真准备，要根据实验目的和实验内容，复习好实验中可能要用到的命令，想好编程的思路，做到胸有成竹，提高上机效率。 (2) 实验过程中积极思考，要深入分析命令、程序的执行结果以及各种屏幕信息的含义、出现的原因并提出解决办法。 (3) 实验后认真总结，要总结本次实验有哪些收获，还存在哪些问题，并写出实验报告。实验报告应包括实验目的、实验内容、流程图（较大程序）、程序（命令）清单、运行结果以及实验的收获与体会等内容。 同学们在上机过程中会碰到各种各样的问题，分析问题和解决问题的过程就是积累经验的过程。只要同学们按照上面3点要求去做，在学完本课程后就一定会有很大的收获。 实验目的： 1. 掌握建立和执行M 文件的方法 2. 掌握利用if 和switch 语句实现选择结构的方法 3. 掌握利用for 语句实现循环结构的方法 4. 掌握绘制二维图形的常用函数 5. 掌握绘制三维图形的常用函数 6. 掌握绘制图形的辅助操作 实验内容： 1. 输入一个百分制成绩，要求输出成绩等级A 、B 、C 、D 、E 。其中90分~100分为A ，80分~89分为B ，70分~79分为C ，60分~69分为D ，60分以下为E 。 要求： a) 分别用if 语句和switch 语句实现。 b) 输入百分制成绩后要判断该成绩的合理性，对不合理的成绩应输出出错信息。 2. 硅谷公司员工的工资计算方法如下： c) 工作时数超过120小时者，超过部分加发15%。 d) 工作时数低于60小时者，扣发700元。 e) 其余按每小时84元计发。 试编程按输入的工号和该员工的工时数，计算应发工资。 3. 设计程序，完成两位数的加、减、乘、除四则运算，即产生两个两位随机整数，再输入一个运算符号，做相应的运算，并显示相应的结果。 4. 建立5×6矩阵，要求输出矩阵第n 行元素。当n 值超过矩阵的行数时，自动转为输出矩阵最后一行元素，并给出出错信息。 5. 根据2 2222 13121116n ++++= π，求π的近似值。当n 分别取100，1000，10000时，结果是多少？ MATLAB 作业二 1、 请将下面给出的矩阵A 和B 输入到MA TLAB 环境中，并将它们转换成符号矩阵。若某 一矩阵为数值矩阵，另以矩阵为符号矩阵，两矩阵相乘是符号矩阵还是数值矩阵。 57651653 550123231001432546 2564206441211 346,3 9636623515 212107600774101 20 172440773473 781248672171107681 5A B ??????????????????? ?==??????????? ?---????????--??? ? 一矩阵为数值矩阵，另以矩阵为符号矩阵，两矩阵相乘是符号矩阵。 2、 利用MA TLAB 语言提供的现成函数对习题1 中给出的两个矩阵进行分析，判定它们是 否为奇异矩阵，得出矩阵的秩、行列式、迹和逆矩阵，检验得出的逆矩阵是否正确。 对A: >> det(A) ans = -3.5432e+04 >> rank(A) ans =7 >> trace(A) ans = 27 对B: >> det(B) ans = -2.6326e-26 >> rank(B) ans =5 >> trace(B) ans = 26 由A和B的行列式的值可以看出，矩阵A为非奇异矩阵，而矩阵B为奇异矩阵，而由norm的结果可以说明对A直接求逆的结果比较准确，而对矩阵B直接求逆的偏差较大，所以对奇异矩阵不适合用inv函数直接求逆。 3、试求出习题1中给出的A和B矩阵的特征多项式、特征值与特征向量，并对它们进行LU 分解。 对A： 对B： MATLAB和Simulink函数Maltlab/simulink 数学函数 sin 正弦asin 反正弦 cos 余弦acos 反余弦 tan 正切atan 反正切 cot 余切acot 反余切 sec 正割asec 反正割 csc 余割acsc 反余割 exp 指数angle 相角 log 自然对数conj 复共轭 log10 常用对数image 复数虚部 sqrt 平方根real 复数实部 命令函数 axis 坐标轴比例尺控制 bode 波特图 c2d 变连续为离散系统 conv 卷积与多项式乘法 d2c 变离散为连续系统 dstep 离散时间阶跃响应 dimpulse 离散时间单位冲激响应 dinitia 离散时间零输入响应 dlsim 离散时间任意输入响应 freqs 拉普拉斯变换频域响应 freqz z变换频域响应 grid 设置网络线 hold on 保留当前图形 impulse 冲激响应 logspace 在对数空间定义角频率范围 lsim 对任意输入的线性系统仿真 margin 求相位裕量、幅值裕量和截止频率 ngrid 为nichols 图画栅格线 nyquist 奈奎斯特图 plot 绘制函数图形 rlocfind 求根轨迹增益 rlocus 根轨迹 roots 多项式的根 semilogx 半对数图 stairs 绘制离散系统时域响应的曲线命令 tf2zp 传递函数到极点零点的转换 tf 求传递函数模型 tilte 图形标题 zp2tf 零极点转换为传递函数 Simulink 函数abs 信号的绝对值 derivative 信号的微分 discrete transfer fcn 离散传递函数display 显示输入信号的值 gain 增益 Integrator 信号的积分 math function 数学函数 out 输出信号 scope 示波器 sine wave 正弦信号 step 阶跃信号 sum 信号的和 transfer fcn 实现线性传递系统transport delay 延迟给定时间 自定义示波器系列之一 自定义示波器系列之一 使用指南 力科示波器MATLAB使用指南 ——力科示波器 —— 美国力科公司 万力劢 基于X-Stream技术的力科示波器既是一个高速信号采集平台，也是一个功能强大的信号处理平台。示波器软件本身提供了丰富的测量和运算功能，能够满足常见的信号处理应用。此外，力科示波器还能借助MATLAB及其丰富的函数、工具箱做更复杂、更个性化的实时信号处理。使用MATLAB实时处理力科示波器采集的数据，可以采取以下三种方式： 遥控（（Remote Control）方式 1.遥控 MATLAB在计算机上运行，通过计算机的GPIB、以太网Ethernet或者串口RS232来控制示波器、读取采集的数据，如图1： 图 1 遥控方式 在这种联机方式下，MATLAB对力科示波器的控制方式又可分为 1）装有MATLAB的计算机作为控制器，通过SCPI程控命令（Standard Commands for Programmable Instruments，即IEEE488.2标准）来控制示波器。这和NI公司的Labview, LabWindows/CVI程控仪器的方式是类似的。 MathWorks公司提供了MATLAB工具箱Instrument Control Toolbox，可以方便地实现SCPI方式的程控。详情请参考该工具箱的帮助文件以及https://www.doczj.com/doc/fa7358908.html,/LeCroy。力科示波器支持的SCPI程控指令请参考LeCroy X-Stream Oscilloscopes Remote Control Manual。 2）力科示波器软件提供了基于COM（Component Object Model）技术的自动化接口（Automation interface）。外部程序通过这种接口可以读取示波器属性、控制示波器和采集数据。这种方式下，力科示波器作为COM Automation server, MATLAB程序作为COM Client。 这种示波器软件与MATLAB程序分别在不同平台上运行的方式，导致两者交互数据的效率受限于GPIB/ENET/RS232接口的数据吞吐率。采集数据量如果比较大，实时性就会很差。如果实际应用必须采取这种联机方式，并且想要实现较高的数据吞吐率，可以选用力科提供的专用仪控总线接口LSIB (LeCroy Serial Interface Bus)，数据传输率可达到325 MB/s。详情请参考https://www.doczj.com/doc/fa7358908.html,/files/pdf/LeCroy_LSIB_Datasheet.pdf 基于MATLAB 的线性代数问题求解 姓名 任淑丽 学号 102093116 班级_10计算科学_ 说明： （1）要求在Notebook 环境下完成； （2）完成后请以姓名(班级学号)实验报告名称命名并存盘; （3）在下一次上课之前由学习委员收齐打包以附件形式后发到 lvximing@https://www.doczj.com/doc/fa7358908.html, ； （5）请自行保存备份,以备后用。 C A B C A B C A B (1)1,2,3.=+=-=*、求 【程序代码】 A=[2,5,-1;2,1,7;1,9,6];B=[5,1,2005;2004,9,2;1,256,178]; C1=A+B C2= A-B C3=A*B C1 = 7 6 2004 2006 10 9 2 265 184 C2 = -3 4 -2006 -2002 -8 5 0 -247 -172 C3 = 10029 -209 3842 2021 1803 5258 18047 1618 3091 A B (2).、求， 【程序代码】 A=[2,5,-1;2,1,7;1,9,6];B=[5,1,2005;2004,9,2;1,256,178]; A B 2515120052172004921961256178-???? ? ?== ? ? ? ?????一、已知矩阵， a=det(A) b=det(B) a = -156 b = 1.0282e+09 C A B 1 22 (3). =? 、求 【程序代码】 A=[2,5,-1;2,1,7;1,9,6];B=[5,1,2005;2004,9,2;1,256,178]; C=A.^2.*sqrt(B) C = 1.0e+03 * 0.0089 0.0250 0.0448 0.1791 0.0030 0.0693 0.0010 1.2960 0.4803 A B AX B XA B (4),. ==、求及的逆矩阵并求解矩阵方程和 【程序代码】 A=[2,5,-1;2,1,7;1,9,6];B=[5,1,2005;2004,9,2;1,256,178]; A1=inv(A) B1=inv(B) A1 = 0.3654 0.2500 -0.2308 0.0321 -0.0833 0.1026 -0.1090 0.0833 0.0513 B1 = 0.0000 0.0005 -0.0000 -0.0003 -0.0000 0.0039 0.0005 -0.0000 -0.0000 A=[2,5,-1;2,1,7;1,9,6];B=[5,1,2005;2004,9,2;1,256,178]; format rat X1=A\B X2=B/A A a abs 绝对值, 模 acos 反余弦 acosh 反双曲余弦 acot 反余切 acoth 反双曲余切 acsc 反余割 acsch 反双曲余割 all 所有元素均非零则为真alpha 透明控制angle 相角 ans 最新表达式的运算结果any 有非零元则为真area 面域图 asec 反正割 asech 反双曲正割 asin 反正弦 asinh 反双曲正弦 atan 反正切 atan2 四象限反正切 atanh 反双曲正切 autumn 红、黄浓淡色 axis 轴的刻度和表现 B b bar 直方图 binocdf 二项分布概率 binopdf 二项分布累积概率binornd 产生二项分布随机数组blanks 空格符号 bode 给出系统的对数频率曲线bone 蓝色调浓淡色阵 box 坐标封闭开关 break 终止最内循环 brighten 控制色彩的明暗butter ButterWorth低通滤波器C c caxis (伪)颜色轴刻度 cd 设置当前工作目录 cdf2rdf 复数对角型转换到实块对角型 ceil 朝正无穷大方向取整cell 创建单元数组char 创建字符串数组或者将其 他类型 变量转化为字符串数组 charfcn Maple函数 Children 图形对象的子对象 clabel 等高线标注 class 判别数据类别 clc 清除指令窗中显示内容 clear 从内存中清除变量和函数 clf 清除当前图形窗图形 close 关闭图形窗 collect 合并同类项 Color 图形对象色彩属性 colorbar 显示色条 colorcube 三浓淡多彩交错色 colordef 定义图形窗色彩 mathlab常用的数学函数 默认分类2010-07-09 08:16:00 阅读30 评论0字号：大中小 abs(x)：纯量的绝对值或向量的长度 angle(z)：复数z的相角(Phase angle) sqrt(x)：开平方 real(z)：复数z的实部 imag(z)：复数z的虚部 conj(z)：复数z的共轭复数 round(x)：四舍五入至最近整数 fix(x)：无论正负，舍去小数至最近整数 floor(x)：地板函数，即舍去正小数至最近整数 ceil(x)：天花板函数，即加入正小数至最近整数 rat(x)：将实数x化为分数表示 rats(x)：将实数x化为多项分数展开 sign(x)：符号函数(Signum function)。 当x<0时，sign(x)=-1； 当x=0时，sign(x)=0; 当x>0时，sign(x)=1。 rem(x,y)：求x除以y的馀数 gcd(x,y)：整数x和y的最大公因数 lcm(x,y)：整数x和y的最小公倍数 exp(x)：自然指数 pow2(x)：2的指数 log(x)：以e为底的对数，即自然对数或log2(x)：以2为底的对数 log10(x)：以10为底的对数 二、MATLAB常用的三角函数 sin(x)：正弦函数 cos(x)：馀弦函数 tan(x)：正切函数 asin(x)：反正弦函数 acos(x)：反馀弦函数 atan(x)：反正切函数 atan2(x,y)：四象限的反正切函数 sinh(x)：超越正弦函数 cosh(x)：超越馀弦函数 tanh(x)：超越正切函数 asinh(x)：反超越正弦函数 acosh(x)：反超越馀弦函数数模-化验结果判别及matlab程序
matlab
mathlab汽车动力匹配程序
【VIP专享】MathcadUniMath(特殊符号输入方法)
mathlab实用题型
mathlab试题 2
mathlab作业二
MATLAB和Simulink函数Maltlabsimulink
力科示波器Matlab使用指南
mathlab的应用
mathlab常用命令
mathlab常用的数学函数
相关主题
文本预览