《数字信号处理》上机实习报告 (10)

《计算机编程与数字信号处理实习》大作业
6月21日
实习任务：从给定的程序（文件包Friday.rar）中，选择一个源程序做详细标注。

（目的：熟悉Matlab程序）
程序名：diyiti
原理：通过对几个程序的详细分析，熟悉掌握了Matlab的一些常用函数的用法和一些简单操作。

包括数组的定义和赋值，如行向量：a=[1 1 1 1],a=（1：3）a=[1:3],列向量：a=[1 1 1 1]等
text,grid on,hold on,print,legend,subplotk,循环语句的用法
等。

6月22日
实习任务：二、能够利用Matlab熟悉地画图，内容包括：X、Y 坐标轴上的label，每幅图上的title，绘画多条曲线时的legend，对图形进行适当的标注等。

（1）在一副图上画出多幅小图；（2）画出一组二维图形；
（3）画出一组三维图形；（4）画出复数的实部与虚部。

（5）完成对一个源程序进行详细注释。

程序名：erweituxing.m
Sanweituxing.m
Shibuxubu.m
原理：在一副图上画出多幅小图是利用subplot（）函数
三维图形用plot3（x,y,z）
复数的实部用real，虚部用imag
6月23-24日
实习任务：三、计算普通褶积与循环褶积，分别使用时间域与频率域两种方法进行正、反演计算，指出循环褶积计算时所存在的边界效应现象；编写一个做相关分析的源程序。

程序名：bianjiexiaoying.m
xiangguan.m
xunhuanzheji.m
zhengfanyan.m
原理：时间域的普通褶积直接用褶积的公式计算，频率域的褶积先对原函数做fft变换，得出的结果是时间域的循环褶积等于频率域的褶积，N的大小与数组的长度相同，如果N小于数组长度则会出现边界效应。

6月25-26日
实习任务：四、设计一个病态（矩阵）系统，分析其病态程度；找出对应的解决方法（提示：添加白噪因子）。

程序名：bingtai.m
6月27日
实习任务：五、设计一个一维滤波处理程序（1、分别做低通、高通、带通、带阻等理想滤波器进行处理；2、窗函数）。

原理：先构造一个函数，然后分别做频率域的低通，高通，带通，
带阻理想滤波器，对原函数做fft变换，最后与滤波器乘积，得出处理结果。

程序名：yiweilvbo.m
6月28-29日
实习任务：六、设计一个二维滤波处理程序（分别做低通、高通等处理）。

程序名：erweilvbo.m
原理：先对一个数组的列向量做fft变换，再对行向量做fft变换，结果为二维fft变换。

按照一维的滤波处理，同时对行和列做低通，高通，带通，带阻滤波处理。

6月30日
实习任务：七、验证时间域的循环褶积对应的是频率域的乘积；线性褶积则不然。

原理：同第三题。

程序名：diqiti.m
7月1-2日
实习任务：八、请用通俗、易懂的语言说明数字信号处理中的一种性质、一条定理或一个算例（顺便利用Matlab对其进行实现）。

程序名：dftxingzhi.m。