(431) 北京理工大学远程教育学院2019-2020年第一学期
《DSP原理及应用》期末大作业
教学站学号姓名成绩
C6000实验1：汇编语言、体系结构和CCS
◆采用simulator配置文件sim6201_simulator.cfg
◆ 1. 新建一个project，把tutorial\sim62xx\hello1的文件添加进去，完成其功能。练
习CCS的基本操作：加载程序、go main、执行程序、设置断点、单步执行、观察变量、观察存储器、测试函数执行的CLK、混合代码显示。
◆ 2. 打开tutorial\sim62xx\volume1的项目文件，完成图形方式观察变量、设置探点、
GEL控制变量、FILE IO、动画显示输入输出的功能。然后单步执行load.asm：如何被C代码代码调用、观察寄存器的变化、测试函数执行时间、如何循环和返回C代码。
◆ 3. 解释如何在C6201上实现32bit int乘32bit int, 结果是32 bit int的过程：在刚才
的hello1中设置3个全局变量int a = 0x10008; int b = 56; int c;在main函数中执行语句 c = a * b; build后加载执行程序, 混合代码显示c = a*b对应的汇编代码，解释之。
◆ 4. Hello1中添加1个C文件sop_c.c ，其中有一个矢量乘的子函数sop(short * a, short
* x, int * y, int n)，然后在main函数中调用。不选择任何优化选项。用混合代码显示，在汇编层次执行sop子函数，观察调用、执行和返回的过程。
C6000实验2：C6000流水线和C运行时环境
◆采用simulator配置文件sim6201_simulator.cfg
◆ 1. 用汇编语言实现两个数组的乘法累加功能：sop_asm.asm中实现函数int
sop_asm(short * a, short * x, int n)；然后，手工优化汇编代码；
◆ 2. C调用汇编函数，汇编函数调用c函数：
在main函数中调用汇编实现的函数addarr3, 在addarr3中调用C子函数addarr2.
addarr3(int * arr1, int * arr2, int * arr3, int * arr4, int n) 是汇编函数，其功能为3个数组的对应位置相加，结果放在arr4[n]中；
汇编函数addarr3中调用C的子函数，它把arr1和arr2相加放到arr3[n]中；
addarr2(int * arr1, int * arr2, int * arr3, int n)
C6000实验3：DSP/BIOS
◆采用simulator配置文件sim6201_simulator.cfg
1. 在volume2的基础上添加一个中断4，并写中断4的服务函数，每发生一次中断把计数器加1。simulator 中，利用pin connect仿真中断；
2. 在实验1的基础上添加一个任务TSK1，优先级为3，其服务函数采用for(){}的结构，该任务每3ms执行一次。
3. 在实验2的基础上添加一个idle线程，观察idle线程的函数每ms执行多少次；
4. 在实验2的基础上添加一个任务TSK2，优先级为4，其服务函数采用for(){}的结构，该任务每2ms执行一次。在第6ms、12ms、18ms….等时刻时，分析TSK1和TSK2的任务切换过程。
C6000实验4：C6000代码优化
1. 采用simulator配置文件sim6201_simulator.cfg
写手工优化的.asm程序：在volume1的load.asm基础上实现16bit数组的乘法累加的函数，并手工优化和采用软件流水优化；要求：
(1)使用LDW和MPY,MPYH指令；对于C64可以使用扩展乘法指令。
(2)画出相关图和模迭代间隔表
2. 采用simulator配置文件sim6701_simulator.cfg .
写一个线形汇编实现的子函数，做复数数组求模，并找出模的最大值的位置和值。复数的实部和虚部为short，求模的结果为unsigned int；
提示：使用INTSP指令把整数转化为浮点数后再开方；
一． 1. 1/e 2. 3 3.1 4.e3 5. ∞ 6. 0 7.∞ 8.0 9.1/2 10.0 11.e2c12.不存在13. 1/12 Matlab实验过程： 1.1/exp(1) syms n; f=(1-1/n)^n; limit(f,n,inf) ans = 1/exp(1) 2.3 syms n; f=(n^3+3^n)^(1/n); limit(f,n,inf) ans = 3 3. 1 syms n; f=(1+sin(2*n))/(1-cos(4*n)); limit(f,n,pi/4) ans = 1 4.e^3 syms x; f=(1+cos(x))^(3*sec(x)); limit(f,x,pi/2) ans = exp(3) 5.inf syms x; f=(x^2)*exp(1/(x^2));
limit(f,x,0) ans = Inf 6.0 syms x; f=(x^2-2*x+1)/(x^3-x); limit(f,x,1) ans = 7.inf syms x; f=((2/pi)*atan(x))^x; limit(f,x,+inf) ans = Inf 8.0 syms x y; f=(1-cos(x^2+y^2))/((x^2+y^2)*exp(x^2+y^2)); limit(limit(f,x,0),y,0) ans = 9.1/2 syms x; f=(1-cos(x))/(x*sin(x)); limit(f,x,0) ans = 1/2 10.0 syms x;
f=atan(x)/(2*x); limit(f,x,inf) ans = 11.exp(2*c) syms c; f=sym('((x+c)/(x-c))^x'); limit(f,'x',inf) ans = exp(2*c) 12.极限不存在 syms x; f=cos(1/x); limit(f,x,0) ans = limit(cos(1/x), x = 0) 13.1/12 syms x; f=1/(x*log(x)^2)-1/(x-1)^2; limit(f,x,1) ans = 1/12 二．观察函数logbx,当b=1/2,1/3,1/4和b=2,3,4时函数的变化特点，总结logbx的图形特点。
学生实验报告 （理工类） 课程名称：信号与线性系统专业班级：M11通信工程 学生学号：1121413017 学生姓名：王金龙 所属院部：龙蟠学院指导教师：杨娟
20 11 ——20 12 学年第 1 学期 金陵科技学院教务处制 实验报告书写要求 实验报告原则上要求学生手写，要求书写工整。若因课程特点需打印的，要遵照以下字体、字号、间距等的具体要求。纸张一律采用A4的纸张。 实验报告书写说明 实验报告中一至四项内容为必填项，包括实验目的和要求；实验仪器和设备；实验内容与过程；实验结果与分析。各院部可根据学科特点和实验具体要求增加项目。 填写注意事项 （1）细致观察，及时、准确、如实记录。 （2）准确说明，层次清晰。 （3）尽量采用专用术语来说明事物。 （4）外文、符号、公式要准确，应使用统一规定的名词和符号。 （5）应独立完成实验报告的书写，严禁抄袭、复印，一经发现，以零分论处。 实验报告批改说明 实验报告的批改要及时、认真、仔细，一律用红色笔批改。实验报告的批改成绩采用百分制，具体评分标准由各院部自行制定。 实验报告装订要求
实验批改完毕后，任课老师将每门课程的每个实验项目的实验报告以自然班为单位、按学号升序排列，装订成册，并附上一份该门课程的实验大纲。
实验项目名称：常用连续信号的表示 实验学时： 2学时 同组学生姓名： 无 实验地点： A207 实验日期： 11.12.6 实验成绩： 批改教师： 杨娟 批改时间： 一、实验目的和要求 熟悉MATLAB 软件；利用MATLAB 软件，绘制出常用的连续时间信号。 二、实验仪器和设备 586以上计算机，装有MATLAB7.0软件 三、实验过程 1. 绘制正弦信号)t Asin t (f 0?ω+=（），其中A=1，πω2=，6/π?=； 2. 绘制指数信号at Ae t (f =），其中A=1，0.4a -=； 3. 绘制矩形脉冲信号，脉冲宽度为2； 4. 绘制三角波脉冲信号，脉冲宽度为4；斜度为0.5； 5. 对上题三角波脉冲信号进行尺度变换，分别得出)2t (f ，)2t 2(f -； 6. 绘制抽样函数Sa （t ）,t 取值在-3π到+3π之间； 7. 绘制周期矩形脉冲信号，参数自定； 8. 绘制周期三角脉冲信号，参数自定。 四、实验结果与分析 1．制正弦信号)t Asin t (f 0?ω+=（），其中A=1，πω2=，6/π?= 实验代码： A=1;
数学实验上机作业整理∈hyd 实验一 1. 计算球体体积（半径r=5） r=5;v=(4/3)*pi*r^3 v =523.5988 2.设矩阵1234567891023416A ?? ? = ? ??? （1）提取A 的第2列赋值给B; A=[1 2 3 4 5;6 7 8 9 10;2 3 4 1 6];B=A(:,2) B = 2 7 3 （2）提取A 的第2行前3个数给C ； A=[1 2 3 4 5;6 7 8 9 10;2 3 4 1 6];C=A(2,[1,2,3]) C = 6 7 8 （3）提取A 第1,3行和2, 4列相交位置元素构成子矩阵D ； A=[1 2 3 4 5;6 7 8 9 10;2 3 4 1 6];D=A([1,3],[2,4]) D = 2 4 3 1 （4）构造矩阵E 使得E 的结构为：132213C E D C ???? ?= ? ?? A=[1 2 3 4 5;6 7 8 9 10;2 3 4 1 6];E=[D [C;C]] E = 2 4 6 7 8 3 1 6 7 8
（5）把A 中间的8换为0； A(2,3)=0;A A = 1 2 3 4 5 6 7 0 9 10 2 3 4 1 6 （6）去掉A 的第2行； A=[1 2 3 4 5;6 7 8 9 10;2 3 4 1 6]; A(2,:)=[] A = 1 2 3 4 5 2 3 4 1 6 3.写出完成下列操作的命令 （1） 建立10阶单位矩阵A; A=eye(10) （2）建立5×6的随机矩阵A ，其元素为[100，200]范围内的随机数； A=rand(5,6)*100+100 （3）将A 对角线元素加30 A+eye(5,6)*30 4.（选做题）设有分块矩阵333223E R A O S ????? =? ??? ，其中E,R,O,S 分别为单位矩阵、随机矩阵、零矩阵和对角矩阵，试通过数值计算验证2 2 E R RS A O S +?? =? ??? 。 S=[1 1;1 1]; E=eye(3);R=rand(3,2); O=zeros(2,3); [E R;O S]^2 [E R+R*S;O S^2] 实验二 1.设矩阵1215346562A -?? ? = ? ?-?? （1）求A 的秩、A 的每个元素3次方; A=[1 2 -1;5 34 6;-5 6 2];
北京邮电大学微机原理硬件实验报告
实验报告一：I/0地址译码和简单并行接口 ——实验一&实验二 一、实验目的 掌握I/O地址译码电路的工作原理；掌握简单并行接口的工作原理及使用方法。 二、实验原理及内容 a) I/0地址译码 1、实验电路如图1-1所示，其中74LS74为D触发器，可直接使用实验台上数 字电路实验区的D触发器,74LS138为地址译码器。译码输出端Y0～Y7在实验台上“I/O地址“输出端引出，每个输出端包含8个地址，Y0：280H～ 287H，Y1：288H～28FH，……当CPU执行I/O指令且地址在280H～2BFH范围内，译码器选中，必有一根译码线输出负脉冲。 例如：执行下面两条指令 MOV DX，2A0H OUT DX，AL（或IN AL，DX） Y4输出一个负脉冲，执行下面两条指令 MOV DX，2A8H OUT DX，AL（或IN AL，DX） Y5输出一个负脉冲。 利用这个负脉冲控制L7闪烁发光（亮、灭、亮、灭、……），时间间隔经过软件延时实现。 2、接线： Y4/IO地址接 CLK/D触发器
Y5/IO地址接 CD/D触发器 D/D触发器接 SD/D触发器接 +5V Q/D触发器接L7（LED灯）或逻辑笔 b) 简单并行接口 1、按下面图4-2-1简单并行输出接口电路图连接线路(74LS273插通 用插座，74LS32用实验台上的“或门”)。74LS273为八D触发器， 8个D输入端分别接数据总线D0～D7，8个Q输出端接LED显示电 路L0～L7。 2、编程从键盘输入一个字符或数字，将其ASCⅡ码经过这 个输出接口输出，根据8个发光二极管发光情况验证正确 性。 3、按下面图4-2-2简单并行输入接口电路图连接电路 (74LS244插通用插座，74LS32用实验台上的“或门”)。 74LS244为八缓冲器，8个数据输入端分别接逻辑电平开关 输出K0～K7，8个数据输出端分别接数据总线D0～D7。 4、用逻辑电平开关预置某个字母的ASCⅡ码，编程输入这 个ASCⅡ码，并将其对应字母在屏幕上显示出来。 5、接线：1）输出 按图4-2-1接线（图中虚线为实验所需接线，74LS32为实验 台逻辑或门） 2）输入 按图4-2-2接线（图中虚线为实验所需接线，74LS32为实 验台逻辑或门） 三、硬件连线图 1、I/O地址译码
练习2﹒1 画出下列常见曲线的图形（其中a=1，b=2,c=3）。 1. 立方抛物线y = 解： x=-4:0.1:4; y=x.^(1/3); plot(x,y) -4 -3-2-101234 0.20.40.60.811.21.4 1.6 2.高斯曲线2 x y e -= 解： fplot('exp(-x^2)',[-4,4])
-4 -3 -2 -1 1 2 3 4 00.10.20.30.40.50.60.70.80.9 1 3、笛卡儿曲线23 3 2 2 33,(3)11at at x y x y axy t t = = +=++ 解：ezplot('x^3+y^3-3*x*y',[-4,4])
-4 -3-2-1 01234 -4-3-2-10123 4x y x 3+y 3-3 x y = 0 或：t=-4:0.1:4; x=3*t./(1+t.^2); y=3*t.^2./(1+t.^2); plot(x,y)
-1.5 -1-0.500.51 1.5 00.5 1 1.5 2 2.5 3 4、蔓叶线233 2 2 2 ,()11at at x x y y t t a x = = = ++- 解：t=-4:0.1:4; x=t.^2./(1+t.^2); y=t.^3,/(1+t.^2); y=t.^3./(1+t.^2); plot(x,y)
00.10.20.30.40.50.60.70.80.91 -4 -3-2-10123 4 或： ezplot('y .^2-x.^3/(1-x)',[-4,4])
本科实验报告 实验名称：实验三使用8251A的串行接口应用实验 课程名称：计算机原理与应用实验实验时间： 任课教师：实验地点： 实验教师： 实验类型：□原理验证■综合设计□自主创新 学生姓名： 学号/班级：组号：学院：同组搭档：专业：成绩：
1. 实验目的 1) 掌握串行通信原理及半双工和全双工的编程方法； 2) 掌握用8251A接口芯片实现微机间的同步和异步通信； 3) 掌握8251A芯片与微机的接口技术和编程方法。 2. 实验原理和内容 8251A是一种可编程的同步/异步串行通信接口芯片，具有独立的接收器和发送器，能实现单工、半双工、双工通信。 1) 8251A内部结构 8251A通过引脚D0~D7和系统数据总线直接接口，用于和CPU传递命令、数据、状态信息。读写控制逻辑用来接收CPU的控制信号、控制数据传送方向。CPU对8251A的读写操作控制表如表3-4所示。 表3-4 CPU对8251A的读写操作控制表 2) 8251A的方式控制字和命令控制字 方式控制字确定8251A的通信方式（同步/异步）、校验方式（奇校/偶校/不校）、字符长度及波特率等，格式如图3-10所示。 命令控制字使8251A处于规定的状态以准备收发数据，格式如图3-11所示。 方式控制字和命令控制字无独立的端口地址，8251A 根据写入的次序来区分。 CPU对8251A初始化时先写方式控制字，后写命令控制字。
3) 状态寄存器 8251状态寄存器用于寄存8251A的状态信息，供CPU查询，定义如图3-12所示。TXRDY位：当数据缓冲器空时置位，而TXRDY引脚只有当条件( 数据缓冲器空?/CTS?TXE)成立时才置位。 溢出错误：CPU没读走前一个字符，下一个字符又接收到，称为溢出错误。
实验1 信号的时域描述与运算 一、实验目的 1. 掌握信号的MATLAB表示及其可视化方法。 2. 掌握信号基本时域运算的MA TLAB实现方法。 3. 利用MA TLAB分析常用信号，加深对信号时域特性的理解。 二、实验原理与方法 1. 连续时间信号的MATLAB表示 连续时间信号指的是在连续时间范围内有定义的信号，即除了若干个不连续点外，在任何时刻信号都有定义。在MATLAB中连续时间信号可以用两种方法来表示，即向量表示法和符号对象表示法。 从严格意义上来说，MATLAB并不能处理连续时间信号，在MATLAB中连续时间信号是用等时间间隔采样后的采样值来近似表示的，当采样间隔足够小时，这些采样值就可以很好地近似表示出连续时间信号，这种表示方法称为向量表示法。表示一个连续时间信号需要使用两个向量，其中一个向量用于表示信号的时间范围，另一个向量表示连续时间信号在该时间范围内的采样值。例如一个正弦信号可以表示如下： >> t=0:0.01:10; >> x=sin(t); 利用plot(t,x)命令可以绘制上述信号的时域波形，如图1所示。 如果连续时间信号可以用表达式来描述，则还可以采用符号表达式來表示信号。例如对于上述正弦信号，可以用符号对象表示如下： >> x=sin(t); >> ezplot(X); 利用ezplot(x)命令可以绘制上述信号的时域波形 Time(seconds) 图1 利用向量表示连续时间信号
t 图 2 利用符号对象表示连续时间信号 sin(t) 2.连续时间信号的时域运算 对连续时间信号的运算包括两信号相加、相乘、微分、积分，以及位移、反转、尺度变换（尺度伸缩）等。 1）相加和相乘 信号相加和相乘指两信号对应时刻的值相加和相乘，对于两个采用向量表示的可以直接使用算术运算的运算符“+”和“*”来计算，此时要求表示两信号的向量时间范围和采样间隔相同。采用符号对象表示的两个信号，可以直接根据符号对象的运算规则运算。 2）微分和积分 对于向量表示法表示的连续时间信号，可以通过数值计算的方法计算信号的微分和积分。这里微分使用差分来近似求取的，由时间向量[N t t t ,,,21?]和采样值向量[N x x x ,,,21?]表示的连续时间信号，其微分可以通过下式求得 1,,2,1,|)('1-?=?-≈ +=N k t x x t x k k t t k 其中t ?表示采样间隔。MA TLAB 中用diff 函数来计算差分 k k x x -+1。 连续时间信号的定积分可以由MATLAB 的qud 函数实现，调用格式为 quad ('function_name',a,b) 其中，function_name 为被积函数名，a 、b 为积分区间。
数学建模作业——实验1 学院：软件学院 姓名： 学号： 班级：软件工程2015级 GCT班 邮箱： 电话： 日期：2016年5月10日
基本实验 1.椅子放平问题 依照1.2.1节中的“椅子问题”的方法，将假设中的“四腿长相同并且四脚连线呈正方形”，改为“四腿长相同并且四脚连线呈长方形”，其余假设不变，问椅子还能放平吗？如果能，请证明；如果不能，请举出相应的例子。 答：能放平，证明如下： 如上图，以椅子的中心点建立坐标，O为原点，A、B、C、D为椅子四脚的初始位置，通过旋转椅子到A’、B’、C’、D’，旋转的角度为α，记A、B两脚，C、D两脚距离地面的距离为f(α)和g(α)，由于椅子的四脚在任何位置至少有3脚着地，且f(α)、g(α)是α的连续函数，则f(α)和g(α)至少有一个的值为0，即f(α)g(α)=0，f(α)≥0，g(α)≥0，若f(0)＞0，g(0)=0，
则一定存在α’∈（0，π），使得 f(α’)=g(α’)=0 令α=π（即椅子旋转180°,AB 边与CD 边互换），则 f(π)=0，g(π)＞0 定义h(α)= f(α)-g(α)，得到 h(0)=f(0)-g(0)＞0 h(π)=f(π)-g(π) ＜0 根据连续函数的零点定理，则存在α’∈（ 0，π），使得 h(α’)= f(α’)-g(α’)=0 结合条件f(α’)g(α’)=0,从而得到 f(α’)=g(α’)=0，即四脚着地，椅子放平。 2. 过河问题 依照1.2.2节中的“商人安全过河”的方法，完成下面的智力游戏：人带着猫、鸡、米过河，船除需要人划之外，至多能载猫、鸡、米之一，而当人不在场时，猫要吃鸡、鸡要吃米，试设计一个安全过河的方案，并使渡河的次数尽量的少。 答： 用i =1,2,3,4分别代表人，猫，鸡，米。1=i x 在此岸，0 =i x 在对岸， ()4321,,,x x x x s =此岸状态，()43211,1,1,1x x x x D ----=对岸状态。安全状态集合为 :
微机原理 实 验 报 告 隋伟 08212013 自动化0801
目录 一、I/O地址译码与交通灯控制实验 (3) 二、可编程定时器/计数器（8253） (6) 三、中断实验（纯DOS） (11) 四、模/数转换器………………………………………………………… 18 五、串行通讯…………………………………………………………… 16 六、课程综合实验(抢答器) …………………………………………… 28 七、自主设计实验——LED显示 (32) 八、参考文献 (35)
一、I/O地址译码与交通灯控制实验 一.实验目的 通过并行接口8255实现十字路口交通灯的模拟控制,进一步掌握对并行口的使用。 二.实验内容 如图5-3，L7、L6、L5作为南北路口的交通灯与PC7、PC6、PC5相连，L2、L1、L0作为东西路口的交通灯与PC2、PC1、PC0相连。编程使六个灯按交通灯变化规律燃灭。 十字路口交通灯的变化规律要求： （1）南北路口的绿灯、东西路口的红灯同时亮3秒左右。 （2）南北路口的黄灯闪烁若干次，同时东西路口的红灯继续亮。 （3）南北路口的红灯、东西路口的绿灯同时亮3秒左右。 （4）南北路口的红灯继续亮、同时东西路口的黄灯亮闪烁若干次。 （5）转（1）重复。 8255动态分配地址: 控制寄存器： 0C40BH A口地址： 0C408H C口地址： 0C40AH
三．程序流程图和程序清单 DATA SEGMENT X DB ? DATA ENDS STACK1 SEGMENT STACK DW 100H DUP(0) STACK1 ENDS CODE SEGMENT ASSUME CS:CODE,DS:DATA,SS:STACK1 START: MOV AX,DATA MOV DS,AX ;---------------INIT---------------- MOV DX,0C40BH ;写控制字
本科实验报告实验名称：随机信号分析实验
实验一 随机序列的产生及数字特征估计 一、实验目的 1、学习和掌握随机数的产生方法。 2、实现随机序列的数字特征估计。 二、实验原理 1、随机数的产生 随机数指的是各种不同分布随机变量的抽样序列（样本值序列）。进行随机信号仿真分析时，需要模拟产生各种分布的随机数。 在计算机仿真时，通常利用数学方法产生随机数，这种随机数称为伪随机数。伪随机数是按照一定的计算公式产生的，这个公式称为随机数发生器。伪随机数本质上不是随机的，而且存在周期性，但是如果计算公式选择适当，所产生的数据看似随机的，与真正的随机数具有相近的统计特性，可以作为随机数使用。 (0，1)均匀分布随机数是最最基本、最简单的随机数。(0，1)均匀分布指的是在[0，1]区间上的均匀分布，即 U(0，1)。实际应用中有许多现成的随机数发生器可以用于产生(0，1)均匀分布随机数，通常采用的方法为线性同余法，公式如下： )(mod ,110N ky y y n n -= N y x n n /= 序列{}n x 为产生的(0，1)均匀分布随机数。 下面给出了上式的3组常用参数： 1、10N 10,k 7==，周期7 510≈?； 2、（IBM 随机数发生器）3116N 2,k 23,==+周期8 510≈?； 3、（ran0）31 5 N 21,k 7,=-=周期9 210≈?； 由均匀分布随机数，可以利用反函数构造出任意分布的随机数。 定理 1.1 若随机变量 X 具有连续分布函数F X (x)，而R 为(0,1)均匀分布随机变量，则有 )(1R F X x -= 由这一定理可知，分布函数为F X (x)的随机数可以由(0,1)均匀分布随机数按上式进行变
实验2 LTI 系统的时域分析 （基础型实验） 一． 实验目的 1. 掌握利用MATLAB 对系统进行时域分析的方法。 2. 掌握连续时间系统零状态响应、冲击响应和阶跃响应的求解方法。 3. 掌握求解离散时间系统响应、单位抽样响应的方法。 4. 加深对卷积积分和卷积和的理解。掌握利用计算机进行卷积积分和卷积和计算的方法。 二． 实验原理与方法 1. 连续时间系统时域分析的MATLAB 实现 1） 连续时间系统的MA TLAB 表示 LTI 连续系统通常可以由系统微分方程描述，设描述系统的微分方程为： (N)(N 1)(M)(M 1)1010(t)(t)...(t)b (t)b (t)...b (t)N N M M a y a y a y x x x ----++=++ 则在MATLAB 中可以建立系统模型如下： 1010[b ,b ,...,b ];a [a ,a ,...,a ];sys tf(b,a); M M N N b --=== 其中，tf 是用于创建系统模型的函数，向量a 和b 的元素是以微分方程求导的降幂次序来排列的，如果有缺项，应用0补齐，例如由微分方程 2''(t)y'(t)3y(t)x(t)y ++= 描述的系统可以表示为： >> b=[1]; >> a=[2 1 3]; >> sys=tf(b,a); 而微分方程由 ''(t)y'(t)y(t)x''(t)x(t)y ++=- 描述的系统则要表示成 >> b=[1 0 -1]; >> a=[1 1 1]; >> sys=tf(b,a); 2） 连续时间系统的零状态响应 零状态响应指系统的初始状态为零，仅由初始信号所引起的响应。MATLAB 提供了一个用于求解零状态响应的函数lism ，其调用格式如下： lism （sys,x,t ）绘出输入信号及响应的波形，x 和t 表示输入信号数值向量及其时间向量。 y= lism （sys,x,t ）这种调用格式不绘出波形，而是返回响应的数值向量。 3） 连续时间系统的冲激响应与阶跃响应
微机原理与接口技术 硬件实验报告 班级: 姓名: 学号: 实验一: 8259中断控制器实验 一、实验目的 1. 掌握8259的工作原理。 2. 掌握编写中断服务程序的方法。 3. 掌握初始化中断向量的方法。 二、实验内容 用单脉冲发生器的输出脉冲为中断源，每按一次产生一次中断申请，点亮或熄灭发光二极管。 三、实验设备 微机实验教学系统实验箱、8086CPU模块 四、连线 ①单脉冲发生器输出P+与8259的IR0相连 ②8259的片选CS8259与CS0相连 ③8259的INT与8086的INT相连 ④8259的INTA与8086的INTA相连 ⑤CS273与CS1相连 ⑥00与LED1相连 其它线均已连好如下图： 五、实验步骤 （1）连线。 （2）编辑程序，编译链接后，调试程序。 （3）调试通过后，在中断服务程序内设置断点，运行程序，当接收到中断请求后，程序停在中断服务程序内的断点处。 （4）撰写实验报告。
六、实验源程序 CODE SEGMENT PUBLIC ASSUME CS:CODE ORG 100H START: MOV DX,4A0H ;写ICW1 MOV AX,13H OUT DX,AX MOV DX,4A2H ;写ICW2 MOV AX,80H ;IR0的中断向量码为80H OUT DX,AX MOV AX,01 OUT DX,AX ;一般嵌套，非缓冲，非自动EOI MOV AX,0 ;写OCW1 OUT DX,AX ;允许中断 ;中断向量存放在（0000H：0200H）开始的四个单元里 MOV AX,0 MOV DS,AX MOV SI,200H ;中断类型号为80H MOV AX,OFFSET HINT ;中断服务程序的入口地址 MOV DS:[SI],AX ADD SI,2 MOV AX,CS MOV DS:[SI],AX STI ;开中断，设置IF=1 JMP $ ;原地跳转 HINT: ;中断服务程序 XOR CX,0FFH ;CX取反 MOV DX,4B0H ;CS273接口的地址，与8个LED灯相连 MOV AX,CX ;输出高低电平控制LED灯的亮灭 OUT DX,AX MOV DX,4A0H ;OCW2的地址 MOV AX,20H ;一般EOI命令，全嵌套方式 OUT DX,AX IRET ;中断返回 CODE ENDS END START 七、实验思考题 1．将P+连线连接到IR1—IR7任意一个；重新编写程序。 将P+接到IR1，在原程序的基础上，把写ICW2的控制字改为81H，再把中断向量的入口地址改为0204H即可。程序如下： CODE SEGMENT PUBLIC ASSUME CS:CODE
实验2 LTI 系统的时域分析 一、实验目的 1.掌握利用MATLAB 对系统进行时域分析的方法。 2.掌握连续时间系统零状态响应、冲激响应和阶跃响应的求解方法。 3.掌握求解离散时间系统响应、单位抽样响应的方法。 4.加深对卷积积分和卷积和的理解。掌握利用计算机进行卷积积分和卷积和计算的方法。 二、实验原理 1、连续时间系统时域分析的MATLAB 实现 1) 连续时间系统的MATLAB 表示 设LTI 因果系统的微分方程一般式为： ) ()()()()()()()(0'1)1(1)(0'1)1(1)(t x b t x b t x b t x b t y a t y a t y a t y a m m m m n n n n ++++=++++---- 则在MATLAB 里，可以建立系统模型如下： b =[b M ,b M?1,…,b 0]; a =[a N ,a N?1,…,a 0] sys =tf (b ,a ); 2）连续时间系统的零状态响应 零状态响应指系统的初始状态为0，仅由输入信号所引起的响应。MATLAB 提供了lsim(sys,x,t)表示求解零状态响应。他不绘出波形，而是返回响应二等地数值向量。 3）连续时间系统的冲激响应与阶跃响应。 MATLAB 提供了impulse 函数来求指定时间范围内，由模型sys 描述的连续时间系统的单位冲击响应。 函数step 用于求解单位阶跃响应。 2、离散时间系统时域分析的MATLAB 实现 1）离散时间系统的MATLAB 表示。 LTI 离散系统通常可以由系统差分方程描述;差分方程为: 0101()(1)...()()(1)...() N M a y n a y n a y n M b x n b x n b x n N +-++-=+-++- 则在MATLAB 里，可以建立系统模型如下： b =[b M ,b M?1,…,b 0]; a =[a N ,a N?1,…,a 0] 2）离散时间系统对任意输入的响应。 Matlab 提供了求LTI 离散系统响应的专用函数filter ，该函数用于求取由差分方程所描述的离散时间系统在指定时间范围内对输入序列所产生的响应。其基本调用格式为： y =filter (b ,a ,x ) 3）离散时间系统的单位抽样响应。 Matlab 提供impz 函数来求解指定时间范围内，由向量b 和a 描述的离散时间系统的单位抽样响应。 3、卷积和与卷积积分 1）离散时间序列的卷积和：
应用统计学 一、单选题共 1某商场2007与2006年相比,商品销售额上涨了16%,销售量增长了18%,则销售价格增减变动的百分比为( )。 B-1.7% 2某公司所属两个分厂,今年与去年相比,由于两个分厂单位新产品成本降低而使公司的总平均成本下降了5%,由于新产品结构的变化使公司总平均成本提高了10%,该公司总平均成本增减变动的百分比为（）。 D4.5% 3某企业的产品产量.产品库存量（）。 D当前是时期数，后者是时点 4变量X与Y之间的负相关是指（）。 C大的X值趋于同小的Y值相关联，小的X值趋于同大的Y值相关联 5物价上涨后,同样多的人民币只能购买原有商品的90%,则物价指数为( )。 B1.11% 6已知从总体中抽取一个容量为10的样本，样本均值的方差等于55，则总体方差等于（）。B550 7用趋势剔除法测定季节变动时，（）。 A不包括长期趋势的影响 8某企业计划规定单位成本降低8%，实际降低了5%，则成本计划完成程度为（）。 C62.5% 9我国1949年末总人口为54167万人,1989年末为111191万人,计算1949-1989年人口平均增长速度开( )次方。 A40 10抽样调查与典型调查的主要区别是（）。 C选取调查单位的方法不同 11对几个特大型商场进行调查，借以了解北京市商业市场商品销售的基本情况。这种调查方式属于（）。 B重点调查 12设产品产量与产品单位成本之间的简单相关系数为-0.78,这说明二者之间存在着( )。 B中度相关
13下列数据属于结构相对数的是（）。 B工业产值占工农业总产值的比重 14某种商品的价格今年比去年上涨了5%，销售额下降了10%，该商品销售量下降的百分比为( ) 。 B14.29% 15某百货公司今年与去年相比,商品销售量增长了10%,零售价格平均下降了10%,则商品零售额（）。 C下降 16某企业产品产量比上年提高了10% ，总成本下降了5% ，则单位成本降低了（）。 C13.64% 17把综合指数变为加权平均数指数形式,是为了（）。 D适应实际资料的要求 18洛纶茨曲线可以用以表示（）。 B累积的次数的分配 19确定连续变量的组限时，相邻组的组限是（）。 B重叠的 20抽样调查抽取样本时，必须遵守的原则是（）。 D随机性 21均值为20，变异系数为0.4，则标准差为（）。 B8 变异系数又称标准差率,是指标准差与平均数的比值,因此,标准差为0.4*20=8 22调查对象与调查单位具有一定的对应关系。如果调查对象是全部商业企业，则调查单位是（）。 A每一个商业企业 23几何平均法计算平均发展速度是（）指标连乘积的N次方根。 A环比发展速度 24F检验主要是用来检验( )。 C回归方程的显著性 25已知一个时间数列的环比增长速度分别为5%、2%、3%，则该时间数列的平均增长速度为（）。 A3.33%
实验三 利用DFT 分析连续信号频谱 一、实验目的 应用离散傅里叶变换(DFT)，分析模拟信号x (t )的频谱。深刻理解利用DFT 分析模拟信号频谱的原理，分析过程中出现的现象及解决方法。 二、 实验原理 连续周期信号相对于离散周期信号，连续非周期信号相对于离散非周期信号，都可以通过时域抽样定理建立相互关系。因此，在离散信号的DFT 分析方法基础上，增加时域抽样的步骤，就可以实现连续信号的DFT 分析。 三、实验内容 1. 利用FFT 分析信号)(e )(2t u t x t -=的频谱。 (1) 确定DFT 计算的各参数（抽样间隔，截短长度，频谱分辨率等）； 答：选取fm=25Hz 为近似的最高频率，则抽样间隔T=)2/(1m f =0.02s 选取6=p T s 进行分析，则截短点数为N==T T p /300 采用矩形窗，确定频域抽样点数为512点。 Matlab 函数如下：%对连续信号x=e(-2t)分析 fsam=50;Tp=6; N=512; T=1/fsam; t=0:T:Tp; x=exp(-2*t); X=T*fft(x,N); subplot(2,1,1);plot(t,x); xlabel('t');title('时域波形 N=512');legend('理论值'); w=(-N/2:N/2-1)*(2*pi/N)*fsam; y=1./(j*w+2); subplot(2,1,2);plot(w,abs(fftshift(X)),w,abs(y),'r-.'); title('幅度谱 N=512');xlabel('w'); legend('理论值','计算值',0); axis([-10,10,0,1.4])
数学实验作业汇总 文稿归稿存档编号：[KKUY-KKIO69-OTM243-OLUI129-G00I-FDQS58-
（1）产生一个5阶魔方矩阵M：M=magic(5) （2）将矩阵M的第3行4列元素赋值给变量t：t=M(3,4) （3）将由矩阵M第2，3，4行第2,5列构成的子矩阵赋给变N：N=M(2:4,2:3:5) （4）将由矩阵M的前3行赋给变量N：N=M(1:3,:) （5）将由矩阵M的后3列赋给变量N：N=M(:,end:-1:end-2) （6）提取M的主对角线元素，并以这些对角线元素构成对角矩阵N：N=diag(diag(M))或N=tril(triu(M)) (7)随机产生1000个100以内的整数赋值给变量t：t=round(rand(1,1000)*100) （8）随机产生100*5个100以内的实数赋值给变量M：M=rand(100,5)*100 （1）删除矩阵M的第7个元素M(7)=[] （2）将含有12个元素的向量t转换成3*4的矩阵：reshape(t,3,4) （3）产生和M同样大小的单位矩阵：eye(size(M)) （4）寻找向量t中非零元素的下标：find(t) （5）逆序显示向量t中的元素：t(end:-1:1) （6）显示向量t偶数位置上的元素：t(2:2:end) （7）利用find函数，将向量t中小于10的整数置为0：t(find(t<10&rem(t,1)==0))=0（8）不用find函数，将向量t中小于10的整数置为0：t(t<10&rem(t,1)==0)=0 （9）将向量t中的0元素用机器0（realmin）来代替：t(find(t=0))=realmin （10）将矩阵M中小于10的整数置为0：M(find(M<10)&rem(M,1)==0)=0
微机原理实验报告 —实验二熟悉汇编程序建立 及其调试方法 姓名：刘莹莹 学号：13024104
一、实验目的 1、熟悉汇编语言源程序的框架结构，学会编制汇程序。 2、熟悉汇编语言上机操作的过程，学会汇编程序调试方法。 二、实验内容 1、学习编写汇编语言源程序的方法，了解数据存放格式。 2、阅读给出的程序，找出程序中的错误。 3、通过调试给出的汇编语言源程序，了解并掌握汇编语言程序的建立、汇编、链接、调试、修改和运行等全过程。 三、实验预习 1、阅读实验指导第一章的内容，了解汇编语言程序建立、汇编、链接、调试的全过程。 2、下面的汇编语言源程序有错误的，试给程序加注释。通过调试手段找出程序中的错误并修改之。写出程序的功能，画出程序流程图。（1）程序修改前： STACKSG：SEGMENT PARA STACK ‘STACK’ ；不能有冒号DB 256 DUP(?) STACKSG ENDS DATASG: SEGMENT PARA ‘DATA’；不能有冒号BLOCK DW 0,-5,8,256,-128,96,100,3,45,6,512 DW 23,56,420,75,0,-1024,-67,39,-2000 COUNT EQU 20 MAX DW ? DATASG ENDS CODESG: SEGMENT ；不能有冒号ASSUME SS:STACKSG , CS:CODESG ASSUME DS:DATASG ORG 100H BEGIN MOV DS, DATASG ；BEGIN: 应改为MOV AX , DATASG MOV DS, AX LEA SI ,BLOCK
实验6 离散时间系统的z 域分析 （综合型实验） 一、实验目的 1） 掌握z 变换及其反变换的定义，并掌握MATLAB 实现方法。 2） 学习和掌握离散时间系统系统函数的定义及z 域分析方法。 3） 掌握系统零极点的定义，加深理解系统零极点分布与系统特性的关系。 二、实验原理与方法 1. z 变换 序列(n)x 的z 变换定义为(z)(n)z n n X x +∞ -=-∞ = ∑ （1） Z 反变换定义为11(n)(z)z 2n r x X dz j π-= ? （2） MATLAB 中可采用符号数学工具箱ztrans 函数和iztrans 函数计算z 变换和z 反变换： Z=ztrans(F)求符号表达式F 的z 变换。 F=iztrans(Z)求符号表达式Z 的z 反变换 2. 离散时间系统的系统函数 离散时间系统的系统函数H(z)定义为单位抽样响应h(n)的z 变换 (z)(n)z n n H h +∞ -=-∞ = ∑ （3） 此外连续时间系统的系统函数还可由系统输入与输出信号z 变换之比得到 (z)(z)/X(z)H Y = （4） 由（4）式描述的离散时间系统的系统时间函数可以表示为 101101...(z)...M M N N b b z b z H a a z a z ----+++=+++ （5） 3. 离散时间系统的零极点分析 MATLAB 中可采用roots 来求系统函数分子多项式和分母多项式的根，从而得到系统的零极点。 此外还可采用MATLAB 中zplane 函数来求解和绘制离散系统的零极点分布图，zplane 函数的调用格式为： zplane(b,a) b 、a 为系统函数分子分母多项式的系数向量（行向量） zplane(z,p) z 、p 为零极点序列（列向量） 系统函数是描述系统的重要物理量，研究系统函数的零极点分布不仅可以了解系统单位抽样响应的变化，还可以了解系统频率特性响应以及判断系统的稳定性； 系统函数的极点位置决定了系统的单位抽样响应的波形，系统函数零点位置只影响冲激响应的幅度和相位，不影响波形。 系统的频率响应取决于系统函数的零极点，根据系统的零极点分布情况，可以通过向量法分析系统的频率响应。
数学实验作业一 对以下问题,编写M文件: (1)用起泡法对10个数由小到大排序. 即将相邻两个数比较,将小的调到前头. 解： 代码如下： zuoye1 clear all;clc; a=[7 2 1 0 9 4 5 -3 8 6]; n=length(a); for ii=1:n-1 if a(ii+1)>=a(ii) t1=a(ii); a(ii)=a(ii+1); a(ii+1)=t1; end for jj=1:n-1 if a(jj+1)>=a(jj) t2=a(jj); a(jj)=a(jj+1); a(jj+1)=t2; end end end a 运行结果显示如下： a = 9 8 7 6 5 4 2 1 0 -3
(2)有一个 矩阵,编程求出其最大值及其所处的位置. 解： 代码如下：zuoye2.m clear; clc; a=[1 2 3 4 5 3 4 5 6 9 6 7 8 8 0 1 2 4 5 6] max=-1; flage1=0; flage2=0 for i=1:4 for j=1:5 if (a(i,j)>max) t=max; max=a(i ,j); a(i,j)=t; flage1=i; flage2=j ； end end end max flage1 flage2 运行结果显示如下： a = 1 2 3 4 5 3 4 5 6 9 6 7 8 8 0 1 2 4 5 6 flage2 = max = 45′
9 flage1 = 2 flage2 = 5 结果： （3）编程求∑=20 1 !n n 。 解： 代码如下：zuoye3.m clear; clc; sum=0; for i=2:11 sum=sum+gamma(i); end sum
4. 问题： 某公司将3种不同含硫量的液体原料（分别记为甲、乙、丙）混合生产两种产品（分别记为A,B ）。按照生产工艺的要求，原料甲、乙必须首先导入混合池中混合，混合后的液体再分别与原料丙混合生产A,B 。一直原料甲、乙、丙的含硫量分别是3%，1%，2%，进货价格分别为6千元/t ，16千元/t ，10千元/t ；产品A,B 的含硫量分别不能超过2.5%，1.5%，售价分别为9千元/t ，15千元/t 。根据市场信息，原料甲、乙、丙的供应量都不能超过500t ；产品A,B 的最大市场需求量分别为100t ，200t 。 (1) 应如何安排生产？ (2) 如果产品A 的最大市场需求量增长为600t ，应如何安排生产？ (3) 如果乙的进货价格下降为13千元/t ，应如何安排生产？分别对(1)、(2)两种情况进 行讨论。 模型： （只考虑问题1，问题2,3只需改变一些约束条件） 设生产时使用原料甲、乙分别为12,x x t ，分别取混合后的液体34,x x t 再加入原料丙 56,x x t 生产产品A,B 。 有质量守恒，可得 1234x x x x +=+ 甲乙混合后的液体的含硫量可表示为 12 12 3%x x x x ++，根据含硫量的要求，可得 12 353512 124646 12 3%*2%* 2.5%*()3%*2%* 1.5%*() x x x x x x x x x x x x x x x x +?+≤+?+?? +?+≤+?+? 根据市场的限制，易得 12563546500 500500100200 x x x x x x x x ≤?? ≤?? +≤??+≤??+≤? 当然还有非负约束 123456,,,,,0x x x x x x ≥ 公司的净利润为（单位：千元）： 35461256123456 9()15()61610()6169155z x x x x x x x x x x x x x x =+++---+=--++-+