rgb2grayx：矩阵x为图像用imread函数读入后得到的返回值，矩阵y为灰度图像的数据矩阵。
本题的第四个要求为采用几种典型的边缘检测算子对灰度图像进行边缘检测。图像边缘包含了图像最重要的信息，图像边缘检测是计算机视觉和数字图像处理等领域研究的重要内容之一。边缘是指图像周围像素灰度有阶跃变化或屋顶变化的那些像素的集合，它存在于目标与背景、区域与区域、基元与基元之间，是图像最基本的特征，以及图像分析与识别的重要环节。边缘检测的实质是采用某种算法来提取出图像中对象与背景的交界线。我们将图像的边缘定义为图像中灰度发生急剧变化的区域边界，图像灰度的变化情况可以用图像灰度分布的梯度来反映，因此可以用局部图像微分来获得边缘检测算子。MATLAB中提供了函数edge可以对灰度图像进行边缘检测。edge函数的调用格式如下：
y1edgeysua
zi：y为灰度图像的数据矩阵，sua
zi是边缘检测算子。几种典型的边缘检测算子包括：sobel算子、prewitt算子、roberts算子、log算子和ca
y算子，返回值y1为边缘检测后图像的数据矩阵。
225第五题：音频信号处理
本题首先要求读入CEG和弦音频信号，然后绘制音频信号的时域波形和频域波形。读入音频信号可以采用wavread函数。wavread函数的调用格式如下：
yfswavreadfile
ame：y就是音频信号，fs是采样频率。绘制出时域波形和频域波形后，观察音频信号的频谱图，我们从中分离出三个能量较为集中的频带。然后设置相应的数字滤波器参数，设计数字滤波器。为了提高滤波器的选择性，我选择了设计椭圆滤波器。MATLAB信号处理工具箱提供了椭圆滤波器的设计函数。本次课程设计中，我用到了ellipord函数和ellip函数，两函数的调用格式如下：NwpoellipordwpwsRpRs：用于计算满足指标的椭圆数字滤波器的最低阶数N和通带边界频率wpo。调用参数wp和ws分别为数字滤波器的通带边界频率和阻带边界频率的归一化值。Rp和Rs分别为通带最大衰减和阻带最小衰减。BAellipNRpRswpo’ftype’：当wpo表示滤波器通带边界频率的标量，而且缺省参数ftype时，该格式返回N阶低通椭圆数字滤波器系统函数的分子和分母多项式系数向量B和A，滤波器通带波纹为RpdB；当ftypehigh时，返回N阶高通椭圆数字滤波器系统函数系数向量B和A。当wpo是表示带通滤波器通带边界频率的二元向量，而且缺省参数ftype时，该格式返回2N阶带通椭圆数字滤波器系统函数的分子和分母多项式系数向量B和A，滤波器通带波纹为RpdB。当ftypestop时，返回r