Matlab作业

MATLAB作业

学号：07082002 姓名：武婷

第一次作业：

1) 熟悉Matlab运行环境;简介一下你都做了什么!

首先，安装了Matlab软件，并启动该软件。

Matlab启动后，即出现Desktop操作桌面，为一个高度集成的工作界面。同时发现有四个界面：指令窗（command window),当前目录浏览器（current

directory），历史指令窗（command history）,工作空间浏览器（workspace brower）。

指令窗（command window)：在该窗内，进行了简单的算术运算，矩阵的不同输入方式的体验，既简单的函数调用等。如在命令窗口中键入一些简单的命令：rand（4，4）生成4*4其元素在0-1之间的随机矩阵，norm(x,1)求向量的1-范数。

再如输入指令：>> 7*(12^2+66)/21，按下回车键即可得到结果： ans = 5.6084e+070

当前目录浏览器：为管理方便，首先建立用户目录并将它设置成当前目录，了解到了有两种方法;交互界面设置法和指令设置法。同时基于对该软件搜索路径知识的学习利用设置路径对话框修改搜索路径，并尝试了两种修改状态：当前

1、 与其他计算机语言相比较，MATLAB语言突出的特点是什么？ 2、 存储在工作空间中的数组能编辑吗？如何操作？

3、 命令历史窗口除了可以观察前面键入的命令外，还有什么用途？ 4、 Matlab提供了两种运行方式，它们是？ 5、 在MATLAB中有几种获得帮助的途径？

6、 有几种建立矩阵的方法？ 7、 在MATLAB中如何建立矩阵??5?4793??1?，并将其赋予变量a？

8、 试用matlab语句生成一个6行5列的矩阵，要求其中的元素全为0 9、 matlab变量命名的规则是？

9、 下列符号哪些是matlab的合法变量名？

3vars、global、help、My_exp、sin、X+Y、_input、 E-4 、AbCd 、 AB_C_

10、若不想显示命令运行的结果，应怎么做？若要注释掉一行，应该怎么做？ 11、请说出下面几个命令的含义：clc，clear, which 12、数组运算和矩阵运算的运算符有什么区别？ 13、求x?4?8i???3?2i3?5i7?6i2?7i9?4i1?4i3?9i7?5i??4?4i?的共轭转置

14、请用冒号语句生成行向量，要求其内容为1到100 ?5?15、有矩阵a?3???73793?

多幅度电平信号差错率仿真设计

一、设计任务

完成多幅度电平信号差错率的仿真，运用MATLAB语言，仿真它们的误码性能，并与理论误码曲线相比较。先产生多组符号序列，再将符号序列映射到对应的幅度电平上，以四电平为例，先产生一个四组符号的序列，将该序列映射到对应的幅度电平上。随机变量信道中的加性噪声可以认为是统计独立，均值为0的高斯分量，分别作出二、四、八、十六电平仿真得到的误码率曲线与和理论上的误码率曲线。

二、物理模型

多幅度电平信号仿真系统模型如下图所示：

出错计数器Am与Am比较^均匀随机数发生器X映射到幅度电平Am+r2检测器(0,?)高斯随机数发生器Am^

图1

三、预期目标

通过以上的物理模型，作出二、四、八、十六电平信号仿真得到的误码率

曲线与和理论上的误码率曲线。 四、系统中所涉及的公式如下：

1、平均比特SNR定义为：

即可得：

1d2 二电平时平均比特SNR为：SNR?

2?25d2 四电平时平均比特SNR为：SNR? 24?21d2 八电平时平均比特SNR为：SNR?

6?2

十六电平时平均比特SNR为:

SNR?85d28?2

2、M电平PAM系统最佳

MATLAB上机作业1

对以下问题,编写M文件：

(1) 用起泡法对10个数由小到大排序。即将相邻两个数比较,将小的调到前头。 function f=qipaofa(x) for j=9:-1:1 for i=1:j

if(x(i)>x(i+1))

t=x(i);x(i)=x(i+1);x(i+1)=t; end end end f=x

x=round(10*rand(1,10)) qipaofa(x);

(2) 有一个4×5矩阵，编程求出其最大值及其所处的位置。 function f=zuidazhi(x) a=1;b=1;c=x(1,1); for i=1:4 for j=1:5

if x(i,j)>c

a=i;b=j;c=x(i,j); end end end

f=[c,a,b]

x=rand(4,5) zuidazhi(x)

20(3) 编程求?n!。

n?1function f=qiuhe(x) b=0; for i=1:x a=prod(1:i); b=b+a; end f=b

qiuhe(20)

(4)一球从100米高度自由落下,每次落地后

梧州学院

MATLAB与系统仿真

专 业：机械设计制造及其自动化班 级：学 号：姓 名：韦邮 箱：指导老师:

09机械5班 0901901012 靖

wei_jing1988@163.com 姚 金 2012年6月27日

1

一、春天是鲜花的季节，水仙花就是其中最迷人的代表，数学上有个水仙花数，它是这样定义的：“水仙花数”是指一个三位数，它的各位数字的立方和等于其本身，比如：153=1^3+5^3+3^3。要求编m文件输出所有存在的水仙花数。（本题10分） 解：

题意中“水仙花数”是指一个三位数说明100—999之间的数，在MATLAB编程M文件，运行可得。 程序如下： for i=100:999; a=floor(i/100); b=floor(i/10)-a*10; c=i-(a*100+b*10); if i==a^3+b^3+c^3; i end end 运行结果：

本程序源程序附件名为：

2

二、带饱和输出特性PID控制系统的被控对象为G0(s)=控制器为Gc(s)=Kp(1?1e?60s，其PID80s?1110.048?TDs)=K

《Matlab编程训练》

作业

专

业

学生姓名 班学

级 号

指导教师 完成日期

实训一 MATLAB语言介绍和数值计算

1．先求下列表达式的值，然后显示MATLAB工作空间的使用情况并保存变量。

2sin850z1?. 21?e

?1234?4??13?1??,B??203?，求下列表达式的值： 347872. 已知 A?????????3657???3?27??(1) A+6*B和A-B+I（其中I为单位矩阵） A+6*B:

A-B+I:

(2) A*B和A.*B A*B程序：

A=[12 34 -4;34 7 87;3 65 7] B=[1 3 -1;2 0 3;3 -2 7] c=A*B

结果：

A.*B程序：

A=[12 34 -4;34 7 87;3 65 7] B=[1 3 -1;2 0 3;3 -2 7] D=A.*B 结果：

(3) A^3和A.^3 A^3程序：

A=[12 34 -4;34 7 87;3 65 7] E=A^3 结果：

A．^3程序：

A=[12 34 -4;34 7 87;3 65 7] C=A.^3

(4) A/B及B\\A A/B程序：

A=[12 34 -

MATLAB

1、 试求出下面线性微分方程的通解。

dy(t)dt55作业5

?13dy(t)dt44?64dy(t)dt33?152dy(t)dt22?176dy(t)dt?80y(t)?e?2t[sin(2t??3)?cos(3t)]&&假设上述微分方程满足已知条件，y(0)?1,y(1)?3,y(?)?2,y(0)?1,y(1)?2试

求出满足该条件的微分方程的解析解。

解： >> syms t y ;

y=dsolve(['D5y+13*D4y+64*D3y+152*D2y+176*Dy+80*y=','exp(-2*t)*(sin(2*t+pi/3)+cos(3*t))'],'y(0)=1','y(1)=3','y(pi)=2','Dy(0)=1','Dy(1)=2'); vpa(y,20) ans =

.20576131687242798354e-2*exp(-2.*t)*cos(3.*t)+.15538705805619602373e-1*exp(-2.*t)*sin(2.*t)+.76830587084294035587e-2*exp(-2.*t)*cos(2.*t)+98.

MATLAB

1、 试求出下面线性微分方程的通解。

dy(t)dt55作业5

?13dy(t)dt44?64dy(t)dt33?152dy(t)dt22?176dy(t)dt?80y(t)?e?2t[sin(2t??3)?cos(3t)]&&假设上述微分方程满足已知条件，y(0)?1,y(1)?3,y(?)?2,y(0)?1,y(1)?2试

求出满足该条件的微分方程的解析解。

解： >> syms t y ;

y=dsolve(['D5y+13*D4y+64*D3y+152*D2y+176*Dy+80*y=','exp(-2*t)*(sin(2*t+pi/3)+cos(3*t))'],'y(0)=1','y(1)=3','y(pi)=2','Dy(0)=1','Dy(1)=2'); vpa(y,20) ans =

.20576131687242798354e-2*exp(-2.*t)*cos(3.*t)+.15538705805619602373e-1*exp(-2.*t)*sin(2.*t)+.76830587084294035587e-2*exp(-2.*t)*cos(2.*t)+98.

基于MATLAB识别图片中的文字 知识点：图像识别 图像处理 计算机视觉

实验目的：以MATLAB为工具，识别图片中的文字。 问题描述：

交通是一个十分重要的问题，红绿灯上的摄像头可以辨别车主的身份。现在，我们是否能够用matlab设计一个图像识别的系统，通过车牌号的来识别相应的数字。 问题分析及模型建立： 第一步：确定车牌号的区域 对于一张图像

在matlab中，在处理图像元素时用（x,y,z）三坐标的形式表达的灰度集时，其中，x和y是空间坐标，z是f 在坐标（x,y）处的值，就可以表达图像在该坐标轴上的点。再将，将f(x,y)的数值简单地显示为一个矩阵，就可以定量地表达了一幅数值图像。

矩阵中的每个元素称为像素。

所以，如果我们想要截取一段车牌号的图像，我们就可以扫描图像每一个点的像素。然后我们进行恢复处理，算子边缘检测，腐蚀，闭操作，删除小对象得到

我们就可以利用RGB的值找到白色区域的边界，然后我们利用边界，截取我们的原图像就可以得到我们所需要的车牌号的区域了。

第二步：进行文字的分割

在进行分割之前，就需要将我们所获得的车牌号的区域的图像进行进一步的处理，突出我们需要的文字部分。 灰度处理

彩色的图片占用的空间比较大，处理会花费很长

第四次 上机作业

1、 从键盘输入一个4位整数，按照如下规则加密后输出。加密规则：每位数字

都加上7，然后用和除以10的余数取代该数字；再把第一位与第三位交换，第二位与第四位交换。 Clear X=ones(1,4);

X（1）=input(’输入第一位：‘)； X（2）=input(’输入第二位：‘)； X（3）=input(’输入第三位：‘)； X（4）=input(’输入第四位：‘)； X=rem(7+x,10);

Y=1000.*x(3)+100.*x(4)+10.*x(1)+x(2)

2、 分别用if 和switch语句实现以下计算，其中a、b、c的值从键盘输入。

??ax2?bx?c,?y??asincb?x,?c?lnb?,x?0.5?x?1.51.5?x?3.5 3.5?x?5.5a=input(‘请输入a：’)； b=input(‘请输入b：’)； c=input(‘请输入c：’)； If(x>=0.5&&x<=1.5) y=a.*x^2+b.*x+c Elseif(x>=1.5&&x<=3.5) y=a.*(sin(b))^c+x

Elseif(x>=3.5&&x<=5.5) y=log(abs(b+c./x)) end

a=input(‘请输入a：’)； b=input(‘请输入b：’)； c=input(‘请输