第四章 MATLAB 仿真系统

通信电路仿真

1 Matlab的基本操作

一、使用函数实现对下列矩阵的左旋和右旋以及反转 已知

答案：如图：矩阵如图

矩阵的左旋如图一 矩阵的右旋如图二 矩阵的左右反转如图三 矩阵的上下反转如图四 二、已知A=[8 9 5 ] B=[-1 3 -2] [36 -7 11] [2 0 3] [21 -8 5] [-3 1 9 ] 计算： 1) A+5*B; 2) A*B和A.*B 3) A^3和A.^3 4) A/B和B\A 5)

[A,B]

通信电路仿真

2 答案：A=[8 9 5；36 -7 11；21 -8 5]; B=[-1 3 -2；2 0 3； -3 1 9]; A+5*B A*B A.*B A^3 A.^3 A/B A\B [A,B] 程序运行结果如下：

通信电路仿真

3

第二章、 Matlab程序设计

一、已知S=1+2+22+23+…+263 求S的值 答案：代码如下： Clear Clc s=0,j=2 for i=1:63 s=s+j^i end s 程序运行结果如下：

二、

通信电路仿真

4 第三章、 Matlab绘图

一、在同一坐标系中画出下列函数的图像：x2，-x2，xsin（x）在[0,2π]上的函数图像 程序代码如下： clear clc x=0:pi/100:2*pi; y1=x.^2; y2=-x.^2; y3=x.*sin(x); title('同一坐标下的函数图像') plot(x,y1,':',x,y2,'h',x,y3,'--') 程序运行结果如下：

二、绘制极坐标图像： 程序代码如下： clear clc x=0:0.01:2*pi; y=2*cos(2*(x-pi/8)); title('极坐标图像') xlabel('x') ylabel('y') polar(x,y) 程序运行结果如下：

通信电路仿真

5

第四章、 Matlab符号运算

matlab讲义

随机过程实验讲义

刘继成

华中科技大学数学与统计学院

前言 (1)

第一章Matlab 简介 (2)

第二章简单分布的模拟 (6)

第三章基本随机过程 (9)

第四章Markov过程 (12)

第五章模拟的应用和例子 (16)

附录各章的原程序 (51)

参考文献 (75)

若想检验数学模型是否反映客观现实，最自然的方法是比较由模型计算的理论概率和由客观试验得到的经验频率。不幸的是，这两件事都往往是费时的、昂贵的、困难的，甚至是不可能的。此时，计算机模拟在这两方面都可以派上用场：提供理论概率的数值估计与接近现实试验的模拟。

模拟的第一步自然是在计算机程序的算法中如何产生随机性。程序语言，甚至计算器，都提供了“随机”生成[0，1]区间内连续数的方法。因为每次运行程序常常生成相同的“随机数”，因此这些数被称为伪随机数。尽管如此，对于多数的具体问题这样的随机数已经够用。我们将假定计算机已经能够生成[0，1]上的均匀随机数。也假定这些数是独立同分布的，尽管它们常常是周期的、相关的、……。

……

本讲义的安排如下，第一章是Matlab简介，从实践动手角度了解并熟悉Matlab环境、命令、帮助等，这将方便于Matlab的初学者。第二章是简单随机变量的模拟，只给出了常用的Matlab 模拟语句，没有堆砌同一种变量的多种模拟方法。对于没有列举的随机变量的模拟，以及有特殊需求的读者应该由这些方法得到启发，或者参考更详细的其他文献资料。第三章是基本随机过程的模拟。主要是简单独立增量过程的模拟，多维的推广是直接的。第四章是Markov过程的模拟。包括服务系统，生灭过程、简单分支过程等。第五章是这些模拟的应用。例如，计算概率、估计积分、模拟现实、误差估计，以及减小方差技术，特别给读者提供了一些经典问题的模拟，通过这些问题的模拟将会更加牢固地掌握实际模拟的步骤。平稳过程的模拟、以及利用平稳过程来预测的内容并没有包含在本讲义之内，但这丝毫不影响该内容的重要性，这也是将会增补进来的主要内容之一。希望读者碰到类似的问题时能够查阅相关资料解决之。

数值分析第四章外推法计算数值微分MATLAB计算实验报

告

数值分析MATLAB计算实验报告

姓名班级学号

⼀、实验名称

⽤MATLAB编程实现数值微分的外推法计算。

⼆、实验⽬的1．掌握数值微分和定义和外推法的计算过程；

2．了解数值微分外推法的计算⽅法并且编写出与其算法对应的MATLAB程序代码；3．体会利⽤MATLAB软件进⾏数值计算。

三、实验内容

⽤外推法计算f(x)=x2e?x在x=0.5的导数。

四、算法描述1.命名函数。

2.如果输⼊未知数少于四个，默认精度10^-3

3.描述T表矩阵坐标

4.依次赋值计算 T表第⼀列

5.根据数值微分计算公式求出T表矩阵的值

6.若达到精度则运算结束，若未达到循环计算

7.输出T表，得出的值就是导数值五、实验结果

六、实验结果分析

此实验通过MATLAB实现外推法数值微分计算，得到相应的数据，⽅便对数据进⾏分析。从结果可以看出，当步长h=0.025时⽤中点微分公式只有3位有效数字，外推⼀次达到5位有效数字，外推两次达到9位有效数字。

七、附录（程序）function g=waituifa(fname,x,h,e)if nargin<4,e=1e-3;

end;

i=1;

j=1;

G(1,1)=(feval(fname,x+h)-feval(fname,x-h))/(2*h);

G(i+1,1)=(feval(fname,x+h/2)-feval(fname,x-h/2))/h;

G(i+1,j+1)=(4^j*G(i+1,j)-G(i,j))/(4^j-1);

while abs(G(i+1,i+1)-G(i+1,i))>e

i=i+1;

G(i+1,1)=(feval(fname,x+h/2^i)-feval(fname,x-h/2^i))/(2*h/2^i);

for j=1:i

G(i+1,j+1)=((4^j)*G(i+1,j)-G(i,j))/(4^j-1);

end

end

G

g=G(i+1,i+1);

第 1 页 （共 5 页） 学习中心/函授站_ 成都学习中心

姓 名 赵洪 学 号 7020140122093

西安电子科技大学网络与继续教育学院

2015学年上学期

《MATLAB与系统仿真》期末考试试题

（综合大作业）

题号 一 二 三 总分

题分 25 12 63

得分

考试说明：

1、大作业于2015年4月3日公布，2015年5月9日前在线提交；

2、考试必须独立完成，如发现抄袭、雷同、拷贝均按零分计。

3、程序设计题（三（8，10））要求写出完整的程序代码，并在matlab软件环境调试并运行通过，连同运行结果一并附上。

一、填空题（1’ ×25=25’）

1、Matlab的全称为 矩阵实验室 。

2、在Matlab编辑器中运行程序的快捷键是：

F5 。

3、Matlab的工作界面主要由以下五个部分组成，它们分别是： 菜单栏

、

工具栏 、 当前工作目录窗口

、

工作空间管理窗口 和

命令窗口

。

4、在Matlab中inf表示： 无穷大

；clc表示： 清空命令窗口中的显示内容 ；more表示：

在命令窗口中控制其后每页的显示内容行数 ；who表示：

查阅Matlad内存变量名 ；whos表示： 列出当前工作空间所有变量

。

5、在Matlab命令窗口中运行命令

Simulink 可以打开Simulink模块库浏览器窗口。

6、求矩阵行列式的函数： det ；求矩阵特征值和特征向量的函数 eig

。

7、Matlab预定义变量ans表示： 没有指定输出变量名 ；eps表示： 系统精度 ；nargin表示： 函数输入参数的个数 。