《信号与系统B》实验报告内容及格式
matlab软件仿真实验（信号与系统）（1）
matlab软件仿真实验（信号与系统）（1）《信号与系统实验报告》学院：信息科学与⼯程学院专业：物联⽹⼯程姓名：学号：⽬录实验⼀、MATLAB 基本应⽤实验⼆信号的时域表⽰实验三、连续信号卷积实验四、典型周期信号的频谱表⽰实验五、傅⽴叶变换性质研究实验六、抽样定理与信号恢复实验⼀MATLAB 基本应⽤⼀、实验⽬的：学习MATLAB的基本⽤法，了解 MATLAB 的⽬录结构和基本功能以及MATLAB在信号与系统中的应⽤。
⼆、实验内容：例⼀已知x的取值范围，画出y=sin(x)的图型。
x=0:0.05:4*pi;y=sin(x);plot(y)例⼆计算y=sin(π/5)+4cos(π/4)例三已知z 取值范围，x=sin(z);y=cos(z);画三维图形。
z=0:pi/50:10*pi;x=sin(z);y=cos(z);plot3(x,y,z)xlabel('x')ylabel('y')zlabel('z')例四已知x的取值范围，⽤subplot函数绘图。
参考程序：x=0:0.05:7;y1=sin(x);y2=1.5*cos(x);y3=sin(2*x);y4=5*cos(2*x);subplot(2,2,1),plot(x,y1),title('sin(x)')subplot(2,2,2),plot(x,y2),title('1.5*cos(x)')subplot(2,2,3),plot(x,y3),title('sin(2*x)')subplot(2,2,4),plot(x,y4),title('5*cos(2*x)')连续信号的MATLAB表⽰1、指数信号：指数信号Ae at在MATLAB中可⽤exp函数表⽰，其调⽤形式为：y=A*exp(a*t) (例取 A=1,a=-0.4)参考程序：A=1;a=-0.4;t=0:0.01:10;ft=A*exp(a*t);plot(t,ft);grid on;2、正弦信号：正弦信号Acos(w0t+?)和Asin(w0t+?)分别由函数cos和sin表⽰，其调⽤形式为：A*cos(w0t+phi) ;A*sin(w0t+phi) (例取A=1，w0=2π，?=π/6) 参考程序：A=1;w0=2*pi; phi=pi/6; t=0:0.001:8;ft=A*sin(w0*t+phi);plot(t,ft);grid on ;3、抽样函数：抽样函数Sa(t)在MATLAB中⽤sinc函数表⽰，其定义为：sinc(t)=sin(πt)/( πt)其调⽤形式为：y=sinc(t)参考程序：t=-3*pi:pi/100:3*pi;ft=sinc(t/pi);plot(t,ft);grid on;4、矩形脉冲信号：在MATLAB中⽤rectpuls函数来表⽰，其调⽤形式为：y=rectpuls(t,width),⽤以产⽣⼀个幅值为1，宽度为width,相对于t=0点左右对称的矩形波信号，该函数的横坐标范围由向量t决定，是以t=0为中⼼向左右各展开width/2的范围，width的默认值为1。
信号与线性系统实验指导书v0.2
《信号与线性系统》实验指导书东华大学信息学院通信与电子信息工程系实验要求及说明一、实验报告内容实验报告包括原理分析、源程序、执行结果分析及实验总结，其中原理分析和实验总结需要手写,其他可打印。
二、实验成绩实验成绩包括出勤（10％）、实验表现（10％）、编程（30％）和实验报告（50％）几部分。
三、其他说明缺席3次及以上取消考试成绩。
目录实验一连续信号的时域分析 1 实验二连续时间系统的时域分析 3 实验三连续信号的频域分析 9 实验四连续系统的频域分析 12 实验五信号采样与重建 15 实验六离散时间信号和系统分析 17 附录 MATLAB主要命令函数表 20实验一连续信号的时域分析一、实验目的1、熟悉MATLAB软件。
2、掌握常用连续信号与离散信号的MA TLAB表示方法。
二、实验设备安装有matlab6.5以上版本的PC机一台。
三、MATLAB使用说明1、在MATLAB可视化绘图中，对于以t为自变量的连续信号，在绘图时统一用plot函数。
例题：绘出t从-10到10的sin(t)的波形。
t=-10:0.05:10;f=sin(t);plot(t,f);title('f(t)=sin(t)');xlabel('t') ;axis([-10,10,-1,1])grid on可得图1所示图形。
图12、此外也可以利用MATLAB的ezplot函数对连续信号画图。
例题：绘出t从-10到10的sin(t)的波形。
clcclear allclose allsyms tf=sin(t)ezplot(f, [-10 10]);xlabel('t');title ('f(t)=sin(t)') ;grid on图2四、实验内容1、用MATLAB表示连续信号：tAeα，cos()A tωϕ+，0sin()A tωϕ+。
2、用MATLAB表示抽样信号(sinc(t))、矩形脉冲信号(rectpuls(t, width))及三角脉冲信号(tripuls(t, width, skew))。
信号与系统实验指导书
信号与系统实验指导手册沈阳工业大学信息科学与工程学院2005年10月前言“信号与系统”是电子工程、通信工程、信息工程、微电子技术、自动化、计算机等电类相关专业的一门重要的专业基础课，为国内、外各高等院校相关专业的主要课程。
由于本课程的理论性、系统性较强，为使抽象的概念和理论形象化、具体化，使学生能够比较深入的理解《信号与系统》课程的基本理论和分析方法，并提高学生分析问题和解决问题的能力。
为此，开设了基于本课程的实验。
《信号与系统》实验指导手册，将《信号与系统》课程的理论知识与“信号与系统”的实验系统设备结合，从内容上对教材起到了一定的补充作用，为学生具体实验进行了指导。
鉴于时间仓促，可能会存在一些不足与错误之处，欢迎大家批评指正，使之完善。
编者2005年10月目录实验一系统的特性测试 (1)实验二信号的采样与恢复 (8)实验三模拟滤波器分析 (14)实验四模拟滤波器的设计 (26)实验一系统的特性测试一、实验目的1、学会利用运算单元，搭建一些简单的实验系统。
2、学会测试系统的频率响应的方法。
3、了解二阶系统的阶跃响应特性。
4、学会对其零状态响应和零输入响应进行分析。
二、实验内容1、根据要求搭建一阶、二阶实验系统。
2、测试一阶、二阶系统的频响特性和阶跃响应。
三、预备知识学习使用波特图测试系统频响的方法。
四、实验仪器1、信号与系统实验箱一台（主板）。
2、线性系统综合设计性模块一块。
3、20MHz双踪示波器一台。
五、实验原理1、基本运算单元（1）比例放大1）反相数乘器由：2211R U R U -= 则有：1122R U R U = 2）同相数乘器 由：54443R R U R U +=则有：()45434R R R U U += （2） 积分微分器1）积分器：由：21211//1R SC U R U -= 则有：()1212121C SR R R U U +-= 2）微分器 由：14131R USC U -= 则有：S C R U U 1134-= （3） 加法器1）反相加法器有：⎪⎪⎭⎫⎝⎛+-=211032R U R U R U2）同相加法器由：()⎪⎪⎩⎪⎪⎨⎧+=⎪⎪⎭⎫ ⎝⎛++=+-+87576434433111R R U R U U R R R R U R U 令643*////R R R R = 则有：()⎪⎪⎭⎫⎝⎛++=4433787*5R U R U R R R R U 2、N 阶系统()()()()()()()()t e E t e dt dE t e dtd E te dt d E t r C t r dt dC t r dtd C t r dt d C m m m m m m n n n n n n ++∧++=++∧++------1111011110根据零状态响应（起始状态为零） ，则对其进行拉氏变换有：()()()()()()()()S e E S Se E S e S E S e S E S r C S Sr C S r S C S r S C m m m m n n n n ++∧++=++∧++----11101110则其传递函数表达式为：()()()nn n n mm m m C S C S C S C E S E S E S E S e S r S H ++∧++++∧++==----111011103、作为一阶系统，一般表达式为：()1010C S C E S E S H ++=一阶系统是构成复杂系统的基本单元，学习一阶的特点有助于对一般系统特性的了解。
卷积信号实验报告
信号与系统上机实验报告一连续时间系统卷积的数值计算140224 班张鑫学号 14071002 一、实验原理计算两个函数的卷积卷积积分的数值运算实际上可以用信号的分段求和来实现，即：如果我们只求当 t = n∆ t1 是r ( t )的值，则由上式可以得到：∆t足够小时，r(t2)就是e(t)和f(t)卷积积分的数值近似值由上面的公式可当1以得到卷积数值计算的方法如下：(1)将信号取值离散化，即以为周期，对信号取值，得到一系列宽度间隔为的矩形脉冲原信号的离散取值点，用所得离散取值点矩形脉冲来表示原来的连续时间信号；(2)将进行卷积的两个信号序列之一反转，与另一信号相乘，并求积分，所得为t=0时的卷积积分的值。
以为单位左右移动反转的信号，与另一信号相乘求积分，求的t<0和t>0时卷积积分的值；(3)将所得卷积积分值与对应的t标在图上，连成一条光滑的曲线，即为所求卷积积分的曲线。
1信号与系统上机实验报告一二、处理流程图三、C程序代码#include"stdafx.h"#include"stdio.h"//#include "stdilb.h"float u(float t){while (t>= 0) return(1);while (t<0) return(0);}float f1(float t){return(u(t+2)-u(t-2));}float f2(float t){return(t*(u(t)-u(t-2))+(4-t)*(u(t-2)-u(t-4)));}int_tmain(int argc, _TCHAR* argv[]){FILE *fp;fp=fopen("juanji.xls","w+");float t,i,j,result=0;for(i=-2;i<=6;i=i+0.1){result=0;for(j=0;j<=4;j=j+0.1)result+=f2(j)*f1(i-j)*0.1;printf("%.1f\t%.2f\t",i,result);fprintf(fp,"%.1f\t%.2f\n",i,result);}printf ("\n");return 0;}四、运行结果五、卷积曲线六、感想与总结卷积是信号与系统时域分析的基本手段，主要用于求解系统的零状态响应。
MATLAB信号与系统实验报告19472[五篇范文]
![MATLAB信号与系统实验报告19472[五篇范文]](https://img.taocdn.com/s3/m/a72999dee109581b6bd97f19227916888486b9f2.png)
MATLAB信号与系统实验报告19472[五篇范文]第一篇：MATLAB信号与系统实验报告19472信号与系统实验陈诉（5）MATLAB 综合实验项目二连续系统的频域阐发目的：周期信号输入连续系统的响应可用傅里叶级数阐发。
由于盘算历程啰嗦，最适适用MATLAB 盘算。
通过编程实现对输入信号、输出信号的频谱和时域响应的盘算，认识盘算机在系统阐发中的作用。
任务：线性连续系统的系统函数为11)(+=ωωjj H，输入信号为周期矩形波如图 1 所示，用MATLAB 阐发系统的输入频谱、输出频谱以及系统的时域响应。
-3-2-1 0 1 2 300.511.52Time(sec)图 1要领：1、确定周期信号 f(t)的频谱nF&。
基波频率Ω。
2、确定系统函数 )(Ω jn H。
3、盘算输出信号的频谱n nF jn H Y&&)(Ω=4、系统的时域响应∑∞-∞=Ω=nt jnn eY t y&)（MATLAB 盘算为y=Y_n*exp(j*w0*n“*t);要求（画出 3 幅图）：1、在一幅图中画输入信号f(t)和输入信号幅度频谱|F(jω)|。
用两个子图画出。
2、画出系统函数的幅度频谱|H(jω)|。
3、在一幅图中画输出信号y(t)和输出信号幅度频谱|Y(jω)|。
用两个子图画出。
解:(1)阐发盘算：输入信号的频谱为(n)输入信号最小周期为=2，脉冲宽度，基波频率Ω=2π/ =π，所以(n)系统函数为因此输出信号的频谱为系统响应为(2)步伐：t=linspace(-3,3,300);tau_T=1/4;%n0=-20;n1=20;n=n0:n1;%盘算谐波次数20F_n=tau_T*Sa(tau_T*pi*n);f=2*(rectpuls(t+1.75,0.5)+rectpuls(t-0.25,0.5)+rectpuls(t-2.25,0.5));figure(1),subplot(2,1,1),line(t,f,”linewidth“,2);%输入信号的波形 axis([-3,3,-0.1,2.1]);grid onxlabel(”Time(sec)“,”fontsize“,8),title(”输入信号“,”fontweight“,”bold“)%设定字体巨细，文本字符的粗细text(-0.4,0.8,”f(t)“)subplot(2,1,2),stem(n,abs(F_n),”.“);%输入信号的幅度频谱xlabel(”n“,”fontsize“,8),title(”输入信号的幅度频谱“,”fontweight“,”bold“)text(-4.0,0.2,”|Fn|“)H_n=1./(i*n*pi+1);figure(2),stem(n,abs(H_n),”.“);%系统函数的幅度频谱xlabel(”n“,”fontsize“,8),title(”系统函数的幅度频谱“,”fontweight“,”bold“)text(-2.5,0.5,”|Hn|“)Y_n=H_n.*F_n;y=Y_n*exp(i*pi*n”*t);figure(3),subplot(2,1,1),line(t,y,“linewidth”,2);%输出信号的波形 axis([-3,3,0,0.5]);grid onxlabel(“Time(sec)”,“fontsize”,8),title(“输出信号”,“fontweight”,“bold”)text(-0.4,0.3,“y(t)”)subplot(2,1,2),stem(n,abs(Y_n),“.”);%输出信号的幅度频谱xlabel(“n”,“fontsize”,8),title(“输出信号的幅度频谱”,“fontweight”,“bold”)text(-4.0,0.2,“|Yn|”)(3)波形：-3-2-1 0 1 2 300.511.52Time(sec)输入信号f(t)-20-15-10-5 0 5 10 15 2000.10.20.30.4n输入信号的幅度频谱|Fn|-20-15-10-5 0 5 10 15 2000.10.20.30.40.50.60.70.80.91n系统函数的幅度频谱|Hn|-3-2-1 0 1 2 300.10.20.30.4Time(sec)输出信号y(t)-20-15-10-5 0 5 10 15 2000.10.20.30.4n输出信号的幅度频谱|Yn| 项目三连续系统的复频域阐发目的：周期信号输入连续系统的响应也可用拉氏变更阐发。
信号与系统实验报告1抽样定理
本科实验报告课程名称：信号与系统实验项目：抽样定理实验地点：北区博学楼机房专业班级：电信1201 学号: ******** 学生姓名：指导教师：***一、实验目的：1、了解电信号的采样方法与过程以及信号恢复的方法。
2、验证抽样定理，加深对抽样定理的认识和理解。
二、原理说明：离散时间信号可以从离散信号源获得，也可以从连续时间信号经抽样而获得。
抽样信号fs（t）可以看成是连续信号f（t）和一组开关函数s（t）的乘积。
即：fs（t）=f（t）×s（t）对抽样信号进行傅里叶分析可知，抽样信号的频谱包含了原连续信号以及无限个经过平移的原信号频谱。
平移后的频率等于抽样频率fs及其各次谐波频率2fs、3fs、4fs、5fs......。
正如测得了足够的实验数据以后，我们可以在坐标纸上把一系列数据点连接起来，得到一条光滑的曲线一样，抽样信号在一定条件下也可以恢复为原信号。
只要用一个截止频率等于原信号频谱中最高频率fmax的低通滤波器，滤除高频分量，经滤波后得到的信号包含了原信号频谱的全部内容，故在低通滤波器的输出可以得到恢复后的原信号。
但原信号得以恢复的条件是fs>2B，其中fs为抽样频率，B为原信号占有的频带宽度。
而fmin=2B为最低的抽样频率，又称为“奈奎斯特抽样率”。
当fs<2B 时，抽样信号的频谱会发生混叠，从发生混叠后的频谱中，我们无法用低通滤波器获得原信号频谱的全部内容。
在实际使用中，仅包含有限频谱的信号是极少的，因此即使fs=2B，恢复后的信号失真还是难免的。
为了实现对连续信号的抽样和抽样信号的复原，可用以下实验原理方案：图1-3 抽样定理实验方框图三、实验内容及步骤：1、方波信号的抽样与恢复。
1）观察方波信号的抽样。
调节函数信号发生器，使其输出频率分别为1KHZ、3KHZ，s（t）的频率分别置3.9KHz、15.6KHz、62.5KHz，观察抽样后的波形，并记录之。
方波原始图62.5KHz的抽样图2）观察恢复后的波形。
信号与系统MATL实验及代码
实验一、MATLAB编程基础及典型实例一、实验目的（1）熟悉MATLAB软件平台的使用；（2）熟悉MATLAB编程方法及常用语句；（3）掌握MATLAB的可视化绘图技术；（4）结合《信号与系统》的特点，编程实现常用信号及其运算。
示例一：在两个信号进行加、减、相乘运算时，参于运算的两个向量要有相同的维数，并且它们的时间变量范围要相同，即要对齐。
编制一个函数型m文件，实现这个功能。
function [f1_new,f2_new,n]=duiqi(f1,n1,f2,n2)a=min(min(n1),min(n2));b=max(max(n1),max(n2));n=a:b;f1_new=zeros(1,length(n));f2_new=zeros(1,length(n));tem1=find((n>=min(n1))&(n<=max(n1))==1);f1_new(tem1)=f1;tem2=find((n>=min(n2))&(n<=max(n2))==1);f2_new(tem2)=f2;四、实验内容与步骤− 2 t （2）绘制信号x(t)= esin( t=0:0.1:30; 23t ) 的曲线，t的范围在0~30s，取样时间间隔为0.1s。
y=exp(-sqrt(2)*t).*sin(2*t/3); plot(t,y);（3）在n=[-10:10]范围产生离散序列：x( n)=⎧2 n ,− 3⎨≤ n ≤ 3,并绘图。
⎩0, Othern=-10:1:10;z1=((n+3)>=0);z2=((n-3)>=0);x=2*n.*(z1-z2);stem(n,x);（4）编程实现如下图所示的波形。
t=-2:0.001:3;f1=((t>=-1)&(t<=1));f2=((t>=-1)&(t<=2));f=f1+f2;plot(t,f);axis([-2,3,0,3]);（5）设序列f1(k)={ 0 , 1 , 2 , 3 , 4 , 5 },f2(k)={ 6 , 5 , 4 , 3 , 2 , 1 , 0 }。
信号实验报告
第一部分正文实验一常用信号观察一、实验目的：1．了解常用波形的输出和特点；2．了解相应信号的参数；3．了解示波器与函数发生器的使用；4．了解常用信号波形的输出与特点。
二、实验原理：描述信号的方法有很多可以是数学表达式（时间的函数），也可以是函数图形（即为信号的波形）。
信号的产生方式有多种，可以是模拟量输出，也可以是数字量输出。
本实验由数字信号发生器产生，是数字量输出，具体原理为数字芯片将数字量通过A/D 转换输出，可以输出广泛频率范围内的正弦波、方波、三角波、锯齿波等等。
示波器可以暂态显示所观察到的信号波形，并具有信号频率、峰值测量等功能。
三、实验内容：1．由数字信号发生器产生正弦波、三角波、方波以及锯齿波并输入示波器观察其波形。
2．使用示波器读取信号的频率与幅值。
四、实验设备：1．函数信号发生器一台2．数字示波器一台。
五、实验步骤：1．接通函数发生器的电源，连接示波器。
2．利用函数发生器产生各种基本信号波形，并将波形结果导入计算机中，保存图像，写出各种信号的数学表达式。
六、实验结果：根据实验测量的数据，绘制各个信号的波形图，并写出相应的数学函数表达式。
该试验包括交流：① 该正弦信号的数学表达式为：)1001sin(4t y π=图1-1输入正弦波（Hz 504，V ±） ② 该方波的数学表达式为： )]02.001.0()02.0([4∑∞-∞=----=k k t u k t u y图1-2 输入方波（Hz 504，V ±） ③ 该三角波的数学表达式为：∑∞-∞=-------+-----=k k t u k t u k t k t u k t u k t y )]}02.002.0()02.001.0()][02.0(02.0[800)]02.001.0()02.0()[02.0(800{图1-3 输入三角波（Hz 504，V ±）④ 该锯齿波的数学表达式为：∑∞-∞=-----=k k t u k t u k t y )]}02.002.0()02.0()[02.0(400{图1-4 输入锯齿波（Hz 504，V ±） 实验的一些问题：数字信号发生器的示值与示波器测量有一定的误差。
信号与系统分析上机实验指导
信号与系统分析上机实验指导吴冬梅石崟闫红梅编西安科技大学2006年3月信号与系统分析上机指导前言“信号与系统分析”课程是电子与通信类专业的主要专业基础课之一，该课程是将学生从电路分析的知识领域引入信号处理与传输领域的关键性课程，对后续专业课起着承上启下的作用。
该课程的基本方法和理论大量应用于计算机信息处理的各个领域，特别是通信、数字语音处理、数字图像处理、数字信号分析等领域，应用更为广泛。
因此，让学生掌握对信号与系统进行分析的基本方法和理论，无论对今后专业课的教学，还是学生毕业后从事专业工作的能力，都具有重要的意义。
长期以来，“信号系统分析”课程一直采用黑板式的单一教学方式，学生仅依靠做习题来巩固和理解教学内容，对课程中大量的应用性较强的内容不能实际动手设计、调试、分析，严重影响和制约了教学效果。
因此，该课程迫切需要进行教学方法和教学手段的改革，即在改进教学方法和调整教学内容的同时，实现在实验环境中，以计算机辅助教学手段，用信号分析的软件帮助学生完成数值计算、信号与系统分析的可视化建模及仿真调试，培养学生主动获取知识和独立解决问题的能力，为学习后继专业课打下坚实的基础。
MATLAB语言在矩阵运算和数值计算等方面具有强大的功能，非常直观、简洁，且效率高，交互性好，因而它在科学与工程计算方面的优势是显而易见的。
采用MATLAB语言作为编程工具，便于学生快速掌握和进行自主创造性的学习。
本上机指导书主要由九个实验组成，内容包括：MATLAB的基本使用；信号的表示及可视化；信号的时域运算、时域变换及MATLAB实现；离散序列卷积和的MATLAB实现；周期信号傅里叶级数；傅立叶变换及性质；利用MATLAB分析系统的频率特性；拉普拉斯变换及其曲面图；离散系统的零极点分析。
本指导书由石崟、闫红梅共同编写，吴冬梅负责全书的修改审订。
在编写过程中得到通信学院教师的大力支持，在此一并表示感谢！由于作者水平有限，本书难免会有错误与不足之处，欢迎读者对本书提出批评和建议。
DSP实验二 语音信号分析与处理,南京理工大学紫金学院实验报告,信号与系统
实验二语音信号分析与处理学号姓名注：1）此次实验作为《数字信号处理》课程实验成绩的重要依据，请同学们认真、独立完成，不得抄袭。
2）请在授课教师规定的时间内完成；3）完成作业后，请以word格式保存，文件名为：学号+姓名4）请通读全文，依据第2及第3 两部分内容，认真填写第4部分所需的实验数据，并给出程序内容。
1. 实验目的(1) 学会MATLAB的使用，掌握MATLAB的程序设计方法(2) 掌握在windows环境下语音信号采集的方法(3) 掌握MATLAB设计FIR和IIR滤波器的方法及应用(4) 学会用MATLAB对语音信号的分析与处理方法2. 实验内容录制一段自己的语音信号，对录制的语音信号进行采样，画出采样后语音信号的时域波形和频谱图，确定语音信号的频带范围；使用MATLAB产生白噪声信号模拟语音信号在处理过程中的加性噪声并与语音信号进行叠加，画出受污染语音信号的时域波形和频谱图；采用双线性法设计出IIR滤波器和窗函数法设计出FIR滤波器，画出滤波器的频响特性图；用自己设计的这两种滤波器分别对受污染的语音信号进行滤波，画出滤波后语音信号的时域波形和频谱图；对滤波前后的语音信号进行时域波形和频谱图的对比，分析信号的变化；回放语音信号，感觉与原始语音的不同。
3. 实验步骤1）语音信号的采集与回放利用windous下的录音机或其他软件录制一段自己的语音（规定：语音内容为自己的名字，以wav格式保存，如wql.wav），时间控制再2秒之内，利用MATLAB提供的函数wavread 对语音信号进行采样，提供sound函数对语音信号进行回放。
[y,fs,nbits]=wavread(file),采样值放在向量y中，fs表示采样频率nbits表示采样位数。
Wavread的更多用法请使用help命令自行查询。
2）语音信号的频谱分析利用fft函数对信号进行频谱分析3）受白噪声干扰的语音信号的产生与频谱分析4）据语音信号的频带情况，设计FIR和IIR两种滤波器5）用滤波器对受污染语音信号进行滤波FIR滤波器fftfilt函数对信号进行滤波，IIR滤波器用filter函数对信号进行滤波6）比较滤波前后信号的波形与频谱7）回放滤波后的语音信号4. 实验数据及实验程序实验数据1）原始语音信号的时域波形和频谱图00.51 1.52 2.53 3.54 4.5x 104-0.2-0.15-0.1-0.0500.050.10.150.20.25声音波形图00.51 1.52 2.53 3.54 4.5x 104100200300400500600声音频谱图2）带限白噪声信号的时域波形和幅频特性00.51 1.52-1-0.500.511.5窄带噪声波形图00.51 1.52-1-0.50.511.5窄带噪声频谱图3）受污染语音信号的时域波形和频谱图0246x 104-0.2-0.15-0.1-0.0500.050.10.150.20.25混合信号波形图0246x 1040100200300400500600混合信号频谱图4）滤波器的频响特性图 FIR 滤波器的频响特性图00.10.20.30.40.50.60.70.80.91-8000-6000-4000-2000Normalized Frequency (⨯π rad/sample)P h a s e (d e g r e e s )0.10.20.30.40.50.60.70.80.91-300-200-1000100Normalized Frequency (⨯π rad/sample)M a g n i t u d e (d B )246x 104-0.3-0.2-0.100.10.20.30.40.50.60.70246x 1040.20.40.60.811.21.4IIR 滤波器的频响特性图0100200300400500600-300-250-200-150-100-5050低通滤波器幅度谱图5）滤波后语音信号的时域波形和频谱图00.51 1.52 2.53 3.54 4.5x 1041002003004005006000246x 104-0.1-0.050.050.10.15恢复信号波形图0246x 104100200300400500600恢复信号频谱图6）滤波前后的语音信号时域波形对比图和频谱对比图0246x 104-0.2-0.15-0.1-0.0500.050.10.150.20.25信号波形图0246x 104-0.1-0.0500.050.10.15恢复信号波形图0246x 104100200300400500600信号频谱图0246x 104100200300400500600恢复信号频谱图实验程序： 1）实验主程序 clc;clear;close%[x,fs,bits]=wavread('录音的名字'); [x,fs,bits]=wavread('录音的名字.wav'); %x=wavread('录音的名字'); sound(x,fs,bits);plot(x);title('声音波形图') figure(2)%y=fft(x,4096); y=fft(x);plot(abs(y));title('声音频谱图')fh=0.9;fl=0.25;n=1;length(x);y1=fh*sinc(fh*(n-5000))-fl*sinc(fl*(n-5000)); figure(5)subplot(1,2,1);plot(y1);title('窄带噪声波形图') y2=fft(y1);subplot(1,2,2);plot(abs(y2));title('窄带噪声频谱图') y3=y1+x; figure(6)subplot(1,2,1);plot(y1+x);title('混合信号波形图') y4=fft(y3);subplot(1,2,2);plot(abs(y4));title('混合信号频谱图')wp=0.5;ws=0.55; detaw=ws-wp; n=ceil(8*pi/detaw); wc=(wp+ws)/2;b1=fir1(n-1,wc/pi,hanning(n));freqz(b1,1,41856)f1=fftfilt(b1,y3);plot(f1)f2=fft(f1);plot(abs(f2))f11=filter (bz,az, y3);figure(8)subplot(1,2,1);plot(f11);title('恢复信号波形图')f22=fft(f11);subplot(1,2,2);plot(abs(f22));title('恢复信号频谱图') sound(f11,fs,bits);figure(9)subplot(1,2,1);plot(x);title('信号波形图')subplot(1,2,2);plot(f11);title('恢复信号波形图') figure(10)subplot(1,2,1);plot(abs(y));title('信号频谱图') subplot(1,2,2);plot(abs(f22));title('恢复信号频谱图') 2）FIR滤波器子程序fh=0.9;f1=0.25;n=1:length(x);h=fh*sinc(fh*(n-5000))-f1*sinc(f1*(n-5000));figure(4)subplot(1,2,1);plot(h);p=fft(h);subplot(1,2,2);plot(abs(p));3）IIR滤波器子程序fs=44100;rp=3;rs=20;wp1=0.5;wss1=0.55;op1=2*fs*tan(wp1/2);os1=2*fs*tan(wss1/2);[N,wc]=buttord(op1,os1,rp,rs,'s')[z,p,k]=buttap(N);[ba,aa]=zp2tf(z,p,k);[b,a]=lp2lp(ba,aa,wc);[bz,az]=bilinear(b,a,fs);H=freqz(bz,az);ma=20*log10(abs(H));figure(7)plot(ma);title('低通滤波器幅度谱图')。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
重 庆 交 通 大 学 实 验 报 告
专 业： 通信工程专业12级 学 号： 631206040 姓 名： 实验所属课程： 信号与系统B 实验室(中心)：通信与交通物联网实验室 指 导 教 师 ： 李益才
实验平均成绩： 实验项目一名称 连续信号的运算实验 实验成绩： 一、实验目的 连续信号的运算包括两个信号的加、乘；还包括一个信号的时间平移。通过本实验，直观上了解连续信号的运算。
二、实验主要内容 1、信号ttf2sin)(1与信号ttf10sin)(1，求)()()(21tftftf，分别画出三个函数的图像。
2、信号ttf2sin)(1与信号ttf10sin)(1，求)()()(21tftftf，分别画出三个函数的图像。
3、信号ttf10sin)(1，请分别画出)2(1tf、)5.0(1tf的图像。 4、信号ttf10sin)(1，请分别画出)2(1tf、)2(1tf的图像。 5、信号ttf10sin)(1，请画出)(1tf的图像。 三、实验步骤及结果 （1）用数值法生成信号f1(t)=sin2πt,f2(t)=sin10πt,f1(t)+f2(t)= sin2πt+ sin10πt,并画出它们的图像。 MATLAB程序： t=linspace(0,10,1000); s1=sin(2*pi*t); s2=sin(10*pi*t); s3=s1+s2; subplot(3,1,1); plot(s1); subplot(3,1,2); plot(s2); subplot(3,1,3); plot(s3); 用数值法生成的信号f1(t),f2(t),f1(t)+f2(t)结果如下：
(2) 用数值法生成信号f1(t)=sin2πt,f2(t)=sin10πt,f1(t)*f2(t)= sin2πt*sin10πt,并画出它们的图像。
MATLAB程序： t=linspace(0,10,1000); s1=sin(2*pi*t); s2=sin(10*pi*t); s3=s1.*s2; subplot(3,1,1); plot(s1); subplot(3,1,2); plot(s2); subplot(3,1,3); plot(s3); 用数值法生成的信号f1(t),f2(t),f1(t)*f2(t)结果如下：
(3)用数值法生成信号f1(2t)=sin20πt,f2(0.5t) = sin5πt,并画出它们的图像。
MATLAB程序： t=linspace(0,10,1000); s1=sin(20*pi*t); s2=sin(5*pi*t); subplot(2,1,1); plot(s1); subplot(2,1,2); plot(s2) 用数值法生成的信号f1(t),f2(t)结果如下：
; (4)用数值法生成信号f1(t)=sin10πt,f2(t+2)=sin10π(t+2),f3(t-2)= sin10π(t-2),并画出它们的图像。
MATLAB程序： t=linspace(0,10,1000); s1=sin(10*pi*t); s2=sin(10*pi*(t+2)); s3=sin(10*pi*(t-2)); subplot(3,1,1); plot(s1); subplot(3,1,2); plot(s2); subplot(3,1,3); plot(s3); 用数值法生成的信号f1(t),f2(t+2), f3(t-2)结果如下： (5)用数值法生成信号f1(t)=sin10πt,f2(-t)=sin10π(-t)并画出它们的图像。
MATLAB程序： t=linspace(0,10,1000); s1=sin(10*pi*t); s2=sin(10*pi*(-t)); subplot(2,1,1); plot(s1); subplot(2,1,2); plot(s2); 用数值法生成的信号f1(t),f2(-t)结果如下： 四、实验结果及分析(包括心得体会，本部分为重点) (1)用数值法生成的信号f1(t),f2(t),f1(t)+f2(t)结果如下：
实验分析：f1(t)=sin2πt,则其输出结果为周期T=1S,幅度为1的正弦函数曲线；f2(t)=sin10πt,则其输出结果为周期
T=0.2S，幅度为1的正弦函数曲线，其图像相当于f1(t)沿横轴压缩为原来的1/5；f1(t)+f2(t)= sin2πt+ sin10πt,则其输出结果为周期T=1S，幅度为1的波动曲线，其图像相当于f1(t)和f2(t)对应各点相加。 (2) 用数值法生成的信号f
1(t),f2(t),f1(t)*f2(t)结果如下：
实验分析：f1(t)=sin2πt,则其输出结果为周期T=1S,幅度为1的正弦函数曲线；f2(t)=sin10πt,则其输出结果为周期
T=0.2S，幅度为1的正弦函数曲线，其图像相当于f1(t)沿横轴压缩为原来的1/5；f1(t)*f2(t)= sin2πt*sin10πt,则其输出结果为周期T=0.5S，幅度为2的波动曲线，其图像相当于f1(t)和f2(t)对应各点相乘。 (3)用数值法生成的信号f1(t),f2 (t)结果如下： 实验分析：f1(t)=sin20πt,则其输出结果为周期T=0.1S,幅度为1的正弦函数曲线；f2(t)=sin0.5πt,则其输出结果为周期
T=4S，幅度为1的正弦函数曲线，其图像相当于f1(t)沿横轴展宽40倍。 (4)用数值法生成的信号f1(t),f2(t+2), f3(t-2)
结果如下：
实验分析：f1(t)=sin10πt,则其输出结果为周期T=0.2S,幅度为1的正弦函数曲线；f2(t)=sin10π(t+2),则其输出结果为周
期T=0.2S，幅度为1的正弦函数曲线，其图像相当于f1(t)沿横轴向左平移两个单位；f3(t)= sin10π(t-2),则其输出结果为周期T=0.2S，幅度为1的正弦函数曲线，其图像相当于f1(t)沿横轴向右平移两个单位。 (5用数值法生成的信号f1(t),f2(-t)结果如下：
实验分析：f1(t)=sin10πt,则其输出结果为周期T=0.2S,幅度为1的正弦函数曲线；f2(t)=sin10π(-t),则其输出结果为周期
T=0.2S，幅度为1的正弦函数曲线，其图像相当于f1(t)沿横轴反折。
心得与体会： 1. 该实验采用的是数值法表示信号，即在区间（0，10）上取1000个点，然后输入要表示的函数， 其次调用subplot产生多个绘图区间，最后调用plot函数产生波形图； 2. 通过该实验，掌握了连续信号的运算包括两个信号的加、乘；还包括一个信号的时间平移； 3. 通过该实验，直观上了解连续信号的运算和特点； 4. 通过该实验，对matlab软件和plot函数有一定的熟悉。 实验项目二名称 观察信号的频谱实验 实验成绩： 一、实验目的 观察周期信号、非周期信号的频谱，了解和掌握周期信号和非周期信号频谱的特点。观察一些常用函数的频谱，掌握它们的频谱特点。
二、实验主要内容 1、画出ttf10sin)(1的频谱图。 2、画出nntgtf)3()(22的频谱图。 3、画出)()()(21tftftf的频谱图。 4、画出常用信号)(t、)(2tg、)(tsign、)(4tet，te4的频谱图。 5、画出)(10sin2ttet的频谱图。 三、实验步骤及结果 （1）用数值法生成信号f1(t)= sin10πt,并画出它们的频谱图像。 MATLAB程序： t=linspace(0,30,1000); s1=sin(10*pi*t); sf1=fftshift(fft(s1)); figure() subplot(1,1,1) plot(abs(sf1)) 用数值法生成的信号f1(t)的频谱图，结果如下：
(2) 用数值法生成信号nntgtf)3()(22，并画出它们的频谱图像。
MATLAB程序： t=linspace(0,30,2250); s1=[ones(1,150),zeros(1,75)]; s2=[ones(1,150),zeros(1,75)]; for i=1:9; s2=[s2,s1]; end sf1=fftshift(fft(s2)); figure() subplot(1,1,1) plot(abs(sf1)) 用数值法生成的信号f2(t)的频谱图，结果如下： （3）用数值法生成信号)()()(21tftftf,并画出它们的频谱图像。 MATLAB程序： t=linspace(0,30,2250); s1=sin(10*pi*t); s21=[ones(1,150),zeros(1,75)]; s2=[ones(1,150),zeros(1,75)]; for i=1:9 s2=[s2,s21]; end s3=s1.*s2; sf3=fftshift(fft(s3)); figure() subplot(1,1,1) plot(abs(sf3)) 4. t= -2: 0.02: 6; x=(t>=0); sf1=fftshift(fft(t)); figure() subplot(1,1,1) plot(abs(sf1)); 用数值法生成的信号f(t)的频谱图，结果如下：