《MATLAB课程设计---MATLAB语音信号的采集与处理》由会员分享，可在线阅读，更多相关《MATLAB课程设计---MATLAB语音信号的采集与处理（27页珍藏版）》请在金锄头文库上搜索。
1、课 程 设 计 报 告设计名称：MATLAB课程设计系别：学生姓名：班级：学号：成绩：指导教师：开课时间: 2021 至 2021 学年第 2 学期 MATLAB语音信号的采集与处理一、实践的目的和要求1.MATLAB软件功能简介 MATLAB是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言如C、Fortra
2、n的编辑模式，代表了当今国际科学计算软件的先进水平。 MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创立用户界面、连接其他编程语言的程序等，主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。 MATLAB的根本数据单位是矩阵，它的指令表达式与数学、工程中常用的形式十分相似，故用MATLAB来解算问题要比用C，FORTRAN等语言完成相同的事情简捷得多，使MATLAB成为一个强大的数学软件。MATLAB特点：1) 高效的数值计算及符号计算功能,能使用户从繁杂的数学运算分析中解脱出来; 2) 具有完备的图形处理功能,实现计算结果和编程的可视化;
3、3) 友好的用户界面及接近数学表达式的自然化语言,使学者易于学习和掌握; 4) 功能丰富的应用工具箱,为用户提供了大量方便实用的处理工具。MATLAB与其他计算机语言相比，它有很多的优点： 友好的工作平台和编程环境 简单易用的程序语言，运算符丰富，语言简洁紧凑，使用方便灵活，库函数极其丰富。 强大的科学计算机数据处理能力，MATLAB既具有结构化的控制语句如for循环、while循环、break语句和if语句，又有面向对象编程的特性。 应用广泛的模块集合工具箱 语法限制不严格，程序设计自由度大。 程序的可移植性很好，根本上不做修改就可以在各种型号的计算机和操作系统上运行。 出色的图形处理功能
4、实用的程序接口和发布平台2.此题目的意义1有利于根底知识的理解通过?MATLAB及在电子信息类课程中的应用?的学习，我掌握了MATLAB的一些根本知识和用法，对我以后的学习带来方便。2有利于逻辑思维的锻炼程序设计是公认的、最能直接有效地训练学生的创新思维，培养分析问题、解决问题能力的学科之一。即使一个简单的程序，从任务分析、确定算法、界面布局、编写代码到调试运行，整个过程学生都需要有条理地构思，这中间有猜想设想、判断推理的抽象思维训练，也有分析问题、解决问题、预测目标等能力的培养。3有利于与其他学科的整合在程序设计中，我们可以解决其它学科有关问题，也利用其它课程的有关知识来解决MATLAB中比
5、拟抽象很难理解的知识。在MATLAB课中整合其它学科的知识，发挥MATLAB的优势。4有利于治学态度的培养。课程设计中，语句的语法和常量变量的定义都有严格的要求，有时输了一个中文标点、打错了一个字母，程序无法正常运行。因此，程序设计初学阶段，我们经常会犯这样的错误，可能要通过几次乃至十屡次的反复修改、才能成功，但这种现象会随着学习的深入而慢慢改观。这当中就有一个严谨治学、一丝不苟的科学精神的培养，又有一个不怕失败、百折不挠品格的锻炼。二、实践内容 1.常用信号的实现 通过编写MATLAB软件仿真程序，掌握MATLAB软件中常用信号的表示方法，绘图函数的使用。内容如下：1单位阶跃信号f(t)=
6、和单位阶跃序列的实现2单位冲激信号和单位脉冲序列的实现3单边衰减指数信号和实指数序列的实现2.信号的时域根本运算 通过调用MATLAB中的相关函数，编写信号的根本运算仿真程序，实现信号的时域根本运算，到达掌握MATLAB相关运算函数的调用格式及功能的目的。 1信号的相加，相乘和移位运算 2序列的卷积运算3.系统分析 通过调用MATLAB相关函数，实现对LTI系统的时域和频域进行分析，熟练掌握LTI系统单位冲激响应和零状态响应得MATLAB求解方法。 1连续LTI系统的时域分析 2连续LTI系统的频域分析4.综合设计 在熟练掌握MATLAB根本的信号处理方法的根底上，按照要求设计特定功能的综合系
7、统，掌握MATLAB信号处理系统的仿真程序设计。综合系统题目1MATLAB模拟滤波器的实现2MATLAB语音信号的采集与处理三、实践仿真程序设计及结果分析1.常用信号的实现 1单位阶跃信号f(t)= 和单位阶跃序列的实现单位阶跃信号f(t)= 的实现程序：t0=0;t1=-1;t2=3;dt=0.01;t=t1:dt:-t0;n=length(t);t3=-t0:dt:t2;n3=length(t3);u=zeros(1,n);u3=ones(1,n3);plot(t,u);hold on;plot(t3,u3);plot(-t0,-t0,0,1);hold off;axis(t1,t2,-0
8、.2,1.5); xlabel(时间(t);ylabel(幅值(f);title(单位阶跃信号); 时域波形：单位阶跃序列的实现程序： k0=0;%单位阶跃出现的位置k1=-3;k2=6; k=k1:k0-1; n=length(k); k3=-k0:k2; n3=length(k3); u=zeros(1,n); u3=ones(1,n3); stem(k,u,filled); hold on; stem(k3,u3,filled); hold off; axis(k1,k2,-0.2,1.5); title(单位阶跃序列); xlabel(时间(k);ylabel(幅值f(k);时域波形：
9、2单位冲激信号和单位脉冲序列的实现单位冲激信号的实现程序：clear;t0=0;t1=-1;t2=5;dt=0.1;t=t1:dt:t2;n=length(t);x=zeros(1,n);x(1,(t0-t1)/dt+1)=1/dt;stairs(t,x);axis(t1,t2,0,1/dt);xlabel(时间(t);ylabel(幅值(f);title(单位冲激信号);时域波形：单位脉冲序列的实现程序： k1=-3; k2=6; k=k1:k2; n=3;%单位脉冲出现的位置 f=(k-n)=0; stem(k,f,filled);title(单位脉冲序列) xlabel(时间(k);yl
10、abel(幅值f(k);时域波形：3单边衰减指数信号和实指数序列的实现单边衰减指数信号的实现程序：clear t1=-1;t2=10;dt=0.1;t=t1: dt: t2;A1=1;a1=0.5;n=A1*exp(-a1*t);plot(t,n);axis(t1,t2,0,1);xlabel(时间(t);ylabel(幅值(f);title(单边衰减指数信号);时域波形：实指数序列的实现程序：clf;k1=-1;k2=10;k=k1:k2;a=-0.6;A=1;f=A*a.k;stem(k,f,filled);title();xlabel(k);ylabel(f(k);时域波形：2.信号的时
11、域根本运算 1信号的相加，相乘和移位运算 实现两个连续信号的相加，即。 程序 clear all; t=0:0.0001:3; b=3; t0=1; u=stepfun(t,t0); n=length(t); for i=1:n u(i)=b*u(i)*(t(i)-t0);end %产生一个斜坡信号 y=sin(2*pi*t); %产生一个正弦信号 f=y+u; %信号相加 plot(t,f); xlabel(时间(t);ylabel(幅值f(t);title(连续信号的相加); 时域波形：实现两个连续信号的相乘，即。程序 clear all;t=0:0.0001:5;b=3;t0=1;u=s
12、tepfun(t,t0);n=length(t);for i=1:n u(i)=b*u(i)*(t(i)-t0);endy=sin(2*pi*t); f=y.*u;plot(t,f);xlabel(时间(t);ylabel(幅值f(t);title(连续信号的相乘) 时域波形：实现两个连续信号的移位，即，常数程序： clear all; t=0:0.0001:2;y=sin(2*pi*(t); yl=sin(2*pi*(t-0.2); plot(t,y,-,t,yl,-); ylabel(幅值f(t);xlabel(时间t);title(信号的位移);时域波形：2序列的卷积运算实现两个连续信号
13、的卷积运算， 程序：%连续函数卷积计算a=1000;t1=-5:1/a:5;f1=stepfun(t1,0);f2=stepfun(t1,-1/a)-stepfun(t1,1/a);subplot(231);plot(t1,f1);axis(-5,5,0,1.2); %xlabel(时间(t);ylabel(f1(t);title(单位阶跃函数);subplot(232);plot(t1,f2);ylabel(f2(t);title(单位冲激函数);y=conv(f1,f2);r=2*length(t1)-1;t=-10:1/a:10;subplot(233);plot(t,y);axis(-5,5,0,1.2);title(f1与f2的卷积);ylabel(y(t);f11=conv(f1,f1);f22=conv(f2,f2);subplot(234);plot(t,f11);title(f1与f1的卷积);ylabel(f11(t);axis(-5,5,0,5000);subplot(235);plot(t,f22);title(f2与f2的卷积);ylabel(f22(t);时域波形：3.系统分析 1连续LTI系统的时域分析 求系统y(2)(t)+6y(1)(t)+8y(t)=3x(1)(t)+9x(t)的冲激响应和阶跃响
