当前位置：文档之家› 西安邮电数字信号处理实验报告

西安邮电数字信号处理实验报告

西安邮电大学

通信与信息工程学院

实验报告

（2016/2017学年第2学期）

实验名称：数字信号处理实验

学生姓名：石龙飞

专业：通信工程

学号：03151307

指导教师：常虹

2017年06月01日

西安邮电大学通院《数字信号处理实验》独立设课实验过程考核表

西安邮电大学通院《数字信号处理实验》独立设课实验成绩鉴定表

1 熟悉MATLAB

1.1实验目的

1、了解MATLAB程序设计语言的基本特点，熟悉MATLAB软件的运行环境。

2、掌握各种信号的建模方式。

3、掌握各种信号的托兴表示方法。

4、掌握变量等有关概念，具备初步的将一般数学模型转化为对应的计机模型并进行处理的能力。

1.2 实验原理

1、图形保持

hold on/off命令控制是保持原有图形还是刷新原有图形，不带参数的hold命令在两种状态之间进行切换。

2、绘图修饰命令

title(图形名称)

xlabel(x轴说明)

ylabel(y轴说明)

text(x,y,图形说明)

legend

1.3 实验内容以及步骤

单位采样序列

代码：

单位阶跃序列代码：

矩形序列

代码：

三角波

代码：

方波

代码：

锯齿波代码：

Sinc函数代码：

1.4 实验结果分析

学习使用MATLAB7.0软件；学习信号的图形表示方法，掌握各种各种信号的的建模方式；实现单位采样序列、单位阶跃序列、举行序列、三角波、方波、锯齿波、Sinc函数。

1.5 实验总结以及体会构成

经过本次试验，我对于用MATLAB来作出一些数字信号中的一些基本波形的图像有了深刻的认识，可以熟练而又准确地以MATLAB为工具做出诸如单位阶跃序列，矩形序列等多种序列的图形，不仅如此，在作图中，我对于MATLAB的一些作图的基本指令，以及需要作出的序列有了深刻的认识与见解；本节课做出这些图形并不是重点，重点是了解MATLAB程序设计语言的基本特点，掌握这一个重要的工具，并且可以将数学模型熟练地转化为计算机模型，并进行处理，这才是本次实验真正的作用与意义，也是作为当代大学生的我们所必须具备的基本能力。

1.6 指导老师评语

2信号的表示

2.1实验目的

1、了解MATLAB 程序设计语言的基本特点，熟悉MATLAB 软件运行

环境。

2、掌握各种信号的建模方式。

3、掌握各种信号的图形表示方法。

4、掌握变量等有关概念，具备初步的将一般数学模型转化为对应的计算机

模型并进行处理的能力

2.2实验原理

常用的 MATLAB 绘图语句有figure、plot、subplot、stem 等，图形修饰语具有title、axis、text 等。

2.2.1 figure 语句

figure 有两种用法。当只有一句figure 命令时，程序会创建一个新的图形窗口，并返回一个整数型的窗口编号。当采用figure(n)时，表示将第n 个图形窗口作为当前的图形窗口，将其显示在所有窗口的最前面。如果该图形窗口不存在，则新建一个窗口，并赋以编号n。

2.2.2 plot 语句

线形绘图函数。用法为plot(x,y,’s’)。参数x 为横轴变量，y 为纵轴变量，s 用以控制图形的基本特征如颜色、粗细等，通常可以省略，常用方法如表1-1 所示。

表 1-1 plot 命令的参数及其含义

2.2.3 subplot 语句

subplot(m,n,i)是分割显示图形窗口命令，它把一个图形窗口分为m 行n 列共m×n 个

小窗口，并指定第i 个小窗口为当前窗口。

2.2.4 二维统计分析图

在MATLAB 中，二维统计分析图形很多，常见的有条形图、阶梯图、杆图和填充图等，所采用的函数分别是：

bar(x,y,选项)

stairs(x,y,选项)

stem(x,y,选项)

fill(x1,y1,选项1,x2,y2,选项2,…)

2.2.5图形保持

hold on/off 命令控制是保持原有图形还是刷新原有图形，不带参数的hold 命令在两种状态之间进行切换。

2.2.6 绘图修饰命令

title(图形名称)

xlabel(x 轴说明)

ylabel(y 轴说明)

text(x,y,图形说明)

legend(图例1,图例2,…)

2.2.7 MATLAB 常用信号生成函数

（1） ZEROS

? 功能：产生全零阵列

? 调用格式：X=ZEROS (N) %产生N 行N 列的全零矩阵

? X=ZEROS(M,N) %产生M 行N 列的全零矩阵

（2） ONES

? 功能：产生全 1 阵列

? 调用格式：X=ONES(N) %产生N 行N 列的全1 矩阵

? X=ONES(M,N) %产生M 行N 列的全1 矩阵

（3） SINC

? 功能：辛格函数

? 调用格式：Y= SINC(X) %

（4） RECTPULS

? 功能：产生矩形脉冲信号

? 调用格式：Y= RECTPULS (T) %产生高度为1、宽度为1、关于T＝0 对称的矩形脉冲? Y= RECTPULS (T，W) %产生高度为1、宽度为W、关于T＝0 对称的矩形脉冲

（5） RAND

? 功能：产生伪随机序列

? 调用格式：Y= RAND (1,N) %产生[0,1]上均匀分布的随机序列

? Y= RANDN (1,N) %产生均值为0，方差为1 的白噪声序列

（6） SAWTOOTH

? 功能：产生周期锯齿波或三角波

? 调用格式：Y= SAWTOOTH (T) %产生幅值为＋1，－1，以2 为周期的方波

? Y=SAWTOOTH(T,WIDTH)% 产生幅值为＋1，－1，以WIDTH *2 为周期的方波

（7） SQUARE

? 功能：产生方波

? 调用格式：Y= SQUARE (T) %产生幅值为＋1，－1，以2 为周期的锯齿波

? Y= SQUARE(T,DUTY) % 产生幅值为＋1，－1，以占空比为DUTY 的方波

? 例： t = 0:.0001:.0625;

? y = SQUARE(2*pi*30*t，80); plot(t,y)％产生一个占空比为80%的方波

（8） FLIPLR

? 功能：序列左右翻转

? 调用格式：Y = FLIPLR(X)

? % X = 1 2 3 翻转后 3 2 1

4 5 66 5 4

（9） CUMSUM 、SUM

? 功能：计算序列累加

? 调用格式：Y = CUMSUM(X) % 向量X 元素累加，记录每一次的累加结果,而SUM只记录最后的结果。

2.3：实验内容以及步骤

1、分别以条形图、阶梯图、杆图和填充图形式绘制曲线y=2sin(x)。

程序如下：

2、采用图形保持，在同一坐标内绘制曲线y1=0.2e-0.5xcos(4πx) 和y2=2e-0.5xcos(πx)。

程序如下：

3、在0≤x≤2π区间内，绘制曲线y1=2e-0.5x 和y2=cos(4πx)，并给图形添加图形标注。

程序如下：

4、实现单位采样序列δ(n)、单位阶跃序列u(n)、矩形序列( ) N R n ，并用图形显示。写出程序及输出图形

5、实现三角波、方波、锯齿波、Sinc 函数，并用图形显示。写出程序及输出图形

2.4：实验结果分析

1、分别以条形图、阶梯图、杆图和填充图形式绘制曲线y=2sin(x)。

2、采用图形保持，在同一坐标内绘制曲线y1=0.2e-0.5xcos(4πx) 和

y2=2e-0.5xcos(πx)。

3、在0≤x≤2π区间内，绘制曲线y1=2e-0.5x 和y2=cos(4πx)，并给图形添加图标注。

4、实现单位采样序列δ (n)、单位阶跃序列u(n)、矩形序列 ( ) N R n ，并用图形显示，写出程序及输出图形。

5、实现三角波、方波、锯齿波、Sinc 函数，并用图形显示。写出程序及输出图形

2.5：实验总结以及体会构成

本次实验主要是熟悉了如何使用MATLAB表示信号，通过MATLAB的函数来实现是主要方法，其次也可以自行编写函数，表示自己需要的特定信号，来适应不同的应用场合。在信号表示中学会画图是非常重要的一件事，几乎所有的信号都可以也需要用图像表示出来，这样对分析和实验都是有着好处的，熟练掌握画图函数，这些便都可以迎刃而解，画图函数有很多种，其参数也是形形色色，记住每个参数的位置意义是非常重要的。

因为是对数字信号的处理，因此免不了会用到很多的理论知识，这时就是理论和时间结合的时候，MATLAB是一个可以仿真信号的工具，可以模拟实际，需要用理论知识来设计解决问题所需要的模型，然后用MATLAB将其仿真出来，看看实际结果是否是我们预期的，这个结果也将是我们在现实中实施该方案所能看到的结果。

2.6：指导老师评语

3 FFT 频谱分析及应用

3.1实验目的

1、通过实验加深对FFT 的理解；

2、熟悉应用FFT 对典型信号进行频谱分析的方法。

3.1实验原理

在各种信号序列中，有限长序列占重要地位。对有限长序列可以利用离散傅立叶变换(DFT)进行分析。DFT 不但可以很好的反映序列的频谱特性，而且易于用快速算法(FFT)在计算机上进行分析。

有限长序列的DFT 是其z 变换在单位圆上的等距离采样，或者说是序列傅立叶的等距离采样，因此可以用于序列的谱分析。FFT 是DFT 的一种快速算法，它是对变换式进行一次次分解，使其成为若干小数据点的组合，从而减少运算量。在MATLAB 信号处理工具箱中的函数fft（x,n）,可以用来实现序列的N 点快速傅立叶变换。

经函数fft求得的序列一般是复序列，通常要求出其幅值和相位。MATLAB中提供了求复数的幅值和相位的函数：abs、angle，这些函数一般和fft同时使用。

3.3.3实验内容以及步骤

被噪声污染得信号，比较难看出所包含得频率分量，如一个由50Hz 和120Hz正弦信号构成的信号，受零均值随机噪声的干扰，数据采样率为1000Hz，试用FFT 函数来分析其信号频率成分，要求：○1 画出时域波形；○2 分析信号功率谱密度。

3.3.4实验结果分析

3.3.5实验总结以及体会构成

本次实验只要是针对有限长序列，我们可以通过MATLAB对其进行离散傅里叶变化，是非常方便的，因为工具内部已经将函数写好，我们只需要将变化的信号交给函数就可以，可以说大大地减小了手工计算量，并且能过直观得将图画出来。

其中需要注意的问题也有很多，比如因为操作失误，将函数名称写错，软件机会报错，这个问题还是比较容易查找的，其次是在表示信号的时候，将一些数值写错，到会最后的变换不是我们需要的，画出的图像有明显的问题，这些都需要细心才能避免，当然还有一些比较头疼的错误，比如函数不会使用，所给的参数错误，或者不符合标准。

3.3.6指导老师评语

4 信号的运算-卷积

4.1实验目的

1、掌握信号的线性卷积运算。

2、掌握信号的循环卷积运算。

3、掌握信号循环卷积计算线性卷积的条件。

4.2实验原理

数字信号卷积定理

4.3 实验内容以及步骤

1、假设卷积下面信号

X(n)=9.0n 0≤n<13

0 else

h(n)= 1 0≤n<12

0 else

选定循环卷积的长度为N = 21。确定1y (n) = x(n)?h(n)的哪些数值与线性卷积2y (n) = x(n) *h(n)结果中的数值相同。编写程序代码并输出图形，并分析错误数据的原因，怎样才能使两者数据相同。

西安交通大学接口技术实验报告

西安交通大学微型计算机接口技术实验报告班级:物联网姓名: 学号：

实验一基本I/O扩展实验一、实验目的 1、了解 TTL 芯片扩展简单 I/O 口的方法，掌握数据输入输出程序编制的方法； 2、对利用单片机进行 I/O 操作有一个初步体会。二、实验内容 74LS244 是一种三态输出的8 总线缓冲驱动器，无锁存功能，当G 为低电平时，Ai 信号传送到Yi，当为高电平时，Yi 处于禁止高阻状态。 74LS273 是一种8D 触发器，当CLR 为高电平且CLK 端电平正跳变时，D0——D7 端数据被锁存到8D 触发器中。实验原理图：三、实验说明利用74LS244 作为输入口，读取开关状态，并将此状态通过74LS273 再驱动发光二极管显示出来，连续运行程序，发光二极管显示开关状态。

四、实验流程图五、实验连线 1、244的cs连接到CPU地址A15，Y7—Y0连接开关K1-K8; 2、273的CS连接到CPU地址A14，Q7-Q0连接到发光二极管L1-L8; 3、该模块的WR,RD连接CPU的WR,RD,数据线AD7-AD0，地址线A7-A0分别与CPU的数据线AD7-AD0，地址线A7-A0相连接。

六、程序源代码（略）七、实验结果通过开关K01 到K08 可以对应依次控制LED 灯的L1 到L8 ，即当将开关Ki 上拨时，对应的Li 被点亮，Ki 下拨时，对应的Li熄灭。此外，如果将开关拨到AAH 时，将会产生LED 灯左移花样显示；如果开关拨到55H 时，将会产生LED 灯右移花样显示。七、实验心得通过本次实验，我了解了TTL 芯片扩展简单I/O 口的方法，同时也对数据输入输出程序编制的方法有一定的了解与掌握，对利用单片机进行I/O 操作有一个初步体会，实验使我对自己在课堂上学的理论知识更加理解，同时也锻炼了我的动手操作能力。

西安邮电大学课程设计报告

一. 实验目的及实验环境实验目的：目前，对于中小超市来说，需要采取一套行之有效的超市管理系统来管理超市的日常进货、销售、员工管理等信息。本次设计的超市管理系统就是想让用户以最少的手工业务处理，来完成超市中的日常工作，增加超市业务处理的效率。一个典型的超市管理系统，除了应该完成超市中必要的销售工作，如收银台，还应该具有后台管理的功能。实验环境：系统的体系结构：C/S结构前台：Java (Eclipse 3.5) 后台数据库：SQL Server 2008 二. 实验内容超市管理系统该系统应具备前台收银操作，还有后台对各种数据的管理。该系统应为多用户使用的，应以用户的权限来进入系统管理模块。友好的界面及方便简单的操作。数据存放在数据库中，系统对数据库可进行增、删、查、改的操作。三．方案设计系统开发的总体任务是实现Java与SQL Server 2008的组合应用，实现用JDBC 方式调用SQL数据库，且有一定的透明性。超市管理系统需要完成的功能主要有： 1、用户分类登录系统及权限管理。 2、超市收银台的销售管理方式。 3、销售商品信息的录入、查询、修改、删除。 4、供应商信息的录入，查询，修改、删除。 5、超市员工的职位基本信息的设定。 6、员工的基本信息的录入，查询，修改和删除。 7、仓库信息的录入与查看。 8、数据间的相互联系：销售商品的数量和库存量的基本联系；商品信息和供应商的相关数据联系；员工与职位间的联系等。 9、数据存放在数据库中，任何数据操作都是对数据库的数据的操作。

系统结构图（功能模块图）各个模块的功能描述： (1)、登录界面：本系统简单的设置有三种职位，用户使用密码登录，成功后给其相应的操作界面。 (2)、仓库管理：是仓库管理员和老板可使用的操作，主要是查看仓库的入库信息和添加入库信息。 (3)、管理商品信息：只有老板可用，主要是查看现有商品的信息、修改、删除及添加新的商品信息。 (4)、管理供应商信息：只有老板可用，主要是查看现有的供应商的信息、修改、删除及添加新供应商的信息。 (5)、管理员工信息：只有老板可用，主要是查看现有的员工信息、修改、删除及增加新员工。 (6)、查看销售纪录。四．测试数据及运行结果 1、Java与数据库SQL 2008的链接： String driverName="com.microsoft.sqlserver.jdbc.SQLServerDriver"; String dbURL="jdbc:sqlserver://localhost:1433;DatabaseName=SuperMarket"; 2、用户登陆(登陆界面，老板界面，仓库管理员界面，销售员界面)：

某厂排气管车间生产计划的优化分析实验报告

————————————————————————————————作者: ————————————————————————————————日期:

西安理工大学实验报告成绩第１页(共页) 课程: 实验日期: 年月日专业班号：组别：交报告日期：年月日姓名：学号:报告退发: （订正、重做）同组者: 教师审批签字: 实验名称：某厂排气管车间生产计划的优化分析实验内容: 1、制定利润最大的生产计划（1）建立线性规划数学模型；车间基本生产工人工时计算如下: S＝27×(22×8 – 2 – 4 ) =4590(小时/月）按产量1%的工废实际可用时间:4590×（1– 1% )=４544（小时/月) 每种设备月实际可利用工时(减去工废后的可用时间）如下表所示: 设备名称月可利用工时(小时／月）设备名称月可利用工时（小时/月) 1.平面铣床139２×( 1– 1% )＝1378 2.卧铣床1４00×（1 – 1% )=１386 2０82×( 1 – 1％）=２061 4．单面铣床７0０×( １–１% )=693 3．组合钻床 5.攻丝床2088×（1– 1% )＝２067 ６．精铣床1396×( 1 – 1%)=1３8２ 1376×（1– 1%）=１３62 8．摇臂钻床1384×( 1 –１% ）=137０7．扩孔钻床９.去毛刺 700×(１– 1% )=693 10.清洗机700×( 1 – 1%)=693 机每种排气管利润（减去料废后的利润）如下表所示：

品项目 1 2 3 4 5 6 7 8 毛坯价格 ① 98１7１０4１02 料废成本(① ×2％) １．96 2.0８ 1.88 2.２4 2.1２ 1.94２.08 2.０4辅料消耗 2 2 2 2 2 2 2 2 动能消耗１10 10 工具等消耗 1 １管理费用１.４５５１.099 1.21 1.44 1.188 1．2２65 1.3０8 1.56 税收15 16 14.8 １7 16.5 14.5 15.6 15．５售价150 １６ 0.1 1４9 １72 166 145．6 15７.8 1５５.8 利润(元/根）11．5８５11.921 13．１１1 ３.32 1 ３.192 １0．９ 335 13.812 １1.7 根据给定的条件及要求,设x1，x2,…ｘ8分别为第一种、第二种……第八种排气管的月产量，月产量单位为1000根/月。将目标函数取为每月的总产品利润,得出关于月计划生产的线性规划模型为: ｍaxＺ= 11５8５ｘ1＋1192１ｘ2+１3110 x3+1３32０x４＋13192 x5+10９33.5x ６+13812 x7＋１１7００x8 ４1.9ｘ1+４3.3ｘ2+４1.７x3+4８．1ｘ4+46．5ｘ５＋40.9 x6+47．8 x7+45.７x8 ≤ ４544 4.0 x1＋4．5 x2+４．8 x3+５.8 x４+５.2ｘ5+４.０ｘ６+４.６ｘ7+５.６x8 ≤1378 3．９x1＋4.5 x２+4．３x3+5．0ｘ4+4.9ｘ5+4.4 x６+5.1 x７+４.8 x８≤１38６ 5．9x1＋5．8ｘ２＋5.7 x3+６．3 x4+６.5 x５+6.0 x６＋6.6ｘ７+6．４ｘ8 ≤2０61 3.５ｘ１+3.0ｘ２+3.7ｘ3＋ 4.０ｘ４+３.８x5+3.0 x6+4.1 x7＋3.4ｘ8≤693 ５.８x1+6.2 x2+５．７ｘ3+6.４x４+6.3 x5+6．0ｘ６+6.５x7＋６.２x8 ≤20６7 5．5 x1＋５.7x２＋4.7ｘ3+6．0 x4+5.９x5+5．２ｘ6＋6.2 x7＋５.6ｘ８≤138２ 3.9ｘ1+3.8x2+ 4.０x3+４.1 x４+3.7 ｘ5+3.5 x６+4．1 x7+3.６x8 ≤13６2

西安邮电大学认识实习报告

西安郵電大学认知实习报告书院系名称：计算机学院学生姓名：张帆帆专业名称：计算机科学与技术班级：1506班实习时间：2016年6月20-216年6月22

认识实习报告 2016年6月20日我们开始了我们为期三天的认识实习，从学校一些重要的实验室到学校的网络中心室，还有西安软件园和中兴手机。虽然说我们从来到这所学校已经快一年了，但是我们对学校的了解还是微乎其微，比如各种各样的实验室，他们有自己不同的发展方向和侧重点，还有各种各样高级的实验设备和经验丰富的指导老师，这些都是一个学校，一个学院甚至一个专业的精华所在，是无数代前人和学长学姐努力的产物，他们把它赋予价值并且发扬光大。相比于对自己学院的核心一无所知，整天待在宿舍，能见识一下自己所在地方的魅力，显得意义非凡。第一天我们参观的是自动化学院的实验室，参观的同时有老师或者学长进行讲解，我们首先参观的是217-221实验室，他们的侧重点是测控和计控原理，其中自控原理的核心思想是闭环控制，老师讲解的时候举了一个空调的例子。一个闭环控制主要部件有控制器、执行器、被控对象和监测系统。即由监测系统检测并反馈给CPU，CPU启动执行器，执行器执行其所控制的被控对象，比如空调会制冷，形成一个闭合的执行系统。211实验室主要是单片机，也叫单板机，是一个装置的核心控制部分。然后还参观了无线传感实验室，主要承担智能传感网、物联网技术等，面向测控技术与仪器。中途老师还带我们参观了航模实验室，虽然核心地区在一个非常小的屋子里，但是里面墙壁，桌面上，地上，摆放着各种各样的飞行器组件及成品，虽然有些成品看起来非常简单，但它却能完成各种各样的飞行任务，保持双翼平衡及旋转角度控制，真的非常神奇，因为完成它们的不是什么传说中的科学家或教授，而是在校大学生，是一些为了梦想不断学习的普普通通的大学生。最后参观的是非常有意思的一个实验室，叫机器人创新实验室，他们主要是做一些有目的性的“机器人”，可以完成特定任务，当然它们的长相不一定是人的形态，是一个比较形象的称呼，主要是由一些传感器和电机以及单片机焊接电路组成，学长还展示了他们的功能，比如避障，摇尾巴，电子秤，会踢足球的机器人，还有门禁系统，我觉得这些和我们东区的嵌入式其实差不多，可是我却并不了解！在这个过程中，老师还展示了一些学长学姐由一个思想的火花把它变成现实

西安交大数字图像处理第二次实验报告

数字图像处理第二次作业

摘要本次报告主要记录第二次作业中的各项任务完成情况。本次作业以Matlab 2013为平台，结合matlab函数编程实现对lena.bmp,elain1.bmp图像文件的相关处理：1.分别得到了lena.bmp 512*512图像灰度级逐级递减8-1显示，2.计算得到lena.bmp图像的均值和方差，3.通过近邻、双线性和双三次插值法将lena.bmp zoom到2048*2048，4. 把lena和elain 图像分别进行水平shear（参数可设置为1.5，或者自行选择）和旋转30度，并采用用近邻、双线性和双三次插值法zoom到2048*2048。以上任务完成后均得到了预期的结果。 1.把lena 512*512图像灰度级逐级递减8-1显示（1）实验原理：给定的lena.bmp是一幅8位灰阶的图像，即有256个灰度色。则K位灰阶图像中某像素的灰度值k(x,y)(以阶色为基准)与原图同像素的灰度值v(x,y)(以256阶色为基准)的对应关系为：式中floor函数为向下取整操作。取一确定k值，对原图进行上式运算即得降阶后的k位灰阶图像矩阵。（2）实验方法首先通过imread()函数读入lena.bmp得到图像的灰度矩阵I，上式对I矩阵进行灰度降阶运算，最后利用imshow()函数输出显示图像。对应源程序为img1.m。（3）处理结果 8灰度级

7灰度级 6灰度级 5灰度级

4灰度级 3灰度级 2灰度级

1灰度级（4）结果讨论：由上图可以看出，在灰度级下降到5之前，肉眼几乎感觉不出降阶后图像发生的变化。但从灰度级4开始，肉眼明显能感觉到图像有稍许的不连续，在灰度缓变区常会出现一些几乎看不出来的非常细的山脊状结构。随着灰度阶数的继续下降，图像开始出现大片的伪轮廓，灰度级数越低，越不能将图像的细节刻画出来，最终的极端情况是退化为只有黑白两色的二值化图像。由此可以得出，图像采样的灰度阶数越高，灰度围越大，细节越丰富，肉眼看去更接近实际情况。 2.计算lena图像的均值方差（1）实验原理对分辨率为M*N的灰度图像，其均值和方差分别为：（2）实验方法首先通过imread（）函数读入图像文件到灰度矩阵I中，然后利用 mean2函数和std2函数计算灰度矩阵（即图像）的均值和标准差，再由标准差平方得到方差。对应源程序：img1.m （3）处理结果均值me =99.0512，标准差st =52.8776，方差sf =2.7960e+03。（4）结果分析图像的均值可反应图像整体的明暗程度，而方差可以反应图像整体的对比度情况，方差越大，图像的对比度越大，可以显示的细节就越多。 3.把lena图像用近邻、双线性和双三次插值法zoom到2048*2048；（1）实验原理图像插值就是利用已知邻近像素点的灰度值来产生未知像素点的灰度值，以便由原始图

西南交通大学限修课数学实验题目及答案四

实验课题四曲面图与统计图第一大题：编程作下列曲面绘图：用平面曲线r=2+cos(t)+sin(t)，t∈(0,π)绘制旋转曲面 t=0:0.02*pi:pi; r=2+cos(t)+sin(t); cylinder(r,30) title('旋转曲面'); shading interp 用直角坐标绘制双曲抛物面曲面网线图，z2=xy (-3

axis off 用直角坐标绘制修饰过的光滑曲面曲面：z 4=sin(x )－cos(y ) x 与y 的取值在(-π,π) [x,y]=meshgrid(-pi:0.02*pi:pi); z4=sin(x)-cos(y); surf(x,y,z4); title('picture 4'); shading interp axis off 用连续函数绘图方法绘制曲面)2 s in (6522x y x z ++=,x ∈[-2pi,2pi], y ∈[-2pi,2pi],并作图形修饰。 ezsurf(@(x,y)(x^2+y^2+6*sin(2*x)),[-2*pi 2*pi -2*pi 2*pi]) title('picture 5'); shading interp axis off 第二大题：按要求作下列问题的统计图： x21是1—10的10维自然数构成的向量，y21是随机产生的10维整数向量，画出条形图。(提示bar(x,y)) x21=1:10; y21=randn(10,1); bar(x21,y21) 随机生成50维向量y22，画出分5组的数据直方图。（提示hist(y,n)）