数字信号处理MATLAB实验-课堂作业PPT文档格式.pptx
《数字信号处理MATLAB实验-课堂作业PPT文档格式.pptx》由会员分享，可在线阅读，更多相关《数字信号处理MATLAB实验-课堂作业PPT文档格式.pptx（36页珍藏版）》请在冰豆网上搜索。
邮件主题邮件主题DSP2012_3班班_XXX_1附件名称附件名称DSP2012_3班班_XXX_1.doc字符、数字均用半角输入第1次实验报告实验报告7%ProgramP2_7clf;
h=321-210-403;
%impulseresponsex=1-23-4321;
%inputsequencey=conv（h,x）;
n=0:
14;
subplot（2,1,1）;
stem（n,y）;
xlabel（Timeindexn）;
ylabel（Amplitude）;
title（OutputObtainedbyConvolution）;
grid;
x1=xzeros（1,8）;
%filters要求输入向量与输出向量的长度相同y1=filter（h,1,x1）;
subplot（2,1,2）;
stem（n,y1）;
title（OutputGeneratedbyFiltering）;
程序案例1：
卷积ok=8;
n1=0:
1:
19;
xa1=sin（2*pi*n1/k）;
subplot（2,2,1）plot（n1,xa1）xlabel（t/T）;
ylabel（x（n）;
xk1=fft（xa1）;
xk1=abs（xk1）;
subplot（2,2,2）stem（n1,xk1）xlabel（k）;
ylabel（X（k）;
8程序案例2：
DFTpn2=0:
15;
xa2=sin（2*pi*n2/k）;
subplot（2,2,3）plot（n2,xa2）xlabel（t/T）;
xk2=fft（xa2）;
xk2=abs（xk2）;
subplot（2,2,4）stem（n2,xk2）xlabel（k）;
第一次课作业要求o实验教程2.1（四）实验内容编制程序产生2种信号，长度可输入确定，函数需要的参数可输入确定，并绘出其图形o实验教程2.2（四）实验内容2编制程序求解下列两个系统的单位冲激响应和阶跃响应，并绘出其图形。
要求分别用filter、conv、impz三种函数完成。
给出理论计算结果和程序计算结果并讨论。
Matlab及其应用9第一次课作业要求N=64，用一个64点的复数FFT程序，一次算出,并绘图。
参考课本5.9.2节10o实验教程2.3（四）实验内容1第一次课作业要求11已知某序列x（n）在单位圆上的N=64等分样点的Z变换为:
用N点IFFT程序计算o实验教程2.3（四）实验内容212p实验教程2.4（四）实验内容求系统的零、极点和幅度频率响应和相位响应。
可能用到的函数:
zplane,tf2zp,freqz,plot,roots,residuez第一次课作业要求13第一次课作业要求14编写一个函数编写一个函数:
functionLconv,Pconv=lconv_pconv（x,y,N）函数用于比较线性卷积和圆周卷积的区别，其中，x,y是两个用于卷积的有限长序列，N为圆周卷积的长度。
Lconv为线性卷积的结果，Pconv为圆周卷积的结果。
要求要求：
用矩阵相乘的形式计算两种卷积，不能使用Matlab中自带的卷积函数。
第一次课作业要求151考虑右图信号，用两种不同方法画出信号在-50Hz50Hz的幅度频谱，a,第一种方法是先写出信号表达式，从傅立叶变换理论推导出幅度频谱的数学表达式，最后用Matlab画出图形b,第二种方法是对信号进行采样，得到离散信号后利用DFT（用Matlab中的函数FFT实现），最后将DFT的结果进行处理得到幅度频谱。
c,用semilogy函数在同一个图形中同时画出两种方法，比较两种方法的区别，解释出现区别的原因。
可能用到的函数可能用到的函数:
abs,fft,fftshift,plot,semilogy,holdon第一次课作业要求第二二次课作业要求o给出一个图像，编程实现对该图像的任意比例的放大及缩小；
编程实现对该图像的任意角度旋转；
解决缩放及旋转时产生的锯齿等不图像不平滑问题。
o编程求解出各自声音信号的基音周期。
至少用两种方法实现语音（元音字母）的基音频率估计；
估计出一个人发不同汉语语音字母发音的基音频率；
估计出一个人对同一个汉语语音字母的不同声调发音的基音频率；
估计出不同人对同一个汉语语音字母发音的基音频率总结出人的发音的一些基本的特点及规律；
16o图像缩小下采样或降采样使得图像符合显示区域的大小；
生成对应图像的缩略图p图像放大上采样或图像插值放大原图像，显示在更高分辨率的显示设备上17图像处理18o缩放算法邻域插值双线性插值双立方插值自适应样条插值及其增强技术小波放大、邻域交换内插和分形放大19图像处理oImread（filename,FTM）20图像处理格式名称格式名称TIFFJPEGGIFBMPPNGXWD描述加标识的图像文件格式联合图像专家图形交换格式Windows位图可移植网络图形XWindow转储可识别扩展名.tif.tiff.jpg.jpeg.gif.bmp.png.xwd图像处理o图像类型亮度图像（Intensityimages）二值图像（Binaryimages）索引图像（Indexedimages）RGB图像（RGBimages）p数据类：
double,unit8,unit16,uint32,int8,int16,int32,single,char,logical2122图像处理器%图像边缘检测clearf=imread（lena512c.bmp,bmp）;
g1=f（:
:
1）;
h=edge（g1,roberts,0.04）;
imshow（h）subplot（1,2,1）imshow（h）subplot（1,2,2）;
imshow（g1）第二次课作业要求23%图像二维DFTClearX=imread（lena512.ras）;
subplot（1,2,1）imshow（uint8（X）Xfft=fft2（X）;
Xfft1=abs（Xfft）;
Xfft1=log10（Xfft1）;
%此处做了一个对数变换。
yclmax=max（max（Xfft1）;
Xfft1=fftshift（Xfft1）;
Xfft2=Xfft1/yclmax*255;
subplot（1,2,2）imshow（uint8（Xfft2）;
图像处理oSizeoImshowoImwriteoImresizeoImrotateofft2Matlab及其应用24oWhosofs=f（1:
2:
end,1:
2,end）基音周期检测o发音器官发音器官：
肺、气管、喉、咽、鼻、口o声带声带：
喉部左右声带，振动形成基本声源o声门声门：
左右声带之间的部位o音调周期音调周期/基音周期基音周期/基音频率基音频率/基频基频声带开启闭合一次的时间基音范围：
低音男性的50Hz到高音女性或儿童的500Hz基音决定了浊音频率的高低p浊音浊音：
声带振动产生的音，所有元音和一些辅音p清音清音：
不由声带振动产生的音，一部分辅音25基音周期检测o声道声道：
喉以上部分，咽、口腔和鼻腔形状变化产生不同声音对发音起决定性作用的器官分布参数系统，有许多自然谐振频率，即共振峰频率共振峰与声道的形状与大小有关语音的频率特性主要由共振峰决定的。
o语音信号处理语音信号处理时域与频域分析语音通信语音识别、语音合成、说话人识别、语音增强26o语音信号的数字化电话语音的频率范围约为60-3400Hz,取样率一般为8kHz高质量语音，频率高端可扩展至7-9KHz,取样率可至15-20kHz语音波形的动态范围达55dB，量化字长应达10bit27基音周期检测o短时性语音信号是非稳态过程，短时间范围内相对稳定（准稳态过程）信号流处理分段（分帧）处理基音范围70-500Hz,则基音周期范围为14ms-2ms帧长一般取10-30ms28基音周期检测o基音检测难点声门激励波形不是完全的周期脉冲声道共振峰的影响噪声的存在基音频率变化范围大29基音周期检测o基音检测方法自相关函数平均幅度差函数小波变换方法倒谱法Hilbert-Huang变换30基音周期检测o自相关函数定义31基音周期检测u序列是周期的，则自相关函数具有相同周期性u偶函数u当k=0时，自相关函数具有极大值，等于确定性信号序列的能量或随机性序列的平均功率32基音周期检测o短时自相关函数o窗口长度N的选择至少大于基音周期的二倍尽可能小，过大将影响短时性33第n2帧3435第三三次课作业要求利用双线性变换法，设计一个0.4dB截止频率为10kHz且在30kHz处有最小阻带衰减为50dB的数字巴特沃兹低通滤波器，其抽样率为100kHz。
利用教材附录A.2节中的（A.9）式求解原型模拟滤波器的阶数，并利用Matlab的M文件buttap设计模拟原型滤波器。
用M文件bilinear将模拟滤波器传输函数变换成期望的数字传输函数。
请用Matlab画出其增益和相位响应，并给出设计的具体步骤。
36第三次课作业要求利用加窗傅里叶级数法，设计一个具有如下指标的线性相位FIR低通滤波器：
通带截止频率在4rad/s处，阻带截止频率在6rad/s处，最大通带衰减为0.2dB，最小阻带衰减为42dB，抽样率为18rad/s。
利用下面的各个窗函数进行设计：
海明窗、汉宁窗和布莱克曼窗。
对于每种情况，给出冲击响应的系数并画出设计的滤波器的增益响应。
分析你的结果。
不要使用M文件fir1。
升级会员 