网站开发的工资,有没有哪种网站推荐一下,阿里云 wordpress搭建网站,广东东莞新增本土确诊0例基于MATLAB的信号频谱分析及实现数 字 信 号 处 理课程设计题目#xff1a; 基于MATLAB的信号频谱分析及实现学院#xff1a; 信息工程专业#xff1a; 通信工程班级#xff1a; 1001学号#xff1a; 2010013448#xff0c; 2010013466姓名#xff1a; 常珍珍 #xff…基于MATLAB的信号频谱分析及实现数 字 信 号 处 理课程设计题目 基于MATLAB的信号频谱分析及实现学院 信息工程专业 通信工程班级 1001学号 2010013448 2010013466姓名 常珍珍 彭婷指导教师 符茂胜目录1、背景32、设计目的33、设计原理34、设计过程45、实验代码及结果45.1 MATLAB源程序45.2 比较结果分析86、实验体会及总结8参考文献91、背景信号的频谱分析就是利用傅里叶分析的方法求出与时域描述相对应的频域描述从中找出信号频谱的变化规律以达到特征提取的目的。实际的待分析信号一般没有解析式直接利用公式进行傅里叶分析非常困难。DFT是一种时域和频域均离散化的傅里叶变换适合数值计算且有快速算法是分析信号的有力工具。本文以连续时间信号为例介绍利用DFT分析信号频谱的基本流程实例列出MATLAB环境下频谱分析的实现程序。2、设计目的1、进一步掌握数字信号处理课程的基本理论、基本方法和基本技术。增进对MATLAB的认识利用MATLAB加深对理论知识的理解。3、了解和掌握使用MATLAB的应用过程和方法为以后的设计打下良好基础。4、通过课程设计培养我们能根据设计要求进行理论知识分析、设计方法总结、典型实例设计等方面的设计综合能力。3、设计原理时域抽样定理时域抽样定理给出了连续信号抽样过程中信号不失真的约束条件对于基带信号信号抽样频率大于等于2倍的信号最高频率即。时域抽样是把连续信号变成适于数字系统处理的离散信号。对连续信号以间隔T抽样则可得到的离散序列为。若则信号与的频谱之间存在其中的频谱为的频谱为。② 离散傅立叶变换(DFT)有限长序列的离散傅立叶变换(DFT)为逆变换为③ 快速傅立叶变换(FFT)离散傅氏变换的快速算法它是根据离散傅氏变换的奇、偶、虚、实等特性对离散傅立叶变换的算法进行改进获得的。 在MATLAB的信号处理工具箱中函数FFT和IFFT用于快速傅立叶变换和逆变换。函数FFT用于序列快速傅立叶变换其调用格式为yfft(x)其中x是序列y是序列的FFTx可以为一向量或矩阵若x为一向量y是x的FFT且和x相同长度若x为一矩阵则y是对矩阵的每一列向量进行FFT。如果x长度是2的幂次方函数fft执行高速基2FFT算法否则fft执行一种混合基的离散傅立叶变换算法计算速度较慢。函数FFT的另一种调用格式为yfft(x,N)式中xy意义同前N为正整数。函数执行N点的FFT若x为向量且长度小于N则函数将x补零至长度N若向量x的长度大于N则函数截短x使之长度为N若x 为矩阵按相同方法对x进行处理。4、设计过程1、利用MATLAB中的wavread命令来读入(采集)语音信号将它赋值给某一向量。再将该向量看作一个普通的信号对其进行FFT变换实现频谱分析。对于波形图与频谱图都可以用 MATLAB画出。我们还可以通过sound命令来对语音信号进行回放以便在听觉上来感受声音变化。2、DFT已知一个语音信号数据采样频率为100Hz试分别绘制N128点DFT的幅频图和N1024点DFT幅频图。3、通过改变系统函数的分子与分母系数比来观察信号滤波前后的幅值与相位。5、5.1 MATLAB源程序(1)编程如下[x,fs,bits]wavread(1.wav,[1024 5120]);sound(x,fs,bits);Xfft(x,4096);magXabs(X);angXangle(X);subplot(2,2,1);plot(x);title(原始信号波形);subplot(2,2,2);plot(X); title(原始信号频谱);subplot(2,2,3);plot(magX);title(原始信号幅值);subplot(2,2,4);plot(angX);title(原始信号相位);程序运行可以听到声音得到的图形为图(1)原始信号波形和频谱图(2)编程如下xwavread(1.wav);sound(x);fs100;N128;yfft(x,N);magyabs(y);f(0:length(y)-1)*fs/length(y);subplot(2,2,1);plot(f,magy);xlabel(频率(Hz));ylabel(幅值);title(N128(
