第二次上机实验报告(共5篇)
1、循环程序设计-1(1)汇编、连接后的截图
TODO:你的截图(必选)
TODO:你的文字解释说明(可选)
说明:masm xun 得到 xun.obj,xun.crf,xun.lst 文件,通过 link xun 得到 xun.map 文件,显示编译成功。
(2).lst 文件的截图,TODO:你的截图(必选)
TODO:你的文字解释说明(可选)
说明:通过 masm 对程序进行编译时生成.lst 文件,通过 notepad++打开.lst 文件,并进行截图(3)反汇编的截图
TODO:你的截图(必选)
TODO:你的文字解释说明(可选)
说明:在 debug 环境下执行 u 指令,显示出反汇编代码。
(4)在完成 DS 赋值后,立即显示各个寄存器的值
TODO:你的截图(必选)
TODO:你的文字解释说明(可选)
说明:按单步 t 之后,下方出现 MOV DS,AX,即下一条即将执行的指令为 MOV DS,AX,再按一次 t,此时 DS 被赋值,此时执行的 r 指令显示的就是 DS 赋值后各个寄存器的值。
(5)在进行计算前,显示数组 M 开始的 n+2 个字的内存值的截图(只能显示这n+2 个字的内存值,多显示、少显示均扣分)
TODO:你的截图(必选)
TODO:你的文字解释说明(可选)
说明:在 debug 环境下,执行 d 指令显示内存,由于不能显示其他字的内存值,所以只能
一行一行截图,可以看到,此时内存值与程序初始定义值相同。
(6)执行完计算后,立即显示各个寄存器的值
TODO:你的截图(必选)
TODO:你的文字解释说明(可选)
说明:先执行 g 指令到执行回到 dos 系统指令,此时显然已经执行完运算,此时执行 r 指令就获得执行完运算后各个寄存器的值。
(7)执行完计算后,显示数组 M 开始的 n+2 个字的内存值的截图(只能显示这n+2 个字的内存值,多显示、少显示均扣分)
TODO:你的截图(必选)
TODO:你的文字解释说明(可选)
说明:执行 d 指令显示出内存值,由于要求是不能显示其他字的值,所以只能一行一行截图,可以看到此时内存值与期望结果相同。
(8)源代码
TODO:你的源代码(必选)
;empty asm file
title I love asm
data segment n
dw 0010h
M
dw
0021h,0084h,0041h,0010h,0027h,0018h,0044h,0011H,0002h,0028h,0024h,0020h,0022h,0014h,0030h,0040h
DATA1 DW ?
ADDRR DW ?
data ends
code segment
assume cs:code, ds:data
main
proc
;assign the data segment base address to DS
mov
ax, data
mov
ds, ax
;TODO...;| add your code between arrows |
;v----------------------------v
MOV CX, n
;循环计数器
LEA DI, M
MOV AX,[DI]
;取第一个数
MOV ADDRR, DI
;记下绝对值最大的数的地址
CMP AX,0
;此数是正数吗?
JNS ZHEN
;是正数,即为绝对值,转去判断下一个数
NEG AX
;不是正数,变为其绝对值
ZHEN:MOV BX,[DI]
CMP BX,0
;此数是正数吗?
JNS COMP
;是正数,即为绝对值,转去比较绝对值大小
NEG BX
;不是正数,变为其绝对值
COMP: CMP AX, BX
;判断绝对值大小
JAE ADDRESS
MOV AX, BX
;(AX)<(BX),使(AX)中为绝对值最大的数
MOV ADDRR, DI
;记下绝对值最大的数的地址
ADDRESS:ADD DI, 2
LOOP ZHEN
MOV
DATA1, AX
;记下此数
;^----------------------------^
;|
The END
|
;method 2: return to dos
mov
ax, 4c00h
int
21h
main
endp code ends end main
2、循环程序设计-2 (1)汇编、连接后的截图
TODO:你的截图(必选)
TODO:你的文字解释说明(可选)
说明:masm xun 得到 xunhuan.obj,xunhuan.crf,xunhuan.lst 文件,通过 link xunhaun 得到xunhaun.map 文件,显示编译成功。
(2).lst 文件的截图
TODO:你的截图(必选)
TODO:你的文字解释说明(可选)
说明:通过 masm 编译 jishu,获得 jishu.lst 文件,要 notepad++打开并进行截图。
(3)反汇编的截图
TODO:你的截图(必选)
TODO:你的文字解释说明(可选)
说明:在 debug 环境下执行 u 指令显示出反汇编代码。
(4)在进行计算前,显示数组 A 中 15 个整数的内存值的截图(只能显示这 15个整数的内存值,多显示、少显示均扣分)
TODO:你的截图(必选)
TODO:你的文字解释说明(可选)
说明:执行 d 指令显示出内存值,由于不能显示其他数内存值,所以一行一行进行截图,可以看到,此时内存值与程序初始定义相同。
(5)在进行计算前,显示数组 B 中 20 个整数的内存值的截图(只能显示这 20个整数的内存值,多显示、少显示均扣分)
TODO:你的截图(必选)
TODO:你的文字解释说明(可选)
说明:执行 d 指令显示出内存值,由于不能显示其他整数内存值,所以一行一行进行截图,可以看到,此时数组 B 中的内存值与程序初始定义相同。
(6)在进行计算前,显示数组 C 中 15 个整数的内存值的截图(只能显示这 15个整数的内存值,多显示、少显示均扣分)
TODO:你的截图(必选)
TODO:你的文字解释说明(可选)
说明:执行 d 指令显示出内存值,由于不能显示其他整数内存值,所以一行一行进行截图,可以看出,此时内存值与程序初始定义相同。
(7)执行完计算后,立即显示各个寄存器的值
TODO:你的截图(必选)
TODO:你的文字解释说明(可选)
说明:先用 g 指令跳到返回 dos 调用指令,此时已经执行完计算,再执行 r 指令,显示的就是执行完运算的各个寄存器的值。
(8)执行完计算后,显示数组 A 中 15 个整数的内存值的截图(只能显示这 15个整数的内存值,多显示、少显示均扣分)
TODO:你的截图(必选)
TODO:你的文字解释说明(可选)
说明:执行 d 指令显示出内存值,由于要求是不能显示其他字的值,所以只能一行一行截图,可以看到此时内存值与期望结果相同。
(9)执行完计算后,显示数组 B 中 20 个整数的内存值的截图(只能显示这 20个整数的内存值,多显示、少显示均扣分)
TODO:你的截图(必选)
TODO:你的文字解释说明(可选)
说明:执行 d 指令显示出内存值,由于要求是不能显示其他字的值,所以只能一行一行截图,可以看到此时内存值与期望结果相同。
(10)执行完计算后,显示数组 C 中 15 个整数的内存值的截图(只能显示这 15个整数的内存值,多显示、少显示均扣分)
TODO:你的截图(必选)
TODO:你的文字解释说明(可选)
说明:执行 d 指令显示出内存值,由于要求是不能显示其他字的值,所以只能一行一行截图,可以看到此时内存值与期望结果相同。
(11)源代码
TODO:你的源代码(必选)
;empty asm file
title I love asm
data segment
A
dw
0027h,0018h,0021h,0028h,0024h,0020h,0022h,0014h,0011h,0010h,0012h,0017h,0084h,0081h,0080H
B
dw
0018h,0021h,0028h,0024h,0020h,0022h,0044H,0041h,0011h,0010h,0012h,0017h,0084h,0081h,0080h,0040h,0042h,0047h,0002h,0007h
C1
DWDUP(?)data ends
code segment
assume cs:code, ds:data
main
proc
;assign the data segment base address to DS
mov
ax, data
mov
ds, ax
;TODO...;| add your code between arrows |
;v----------------------------v
mov si,0
MOV DI,0
mov bx,A
jiao:MOV CX,B[SI]
;想法是 di 指 A 中的元素,di 递增,拿 A 中的一个元素每次与 B 中的所有元素作比较
CMP A[DI],CX
;直到遇到相同为止,如果没有遇到相同,有设置对 si 值比较指令,b 中所有元素相同时
JZ JI
;会跳转到 JI
INC SI
INC SI
CMP SI,42
JZ JI
JMP jiao
JI:mov ax,B[SI]
mov C1[DI],AX
MOV SI,0
INC DI
INC DI
CMP DI,32
JZ EXIT
JMP jiao
;^----------------------------^
;|
The END
|
;method 2: return to dos
EXIT: mov
ax, 4c00h
int
21h
main
endp code ends end main
3、分支程序设计 (1)汇编、连接后的截图
TODO:你的截图(必选)
TODO:你的文字解释说明(可选)
说明:masm xun 得到 jishu.obj,jishu.crf,jishu.lst 文件,通过 link xun 得到 jishu.map 文件,显示编译成功。
(2).lst 文件的截图
TODO:你的截图(必选)
TODO:你的文字解释说明(可选)
说明:通过 masm 形成 jishu.lst 文件,通过 notepad++打开,并进行截图。
(3)反汇编的截图
TODO:你的截图(必选)
TODO:你的文字解释说明(可选)
说明:在 debug 环境下执行 u 指令,显示出反汇编代码。
(4)在完成 DS 赋值后,立即显示各个寄存器的值
TODO:你的截图(必选)
TODO:你的文字解释说明(可选)
说明:第一步 t 指令执行底下为 MOV DS,AX,说明接下来执行的指令是 MOV DS,AX,所以再执行一次 t 指令,然后用 r 指令显示寄存器当前的值。
(5)在进行计算前,显示在数据段中定义的学号字符串的内存值的截图(只能显示完整的学号字符串,多显示、少显示均扣分)
TODO:你的截图(必选)
TODO:你的文字解释说明(可选)
说明:执行 d 指令显示出内存值,由于不能显示其他内存值,所以一行一行进行截图,可以看出,此时内存值与程序初始定义相同。
(6)在进行计算前,显示在数据段中定义的 COUNT 数组的内存值的截图(只能显示完整的 COUNT 数组内容,多显示、少显示均扣分)
TODO:你的截图(必选)
TODO:你的文字解释说明(可选)
说明:执行 d 指令显示出内存值,由于不能显示其他内存值,所以一行一行进行截图,可以看出,此时内存值与程序初始定义相同。
(7)执行完计算后,立即显示各个寄存器的值
TODO:你的截图(必选)
TODO:你的文字解释说明(可选)
说明:通过执行 g 指令跳到调用 dos 系统的指令,此时已经执行完运算,再执行 r 指令,显示出的就是执行完计算后各个寄存器的值。
(8)执行完计算后,显示在数据段中定义的学号字符串的内存值的截图(只能显示完整的学号字符串,多显示、少显示均扣分)
TODO:你的截图(必选)
TODO:你的文字解释说明(可选)
说明:执行 d 指令显示出内存值,由于不能显示其他内存值,所以一行一行进行截图,可以看出,此时内存值与程序期望相同。
(9)执行完计算后,显示在数据段中定义的 COUNT 数组的内存值的截图(只能显示完整的 COUNT 数组内容,多显示、少显示均扣分)
TODO:你的截图(必选)
TODO:你的文字解释说明(可选)
说明:说明:执行 d 指令显示出内存值,由于不能显示其他内存值,所以一行一行进行截图。学号为 2184411027,0-9 个数分别为 1,3,2,0,2,0,0,1,1,0,结果正确。
(10)程序在 DOSBox 下直接运行的截图
TODO:你的截图(必选)
TODO:你的文字解释说明(可选)
(11)源代码
TODO:你的源代码(必选)
;empty asm file
title I love asm
data segment
A
dw
0027h,0018h,0021h,0028h,0024h,0020h,0022h,0014h,0011h,0010h,0012h,0017h,0084h,0081h,0080H
B
dw
0018h,0021h,0028h,0024h,0020h,0022h,0044H,0041h,0011h,0010h,0012h,0017h,0084h,0081h,0080h,0040h,0042h,0047h,0002h,0007h
C1
DWDUP(?)data ends
code segment
assume cs:code, ds:data
main
proc
;assign the data segment base address to DS
mov
ax, data
mov
ds, ax
;TODO...;| add your code between arrows |
;v----------------------------v
mov si,0
MOV DI,0
mov bx,A
jiao:MOV CX,B[SI]
;想法是 di 指 A 中的元素,di 递增,拿 A 中的一个元素每次与B 中的所有元素作比较
CMP A[DI],CX
;直到遇到相同为止,如果没有遇到相同,有设置对 si 值比较指令,b 中所有元素相同时
JZ JI
;会跳转到 JI
INC SI
INC SI
CMP SI,42
JZ JI
JMP jiao
JI:mov ax,B[SI]
mov C1[DI],AX
MOV SI,0
INC DI
INC DI
CMP DI,32
JZ EXIT
JMP jiao
;^----------------------------^
;|
The END
|
;method 2: return to dos
EXIT: mov
ax, 4c00h
int
21h
main
关键词:上机实验 同步协同上机实验系统 实验评价 学习评价
中图分类号:G434文献标识码:A文章编号:1673-8454(2007)11-0017-03
一、引言
上机实验是高校实践教学的重要环节。以同步协同上机实验系统作为支撑环境的上机实验打破了教师规定实验题目和步骤,学生按部就班进行实验的传统模式。基于同步协同上机实验系统的上机实验完全以学生为主体,教师只起组织和引导作用;实验内容完全开放;学生在交流和协作中,进行研究性自主学习;实验学习成绩的评定不再单纯根据实验结果或实验报告,实现了实验学习评价的多元化和智能化。因此,研究基于同步协同上机实验系统的实验学习评价对于提高上机实验的教学质量非常重要。
二、同步协同上机实验系统简介
同步协同上机实验系统的目标是适应不同层次学生的需要,实现教学内容的开放性;促进学生之间交流和协作,以及教学双方的有效交流与互动;实现实验学习评价的多元化和智能化。系统的功能和模块结构划分如图1所示,包括协同上机实验、实验成绩管理、实验教学管理和系统管理四个模块。
采用同步协同上机实验系统作为支撑环境的上机实验过程一般是:(1)学生根据课程的进度预约实验;预约成功后,得到实验的时间、地点、目标要求和可选的实验任务;然后收集资料,自主学习准备实验。(2)在上机实验时,学生首先登录同步协同上机实验系统,协同完成实验任务,评价其他同学的实验方案、算法或运行结果;遇到问题可到讨论区求助;教师引导、辅助学生完成实验任务,评价实验任务。(3)在上机实验结束时,学生提交自我评价,系统辅助生成实验报告,并结合同学评价和教师评价成绩得出学生的上机实验成绩。
这种上机实验确定了具体的实验教学目标,并分解出了用于实验考核的基本实验教学目标。另外,根据实验目标设计了多项实验任务,供学生选择和修改;每项实验任务能够达到一项或多项实验目标。
三、上机实验学习评价的目的、功能和特点
上机实验学习评价是实验教学环节的重要组成部分,其目的在于促进学生学习。因此,上机实验学习评价旨在增强学生对于操作技能和理论知识的掌握和应用水平的了解;激励学生上机实验的积极性和主动性;端正学生的实验态度。基于同步协同上机实验系统的上机实验学习评价具有以下特点:
1.具有人文特性
基于同步协同上机实验系统的实验学习评价是体现“以人为本”的主体性评价,即“人文性”评价。因此,能最大限度地激发学生的内部动力,促进学生的自我发展。
2.形成性评价与总结性评价相结合
对学习结果的评价有助于帮助学生总结实验过程中的多种收获与体验,强化学生通过实验掌握的知识;而对学习过程的评价,可在不同程度上激发学生上机实验的积极性和主动性,帮助学生不断调整和改善实验过程,保证上机实验的学习质量。因此,不仅要重视学生的上机实验结果,而且要注意学生上机实验的整个过程。另外形成性评价应贯穿于上机实验的全过程,注重成绩背后的动机取向和努力程度,反映学生的操作能力、知识应用能力、协作能力、组织能力和创新能力等,以促进各种能力的培养。在不同课程的上机实验中,还应注重形成性评价与总结性评价各自权重的探索。
3.评价主体的多元性
评价是一种价值判断过程,这种价值是多元的。在评价情境中,不论评价者还是被评价者,不论教师还是学生,都是平等的主体。评价的主客体也应是多元的。[1] 因此基于同步协同上机实验系统的实验学习评价包括互相评价、自我评价、小组评价和教师评价四个部分。
互相评价是指学生在上机实验过程中,参考或验证其他同学的实验方案、算法和运行结果后给出的评价。其目的在于发现其他同学的优点及促进同学间的交流。
自我评价包括在实验准备阶段和在上机实验过程中的自我评价,以及在上机实验结束时提交的自我评分。
小组评价是指组长对于组内成员的评价。协同任务组是学生自发形成的;任务组的学生可以相互交流与协作共同完成一个任务;一个学生可以参加多个任务组。
教师评价指辅导教师在实验过程中,对于实验任务或学生的评价。对于每个学生和实验任务,教师评价不是必需的,其目的是通过评价引导上机实验的进行。
4.评价手段的数字化
在基于同步协同上机实验系统的上机实验中,学生通过同步协同上机实验系统进行交流与协作。辅导教师可以面对面辅导,也可以通过同步协同上机实验系统完成应答提问、删除与实验无关的讨论、删除不符合实验要求的任务和发布公开消息等辅导任务。因此大部分实验教学活动都基于同步协同上机实验系统完成,数字化信息是上机实验的重要载体。评价的表达以及数据的收集和整理呈现出数字化特征。
5.评价方式的开放性
基于同步协同上机实验系统的上机实验内容是开放的;不仅可以选择和修改,还可以自己设计实验任务。实验学习评价的内容也具有适当的开放性,而且评价既要适合于课程,又不局限于课程,能够体现学生在上机实验中的能力培养。通过评价不仅促进学生的学习,还要促进学生学会如何学习。
6.评价形式和标准的多样化
相对于一元智力理论和智力评价理论,多元智能理论能够更为准确地描绘和评价人类能力。[2] 实验教学在培养学生的动手能力和独立分析问题、解决问题的能力等方面起着决定性的作用。[3] 这决定了上机实验评价必须对学生进行多层面能力的考核,因此,评价形式和标准应该是多元的。在基于同步协同上机实验系统的上机实验中,可以针对实验任务的设计、实验方案的设计、程序运行情况、提出的问题、以及解决的问题进行多种形式的评价,各种评价贯穿于整个实验过程中。实验学习评价改变了单纯用实验数据结果和实验报告来衡量学生的成绩的方式,注重学生个体的实际情况和差异性。
四、上机实验评价的实施
同步协同上机实验系统作为上机实验的支撑环境,对上机实验的评价提供全方位的支持。基于同步协同上机实验系统的实验学习评价实施,一般要经过以下五个环节:
1.确定实验目标和可选的实验任务
实验目标作为实验评价的基础,根据实验教学大纲确定。为便于评价,实验目标必须详细、明确而具体。在此基础上,从中分解出基本目标。基本目标指学生通过实验最低限度掌握的知识点和基本技能或具备的基本能力。
实验教师根据实验目标设计多项实验任务,每项任务对应一个或多个实验目标。学生可以从中选择或修改实验任务,也可以自己设计实验任务;但是实验任务与实验目标的对应关系必须合理,否则会影响评价结果,甚至被辅导教师删除。
2.构建实验学习评价的评价模型
评价模型作为评价的依据和前提,应围绕实验目标的完成和多层次能力的培养来构建,旨在促进学生的发展。评价模型主要包括评价指标、质量等级和评分方法三个部分。评价指标反映达到实验目标和具备能力的各个方面;质量等级描述评价指标在质量上从好到差的档次,一般为优、良、及格和不及格;评分方法指整体评分和分项评分的计算方法,包括各个评价指标的标准分值。评价模型由实验教师根据学科特点和上机实验经验,与学生共同商定。在与学生共同商量的过程中,让学生掌握每项指标的含义,以及互相评价和自我评价的方法。例如基于同步协同上机实验系统的《VB程序设计》上机实验的自我评价表见表1,其中总分S=1.0a+0.8b+0.6c+0.4d(其中a、b、c、d分别是等级A、B、C、D对应指标的标准分值之和)。对于表1的评价结果,系统自动计算总分为:
S=1.0×(30+30+2)+0.8×(6+3+2+4)+0.6×(6+6+2)+0.4×9=86
3.在上机实验中不断进行评价
在上机实验过程中,学生要参考其他同学的实验内容或查阅问题的回答,系统强制必须在继续其它工作前进行评价;学生也可以自觉对其他同学的提问、问题的回答、任务设计、实验方案、算法和运行情况等进行评价。为了保证评价的公正性,系统仅显示每个同学的网名;为了激励学生的实验,学生可以随时查看自己当前的总评结果。辅导教师也可以对提问、问题的回答、任务设计、实验方案、算法和运行情况等进行评价,但不是必需的。学生在实验结束提交时,系统强制学生提交自我评价信息。
4.系统自动进行公正的评定与分析
在全部学生提交实验后,系统自动根据评价模型计算评价结果,并对评价结果进行分析,为每个学生提供改进性建议。
5.对评价结果进行反思
实验结束后,学生可以随时登录系统,查阅自己的实验成绩和改进建议,反思自己的行为。通过反思了解自己有哪些不足,哪些方面需要继续努力,促进实验学习。
五、结论
文章探讨了基于同步协同上机实验系统的实验学习评价,指出了评价的目标、应具备的特点和实施过程。这种实验学习评价能够激发学生在上机实验中进行创新性自主学习的积极性和主动性;有利于培养学生的多种能力、协作和创新精神;实现了实验成绩评定的多元化和智能化,降低了辅导教师的工作量;适合进行实验教学的多种改革研究。[4][5][6] 实践证实,基于同步协同上机实验系统的实验学习评价可操作性强,节约了辅导教师和学生的大量时间。
参考文献:
[1]赵蓉. 改进课堂学习评价的对策[J]. 思想政治课教学, 2006.1:81-83.
[2]任一波. 由多元智能理论引发对高职生学习评价的思考[J]. 浙江工商职业技术学院学报, 2006.1(5):91-93.
[3]原丕业,杜秀云,王军英. 高校实验教学改革初探[J]. 中国成人教育,2004.8:52-53.
[4]郑胜林,潘保昌, 徐杜等. 互动式实验教学模式之研究[J].高教探索,2004.1:80-81.
[5]贾雪枫. 创造性研究型教学模式初探[J]. 教育探索,2004.1:18-21.
一、实验目得 初步熟悉 MATLAB 工作环境,熟悉命令窗口,学会使用帮助窗口查找帮助信息。
命令窗口 二、实验内容(1)熟悉MATLAB平台得工作环境.(2)
熟悉MATLAB 得5 个工作窗口。
(3)MATLAB 得优先搜索顺序.三、实验步骤 1、熟悉MATLAB 得5 个基本窗口 ① mand Window(命令窗口)
② Workspace(工作空间窗口)
③ mand History(命令历史记录窗口)
④ Current Directory(当前目录窗口)
⑤ Help Window(帮助窗口)
(1)命令窗口(mand Window)。
在命令窗口中依次输入以下命令:
>〉x=1 〉> y=[1 2 3
4 5 68 9]; 〉> z1=[1:10],z2=[1:2:5]; >> w=linspace(1,10,10); >〉 t1=ones(3),t2=ones(1,3),t3=ones(3,1)〉〉 t4=ones(3),t4=eye(4)x =1 z1 =
3
9t1 =
1
1
t2 =
1
t3 =
1 t4 =
1
1
1
t4 =
1
0
0
0
0
0
0
0
0
0
0
0
0
思考题: ①
变量如何声明, , 变量名须遵守什么规则、就是否区分大小写。
答:(1)变量声明
1、局部变量 每个函数都有自己得局部变量,这些变量只能在定义它得函数内部使用。当函数运行时,局部变量保存在函数得工作空间中,一旦函数退出,这些局部变量将不复存在。
脚本(没有输入输出参数,由一系列MATLAB命令组成得M文件)没有单独得工作空间,只能共享调用者得工作空间。当从命令行调用,脚本变量存在基本工作空间中;当从函数调用,脚本变量存在函数空间中. 2、全局变量 在函数或基本工作空间内,用global声明得变量为全局变量。例如声明a为全局变量:
>〉global a 声明了全局变量得函数或基本工作空间,共享该全局变量,都可以给它曲赋值。
如果函数得子函数也要使用全局变量,也必须用global声明.3、永久变量 永久变量用persistent声明,只能在M文件函数中定义与使用,只允许声明它得函数存取.当声明它得函数退出时,MATLAB不会从内存中清除它,例如声明a为永久变量: 〉〉persistent a(2)变量命名规则如下: :
始于字母,由字母、数字或下划线组成;区分大小写;可任意长,但使用前N个字符.N与硬件有关,由函数namelengthmax返回,一般N=63; 不能使用关键字作为变量名(关键字在后面给出); 避免使用函数名作为变量名.如果变量采用函数名,该函数失效.(3)区分大小写
②
试说明分号、逗号、冒号得用法。
答:分号:加上分号“;“其作用就是将计算机结果存入内存,但不显示在屏幕上,反之,语句结尾若不加“;”,则表示在语句执行后,在将计算结果存入内存得同时,还将运算结果显示出来.逗号:分割列表
冒号:从什么到什么,比如1:10意思就是——1到10 ③ ③
l inspace())
称为“线性等分”函数, , 说明它得用法。可使用help命令,格式如下:
>>help linspace
④
数 说明函数 ones())、、zero s()、eye()得用法。
答;(1)ones()函数:全部元素都为 1 得常数矩阵;(2)zeros()函数:全部元素都为 0 得矩阵;(3)eye()函数:单位矩阵;
(4)linspace()函数:如 a=linspace(n1,n2,n3),表示在线性空
间上,行矢量得值从 n1 到 n2(2)工作空间窗口(Workspace). 单击工作空间窗口右上角得按钮,将其从MATLAB 主界面分离出来。
① 在工作空间查瞧各个变量,或在命令窗口用who, whos(注意大小写)查瞧各个 变量。
② 在工作空间双击变量,弹出Array Editor 窗口(数组编辑器窗口),即可修改变量。
③ 使用save 命令把工作空间得全部变量保存为my_var、mat 文件。
〉〉save my_var、mat
④ 输入下列命令: 〉〉clear all %清除工作空间得所有变量
观察工作空间得变量就是否被清空。使用load 命令把刚才保存得变量载入工作空间。
>〉load my_var、mat
⑤ 清除命令窗口命令: >〉clc(3)历史命令窗口(mand History)。
打开历史命令窗口,可以瞧到每次运行MATLAB 得时间与曾在命令窗口输入过得命 令,练习以下几种利用历史命令窗口重复执行输入过得命令得方法。
① 在历史命令窗口中选中要重复执行得一行或几行命令,右击,出现快捷菜单,选择 Copy,然后再Paste 到命令窗口。
② 在历史命令窗口中双击要执行得一行命令,或者选中要重复执行得一行或几行命令 后,用鼠标将其拖动到命令窗口中执行。
③ 在历史命令窗口中选中要重复执行得一行或几行命令,右击,出现快捷菜单,选择 Evaluate Selection,也可以执行。
④ 或者在命令窗口使用方向键得上下键得到以前输入得命令。例如,按方向键“↑” 一次,就重新将用户最后一次输入得命令调到MATLAB 提示符下。重复地按方向上键 “↑”,就会在每次按下得时候调用再往前一次输入得命令。类似地,按方向键“↓”得时 候,就往后调用一次输入得命令.按方向键“←“或者方向键“→”就会在提示符得命令 中左右移动光标,这样用户就可以用类似于在字处理软件中编辑文本得方法编辑这些命令。
(4)
当前目录命令窗口(Current Directory).MATLAB 得当前目录即就是系统默认得实施打开、装载、编辑与保存文件等操作时得文 件夹。打开当前目录窗口后,可以瞧到用“save”命令所保存得my_var、mat 文件就是保存在 目录C:MATLAB6p5work 下。
(5)帮助窗口(Help Window)。
单击工具栏得图标,或选择菜单View|Help,或选择菜单Help|MATLAB Help 都能 启动帮助窗口。
① 通过Index 选项卡查找log2()函数得用法,在Search index for 栏中输入需要查找得 词汇“log2”,在左下侧就列出与之最匹配得词汇条目,选择“log2[1]”,右侧得窗口就 会显示相应得内容。
② 也可以通过Search 选项卡查找log2()函数得用法.Search 选项卡与Index 选项卡 不同,Index 只在专用术语表中查找,而Search 搜索得就是整个HTML 帮助文件。
2、MATLAB 得数值显示格式设置 屏幕显示方式有紧凑(pact)与松散(Loose)两种,其中Loose 为默认方式.>〉a=ones(1,30)
〉>format pact
〉〉a 数字显示格式有short、long、short e、long e 等,请参照教材得列表练习一遍。
〉〉format long >>pi
〉〉format short >〉pi
〉>format long 〉〉pi
>>format + 〉〉pi
〉>—pi
3、变量得搜索顺序 在命令窗口中输入以下指令:
>>pi
〉>sin(pi);
>>exist(’pi”)
〉>pi=0;>〉exist(’pi“)
>〉pi
〉〉clear pi 〉>exist('pi’)>>pi
思考题:① 3 次执行exist(’pi’)得结果一样吗?如果不一样,试解释为什么? 答:不一样,pi原来就是库存函数,但就是如果被赋值则系统默认被赋予得值为pi后来得值,但就是当执行clear pi之后所赋得值被清空,因此pi得值又成为3、1416 ② 圆周率 pi 就是系统得默认常量,为什么会被改变为 0? 答:pi原来就是库存函数,但就是如果被赋值则系统默认被赋予得值为 pi 后来得值,但就是当执行 clear pi之后所赋得值被清空,因此 pi 得值又成为3、1416
实验二
MA TLAB语言基础 一、实验目得 基本掌握 MATLAB 向量、矩阵、数组得生成及其基本运算(区分数组运算与矩阵运算)、常用得数学函数。交接字符串得操作.二、实验内容(1)向量得生成与运算。
(2)矩阵得创建、引用与运算。
(3)
多维数组得创建及运算。
(4)字符串得操作。
三、实验步骤 1、向量得生成与运算 1)
向量得生成 直接输入法: A =
>〉 B=[1;2;3;4;5]
B =
1
2
5 冒号生成发:
〉〉 A=1:2:10 ,B=1:10,C=10:—1:1
A =
3
5
7
B =
5
C =
10
8
4
1 函数法:
Linspace()就是线性等分函数,logspace()就是对数等分函数。
>〉 A=linspace(1,10),B=linspace(1,30,10)A =
Columns 1 through 91、0000
1、0909
1、1818
1、2727
1、3636
1、4545
1、5455
1、6364
1、7273
Columns 10 through 181、8182
1、9091
2、0000
2、0909
2、1818
2、2727
2、3636
2、4545
2、5455
Columns 19 through 272、6364
2、7273
2、8182
2、9091
3、0000
3、0909
3、1818
3、2727
3、3636
Columns 28 through 36
3、4545
3、5455
3、6364
3、7273
3、8182
3、9091
4、0000
4、0909
4、1818
Columns 37 through 454、2727
4、3636
4、4545
4、5455
4、6364
4、7273
4、8182
4、9091
5、0000
Columns 46 through 545、0909
5、1818
5、2727
5、3636
5、4545
5、5455
5、6364
5、7273
5、8182
Columns 55 through 63
5、9091
6、0000
6、0909
6、1818
6、2727
6、3636
6、4545
6、5455
6、6364
Columns 64 through 72
6、7273
6、8182
6、9091
7、0000
7、0909
7、1818
7、2727
7、3636
7、4545
Columns 73 through 81
7、5455
7、6364
7、7273
7、8182
7、9091
8、0000
8、0909
8、1818
8、2727
Columns 82 through 908、3636
8、4545
8、5455
8、6364
8、7273
8、8182
8、9091
9、0000
9、0909
Columns 91 through 999、1818
9、2727
9、3636
9、4545
9、5455
9、6364
9、7273
9、8182
9、9091
Column 100
10、0000
B =
Columns 1 through 91、0000
4、2222
7、4444
10、666713、8889
17、1111
20、333323、555626、7778
Column 1030、0000 >> A=logspace(0,4,5)A =
1
10
100
1000
10000 练习:
使用 logspace()创建1—4得有 10 个元素得行向量。
答案:
>〉 A=logspace(1,10,4*pi)
A =1、0e+010 *
Columns 1 through 9
0、0000
0、0000
0、0000
0、0000
0、0000
0、0000
0、0001
0、0005
0、0035
Columns 10 through 120、0231
0、1520
1、0000 2)向量得运算 维数相同得行向量只见可以相加减,维数相同得列向量也可以相加减,标量可以与向量直接相乘除。
〉〉 A=[1 2 3 4 5],B=3:7,A =
4
B =
7
>〉 AT=A”,BT=B“, AT =
BT =
6
>〉 E1=A+B, E2=A-B
E1 =
6
10E2 =
—2
—2
—2
〉〉 F=AT—BT, F =
-2
—2
〉〉 G1=3*A, G2=B/3, G1 =
12
15
G2 =1、0000
1、3333
1、6667
2、0000
2、3333 向量得点积与叉积运算。
〉〉 A=ones(1,10);B=(1:10); BT=B”;>> E1=dot(A,B)
E1 =
55
〉〉 E2=A*BT
E2 =
〉〉 clear >> A=1:3,B=3:5,A =
3
B =
〉〉 E=cross(A,B)
E = -2
—2 2.矩阵得创建、引用与运算
1)矩阵得创建与引用
矩阵就是由元素构成得矩阵结构,行向量与列向量就是矩阵得特殊形式。
直接输入法:
〉〉 A=[1 2 3;4 5 6]
A =
1
3
4
6
〉〉 B=[ 1 4 75 86 9 ] B =
1
7
8
〉> A(1)
ans =〉〉 A(4:end)
ans =
6 〉〉 B(:,1)ans =
2〉> B(:)ans =
1
8>〉 B(5)
ans =抽取法 〉> clear 〉> A=[1 2 3 4;5 6 7 8;9 10 11 12;13 14 15 16] A =
3
4
10
12
13
16 >〉 B=A(1:3,2:3)
B =
6
10
11 〉> C=A([1 3],[2 4])
C =
10
12 〉〉 A([1 3;2 4])ans =
9
5 函数法: 〉〉 A=ones(3,4)A =
1
1
>〉 B=zero(3)
??? Undefined function or method “zero’ for input arguments of type ”double’、>> B=zeros(3)B =
0
0
0
0
0
0
0
0
0 >> C=eyes(3,2)??? Undefined function or method “eyes’ for input arguments of type ”double’、〉> C=eye(3,2)C =
0
0
0
0 >> D=magic(3)D =
6
拼接法 〉〉 clear 〉> A=ones(3,4)A =
1
1
1
1
>〉 B=zeros(3)
B =
0
0
0
0
0
0
0
0
0 >〉 C=eye(4)C =
0
0
0
0
0
0
0
0
0
0
0
0
>〉 D=[A B] D =
1
0
0
0
0
0
0
0
0
0 >〉 F=[A;C] F =
1
1
1
0
0
0
0
1
0
0
0
0
0
0
0
0
拼接函数与变形函数法:
〉〉 clear 〉> A=[0 1;1 1] A =
0
1
1
〉> B=2*ones(2)
B =
2
>> cat(1,A,B,A)ans =
0
1
1
2
2
0
>〉 cat(2 A,B,A)
??? cat(2 A,B,A)Error: Unexpected MATLAB expression、>〉 cat(2,A,B,A)ans =
0
1
0
1
2
2
1
1 >〉 repmat(A,2,2)ans =
0
1
0
1
1
1
0
0
1
1
1
〉> repmat(A,2)ans =
0
0
1
0
1
0
1
练习:使用函数法、拼接法、拼接函数法与变形函数法,按照要求创建以下矩阵:A 为
得全1矩阵、B 为得 0 矩阵/C 为得单位矩阵、D 为得魔方阵、E 由 C
与D纵向拼接而成,F 抽取 E 得2—--5 行元素生成、G 由 F 经变形为得矩阵而得、以 G 为子矩阵用复制函数(repmat)生成得大矩阵H。
答案 :
>〉 A=ones(3,4)
A =
1
1
1
1 〉〉 B=zeros(3,3)B =
0
0
0
0
0
0
0
0
0 >〉 C=eye(3)
C =
1
0
0
0
0
0
0
>〉 D=magic(3)D =
5
4
>> E=[C;D] E =
1
0
0
0
0
0
0
8
〉〉 F=(2:5,:)??? F=(2:5,:)Error: Expression or statement is incorrect—-possibly unbalanced(, {, or [、>〉 F=E(2:5,:)F =
0
0
0
0
〉〉 G=respace(E,3,4)
??? Undefined mand/function “respace’、〉>
G=respace(F,3,4)
??? Undefined mand/function ”respace“、〉> G=reshape(F,3,4)
G =
0
0
1
6
0
0
〉〉 H=repmat(G,2,2)H =
0
0
0
5
0
8
0
0
7
0
0
0
3
1
0
1
0
0
1
0
0
0
0
2)矩阵得运算 矩阵得加减、数乘与乘法
已知矩阵: >> A=[1 2
-1], A =
-1 >〉 B=[-1 0
2] B =
0
>〉 A+B ans =
0
4
1 〉〉 2*A ans =
4
-2 >〉 2*A-3*B ans =
-8 〉〉 A*B ans =
4-4
—2 矩阵得逆矩阵 〉〉 format rat;A=[1 0 1;2 1 2;0 4 6] A =
1
0
0
>> A1=inv(A)A1 =
—1/3
2/3
-1/6
-2
1
0
4/3
-2/3
1/6
>〉 A*A1
ans =
0
0
0
0
0
0
矩阵得除法 >〉 a=[1 2 1;3 1 4;2 2 1],b=[1 1 2],d=b’ a =
1
3
2
b =
1
d =
1
1
>> c1=b*inv(a),c2=b/a c1 =
6/7
3/7
-4/7
c2 =
6/7
3/7
-4/7
>> c3=inv(a)*d, c4=ab c3 =
2/7
-4/7
??? Error using ==〉 mldivide Matrix dimensions must agree、>〉 c3=inv(a)*d,c4=ad c3 =
2/7
-4/7
c4 =
2/7
-4/7
练习: 按下列要求求出各种得矩阵运算得值
求矩阵得秩、特征值与特征向量、矩阵得乘幂与开方;
矩阵得指数与对数
矩阵得提取与翻转 答案:
〉〉 A=[6 3
4 3
—2 5 7 —4
8-1 —3 —7] A =
6
-4
-1
—3
—7 〉> B=rank(A)B =〉> rb=rank(A)
rb =>> [X,Lamda]=eig(A)
??? Error using ==〉 eig Matrix must be square、〉〉 [X,Lamda]=eigs(A)??? Error using ==〉 eigs A must be a square matrix or a function which putes A*x、〉〉 C=[6 3 4
-2 5 7
—1-3] C =
4
—2
8
—1
—3 >> [X,Lamda]=eigs(C)
X =0、8013
-0、1094
-0、1606
0、3638
—0、6564
0、86690、4749
0、7464
-0、4719 Lamda =9、7326
0
0
0
-3、2928
0
0
0
1、5602 〉〉 [X,Lamda]=eig(C)
X =
0、8013
—0、1094
—0、16060、3638
—0、6564
0、86690、4749
0、7464
-0、4719 Lamda =9、7326
0
0
0
-3、2928
0
0
0
1、5602 〉> [X,Lamda]=eig(C)
X =
0、8013
—0、1094
-0、1606
0、3638
—0、6564
0、8669
0、4749
0、7464
—0、4719 Lamda =9、7326
0
0
0
—3、2928
0
0
0
1、5602 〉> D=A^2 ??? Error using ==> mpower Matrix must be square、〉> D=C^2 D =
33
34
〉〉 E=sqrtm(C)E =
2、2447 + 0、2706i0、6974 — 0、1400i0、9422-0、3494i
-0、5815 + 1、6244i
2、1005-0、8405i
1、7620-2、0970i
1、9719 - 1、8471i
-0、3017 + 0、9557i
0、0236 + 2、3845i >> F=expm(C)F =
1、0e+004 *
1、0653
0、5415
0、6323
0、4830
0、2465
0、2876
0、6316
0、3206
0、3745 〉〉 G=logm(C)Warning: Principal matrix logarithm is not defined for A with
nonpositive real eigenvalues、A non-principal matrix
logarithm is returned、〉 In funm at 153
In logm at 27 G =
1、7129 + 0、4686i0、5305-0、2425i0、5429-0、6049i
1、1938 + 2、8123i0、3658 — 1、4552i
—0、5514-3、6305i
-0、0748 - 3、1978i
0、7419 + 1、6546i1、8333 + 4、1282i >> H=fliplr(C)H =
5
-2
—1
>> I=triu(C)I =
0
0
0
-3 〉〉 J=tril(C)
J =
0
0
-2
0
-3 >> K=diag(C)K =
6-3 3.多维数组得创建及运算 1)多维数组得创建 >〉 A1=[1,2,3;4 5 6;7,8,9];A2=reshape([10:18],3,3)A2 =
10
13
16
14
17
12
15〉〉 T1(:,:,1)=ones(3);T1(:,:,2)=zeros(3)T1(:,:,1)
=
1
1
1
T1(:,:,2)=
0
0
0
0
0
0
0
0
0 〉> T2=ones(3,3,2)
T2(:,:,1)=
1
1
1 T2(:,:,2)=
1
1
1 >> T3=cat(3,A1,A2),T4=repmat(A1,[1,1,2])
T3(:,:,1)=
T3(:,:,2)=
10
13
16
17
15T4(:,:,1)
=
2
6
8
T4(:,:,2)=
7
8
2)多维数组得创建
数组运算用小圆点加在运算符得前面表示,以区分矩阵得运算。特点就是两个数组相对应得元素进行运算。
〉〉 A=[1:6];B=ones(1,6);>> C1=A+B,C2=A—B C1 =
2
C2 =
0
3
>> C3=A、*B,C4=B、/A,C5=A、B C3 =
6 C4 =
1、0000
0、5000
0、3333
0、2500
0、2000
0、1667 C5 = 1、0000
0、5000
0、3333
0、2500
0、2000
0、1667
关系运算或逻辑运算得结果都就是逻辑值.>〉 I=A〉3,C6=A(I)
I =
0
0
0
1
1 C6 =
4
6 〉〉 A1=A-3,I2=A1&A A1 =
—2
—1
0
2
I2 =
1
1
0
1
〉〉 I3=~I I3 =
1
0
0
0 4.字符串得操作 1)字符串得创建
>〉 S1=”Ilike MATLAB’ S1 = Ilike MATLAB >> S2=“I'’m a stuent、” S2 = I“m a stuent、〉> S3=[S2,”and’,S1] S3 = I“m a stuent、andIlike MATLAB 2)求字符串长度 〉〉 length(S1)ans =>> size(S1)ans =
13)字符串与一维数值数组得相互转换 >> CS1=abs(S1)CS1 =
101
32
65
84
66 〉> CS2=double(S1)CS2 =
108
105
101
32
77
65
76
>> char(CS2)
ans = Ilike MATLAB >> setstr(CS2)
ans = Ilike MATLAB 练习:用char()与向量生成得方法创建如下字符串AaBbCc、、、、、、XxYyZz、〉> S1=65:90;S2=97:122;〉〉 C=[S1;S2];>>
C=C(:)’;〉>
S3=double(C);char(S3)ans = AaBbCcDdEeFfGgHhIiJjKkLlMmNnOoPpQqRrSsTtUuVvWwXxYyZz 实验三
MAL TAB B 数值运算 一、实验目得
掌握 MATLAB 得数值运算及其运算中所用到得函数,掌握结构数组与细胞数组得操作。
二、实验内容
1)
多项式运算。
2)
多项式插值与拟合。
3)
数值为积分。
4)
结构数组与细胞数组。
三、实验步骤
1、多项式运算
1)
多项式表示。在 MATLAB 中,多项式表示成向量得形式.如: 在 MATLAB 中表示为 >> s=[1 3-5 0 9] 2)
多项式得加减法相当于向量得加减法,但必须注意阶次要相同。如不同,低阶次得要补 0。如多项式与多项式相加。
〉> s1=[0 0 2 3 11] 〉〉 s2=[1 2 —5 4 7] 〉〉 s3=s1+s2 答;s1 =
0
0
3s2 =
—5
4
s3 =
—3
18 3)
多项式得乘、除法分别用函数 conv与得 deconv实现。
〉> s1=[2 3 11] >> s2=[1 3 —5 4 7] >〉 s3=conv(s1,s2)
>> s4=deconv(s3,s1)答;s1 =
11 s2 =
1
3
4
s3 =
—29
65
s4 =
1
7 4)
多项式求根用函数roots。
>〉 s1=[2 4 2] >> roots(s1)
答;s1 =
ans =
—1
5)
多项式求值用函数 polyval >> s1=[2 4 1-3] 〉〉 polyval(s1,3)
〉> x=1:10 〉〉 y=polyval(s1,x)答;s1 =
2
-3 ans =
90 x =
1
6
9
10 y =
Columns 1 through 8
193
352
579
886
1285
Columns 9 through 10 练习:求得“商”及余数.〉> s1=[1 0 1];s2=[1 3];s3=[1 1]; 〉> s4=[1 0 2 1]; >> [q,r]=deconv(conv(conv(s1,s2),s3),s4)答;q =
1
4 r =
0
0
-1 2、多项式插值与拟合 有一组实验数据如附表 1—1 所示。
附表 1-1 X 1 2 3 4 5 6 7 8 9 10 Y 16 32 7
1
分别用拟合(二阶至三阶)与插值(线性与三次样条)得方法来估算X=9、5 时 Y 得值。以下就是实现一阶拟合得语句。
〉> x=1:10 >〉 y=[16 32 70 142 260 436 682 1010 1342 1960] 〉> p1=polyfit(x,y,1)
%一阶拟合 >> y1=polyval(p1,9、5)
%计算多项式 P1 在x=9、5 得值 答;x =
2
9y =
Columns 1 through 8
32
70
142
260
436
682
1010
Columns 9 through 10
1342
1960 p1 =
200、9818 -510、4000 y1 =1、3989e+003 3、数值微积分
1)
差分使用 diff 函数实现。
>〉 x=1:2:9 〉〉 diff(x)答;x =
5
ans =
2 2)
可以用因变量与自变量差分得结果相处得到数值微分。
〉> x=linspace(0,2*pi,100);>〉 y=sin(x); >〉 plot(x,y)>〉 y1=diff(y)、/diff(x); 〉〉 plot(x(1:end-1),y1)答;
3)
cumsum 函数求累计积分,trapz 函数用梯形法求定积分,即曲线得面积。
〉〉 x=ones(1,10)
>> cumsum(x)>> x=linspace(0,pi,100); >> y=sin(x); 〉〉 trapz(x,y)〉> p=cumsum(y); >> p(100)*pi/(100—1)答;x =
1
1
ans =
9ans =
1、9998 ans =
1、9998 练习:图 A1就是瑞士地图,为了算出其国土面积,首先对地图作如下测量:以由西向东方向为 X 轴,由南向北方向为 Y 轴,选择方便得原点,并将从最西边界点到最东边界点在 X 轴上得区间适当划分为若干段,在每个分点得Y方向测出南边界点与北边界点得 Y 坐标Y1与 Y2,这样就得到表 1,根据地图比例尺回到18mm相当于 40Km,试由测量数据计算瑞士国土近似面积,与其精确值 41228 比较。地图得数据见附表 1—2(单位mm).附表 1—2 X 7 10、5 13 17、5 34 40、5 44、5 48 56 61 68、5 76、5 80、5 91 Y1 44 45 47 5
34 41 45 46 Y2 44 59 7
18 1 18 续表 X 96 101 104 106、5 111、5 118 123、5 1 36、5 142 146 150 157 158 Y1 43 37 33 28 32 65 55 54 52 50 66 66 68 Y2 121 124 1 21 1 21 121 116 1 22 83 81 82 86 85 68 提示:由高等数学得知识,一条曲线得定积分就是它与 x 轴所围成得面积,那么两条曲线所围成得面积可由两条曲线得定积分相减得到。、结构数组与细胞数组 1)
机构数组得创建. 〉> student、number=’20050731001’;〉〉 student、name=’Jack”;〉〉 student(2)、number=’20050731002’;〉〉 student(2)、name=’Lucy“; >> student 或者用 struct 函数创建。
>> student=struct(”number“,{’001”,“002’},”name',{'Jack’,“Lucy’});答;student =
1x2 struct array with fields:
number
name 2)
机构数组得操作。
〉> student(1)、subject=[]
%添加 subject 域并赋予空值 〉〉 student(1)、score=[] >〉(studeng)〉>fieldnames(studeng)
〉>fieldnames(student)
〉〉 getfield(student,{2},'name')〉〉 student=rmfield(student,”subject’)
%删除 subject 域 〉〉 student=setfield(student,{1},“score’,90);>> student(2)、score=88;
%比较与上一条语句就是否效果一样 答;student =
1x2 struct array with fields:
number
name
subject student =
1x2 struct array with fields:
number
name
subject
score ??? Undefined function or variable ’studeng’、练习:创建一结构数组stusorce,其域为:No,Name,English,Math,Chinese,Total,Average。结构数组得大小为 2×2。
3)
细胞数组得创建。
〉> A={’How are you!”,ones(3);[1 2;3 4],{“cell’}};
%直接创建
〉〉 B(1,1)={’Hello world”};
%由各个细胞元素创建
>> B(1,2)={magic(3)};
>〉 B(2,1)={[1 2 3 4]};答
或者用cell 函数先创建空得细胞数组,然后再给各个元素赋值 c=cell(1,2); >> c(1,1)={’Hello world’};>> c(1,2)={magic(3)}; 〉〉 c(1,3)={[1 2 3 4]}; 4)
细胞数组得操作。
〉> ans1=A(1,1)>> ans2=A(1,1)>〉 whos ans1 ans2 〉〉 celldisp(A)〉> a1=A{2,1}(1,2)>〉 [a2 a3]=deal(A{1:2})
答;ans1 =
’How are you!’ ans2 =
’How are you!’
Name
Size
Bytes
Class
Attributes
ans1
1x1
84
cell
ans2
1x1
84
cell
A{1,1} = How are you!
A{2,1} =
A{1,2} =
1
A{2,2}{1} = cell a1 =a2 = How are you!a3 =
1
2
实验四
MA LT AB B 符号运算
一、实验目得
掌握符号变量与符号表达式得创建,掌握MALTAB得symbol工具箱得一些基本运用。
二、实验内容
1)
符号变量、表达式、方程及函数得表示。
2)
符号微积分运算.
3)
符号表达式得操作与转换.4)
符号微分方程求解.三、实验步骤1、符号运算得引入
在数值运算中如果求,则可以不断让得让x趋近0,一球得表达式趋近什么数,但终究不能令 x=0,因为在数值运算中 0 不就是能作除数得。MATLAB得符号运算能解决这内问题。输入如下命令: 〉> f=sym(’sin(pi*x)/x')
>> limit(f,’x“,0)答;f = sin(pi*x)/x ans = pi 2 2、符号常量、符号变量、符号表达式得创建
1)
使用 sym()创建 输入以下命令,观察 Workspace 中 A、B、f就是什么内性得数据,占用多少字节得内存空间。
>〉 A=sym(”1’)
%符号常量 〉〉 B=sym(’x’)
%符号变量 >> f=sym(’2*x^2+3*x-1’)
%符号表达式 >> clear >> f1=sym(“1+2”)
%有单引号,表示字符串 >〉 f2=sym(1+2)
%无单引号 >〉 f2=sym(1+2)
〉> f4=sym(“2*x+3’)
%为什么出错 〉〉 x=1 〉〉 f4=sym(2*x+3)
答;A = 1 B = x f = 2*x^2+3*x-1 f1 = 1+2 f2 = 3 f2 = 3 f4 = 2*x+3 x =f4 = 5 通过瞧 MATLAB 得帮助可知,sym()得参数可以使字符串或就是数值类型,无论就是哪种类型都会生成符号类型数据。
2)
使用 syms创建 〉> clear >〉 syms x y z 〉> x,y,z 〉〉 f1=x^2+2*x+1 >〉 f2=exp(y)+exp(z)^2 >> f3=f1+f2 答;x = x y = y z = z f1 = x^2+2*x+1 f2 = exp(y)+exp(z)^2 f3 = x^2+2*x+1+exp(y)+exp(z)^2 3 3、符号矩阵创建
>〉 syms a1 a2 a3 a4 〉〉 A=[a1 a2;a3 a4] 〉> A(1),A(3)答;A = [ a1,a2] [ a3, a4] ans = a1 ans = a2 4 4、符号算术运算
1)
符号向量相乘、相除 符号量相成与数值量相乘一样,分成矩阵乘与数组乘。
〉> a=sym(5);b=sym(7); 〉> c1=a*b >〉 c2=a/b >> a=sym(5);B=sym([3 4 5]);>〉 C1=a*B,C2=aB >> syms a b >> A=[5 a;b 3];B=[2*a b;2*b a];〉> C1=A*B,C2=A、*B >〉 C3=AB,C4=A、/B 答;c1 =
35 c2 = 5/7 C1 = [ 15, 20, 25] C2 = [ 3/5, 4/5,1] C1 = [ 10*a+2*a*b,5*b+a^2] [
2*a*b+6*b,b^2+3*a] C2 = [
10*a,a*b] [ 2*b^2,3*a] C3 = [
2*a*(b—3)/(—15+a*b),(a^2-3*b)/(—15+a*b)] [
2*b*(a-5)/(-15+a*b), —(5*a-b^2)/(-15+a*b)] C4 = [ 5/2/a,a/b] [
1/2, 3/a] 2)
符号数值任意精度控制与运算 任意精度得 VPA 运算可以使用命令 digits(设定默认得精度)与vpa(对指定对象以新得精度进行计算)来实现。
>〉 a1=sym(’2*sqrt(5)+pi')
〉> a=sym(’2*sqrt(5)+pi’)〉〉 b=sym(2*sqrt(5)+pi)
>> digits >> vpa(a)>> digits(15)>〉 vpa(a)〉〉 c1=vpa(a,56)
〉〉 c2=vpa(b,56)
答 a1 = 2*sqrt(5)+pi a = 2*sqrt(5)+pi b = 8572296331135796*2^(-50)
Digits = 32 ans = 7、637263128 ans = 7、637 c1 = 7、6372631285535581572696
c2 = 7、6372726175781250000000 注意观察c1 与c2 得数值类型,c1 与 c2 就是否相等。
3)
符号类型与数值类型得转换 使用命令sym可以把数值型对象转换成有理数性符号对象,命令vpa可以讲数值型对象转换为任意精度得 VPA 型符号对象.使用 double,numeric 函数可以将有理数型与 VPA 型符号对象转换成数值对象、〉> clear >〉 a1=sym('2*sqrt(5)+pi’)
>〉 b1=double(a1)
%符号转数值 >> b2=isnumeric(b1)
%判断就是否转换成了数值 >> a2=vpa(a1,70)
%数值转符号 答;a1 = 2*sqrt(5)+pi b1 =7、6137 b2 =
1 a2 = 7、6372631285535589083128858 5 5、、符号表达式得操作与转换
1)独立变量得确定原则 独立变量得确定原则:在符号表达式中默认变量就是惟一得.MATLAB 会对单个英文小写字母(除 i、j外)进行搜索,且以 x 为首选独立变量。如果表达式中字母不唯一,且无 x,就选在字母表最接近x 得字母.如果有相连得字母,则选择在字母表中较后得那一个。例如:中,y 就是默认独立变量。,t 就是默认独立变量。
输入以下命令,观察并分析结果。
>> clear >〉 f=sym(”a+b+i+j+x+y+xz’)
>〉 findsym(f)
〉> findsym(f,1)>〉 findsym(f,2)>> findsym(f,3)
>〉 findsym(f,4)〉〉 findsym(f,5)>〉 findsym(f,6)答;f = a+b+i+j+x+y+xz ans = a, b, j,x,xz, y ans = x ans = x,xz ans =
x,xz,y ans = x,xz,y,j ans = x,xz,y,j,b ans = x,xz,y,j,b,a 2)符号表达式得化简 符号表达式化简主要包括表达式美化(pretty)、合并同类项(collcet)、多项式展开(expand)、因式分解(factor)、化简(simple或 simplify)等函数。
①合并同类项(collect)。分别按 x 得同幂项与e指数同幂项合并表达式: . 〉〉 syms x t; >〉 f=(x^2+x^exp(—t)+1)*(x+exp(-t));〉〉 f1=collect(f)〉〉 f2=collect(f,’exp(-t)“)答;f1 = x^3+exp(-t)*x^2+(x^exp(-t)+1)*x+(x^exp(-t)+1)*exp(-t)
f2 =(x^2+x^exp(-t)+1)*exp(-t)+(x^2+x^exp(—t)+1)*x ②对显示格式加以美化(pretty)。针对上例,用格式美化函数可以使显示得格式更符合数学书写习惯。
>〉 pretty(f1)>〉 pretty(f2)
答;f1 = x^3+exp(-t)*x^2+(x^exp(—t)+1)*x+(x^exp(—t)+1)*exp(-t)
f2 =(x^2+x^exp(-t)+1)*exp(-t)+(x^2+x^exp(—t)+1)*x >〉
pretty(f1)pretty(f2)
3
2
exp(-t)
exp(-t)
x
+ exp(—t)
x
+(x
+ 1)
x +(x
+ 1)exp(—t)
2
exp(—t)
exp(-t)
(x
+ x
+ 1)exp(-t)+(x
+ x
+ 1)x 注意与直接输出得 f1 与 f2 对比。
③多项式展开(expand)。展开成 x 不同次幂得多项式、〉> syms x 〉> f=(x—1)^12;〉> expand(f)〉〉 pretty(expand(f))
答;ans = 1+x^12-12*x^11+66*x^10—220*x^9+495*x^8—792*x^7+924*x^6-792*x^5+495*x^4-220*x^3+66*x^2—12*x
7
6
5
+ x
- 12 x
+ 66 x
— 220 x
+ 495 x
x
+ 924 x
- 792 x
+ 495 x
- 220 x
+ 66 x
- 12 x ④ 因式分解(factor)。将表达式做因式分解。
>> syms x;f=x^12—1; 〉> pretty(factor(f))
答;ans = 1+x^12—12*x^11+66*x^10—220*x^9+495*x^8-792*x^7+924*x^6—792*x^5+495*x^4—220*x^3+66*x^2-12*x
12
11
8
6
5
+ x
— 12 x
+ 66 x
— 220 x
+ 495 x
x
+ 924 x
— 792 x
3
+ 495 x
— 220 x
+ 66 x
- 12 x 〉〉 syms x;f=x^12—1;pretty(factor(f))
2
2
(x-1)(1 + x
+ x)(1 + x)
(1 — x + x)(1 + x)(x
— x
+ 1)⑤化简(simple或 simplify)。
将函数化简.〉〉 clear 〉〉 syms x;f=(1/x^3+6/x^2+12/x+8)^(1/3); 〉> g1=simple(f)>〉 g2=simplify(f)
答;g1 =(2*x+1)/x g2 =((2*x+1)^3/x^3)^(1/3)6 6、符号表达式得变量替换
subs 函数可以对符号表达式中得符号变量进行替换 >> clear >> f=sym(’(x+y)^2+4*x+10’)〉〉 f1=subs(f,’x',’s')
%使用 s 替换x >〉 f2=subs(f,”x+y’,“z’)答;f =(x+y)^2+4*x+10 f1 =(s+y)^2+4*s+10 f2 = z^2+4*x+10 7、符号极限、符号积分与微分
1)
求极限函数得调用格式 limit(F,x,a)
%返回符号对象 F 当 x→a 时得极限 limit(F,a)
%返回符号对象 F 当独立变量*→a 时得极限
limit(F)
%返回符号对象 F 当独立变量→0(a=0)时得极限 limit(F,x,a,’right’)
%返回符号对象 F 当 x→a 时得右极限 limit(F,x,a,’left”)
%返回符号对象F当 x→a 时得左极限 例一:
〉〉 clear >> f=sym(“sin(x)/x+a*x”)
>〉 limit(f,“x’,0)
%以 x 为自变量求极限 >> limit(f,”a’,0)
%以 a 为自变量求极限 >> limit(f)
%在默认情况下以 x 为自变量求极限 >> findsym(f)
%得到变量并且按字母表顺序排列 答‘f = sin(x)/x+a*x ans = 1 ans = sin(x)/x ans = 1 ans = a,x 例二: >> clear >〉 f=sym('sqrt(1+1/n));
〉〉 limit(f,n,inf)
%求 n 趋于无穷大时得极限
2)
求积分函数得调用格式 int(F)
%求符号对象 F 关于默认变量得不定积分 int(F,v)
%求符号对象 F 关于指定变量 v 得不定积分 int(F,a,b)
%求符号对象 F 关于默认变量得从 a 到 b 得定积分 int(F,v,a,b)
%求符号对象 F 关于指定变量得从 a 到b得定积分 3)
求微分方程得调用格式 diff(F)
%求符号对象 F 关于默认变量得微分 diff(F,v)
%求符号对象 F 关于指定变量 v 得微分 diff(F,n)
%求符号对象 F 关于默认变量得 n 阶微分,n 为自然数 1、2、3…… diff(F,v,n)%求符号对象 F 关于指定变量 v 得 n 阶微分 8 8、符号方程求解
1)常规方程求解函数得调用格式 g=solve(eq)
%求方程(或表达式或字串)eq 关于默认变量得解
g=solve(eq,var)
%求方程(或表达式或字串)eq 关于指定变量 var 得解 g=solve(eq1,eq2,…、,eqn,var1,var2,…,varn)
%求方程(或表达式或字串)eq1,eq2,eq3,……eqn 关于指定变量组var1,var2,……,varn)得解 求一元二次方程得解.其求解方法有多种形式:
① seq=solve(“a*x^2+b*x+c')
② seq=solve(’a*x^2+b*x+c=0”)③ eq=’a*x^2+b*x+c“;④ eq=”a*x^2+b*x+c=0’; seq=solve(eq)
⑤sym x a b c
eq=a*x^2+b*x+c seq=solve(eq)
2)常微分方程求解 求解常微分方程得函数就是 dsolve。应用此函数可以求得常微分方程(组)得通解,以及给定边界条件(或初始条件)后得特解。
常微分方程求解函数得调用格式:
r=dsolve(“eq1,eq2,…’,’cond1,cond2,…’,’v’)r=dsolve(’eq1’,’eq2”,…,“cond1’,’cond2’,…,’v”)说明: ① 以上两式均可给出方程 eq1,、qeq2 对应初始条件 cond1、cond2 之下得一 v 作为解变量得各微分方程得解。
② 常微分方程解得默认变量为 t。
③ 第二式中最多可接受得输入式就是 12 个。
④ 微分方程得表达方法。
在用 MATLAB 求解常微分方程时,用大写字母 Dy表示,用D2y表示,依此类推。
边界条件以类似于 y(a)=b 给出。其中 y 为因变量,a、b 为常数、如果初始条件给得不够,求出得解为含有 C1、C2 等待定常数得通解。
例一 求微分方程得通解、练习:(1)求。
(2)求函数得积分;求函数得导数(3)计算定积分(4)求下列线性方程组得解
(5)求解但 y(0)=2,在 z(0)=7 时,微分方程组得解。
实验五
MATLA B 程序设计
一、实验目得
掌握 MATLAB 程序设计得主要方法,熟练编写 MATLAB函数、二、实验内容
(1)M 文件得编辑。
(2)程序流程控制结构。
(3)子函数调用与参数传递。
(4)局部变量与全局变量。
三、实验步骤1、M M 文件得编辑
选择MATLAB 得菜单,打开新得M文件进行编辑,然后输入以下内容,并保存文件名为 expl、m。
s=0;for n=1:100
s=s+n;end s 答;s =
5050 保存好文件后,在命令窗口输入expl 即可运行该脚本文件,主义观察变量空间。紧接着创建 M 函数文件,然后输入以下内容,并保存文件名为 expl2、m。
function s=expl2(x)
s=0;for n=1:x
s=s+n; end
保存好文件后,在命令窗口输入
>〉 clear 〉> s=expl2(100)以 open 命令可以打开 M 文件进行修改。
〉〉open conv
%打开 conv 函数 2 2、程序流程控制结构
1)
for 循环结构
>〉 for n=1:10 n end 答;n =n =n =n =n =n =
6 n =n =
n =n =
另一种形式得 for 循环: 〉> n=10:—1:5;>〉 for i=n
%循环得次数为向量 n 得列数 i end 答;i =
10 i =i =i =i =i =2)
while 循环结构 在命令窗口输入:
clear x=1;while 1 x=x*2 end 将会瞧到 MATLAB 进入死循环因为 while 判断得值恒为真,这时须按 Ctrl+C键来中断运行,并且可瞧到 x 得值为无穷大。
练习:
(1)请把 exp2、m函数文件用 while 循环改写。
(2)用公式求得近似值,直到最后一项得绝对值小于为止,试编写M脚本文件、3)
if-else—end 分支结构
if—else—end 分支结构有如下 3 种形式。
(a)
if
表达式 语句组 1 end(b)
if
表达式 语句组 1 else 语句组2 end(c)
if
表达式 A
语句组 1 else if 表达式 B 语句组 2 else if 语句组 3 …… else
语句组 n end 4)
switch—case结构 创建 M 脚本文件 exp3、m,输入以下内容并在命令窗口中运行。
%功能:判断键盘输入得数就是奇数还就是偶数 n=input(’n=');if isempty(n)
error(’please input n');
n=input(“n=’);end switch mod(n,2)
case 1
A=’奇数'
case 0
A=’偶数” End 答;n=input('n=’); if isempty(n)
error('please input n’);
n=input(“n=’); end switch mod(n,2)
case 1
A=”奇数“
case 0
A=”偶数' end n=1 A = 奇数 3、子函数与参数传递
有一个函数,试编写实现该函数得函数文件.function g=expl4(x)
%主函数 g=0;for n=1:x
g=g+fact(n);
%调用子函数 end
function y=fact(k)
%子函数 y=1; for n=1:k
y=y*n; end 输入参数可以有函数 nargin 计算,下面得例子 sinplot(),当知输入一个参数 w 时,sinplot()函数会给p赋予默认值 0。4、局部变量与全局变量
一、实验目的与要求
1、掌握Java中类和接口的基本知识;
2、继承的基本知识;
3、集合的基本操作
二、实验内容 1.PIMCollection 类
创建一个可以管理PIMEntities实体的集合,该类实现了Collection接口,且必须实现了以下方法:(1).getNotes();返回PIMCollection中所有的PIMNote项,如果不包含PIMNote项,则返回一个空集合
(2).getTodos();返回集合PIMCollection中所有的PIMTodo项,如果不包含PIMTodo项,则返回一个空集合
(3).getAppointment();返回集合PIMCollection中所有的PIMAppointment项,如果不包含PIMAppointment项。则返回一个空集合(4).getContacts();返回结合PIMCollection中所有的PIMContact项,如果不包含PIMContact项,则返回一个空集合
(5).getItemsForDate(Date d);返回集合PIMCollection中所有与d匹配的PIMEntities,如果没有匹配d的项,则返回一个空集合
(只有PIMTodo和PIMAppointment项,因为PIMNote和PIMContact跟日期没有关系)。
2.TestPIMCollection类(自己的测试类)
向PIMCollection集合类中加入了一些PIMEntity实体类,并用PIMCollection类中实现的方法分别筛选出了其中的PIMNote、PIMTodo、PIMAppointment、PIMContact、符合指定日期d的集合并打印输出。
三、实验器材
计算机+windows操作系统+eclipse
四、实验步骤
1.程序源代码(压缩在文件夹里)2.编译执行程序 3.记录执行结果
五、实验结果
1.TestPIMCollection的测试代码:
2.程序的执行结果:
如图所示:程序分别输出了Todo、Note、Appointment、Contact、以及匹配指定日期的结合的结果。
六、实验小结
1、熟悉了Java集合的基本内容和操作,也包括了泛型的一些知识。
2、通过这个实验巩固了java的基础知识:类和接口的有关知识,继承的应用。
3、掌握了在eclipse下编译执行Java程序的知识。
4、通过几次的编程,熟悉了java的编程规则。
---中海2000物流信息系统
班级:XXX班
姓名:XXX
学号:XXX号
物流信息系统上机实验报告
-----中海2000物流信息系统
班级:XX班学号:XX3号姓名:XXX
通过两天的上机实践学习,我了解并熟悉了中海物流信息系统的各个流程及其管理。中海2000物流管理信息系统是在总结中海物流集团多年从事第三方物流服务的基础上,经过多种物流业务模式考验的现代物流企业综合物流业务管理信息系统。该信息系统共有基础数据、客户管理、配送管理、仓储管理、调度管理、费用产生、统计调查、客户分析几大功能模块。接下来,我就结合这两天的实践学习,阐述一下我对该信息系统的认识及对物流信息系统管理的体会。首先,我简单介绍一下本系统的几大功能模块:
一、基础数据
系统的运行要根据使用的具体环境对系统进行具体设置。基础数据部分即是对操作人员及密码,系统权限的设置。
每一个操作员可以灵活的设置自己的密码,保证用户的合法性。对于系统的使用人员给予一定的权限,明确每个操作人员的权限和职责,避免与业务无关的人员对系统进行非法操作,严格控制系统用户的合法性,保证系统的安全。
二、客户管理
该模块主要管理客户的各种相关信息及为其提供相应的查询服务。它包括客户资料、合同管理、工厂资料、运输资料四部分。
客户资料分为两个等级(1)客户本身的资料:客户名称、联系方式、收费状态、报价资料等,由商务部门维护(2)客户下的工厂资料:供应商、购买商、供应购买商,由操作部门维护。
三、配送管理
该模块分为:部件维护、入仓、出仓、查询。
(一)入仓
1、部件维护为客户输入部件信息,并进行管理。
2、入仓资料
(1)入仓资料客户的供应商要求部件入仓,输入入仓件数及具体信息。
(2)入仓配车为客户选择运输方式及线路。
(3)入仓报关为客户确认入仓部件的报关要求。
(4)入仓卸车为客户选择卸车方式并卸车。
(5)入仓验货经验收,确保货物全部符合客户要求。
(6)分配仓位仓管员把部件存放在其包租的仓位中。
(7)报关确认
(8)入仓确认操作员核对所有入仓资料无误后,进行入仓确认。
3、费用处理客户部件入仓时产生处理费用,加班费用及其他费用等。
(二)出仓
1、出仓选货客户选择要购买的部件,操作员提出选货单给仓管员进行选货操作。
2、出仓资料
(1)出仓资料 操作员准备出仓资料,作为仓管员出货检查的依据。
(2)出仓配车 为客户选择运输方式及线路。
(3)出仓报关为客户确认出仓部件的报关要求。
(4)出仓装卸 选择卸货方式
(5)报关确认
(6)出仓确认 核对所有出仓资料,无误后进行出仓确认。
3、费用处理出仓时产生处理费用,加班费用,及其他费用等。
四、仓储管理
分别进行入仓部件查询,出仓部件查询,货物编号查询,货物名称查询,呆滞货物查询,进出库查询,仓位库存查询。
五、调度管理
根据部件情况完成相应的调度配载工作。
六、费用产生
物流计费项目主要包括:仓储保管费,流通加工费,运输配送费。仓储保管费按体积、重量、箱数计收;流通加工费以所耗人力工时及材料为基准计收,而运输配送费则以紧急度,车辆,货品数,装卸次数,运输配送区域,运行管理等为计费依据。计费标准的确定必须考虑市场的接受度和运营成本。
其次,我总结一下该系统的特点如下:
中海物流信息系统属于第三方物流企业的物流管理信息系统。第三方物流企业服务于生产企业与流通企业以及消费者,以提供第三方物流服务为主。主要有可得性,及时性,精确性,灵活性等特点。
再次,是我通过上机实践获得的体会。有以上各个方面的了解,不难发现好的物流信息系统可以降低成本,降低风险,提高服务水平,为决策提供科学依据。使用特定的信息管理系统,进行信息的收集,加工,传递,反馈,从而极大地提高劳动效率,提高经济效益。
物流信息系统首先是为物流活动提供了一个作业平台,所有的作业活动成果都要通过该系统来反映,为其他作业和决策提供依据;同时由于该系统是一个定向开放的系统,所以作业的信息会在一定的范围内进行传递和沟通。物流信息系统也是一个决策支持系统,它具有数据汇总查询功能,以便决策者及时掌握各种物流数据信息,为决策提供依据。物流信息系统一般会记录每一订单的订单号,货物名称,下单日期等各种有用的信息,也可以提高信息的跟踪能力。最后,得出结论如下:
【第二次上机实验报告】推荐阅读:
上机实验报告8函数11-14
数据处理上机实验报告06-25
电算化上机实验报告11-15
工程热力学上机实验报告07-15
第二次上机实习作业04-24
《c语言程序设计》上机实验报告要求05-03
java上机实验考试12-23
数据库上机实验题06-30
网络数据库上机实验01-26
数据结构上机实验题02-15