matlab实验报告

《通信系统仿真》

实验报告

姓名：杨仕浩

专业：通信工程

学号： 2014111347

实验一 MATLAB环境的熟悉与基本运算

一、实验目的：

1．熟悉MATLAB开发环境

2．掌握矩阵、变量、表达式的各种基本运算

二、实验基本知识：

1.熟悉MATLAB环境:

MATLAB桌面和命令窗口、命令历史窗口、帮助信息浏览器、工作空间浏览器文件和搜索路径浏览器。

变量命名规则如下：

（1）变量名可以由英语字母、数字和下划线组成

（2）变量名应以英文字母开头

（3）长度不大于31个

（4）区分大小写

MATLAB中设置了一些特殊的变量与常量，列于下表。

三、实验内容

1、学习安装MATLAB软件。

2、学习使用help命令，例如在命令窗口输入help eye，然后根据帮助说明，

学习使用指令eye（其它不会用的指令，依照此方法类推）

3、学习使用clc、clear，观察command window、command history和workspace

等窗口的变化结果。

4、初步程序的编写练习，新建M-file，保存（自己设定文件名，例如exerc1、

exerc2、 exerc3……），学习使用MATLAB的基本运算符、数组寻访指令、

标准数组生成函数和数组操作函数。

注意：每一次M-file的修改后，都要存盘。

练习A：

（1） help rand，然后随机生成一个2×6的数组，观察command window、command history和workspace等窗口的变化结果。

（2）学习使用clc、clear，了解其功能和作用。

（3）输入C=1:2:20，则C（i）表示什么？其中i=1,2,3, (10)

（4）输入A=[7 1 5；2 5 6；3 1 5]，B=[1 1 1; 2 2 2; 3 3 3]，在命令窗口中执行下列表达式，掌握其含义：

A(2, 3) A(:,2) A(3,:) A(:,1:2:3) A(:,3).*B(:,2)

A(:,3)*B(2,:) A*B A.*B A^2 A.^2 B/A B./A （5）二维数组的创建和寻访，创建一个二维数组（4×8）A，查询数组A 第2行、第3列的元素，查询数组A第2行的所有元素，查询数组A

第6列的所有元素。

（6）两种运算指令形式和实质内涵的比较。设有3个二维数组A

2×4，B

2

×4，C

2×2

，写出所有由2个数组参与的合法的数组运算和矩阵指令。

（7）学习使用表8列的常用函数（通过help方法）

（8）学习使用表9数组操作函数。

练习B

（9）用reshape指令生成下列矩阵（不超过3次操作），并取出方框内的数组元素。

四、试验结果

（1）

（2）clc:清除命令窗口所有内容，数值不变；clear:初始化变量的值。（3）

C(i)表示C中的第i个的数值。

（4）

（5）（6）

（9）

实验二 MATLAB数值及符号运算

一、实验目的：

1、掌握矩阵的基本运算

2、掌握矩阵的数组运算

3、掌握多项式的基本运算

4、会求解代数方程

5、掌握创建符号表达式和矩阵的方法

6、掌握符号表达式的微分和积分运算

二、实验基本知识：

1、多项式运算

poly ——产生特征多项式系数向量

roots ——求多项式的根

p=poly2str(c,‘x’)（以习惯方式显示多项式）

conv,convs多项式乘运算

deconv多项式除运算

多项式微分

polyder(p): 求p的微分

polyder(a,b): 求多项式a,b乘积的微分

2、代数方程组求解

1）.恰定方程组的解

方程ax=b(a为非奇异)两种求解方法： x=inv(a) b —采用求逆运算解方程

x=a\b —采用左除运算解方程

2）超定方程组的解

方程 ax=b ,m

3）欠定方程组的解

matlab可求出两个解：用除法求的解x是具有最多零元素的解

是具有最小长度或范数的解，这个解是基于伪逆pinv求

得的。

3、符号矩阵的创建

1）用matlab函数sym创建矩阵（symbolic 的缩写)

命令格式：A=sym('[ ]')

2）用字符串直接创建矩阵

4、符号矩阵与数值矩阵的转换

将数值矩阵转化为符号矩阵，函数调用格式：double(A)

将符号矩阵转化为数值矩阵，函数调用格式： numeric(A)

5、符号微积分与积分变换

diff(f) —对缺省变量求微分

diff(f,v) —对指定变量v求微分

diff(f,v,n) —对指定变量v求n阶微分

int(f) —对f表达式的缺省变量求积分

int(f,v) —对f表达式的v变量求积分

int(f,v,a,b) —对f表达式的v变量在（a,b)区间求定积分

6、符号代数方程求解

solve(f) —— 求一个方程的解

solve(f1,f2, …fn) —— 求n 个方程的解 7、符号微分方程求解指令：dsolve

命令格式：dsolve(f,g)： f —— 微分方程，可多至12个微分方程的求

解；g 为初始条件，默认自变量为 'x',可任意指定自变量't', 'u'等，微分方程的各阶导数项以大写字母D 表示

三、实验内容：

1、已知a=[1 2 3;4 5 6;7 8 0]，求其特征多项式并求其根。

2、已知多项式a(x)=x 2+2x+3，b(x)=4x 2+5x+6，求a ，b 的积并微分。

3、求解方程 1）???=+=+ 133x 28

221

21x x x

2）???

??=+=+=+3

43 23x 21221

2121x x x x x

3）???=++=++ 243x 21

32321

321x x x x x

4、用两种方法创建符号矩阵，A =[ a, 2*b]

[3*a, 0]

5、计算二重不定积分

6、对符号方程f = ax2+bx+c 求解。1）对x 求解，2）对a 求解。

7、求解微分方程0)0(,1)0(,02222===++dx dy

y y dx dy dx

y d 。

8、已知系统冲激响应为h(n)=[0.05,0.24,0.40,0.24,0.15,-0.1,0.1] ，系统输入u(n)由指令randn('state',1);u=2*(randn(1,100)>0.5)-1产生，该输入信号

dxdy

xe xy

??-

的起始作用时刻为0。试用直杆图（提示：用stem指令）画出分别显示该系统输入、输出信号的两张子图

四、实验结果

（1）（2）

（3）（4）（5）（7）

（6）

实验三 MATLAB语言的程序设计

一、实验目的：

1、熟悉MATLAB程序编辑与设计环境

2、掌握各种编程语句语法规则及程序设计方法

3、函数文件的编写和设计

4、了解和熟悉跨空间变量传递和赋值

二、实验基本知识：

1、程序流程控制语句

for循环结构

语法：for i＝初值：增量：终值

语句1

……

语句n

end

说明：1．i＝初值：终值，则增量为1。

2．初值、增量、终值可正可负，可以是整数，也可以是小数，只须符合数学逻辑。

while 循环结构

语法：while 逻辑表达式

循环体语句

end

说明：1、whiIe结构依据逻辑表达式的值判断是否执行循环体语勾。若表达式的值为真，执行循环体语句一次、在反复执行时，每次都要进行判断。若表

达式的值为假，则程序执行end之后的语句。

2、为了避免因逻辑上的失误，而陷入死循环，建议在循环体语句的适当位

置加break语句、以便程序能正常执行。（执行循环体的次数不确定；

每一次执行循环体后，一定会改变while后面所跟关系式的值。）

3、while循环也可以嵌套、其结构如下：

while逻辑表达式1

循环体语句1

while逻辑表达式2

循环体语句2

end

循环体语句3

end

if-else-end分支结构

if 表达式 1

语句1

else if 表达式 2(可选)

语句2

else(可选)

语句3

end

end

说明：1.if结构是一个条件分支语句，若满足表达式的条件，则往下执行；若不满足，则跳出if结构。

2．else if表达式2与else为可选项，这两条语句可依据具体情况取舍。

3．注意：每一个if都对应一个end，即有几个if，记就应有几个end。switch-case结构

语法：switch表达式

case常量表达式1

语句组1

case常量表达式2

语句组2

……

otherwise

语句组n

end

说明：1．switch后面的表达式可以是任何类型，如数字、字符串等。

2．当表达式的值与case后面常量表达式的值相等时，就执行这个case后面的语句组如果所有的常量表达式的值都与这个表达式的值不相等时，则执行otherwise后的执行语句。

程序流程控制指令

break、return、pause

2、函数文件的结构和编写方法

函数只能在函数体内对变量进行操作，也就是只能访问函数本身工作空间中的变量。

M函数文件的结构

（1）函数定义行（function）

（2） H1行（函数帮助文本的第一行）

（3）函数帮助文本

（4）函数体

（5）注释

函数文件编写后，保存时，其文件名必须与函数名相同。注意：函数名不要

与MATLAB 自身的函数命令相同。

三、实验内容：

练习A

1、熟悉MATLAB 程序编辑与设计环境

2、用for 循环语句实现求1～100的和。

3、用for 循环语句实现编写一个求n 阶乘的函数文件。

4、一行100列的Fibonacci 数组a ，元素满足Fibonacci 规则：a(1)=a(2)=1， a(i)=a(i-1)+a(i-2)(i=3,4,…)，求该数组中第一个大于10000的元素。 练习B 思考题

4、已知一维数组A=[2,4,5,8,10],B=[4,6,9,3,4]，用for 循环语句实现

∑=+-n

i i n i

B

A 1

1

。求和函数可用sum()。.

5、利用for 循环求1！+2！+3！+ … +5！的值。.

6、编写M 函数表示曲线y2=e -t/3sin(3t)以及它的包络线y1=e -t/3，并从命令窗口输入命令语句绘制曲线，t 的取值范围是[0,4π]。

7、设01

.0)3(1

1.0)2(1)(4

2+-++-=

x x x f ，编写一个M 函数文件，使得调用f(x)

时，x可用矩阵代入，得出的f(x) 为同阶矩阵。

四、实验结果

（1）

（2）

（3）（4）

（5）

（6）

（7）

实验四采用SIMULINK的系统仿真

一、实验目的：

1、熟悉SIMULINK 工作环境及特点

2、掌握线性系统仿真常用基本模块的用法

3、掌握SIMULINK 的建模与仿真方法

二、实验基本知识：

1．了解SIMULINK模块库中各子模块基本功能

MATLAB实验报告50059

实验一MATLAB操作基础 实验目的和要求： 1、熟悉MATLAB的操作环境及基本操作方法。 2、掌握MATLAB的搜索路径及设置方法。 3、熟悉MATLAB帮助信息的查阅方法 实验内容： 1、建立自己的工作目录，再设置自己的工作目录设置到MA TLAB搜索路径下，再试 验用help命令能否查询到自己的工作目录。 2、在MA TLAB的操作环境下验证课本；例1-1至例1-4，总结MATLAB的特点。 例1-1

例1-2 例1-3 例1-4

3、利用帮助功能查询inv、plot、max、round等函数的功能。 4、完成下列操作： （1）在matlab命令窗口输入以下命令： x=0:pi/10:2*pi; y=sin(x); （2）在工作空间窗口选择变量y，再在工作空间窗口选择回绘图菜单命令或在工具栏中单击绘图命令按钮，绘制变量y的图形，并分析图形的含义。

5、访问mathworks公司的主页，查询有关MATLAB的产品信息。 主要教学环节的组织： 教师讲授实验目的、开发环境界面、演示实验过程，然后同学上机练习。 思考题： 1、如何启动与退出MA TLAB集成环境？ 启动： （1）在windows桌面，单击任务栏上的开始按钮，选择‘所有程序’菜单项，然后选择MA TLAB程序组中的MA TLABR2008b程序选项，即可启动 MATLAB系统。 （2）在MA TLAB的安装路径中找到MA TLAB系统启动程序matlab.exe，然后运行它。 （3）在桌面上建立快捷方式后。双击快捷方式图标，启动MA TLAB。 退出： （1）在MA TLAB主窗口file菜单中选择exitMATLAB命令。 （2）在MA TLAB命令窗口中输入exit或quit命令。 （3）单击MATLAB主窗口的关闭按钮。 2、简述MATLAB的主要功能。 MATLAB是一种应用于科学计算领域的数学软件，它主要包括数值计算和符 号计算功能、绘图功能、编程语言功能以及应用工具箱的扩展功能。 3、如果一个MATLAB命令包含的字符很多，需要分成多行输入，该如何处理？

《MATLAB与数值分析》第一次上机实验报告

电子科技大学电子工程学院标准实验报告（实验）课程名称MATLAB与数值分析 学生姓名：李培睿 学号：2013020904026 指导教师：程建

一、实验名称 《MATLAB与数值分析》第一次上机实验 二、实验目的 1. 熟练掌握矩阵的生成、加、减、乘、除、转置、行列式、逆、范数等运算 操作。（用.m文件和Matlab函数编写一个对给定矩阵进行运算操作的程序） 2. 熟练掌握算术符号操作和基本运算操作，包括矩阵合并、向量合并、符号 转换、展开符号表达式、符号因式分解、符号表达式的化简、代数方程的符号解析解、特征多项式、函数的反函数、函数计算器、微积分、常微分方程的符号解、符号函数的画图等。（用.m文件编写进行符号因式分解和函数求反的程序） 3. 掌握Matlab函数的编写规范。 4、掌握Matlab常用的绘图处理操作，包括：基本平面图、图形注释命令、 三维曲线和面的填充、三维等高线等。（用.m文件编写在一个图形窗口上绘制正弦和余弦函数的图形，并给出充分的图形注释） 5. 熟练操作MATLAB软件平台，能利用M文件完成MATLAB的程序设计。 三、实验内容 1. 编程实现以下数列的图像，用户能输入不同的初始值以及系数。并以x， y为坐标显示图像 x(n+1) = a*x(n)-b*(y(n)-x(n)^2); y(n+1) = b*x(n)+a*(y(n)-x(n)^2) 2. 编程实现奥运5环图，允许用户输入环的直径。 3. 实现对输入任意长度向量元素的冒泡排序的升序排列。不允许使用sort 函数。 四、实验数据及结果分析 题目一： ①在Editor窗口编写函数代码如下：

(完整word版)基于单片机仿真软件proteus的流水灯实验报告

单片机原理》实验报告 题目：流水灯设计 姓名：刘伟 学号：1042157110 专业：10 计科特色班院系：信息工程学院指导老师：史先桂 完成时间：2014 年 4 月19 日 安徽新华学院教务处制

一、实验目的 1、熟练掌握单片机仿真软件proteus 使用方法和注意事项。 2、了解简单单片机应用系统的设计方法。 3、帮助学生养成良好实验习惯。 二、实验内容单片机仿真软件proteus 上实现8个发光LED “流水”的现象，实现两个流水灯情况：1、先奇数 灯亮，再偶数灯亮；2、实现流水灯从两边向中间亮,再从中间到两边亮。并通过编写程序控制流水现象。 三、实验说明依照实验的硬件电路原理，在单片机仿真软件proteus 上进行硬件电路的模拟，然后进行实验。在发光二极管两次点亮的间隔中加延时程序，让每次点亮停留一段时间，像这样人眼就可以看到“流水” 的现象。 四、实验环境 硬件：pc 机；软件：单片机仿真软件proteus。 五、实验原理图 六、实验参考程序 #include #include void delay_ms(int n) { int i,j; for(i=0;i

Matlab数学实验报告一

数学软件课程设计 题目非线性方程求解 班级数学081 姓名曹曼伦

实验目的：用二分法与Newton迭代法求解非线性方程的根； 用Matlab函数solve、fzero、fsolve求解非线性方程（组）的解。 编程实现二分法及Newton迭代法； 学会使用Matlab函数solve、fzero、fsolve求解非线性方程（组）的解。 通过实例分别用二分法及迭代法解非线性方程组并观察收敛速度。 实验内容： 比较求exp(x)+10*x-2的根的计算量。(要求误差不超过十的五次方） (1)在区间（0，1)内用二分法； (2)用迭代法x=(2-exp(x))/10，取初值x=0 。 试验程序 (1)二分法: format long syms x s=exp(x)+10*x-2 a=0; b=1; A=subs(s,a) B=subs(s,b) f=A*B %若f<0，则为由根区间 n=0; stop=1.0e-5; while f<0&abs(a-b)>=stop&n<=100; Xk=(a+b)/2; %二分 M= subs(s, Xk); if M* A<0 symbol=1 %若M= subs(s, Xk)为正，则与a二分 b= Xk else symbol=0 % 若M= subs(s, Xk)为负，则与b二分 a= Xk end n=n+1 end Xk n (2)牛顿迭代法； format long

syms x s= (2-exp(x))/10; %迭代公式 f=diff(s); x=0; %迭代初值 a=subs(f,x); %判断收敛性(a是否小于1) s=(2-exp(x))/10; stop=1.0e-5; %迭代的精度 n=0; while a<1&abs(s-x)>=stop&n<=100; x=s %迭代 s=(2-exp(x))/10; n=n+1 end 实验结果： （1）二分法： symbol =1 b =0.50000000000000 n =1 symbol =1 b =0.25000000000000 n =2 symbol =1 b =0.12500000000000 n =3 symbol =0 a =0.06250000000000 n =4 symbol =1 b =0.09375000000000 n =5 symbol =0 a =0.07812500000000 n =6 symbol =1 b =0.09054565429688 n =15 symbol =1 b =0.09053039550781 n =16 symbol =0 a =0.09052276611328 n =17 Xk =0.09052276611328 n =17 （2）迭代法 由x =0.10000000000000 n =1 x =0.08948290819244 n =2 x =0.09063913585958 n =3 x =0.09051261667437 n =4 x =0.09052646805264 n =5 试验结果可见用二分法需要算17次，而用迭代法求得同样精度的解仅用5次，但由于迭代法一般只具有局部收敛性，因此通常不用二分法来求得非线性方程的精确解，而只用它求得根的一个近似解，再用收敛速度较快的迭代法求得其精确解。

Matlab实验报告3

实验三函数的可视化与Matlab作图 一、按要求绘制如下曲线（面）： 1. 在[0，4pi]上画sin(x),cos(x)在同一图像中，其中cos(x)图像用红色小圆圈，并在函数图上标注“y=sin(x)”，”y=cos(x)”,X轴，Y轴，标题为“正弦余弦函数图像。”答：>> clear >> clf, x=linspace(0,4*pi,200);y1=sin(x);y2=cos(x); plot(x,y1,'k-',x,y2,'ro') >> title('正弦余弦函数图像。') >> legend('y=sin(x)','y=cos(x)') >> ylabel('\it{Y轴}'); >> xlabel('\it{X轴}'); 2.任意绘制彗星曲线图。 答：>> clf; >> x=[1:10]; y=[5 6 3 4 8 1 10 3 5 6]; >> z=0:0.1:100; x=sin(z);y=cos(z).*10; >> %三维彗星图 comet3(x,y,z) >> %二维彗星图

t = -pi:pi/200:pi; comet(t,tan(sin(t))-sin(tan(t)))

3.在多窗口中绘制y=sin(t)*sin(t);y1=sin(3*t+2.5);y2=sin(5*t+5)并加以标注。答：>> clf; t=0:0.1:4*pi; subplot(3,1,1),plot(sin(t).*sin(t)),legend('y=sin(t)*sin(t)') subplot(3,1,2),plot(sin(3*t+2.5)),legend('y1=sin(3*t+2.5)') subplot(3,1,3),plot(sin(5*t+5)),legend('y2=sin(5*t+5)') 4.自拟题目绘制三维线图。 绘制以下方程y1=sin(t),y2=cos(t),x=t在t=[0,2π] 对应的三维曲线。 >> clf; >> t=0:pi/10:2*pi; >> y1=sin(t);y2=cos(t); >> plot3(y1,y2,t);grid on; >> xlabel('Dependent Variable Y1'); >> ylabel('Dependent Variable Y2'); >> zlabel('Dependent Variable X'); >> title('Sin and Cos Curve');

实验一proteus仿真软件使用方法

实验一 Proteus仿真软件使用方法 一．实验目的： （1）了解Proteus仿真软件的使用方法。 （2）了解51单片机编程器Keil与Proteus仿真软件的联用方法。 二．实验要求： 通过讲授和操作练习，学会正确使用Proteus仿真软件及Keil编程及其联合调试。 三．实验内容： （1）Proteus 仿真软件介绍 Proteus 软件是由英国LabCenter Electronics 公司开发的EDA工具软件，由ISIS 和ARES两个软件构成，其中ISIS是一款便捷的电子系统仿真平台软件，ARES是一款高级的布线编辑软件。它集成了高级原理布线图、混合模式SPICE电路仿真、PCB设计以及自动布线来实现一个完整的电子设计。 通过Proteus ISIS软件的VSM(虚拟仿真技术), 用户可以对模拟电路、数字电路、模数混合电路，以及基于微控制器的系统连同所有外围接口电子元器件一起仿真。 图1是Proteus ISIS的编辑窗口： 图1 ISIS的编辑界面 图中最顶端一栏是“标题栏”，其下的“File View Edit ……”是“菜单栏”，再下面的一栏是“命令工具栏”，最左边的一栏是“模式选择工具栏”；左上角的小方框是“预览窗口”，其下的长方框是“对象选择窗口”，其右侧的大方框是“原理图编辑窗口”。 选择左侧“模式选择工具栏”中的图标，并选择“对象选择窗口”中的P按钮，就会出现如图2的元器件选择界面：

图2 元器件库选择界面 在元器件列表框中点击你需要的器件类型（例如：电阻-Resistors,单片机芯片-MicroprocessorICs, LED-Optoelectronics）或在左上角的关键字（Keywords）框中输入你需要的器件名称的关键字（如：信号源 - Clock, 运放 - CA3140等），就会在图2中间的大空白框列出你所需的一系列相关的元件。此时，你可用鼠标选中你要的元件，则图2右上角的预览框会显示你所要元件的示意图，若就是你要的元器件，则点击OK按钮，该元器件的名称就会列入位于图1左侧的“对象选择窗口”中（参见图1左侧下方框）。 所需元器件选择好后，在“对象选择窗口”选择某器件，就可以将它放到图1中的“原理图编辑窗口”中（若器件的方向不合适，你可以利用图1左下角的旋转按钮来改变它）。将所要的元器件都选好后，将它们安放到合适的位置，就可以用连接线把电路连接好，结果存盘（请按规定的目录存盘，并记住其路径/目录/文件名[学号-实验序号]）。 （2）51单片机编程器– Keil V3的使用 Keil编程器可用于MCS-51单片机软件编程与调试，它的工作界面如图3所示： Keil编程器是Keil Software Inc/Keil Electronic GmbH 开发的基于80C51内核的微处理器软件开发平台，可以完成从工程（Project）的建立和管理、程序的编译和连接、目标代码的形成、软件仿真等一套完整的软件开发流程。它与Proteus挂接，可以进行单片机应用系统的硬件仿真。 汇编语言编程方法： ①打开“File”菜单→选择新建“New...”→在弹出的文本框（Text1）中编写所需的汇编语言程序→程序写好后，保存（从File→Save As..→选择某目录，文件名.ASM, 存盘）； ②打开“Project”菜单→选择新建工程“New Project...”→在弹出的窗口填写：工程名→保存（文件名的后缀是 .uv2 。此时图3的工程窗口中将建立Target1

数学实验“几种常见的求积分近似解的方法”实验报告(内含matlab程序)

西京学院数学软件实验任务书

实验二十一实验报告 一、实验名称：Romberg 积分法，Gauss 型积分法，高斯-勒让德积分法，高斯-切比雪夫积分法，高斯-拉盖尔积分法，高斯-埃尔米特积分法。 二、实验目的：进一步熟悉Romberg 积分法，Gauss 型积分法，高斯-勒让德积分法，高斯-切比雪夫积分法，高斯-拉盖尔积分法，高斯-埃尔米特积分法。 三、实验要求：运用Matlab/C/C++/Java/Maple/Mathematica 等其中一种语言完成程序设计。 四、实验原理： 1．Romberg 积分法： 龙贝格积分法是用里查森外推算法来加快复合梯形求积公式的收敛速度，它的算法如下，其中()i m T 是通过一系列逼近原定积分的龙贝格分值. 计算(0)1[()()]2 b a T f a f b -= + 对1,2,3,k n = ，计算下列各步： 21()(1)1 111 1(21)()[()]222k k k k k j b a j b a T T f a ---=---=++∑

对1,2,,m k = 和,1,2,,1i k k k =-- ，计算111 441 m i i i m m m m T T T --+-=- 随着计算的步骤的增加，()i m T 越来越逼近积分()b a f x dx ?。 2．Gauss 型积分法： 高斯积分公式的思想是用n 个不等距的节点123,,,n x x x x 对被积函数进行插值，然后对插值后的函数进行积分，其积分公式为： 1 1 1 ()()n k k k f x dx A f x -=≈∑? 如果积分区间不是[1,1]-，则需转换到此区间： 11()()222 b a b a b a b a f x dx f t dt ---+= +? ? 其中系数k A 、节点k x 与n 的关系如下表所示： 3．高斯-切比雪夫积分法： 第一类切比雪夫积分形式为： 1 1 ()()n k k k f x dx A f x -=≈∑? 其中k A n π= ，21cos 2k k x n π-= 4．高斯-拉盖尔积分法： 高斯-拉盖尔公式有两种形式： 1 ()()n x k k k e f x dx A f x +∞ -=≈∑?

matlab实验报告3详解

实验四、LTI系统的响应 课程名称： MATLAB应用技术专业班级：通信1422 学生学号： 1430119231 学生姓名：周妍智 所属院部：电子信息工程系指导教师：徐树梅 2015 —— 2016 学年第二学期

实验项目名称： LTI 系统的响应 实验学时： 16 学生姓名： 周妍智 实验地点： 微机11 实验日期： 2016.4.17 实验成绩： 批改教师： 徐树梅 批改时间： 一、 实验目的 1. 熟悉连续时间系统的单位冲激响应、阶跃响应的意义及求解方法 2. 熟悉连续（离散）时间系统在任意信号激励下响应的求解方法 3. 熟悉应用MATLAB 实现求解系统响应的方法 二、 实验原理 1.连续时间系统 对于连续的LTI 系统，当系统输入为f (t )，输出为y (t )，则输入与输出之间满足如下的线性常系数微分方程： () ()0 ()()n m i j i j i j a y t b f t ===∑∑，当系统输入为单位冲激信号δ(t )时产生 的零状态响应称为系统的单位冲激响应，用h(t)表示。若输入为单位阶跃信号ε(t )时，系统产生的零状态响应则称为系统的单位阶跃响应，记为g(t)，如下图所示。 系统的单位冲激响应h (t )包含了系统的固有特性，它是由系统本身的结构及参数所决定的，与系统的输入无关。我们只要知道了系统的冲激响应，即可求得系统在不同激励下产生的响应。因此，求解系统的冲激响应h(t )对我们进行连续系统的分析具有非常重要的意义。 在MATLAB 中有专门用于求解连续系统冲激响应和阶跃响应， 并绘制其时域波形的函数impulse( ) 和step( )。如果系统输入为f (t )，冲激响应为h(t)，系统的零状态响应为y (t )，则有：()()()y t h t f t =*。 若已知系统的输入信号及初始状态，我们便可以用微分方程的经典时域求解方法，求出系统的响应。但是对于高阶系统，手工计算这一问题的过程非常困难和繁琐。 在MATLAB 中，应用lsim( )函数很容易就能对上述微分方程所描述的系统的响应进行仿真，求出系统在任意激励信号作用下的响应。lsim( )函数不仅能够求出连续系统在指定的任意时间范围内系统响应的数值解，而且还能同时绘制出系统响应的时域波形图。 以上各函数的调用格式如下： ⑴ impulse( ) 函数 函数impulse( )将绘制出由向量a 和b 所表示的连续系统在指定时间范围内的单位冲激响应h (t )的时域波形图，并能求出指定时间范围内冲激响应的数值解。

基于单片机仿真软件proteus的流水灯实验报告

《单片机原理》实验报告题目：流水灯设计 姓名：刘伟 学号：1042157110 专业：10计科特色班 院系：信息工程学院 指导老师：史先桂 完成时间：2014 年 4 月19日安徽新华学院教务处制

一、实验目的 1、熟练掌握单片机仿真软件proteus使用方法和注意事项。 2、了解简单单片机应用系统的设计方法。 3、帮助学生养成良好实验习惯。 二、实验内容 单片机仿真软件proteus上实现8个发光LED“流水”的现象，实现两个流水灯情况：1、先奇数灯亮，再偶数灯亮；2、实现流水灯从两边向中间亮,再从中间到两边亮。并通过编写程序控制流水现象。 三、实验说明 依照实验的硬件电路原理，在单片机仿真软件proteus上进行硬件电路的模拟，然后进行实验。在发光二极管两次点亮的间隔中加延时程序，让每次点亮停留一段时间，像这样人眼就可以看到“流水”的现象。 四、实验环境 硬件：pc机；软件：单片机仿真软件proteus。 五、实验原理图 六、实验参考程序 #include #include void delay_ms(int n) //延时n毫秒 { int i,j; for(i=0;i

void main() { int i, num; unsigned char p1; unsigned char p0, p; while(1) { for(i=0; i<3; i++) //间隔500ms先奇数亮再偶数亮，循环三次 { P1=0xaa; P0=p1; delay_ms(500); P1=0x55; P0=p1; delay_ms(500); } p1= 0xfe; num=3; while(num-- > 0) //一个灯上下循环三? { for(i=0; i<8; i++) { P0=p1; delay_ms(100); p1=_cror_(p1,1); } } //两个分别从两边往中间流动三? p1= 0xfe; p0= 0x7f; num=3; while(num-- > 0) { for(i=0; i<4; i++) //让LED灯从两边向中间点亮 { p1=_crol_(p1,1); p0=_cror_(p0,1); p=p1&p0; delay_ms(100); P1=p;

matlab数学实验报告5

数学实验报告 制作成员班级学号 2011年6月12日

培养容器温度变化率模型 一、实验目的 利用matlab软件估测培养容器温度变化率 二、实验问题 现在大棚技术越来越好，能够将温度控制在一定温度范围内。为利用这种优势，实验室现在需要培植某种适于在8.16℃到10.74℃下能够快速长大的甜菜品种。为达到实验所需温度，又尽可能地节约成本，研究所决定使用如下方式控制培养容器的温度：1，每天加热一次或两次，每次约两小时； 2，当温度降至8.16℃时，加热装置开始工作；当温度达到10.74℃时，加热装置停止工作。 已知实验的时间是冬天，实验室为了其它实验的需要已经将实验室的温度大致稳定在0℃。下表记录的是该培养容器某一天的温度 时间(h)温度（℃）时间（h）温度（℃）09.68 1.849.31 0.929.45 2.959.13 3.878.981 4.989.65 4.988.811 5.909.41 5.908.691 6.839.18 7.008.5217.938.92 7.938.3919.048.66 8.978.2219.968.43 9.89加热装置工作20.848.22 10.93加热装置工作22.02加热装置工作10.9510.8222.96加热装置工作12.0310.5023.8810.59 12.9510.2124.9910.35 13.889.9425.9110.18 三、建立数学模型 1，分析：由物理学中的傅利叶传热定律知温度变化率只取决于温度

差，与温度本身无关。因为培养容器最低温度和最高温度分别是：8.16℃和10.74℃；即最低温度差和最高温度差分别是：8.16℃和10.74℃。而且，16.8/74.10≈1.1467，约为1，故可以忽略温度对温度变化率的影响2， 将温度变化率看成是时间的连续函数，为计算简单，不妨将温度变化率定义成单位时间温度变化的多少，即温度对时间连续变化的绝对值（温度是下降的），得到结果后再乘以一系数即可。 四、问题求解和程序设计流程1）温度变化率的估计方法 根据上表的数据，利用matlab 做出温度-时间散点图如下： 下面计算温度变化率与时间的关系。由图选择将数据分三段，然后对每一段数据做如下处理：设某段数据为{(0x ,0y ),(1x ,1y ),(2x , 2y ),…,(n x ,n y )},相邻数据中点的平均温度变化率采取公式： 温度变化率=（左端点的温度-右端点的温度）/区间长度算得即：v( 2 1i i x x ++)=(1+-i i y y )/(i i x x - +1). 每段首尾点的温度变化率采用下面的公式计算：v(0x )=(30y -41y +2y )/(2x -0x )v(n x )=(3n y -41+n y +2+n y )/(n x -2-n x )

MATLAB实验报告(1-4)

信号与系统MATLAB第一次实验报告 一、实验目的 1.熟悉MATLAB软件并会简单的使用运算和简单二维图的绘制。 2.学会运用MATLAB表示常用连续时间信号的方法 3.观察并熟悉一些信号的波形和特性。 4.学会运用MATLAB进行连续信号时移、反折和尺度变换。 5.学会运用MATLAB进行连续时间微分、积分运算。 6.学会运用MATLAB进行连续信号相加、相乘运算。 7.学会运用MATLAB进行连续信号的奇偶分解。 二、实验任务 将实验书中的例题和解析看懂，并在MATLAB软件中练习例题，最终将作业完成。 三、实验内容 1.MATLAB软件基本运算入门。 1). MATLAB软件的数值计算： 算数运算 向量运算：1.向量元素要用”[ ]”括起来，元素之间可用空格、逗号分隔生成行向量，用分号分隔生成列向量。2.x=x0:step:xn.其中x0位初始值，step表示步长或者增量，xn 为结束值。 矩阵运算：1.矩阵”[ ]”括起来；矩阵每一行的各个元素必须用”,”或者空格分开； 矩阵的不同行之间必须用分号”;”或者ENTER分开。2.矩阵的加法或者减法运算是将矩阵的对应元素分别进行加法或者减法的运算。3.常用的点运算包括”.*”、”./”、”.\”、”.^”等等。 举例：计算一个函数并绘制出在对应区间上对应的值。

2).MATLAB软件的符号运算：定义符号变量的语句格式为”syms 变量名” 2.MATLAB软件简单二维图形绘制 1).函数y=f(x)关于变量x的曲线绘制用语：>>plot(x,y) 2).输出多个图像表顺序：例如m和n表示在一个窗口中显示m行n列个图像，p表 示第p个区域，表达为subplot(mnp)或者subplot(m,n,p) 3).表示输出表格横轴纵轴表达范围：axis([xmax,xmin,ymax,ymin]) 4).标上横轴纵轴的字母：xlabel(‘x’),ylabel(‘y’) 5).命名图像就在subplot写在同一行或者在下一个subplot前：title(‘……’) 6).输出：grid on 举例1： 举例2：

MATLAB实验报告

数字信号处理及MATLAB 实验报告 班级： 学号： 姓名：

4.7.2 例4,2 设x(n)是由两个正弦信号及白噪声的叠加，试用FFT文件对其作频谱分析。程序清单 %产生两个正弦加白噪声 N=256; f1=.1;f2=.2;fs=1; a1=5;a2=3; w=2*pi/fs; x=a1*sin(w*f1*(0:N-1))+a2*sin(w*f2*(0:N-1))+randn(1,N); %应用FFT求频谱 subplot(2,2,1); plot(x(1:N/4)); title('原始信号'); f=-0.5:1/N:0.5-1/N; x=fft(x); y=ifft(x); subplot(2,2,2); plot(f,fftshift(abs(x))); title('频域信号'); subplot(2,2,3); plot(real(x(1:N/4))); title('时域信号');

例4.3 设x（n）为长度N=6的矩形序列，用MATLAB程序分析FFT取不同长度时x（n）频谱的变化。N=8,32,64,时x（n）的FFT MATLAB实现程序如下。 x=[1,1,1,1,1,1]; N=8; y1=fft(x,N); n=0:N-1; subplot(3,1,1);stem(n,abs(y1),'.k');axis([0,9,0,6]); N=32; y2=fft(x,N); n=0:N-1; subplot(3,1,2);stem(n,abs(y2),'.k');axis([0,40,0,6]); N=64; y3=fft(x,N); subplot(3,1,3);stem(n,abs(y3),'.k');axis([0,80,0,6]);

浅析Matlab数学实验报告

数学实验报告 姓名: 班级: 学号: 第一次实验任务 过程: a=1+3i; b=2-i; 结果: a+b =3.0000 + 2.0000i a-b =-1.0000 + 4.0000i a*b = 5.0000 + 5.0000i a/b = -0.2000 + 1.4000i 过程: x=-4.5*pi/180; y=7.6*pi/180; 结果: sin(abs(x)+y)/sqrt(cos(abs(x+y))) =0.2098 心得:对于matlab 中的角度计算应转为弧度。 (1)过程: x=0:0.01:2*pi; y1=sin(x); y2=cos(x); y3=exp(x); y4=log(x); plot(x,y1,x,y2,x,y3,x,y4) plot(x,y1,x,y2,x,y3,x,y4) 结果: (2)过程:>> subplot(2,2,1) >> plot(x,y1) >> subplot(2,2,2) >> plot(x,y2) ./,,,,2,311b a b a b a b a i b i a ?-+-=+=计算、设有两个复数 6,7,5.4)

cos()sin(2=-=++y x y x y x ,其中、计算的图形。 下分别绘制)同一页面四个坐标系)同一坐标系下(、在( x y e y x y x y x ln ,,cos ,sin 213==== >> subplot(2,2,3) >> plot(x,y3) >> subplot(2.2.4) >> subplot(2,2,4) >> plot(x,y4) 结果: 心得:在matlab中,用subplot能够实现在同一页面输出多个坐标系的图像,应注意将它与hold on进行区别,后者为在同一坐标系中划出多条曲线。 5、随机生成一个3x3矩阵A及3x2矩阵B,计算(1)AB,(2)对B中每个元素平方后得到的矩阵C,(3)sinB,(4)A的行列式,(5)判断A是否可逆,若可逆,计算A的逆矩阵,(6)解矩阵方程AX=B,(7)矩阵A中第二行元素加1,其余元素不变,得到矩阵D,计算D。 过程:A=fix(rand(3,3).*10) ; B=fix(rand(3,3).*10);

计算机组成原理实验proteus的应用

南京理工大学紫金学院 计算机组成原理 实验报告

实验报告一：加法器实验操作 一、实验目的： 1.掌握 proteus 软件常用命令的使用方法 2.掌握加法器的基本使 用二、实验内容： 1． proteus 软件命令使用 2．串联加法器的连接 3．并联加法器的连接 三、实验步骤： (一 )proteus 软件界面的基本使用 ①通过 File 界面创建新文件或打开已有文件 ②所用的原件模型必须在蓝方框内放置

③Proteus 工具栏功能如上图 ④Proteus 中共有 36 种大的类别元件库，及超过 8000 种以上的具体元 件库文件添加 Proteus 中的元件并不是很全，有时需要添加第三方库文件才可进行仿真，可以通过 以下俩种方式进行添加 1、将第三方库文件拷贝至Proteus 程序目录下的LIBRARY 目录下，相应的元件模型文件也要

拷贝到 MODELS 目录下 2、将第三方库文件统一放至一个文件夹中，同时元件模型文件也要统一放至一个文件夹中，打 开Proteus 菜单 SYSTEM 下的 SET PATH?在弹出的 Path Configuration 对话框的 Library folders 中添加库文件目录,在 Simulation and folders 中添加元件模型文件目录 ⑤点击左侧工具栏按钮进入元件模式，再次点按钮，即可调出元件库 库元件分类说明： Analog Ics 模拟电路集成库 Capacitors 电容库 CMOS 4000 series CMOS 4000 库 Connectors 插座，插针，等电路接口连接库 Data Converters ADC,DAC 数 /模、模 /数库 Debugging Tools 调试工具 Diodes 二极管库 ECL 10000 Series ECL 10000 库 Electromechanica l 电机库 Inductors 电感 库 Laplace Primitives 拉普拉斯变换库 Memory ICs 存储元件库 Microprocessor ICs CPU 库 Miscellaneous 元件混合类型库 Modeling Primitives 简单模式库 Operational Amplifiers 运放库 Optoelectroni cs 光电元件库 PLDs & FPGAs 可编程逻辑器件 Resistors 电阻 库

Matlab程序设计实验报告

实验七Matlab 程序设计 实验目的： 1、掌握建立和执行M 文件的方法； 2、掌握实现选择结构的方法； 3、掌握实现循环结构的方法。 实验内容： 1. 编写用 5 次多项式拟合函数y=sin(x), x [0, 2 ]的脚本M 文件，要求绘图观察拟合的效果。 function shiyan1 x=0:0.5:2*pi y=sin(x) p=polyfit(x,y,5) x1=0:0.2:2*pi y1=polyval(p,x1) plot(x,y, 'b' ,x1,y1, '*r' x =

Columns 1 through 9 0 0.5000 1.0000 1.5000 2.0000 2.5000 3.0000 3.5000 4.0000 Columns 10 through 13 4.5000 5.0000 5.5000 6.0000 y = Columns 1 through 9 0 0.4794 0.8415 0.9975 0.9093 0.5985 0.1411 -0.3508 -0.7568 Columns 10 through 13 -0.9775 -0.9589 -0.7055 -0.2794 p = -0.0056 0.0881 -0.3967 0.2671 0.8902 0.0029 x1 = Columns 1 through 10 0 0.2000 0.4000 0.6000 0.8000 1.0000 1.2000 1.4000 1.6000 1.8000 Columns 11 through 20

2. 2.2000 2.4000 2.6000 2.8000 3.0000 3.2000 3.4000 3.6000 1.8001 Columns 21 through 30 4.0 4.2000 4.4000 4.6000 4.8000 5.0000 5.2000 5.4000 5.6000 5.8000 Columns 31 through 32 6.0 6.2000 y1 = Columns 1 through 10 0.29 0.1886 0.3786 0.5585 0.7172 0.8461 0.9391 0.9926 1.0048 0.9761 Columns 11 through 20 0.9083 0.8048 0.6701 0.5098 0.3301 0.1381 -0.0590 -0.2538 -0.4389 -0.6073 Columns 21 through 30 -0.7524 -0.8685 -0.9505 -0.9949 -0.9991 -0.9626 -0.8863 -0.7732 -0.6288 -0.4606 Columns 31 through 32

Proteus仿真软件使用方法

实验八 Proteus仿真软件使用方法 1.实验目的： （1）了解Proteus仿真软件的使用方法。 （2）了解51单片机编程器Keil与Proteus仿真软件的联用方法。 2.实验要求： 通过讲授和操作练习，学会正确使用Proteus仿真软件及Keil编程及其联合调试。 3.实验内容： （1）Proteus 仿真软件介绍 Proteus 软件是由英国LabCenter Electronics 公司开发的EDA工具软件，由ISIS和ARES两个软件构成，其中ISIS是一款便捷的电子系统仿真平台软件，ARES是一款高级的布线编辑软件。它集成了高级原理布线图、混合模式SPICE电路仿真、PCB设计以及自动布线来实现一个完整的电子设计。 通过Proteus ISIS软件的VSM(虚拟仿真技术), 用户可以对模拟电路、数字电路、模数混合电路，以及基于微控制器的系统连同所有外围接口电子元器件一起仿真。 图8-1是Proteus ISIS的编辑窗口： 图8-1 ISIS的编辑界面 图中最顶端一栏是“标题栏”，其下的“File View Edit ……”是“菜单栏”，再下面的一栏是“命令工具栏”，最左边的一栏是“模式选择工具栏”；左上角的小方框是“预览窗口”，其下的长方框是“对象选择窗口”，其右侧的大方框是“原理图编辑窗口”。 选择左侧“模式选择工具栏”中的图标，并选择“对象选择窗口”中的P按钮，就会出现如图8-2的元器件选择界面：

图8-2 元器件库选择界面 在元器件列表框中点击你需要的器件类型（例如：电阻-Resistors,单片机芯片-MicroprocessorICs, LED-Optoelectronics）或在左上角的关键字（Keywords）框中输入你需要的器件名称的关键字（如：信号源 - Clock, 运放 - CA3140等），就会在图8-2中间的大空白框列出你所需的一系列相关的元件。此时，你可用鼠标选中你要的元件，则图8-2右上角的预览框会显示你所要元件的示意图，若就是你要的元器件，则点击OK按钮，该元器件的名称就会列入位于图8-1左侧的“对象选择窗口”中（参见图1左侧下方框）。 所需元器件选择好后，在“对象选择窗口”选择某器件，就可以将它放到图8-1中的“原理图编辑窗口”中（若器件的方向不合适，你可以利用图1左下角的旋转按钮来改变它）。将所要的元器件都选好后，将它们安放到合适的位置，就可以用连接线把电路连接好，结果存盘（请按规定的目录存盘，并记住其路径/目录/文件名[学号-实验序号]）。 （2）51单片机编程器– Keil V3的使用 Keil编程器可用于MCS-51单片机软件编程与调试，它的工作界面如图8-3所示： Keil编程器是Keil Software Inc/Keil Electronic GmbH 开发的基于80C51内核的微处理器软件开发平台，可以完成从工程（Project）的建立和管理、程序的编译和连接、目标代码的形成、软件仿真等一套完整的软件开发流程。它与Proteus挂接，可以进行单片机应用系统的硬件仿真。 汇编语言编程方法： ①打开“File”菜单→选择新建“New...”→在弹出的文本框（Text1）中编写所需的汇编语言程序→程序写好后，保存（从File→Save As..→选择某目录，文件名.ASM, 存盘）； ②打开“Project”菜单→选择新建工程“New Project...”→在弹出的窗口填写：工程名→保存（文件名的后缀是 .uv2 。此时图8-3的工程窗口中将建立Target1及 Source Group 1）；

matlab实验报告

Matlab实验报告 ——定积分的近似计算 学生姓名： 学号： 专业：数学与应用数学专业

数学实验报告 实验序号:1001114030 日期：2012年10月20日 班级应一姓名陈璐学号1001114030 实验名称：定积分的近似运算 问题背景描述： 利用牛顿—莱布尼茨公式虽然可以精确地计算定积分的值，但它仅适合于被积分函数的原函数能用初等函数表达出来的情形。如果这点办不到或不容易办到， 这就有必要考虑近似计算的方法。在定积分的很多应用问题中，被积函数甚至没 有解析表达式，可能只是一条实验记录曲线，或者是一组离散的采样值，这时只 能应用近似方法去计算相应的定积分。 实验目的： 本实验将主要研究定积分的三种近似计算算法：矩形法、梯形法、抛物线发。对于定积分的近似数值计算，Matlab有专门函数可用。 实验原理与数学模型： 1.sum（a）：求数组a的和。 2.format long：长格式，即屏幕显示15位有效数字。 3.double（）：若输入的是字符则转化为相应的ASCII码；若输入的是整型数之则转化为 相应的实型数值。 4.quad（）：抛物线法求数值积分。格式：quad（fun，a，b）。此处的fun是函数，并且

为数值形式，所以使用*、/、^等运算时要在其前加上小数点。 5.trapz():梯形法求数值积分。格式：trapz（x，y）。其中x为带有步长的积分区间;y为数 值形式的运算。 6.fprintf（文件地址，格式，写入的变量）：把数据写入指定文件。 7.syms 变量1变量2……：定义变量为符号。 8.sym（'表达式'）：将表达式定义为符号。 9.int（f，v，a，b）：求f关于v积分，积分区间由a到b。 10.subs（f，'x'，a）：将a的值赋给符号表达式f中的x，并计算出值。若简单地使用subs （f），则将f的所有符号变量用可能的数值代入，并计算出值。 实验所用软件及版本：Matlab 7.0.1

金融MATLAB实验报告三答案详解

安徽财经大学金融证券实验室实验报告 实验课程名称《金融》TLABMA 金融学院部课系开

级班 学号 姓名 师导指教日年月 1．

2 一、期权定价分析 1.black-scholes方程求解

例1：假设欧式股票期权，六个月后到期，执行价格90元，现价为102元，无股利支付， 股价年化波动率为55%，无风险利率为8%，计算期权价格。 解：clear Price=102; >>Strike=90; >>Rate=0.08; >>Time=6/12; >>V olatility=0.55; [CallDelta,PutDelta]=blsprice(Price,Strike,Rate,Time,V olatility) 计算结果： CallDelta= 23.5648 PutDelta= 8.0358 2.期权价格与波动率关系分析 Price=102; >>Strike=90; >>Rate=0.08; >>Time=6/12; V olatility=0.08:0.01:0.5; >>N=length(V olatility) Call=zeros(1,N); Put=zeros(1,N); for i=1:N [Call(i),Put(i)]=blsprice(Price,Strike,Rate,Time,V olatility(i)); N= 43 end plot(Call,'b--'); hold on plot(Put,'b'); xlabel('V olatility') ylabel('price') legend('Call','Put')