南昌大学信号与系统实验报告
南昌大学信号与系统MATLAB软件实验报告
信号与系统软件实验报告班级:姓名:学号:指导教师:实验1一．实验要求：绘出门函数()(2)(2)f t t t εε=+--的波形二．实验原理：在MATLAB 中，有一个专门用于表示单位阶跃信号的函数，即stepfun( )函数，它是用数值计算法表示的单位阶跃函数()t ε。
其调用格式为：stepfun(t,t0),其中，t 是以向量形式表示的变量，t0表示信号发生突变的时刻，在t0以前，函数值小于零，t0以后函数值大于零。
有趣的是它同时还可以表示单位阶跃序列()k ε，这只要将自变量以及取样间隔设定为整数即可。
有关单位阶跃序列()k ε的表示方法，我们后面有专门论述，下面通过一个例子来说明如何调用stepfun( )函数来表示单位阶跃函数。
三．Matlab 程序t=-4:0.01:4; %定义时间样本向量t1=-2; %指定信号发生突变的时刻 u1=stepfun(t,t1); %产生左移位的阶跃信号(t+2) t2=2; %指定信号发生突变的时刻 u2=stepfun(t,t2); %产生右移位的阶跃信号(t-2) g=u1-u2; %表示门函数plot(t,g) %绘制门函数的波形axis([-4,4,-0.5,1.5]) %设定坐标轴范围-4<x<4 ，-0.5<y<1.5四．实验结果实验2一．实验要求：若某连续系统的输入为e (t )，输出为r (t )，系统的微分方程为：''()5'()6()3'()2()y t y t y t f t f t ++=+①求该系统的单位冲激响应h (t )及其单位阶跃响应g (t )。
②若2()()t f t e t ε-= 求出系统的零状态响应y(t )二．实验原理：对于连续的LTI 系统，当系统输入为f (t )，输出为y (t )，则输入与输出之间满足如下的线性常系数微分方程：()()0()()nmi j i j i j a y t b f t ===∑∑，当系统输入为单位冲激信号δ(t )时产生的零状态响应称为系统的单位冲激响应，用h(t)表示。
南昌大学信号与系统实验课程3周期信号的频谱测试
南昌大学实验报告学生姓名：学号：班级：实验类型：□验证□综合■设计□创新实验日期：2011-04-15 实验成绩：硬件实验一：周期信号的频谱测试(一)实验目的1，掌握周期信号频谱的测试方式；2，了解典型信号频谱的特点。
(二)实验内容1，测试正弦波的幅度频谱：将信号源、示波器按图连接好，信号源CH1的输出波形调为正弦波，输出频率自选，输出信号幅度自选，并记录幅度与频率的参数。
测出前五次谐波分量，将其数据填入表中。
2，测试三角波的幅度频谱：在实验步骤一的基础上将信号源CH1的输出波形调为三角波，频率自选，幅度自选，并记录幅度和周期的参数。
测出前五次谐波分量。
将测量数据填入表中。
3，测量周期矩形脉冲的幅度频谱：将信号源的输出线接“脉冲”输出端，信号频率，幅度和脉宽自选，测出信号前五次谐波分量，填入表中。
4，测量人体的幅度频谱：用手连接到示波器信号线正极和地线，测得人体频谱分量，测出信号前五次谐波分量，填入表中。
人体的基准频率为50Hz。
(三)实验过程1，测试正弦波的幅度频谱：将信号源、示波器按图连接好，信号源输出波形调为正弦波，设置正弦波的基准频率为100Hz，任意调节输出信号幅度，记录幅度与频率的参数。
测得前南昌大学实验报告学生姓名：学号：班级：实验类型：□验证□综合■设计□创新实验日期：2011-04-15 实验成绩：五次谐波分量，将其数据填入下表中。
f(Hz) 100 200 300 390 490Cn(dB) -5.55 -44.7 -48.7 -51.1 -51.1 Cn(V) 0.394 4.79e-5 1.908e-6 4.37e-6 4.37e-61 0 0 0 0计算值Cn根据上述测量值和理论值画出相应图形如下：上述为Han窗函数下的情况表格，下再附Flattop窗函数下的表格：f(Hz) 100 210 300 370 480Cn(dB) -11.9 -50.3 -50.3 -50.7 -51.9下述图形进行更形象的比较：南昌大学实验报告学生姓名：学号：班级：实验类型：□验证□综合■设计□创新实验日期：2011-04-15 实验成绩：2，测试三角波的幅度频谱：在实验步骤一的基础上将信号源输出波形调为三角波，设置基准频率数值为10HZ，任意幅度调节，并记录幅度和频率的参数。
《信号与系统》课程实验报告
《信号与系统》课程实验报告《信号与系统》课程实验报告一图1-1 向量表示法仿真图形2．符号运算表示法若一个连续时间信号可用一个符号表达式来表示，则可用ezplot命令来画出该信号的时域波形。
上例可用下面的命令来实现（在命令窗口中输入，每行结束按回车键）。
t=-10:0.5:10;f=sym('sin((pi/4)*t)');ezplot(f,[-16,16]);仿真图形如下：图1-2 符号运算表示法仿真图形三、实验内容利用MATLAB实现信号的时域表示。
三、实验步骤该仿真提供了7种典型连续时间信号。
用鼠标点击图0-3目录界面中的“仿真一”按钮，进入图1-3。
图1-3 “信号的时域表示”仿真界面图1-3所示的是“信号的时域表示”仿真界面。
界面的主体分为两部分：1) 两个轴组成的坐标平面（横轴是时间，纵轴是信号值）；2) 界面右侧的控制框。
控制框里主要有波形选择按钮和“返回目录”按钮，点击各波形选择按钮可选择波形，点击“返回目录”按钮可直接回到目录界面。
图1-4 峰值为8V，频率为0.5Hz，相位为180°的正弦信号图1-4所示的是正弦波的参数设置及显示界面。
在这个界面内提供了三个滑动条，改变滑块的位置，滑块上方实时显示滑块位置代表的数值，对应正弦波的三个参数：幅度、频率、相位；坐标平面内实时地显示随参数变化后的波形。
在七种信号中，除抽样函数信号外，对其它六种波形均提供了参数设置。
矩形波信号、指数函数信号、斜坡信号、阶跃信号、锯齿波信号和抽样函数信号的波形分别如图1-5～图1-10所示。
图1-5 峰值为8V，频率为1Hz，占空比为50%的矩形波信号图1-6 衰减指数为2的指数函数信号图1-7 斜率=1的斜坡信号图1-8 幅度为5V，滞后时间为5秒的阶跃信号图1-9 峰值为8V，频率为0.5Hz的锯齿波信号图1-10 抽样函数信号仿真途中，通过对滑动块的控制修改信号的幅度、频率、相位，观察波形的变化。
南昌大学信号与系统实验课程4傅里叶变换
南昌大学实验报告学生姓名： 学 号： 班级： 实验类型：□ 验证 □ 综合 ■ 设计 □ 创新 实验日期： 2011-04-23 实验成绩：傅里叶变换(一)实验目的1，请绘制矩形脉冲⎪⎩⎪⎨⎧<=其他021||1)(t t f 的波形（t ∈[-1，1]）和频谱F(ω)（ω∈[-8pi ，8pi]），并利用你计算得到的频谱恢复时域信号fs(t)，比较和原函数f(t)的差别。
2，绘制周期T1=1、幅度E=1的对称方波的前10项傅里叶的系数（三角函数形式），并用前5项恢复原信号。
3，根据锯齿波信号，分别取一个周期的抽样数据x1(w)，0<=t<5，计算其傅里叶变换x1(w)和x(w)，比较有何不同并解释原因。
(二)实验内容1，根据函数编写matlab 程序，从而实现其功能，程序函数及说明如下： T=2; %定义时域抽样区间长度 N=200; %定义时域抽样点数t=linspace(-T/2,T/2-T/N,N); %定义时域抽样点f=0*t; %初始化时域信号f(t>-1/2&t<1/2)=1; %为时域信号赋值OMG=16*pi; %定义频域抽样区间长度 K=100; %定义频域抽样点数omg=linspace(-OMG/2,OMG/2-OMG/K,K); %定义频域抽样点F=0*omg; %初始化频谱for k=1:K %循环计算每个频谱抽样点的频谱 for n=1:N %循环实现求和运算功能 F(k)=F(k)+T/N*f(n)*exp(-j*omg(k)*t(n));endendfs=0*t; %初始化合成信号for n=1:N %循环计算每个时域抽样点的合成信号 for k=1:K %循环实现求和运算功能 fs(n)=fs(n)+OMG/2/pi/K*F(k)*exp(j*omg(k)*t(n)); End南昌大学实验报告学生姓名：学号：班级：实验类型：□验证□综合■设计□创新实验日期：2011-04-23 实验成绩：endsubplot(1,2,1) %绘制输出图形plot(t,f,t,fs)subplot(1,2,2)plot(omg,F)根据上述函数得到之后的波形图及对应的频谱图：2，根据函数编写matlab程序，从而实现其功能，程序函数及说明如下：E=1;%定义方波幅度T1=1;%定义方波周期omg1=2*pi/T1;%定义基波N=1000;%定义时域抽样点数t=linspace(-T1/2,T1/2-T1/N,N)';%生成时域抽样点f=0*t;%初始化时域信号f(:)=-E/2;f(t>-T1/4&t<T1/4)=E/2;k1=-10; %确定系数的起始下标k2=10;k=[k1:k2]'; %生成系数下标系列F=1/N*exp(-j*kron(k*omg1,t.'))*f; %求指数形式傅里叶级数的系数a0=F(11); %转换到三角函数形式的系数ak=F(12:21)+F(10:-1:1);f1=cos(kron(t,[0:1]*omg1))*[a0;ak(1:1)];南昌大学实验报告学生姓名：学号：班级：实验类型：□验证□综合■设计□创新实验日期：2011-04-23 实验成绩：f3=cos(kron(t,[0:3]*omg1))*[a0;ak(1:3)];fs=cos(kron(t,[0:5]*omg1))*[a0;ak(1:5)]; %用前五个系数合成原函数figure;subplot(1,2,1);plot(t,f,t,f1,t,f3,t,fs); %绘制图形subplot(1,2,2);stem(ak,'k') %绘制杆图形式曲线根据上述函数得到之后的波形图及对应的频谱图：3，根据函数编写matlab程序，计算和比较不同周期下的傅里叶变换差别。
信号与系统实验实验报告
信号与系统实验实验报告一、实验目的本次信号与系统实验的主要目的是通过实际操作和观察，深入理解信号与系统的基本概念、原理和分析方法。
具体而言，包括以下几个方面：1、掌握常见信号的产生和表示方法，如正弦信号、方波信号、脉冲信号等。
2、熟悉线性时不变系统的特性，如叠加性、时不变性等，并通过实验进行验证。
3、学会使用基本的信号处理工具和仪器，如示波器、信号发生器等，进行信号的观测和分析。
4、理解卷积运算在信号处理中的作用，并通过实验计算和观察卷积结果。
二、实验设备1、信号发生器：用于产生各种类型的信号，如正弦波、方波、脉冲等。
2、示波器：用于观测输入和输出信号的波形、幅度、频率等参数。
3、计算机及相关软件：用于进行数据处理和分析。
三、实验原理1、信号的分类信号可以分为连续时间信号和离散时间信号。
连续时间信号在时间上是连续的，其数学表示通常为函数形式；离散时间信号在时间上是离散的，通常用序列来表示。
常见的信号类型包括正弦信号、方波信号、脉冲信号等。
2、线性时不变系统线性时不变系统具有叠加性和时不变性。
叠加性意味着多个输入信号的线性组合产生的输出等于各个输入单独作用产生的输出的线性组合；时不变性表示系统的特性不随时间变化，即输入信号的时移对应输出信号的相同时移。
3、卷积运算卷积是信号处理中一种重要的运算，用于描述线性时不变系统对输入信号的作用。
对于两个信号 f(t) 和 g(t)，它们的卷积定义为：＼（f g)（t) ＝＼int_{＼infty}＾｛＼infty} f(＼tau) g(t ＼tau) d\tau ＼在离散时间情况下，卷积运算为：＼（f g)n ＝＼sum_{m ＝＼infty}＾｛＼infty} fm gn m ＼四、实验内容及步骤实验一：常见信号的产生与观测1、连接信号发生器和示波器。
2、设置信号发生器分别产生正弦波、方波和脉冲信号，调整频率、幅度和占空比等参数。
3、在示波器上观察并记录不同信号的波形、频率和幅度。
z变换实验报告
南昌大学实验报告(信号与系统)学生姓名：肖江学号：6100210030 专业班级：电子103班实验类型：□验证□综合□设计□创新实验日期：2012/6/1 实验成绩：Z变换、离散时间系统的Z域分析一、实验目的1、学会用matlab求解z变换与逆z变换。
2、学会离散系统零极点分布图的绘制，理解离散系统零极点分布图的含义。
3、求解离散系统的频率响应特性。
二、实验说明1、一离散系统的差分方程为y(n)-by(n-1)=x(n),若激励为x(n)=a n u(n)，起始值y(-1)=0,求响应y(n)。
2、当H(s)极点位于z平面中各方框附近的位置，画出对应的h(n)波形填入方框中。
3、求系统差分方程为y(n)-1.1y(n-1)+0.7y(n-2)=x(n-1),的系统的频率响应特性。
三、实验内容1、syms n a b z%定义符号n a b zx=a^n; %定义激励信号X=ztrans(x); %计算激励信号的变换H=1/(1-b*z^(-1)); %写出系统z变换式Y=H*X; %计算输出的变换式y1=iztrans(Y); %计算输出时域表达式y=simplify(y1) %化简表达式2、pos=[26,19,18,17,24,27,13,11,9,23,28,7,4,1,22];figure,id=1; %生成新图框，子图id初始化为1for r=0.8:0.2:1.2 %极点的幅度依次为0.8,1.0,1.2for theta=0:pi/4:pi %极点的弧度依次为0，Π/4，Π/2，3Π/4,Πp=r*exp(j*theta);if theta~=0&theta~=pip=[p;p']; %如果极点不在实轴上添加一个共轭极点end[b a]=zp2tf([],p,1); %由零极点得到传递函数subplot(4,7,pos(id));[h,t]=impz(b,a,20); %计算20个点的单位样值响应stem(t,h,'k-','MarkerSize',5);%绘制单位样值响应id=id+1; %子图序号加1end%退出弧角循环end%退出幅度循环3、a=[1,-1.1,0.7];b=[0,1];subplot(2,1,1),zplane(b,a); %绘制零极点分布图subplot(2,1,2),impz(b,a); %绘制单位样值响应figure,freqz(b,a) %绘制频率特性4、a=[1,-1.1,0.6];b=[0.6,-1.1,1];subplot(2,1,1),zplane(b,a); %绘制零极点分布图subplot(2,1,2),impz(b,a); %绘制单位样值响应figure,freqz(b,a); %绘制频率响应n=[0:40]'; %生成时间点x1=sin(0.1*pi*n); %生成单频信号x2=0*n; %准备方波信号x2(mod(n,10)<5)=1; %生成周期为10的方波信号y1=filter(b,a,x1); %分别对两个信号滤波y2=filter(b,a,x2);figuresubplot(2,1,1),stem(n,x1); %绘制单频信号及其输出波形subplot(2,1,2),stem(n,y1);figuresubplot(2,1,1),stem(n,x2); %绘制方波信号及其输出波形subplot(2,1,2),stem(n,y2);四、实验结果1、y =(a^(1+n)-b^(1+n))/(a-b)2、输出波形如下3、输出波形如下：4、输出波形如下：五、实验总结通过本次实验的学习，对离散系统有了更多的了解，通过用matlab画出离散系统的零极点分布图，使我对离散系统的零极点分布与其对用的频响特性有了深刻的了解；同时对全通网络的相频失真有了进一步了解，幅度没有失真，但对不同的频率信号的相移不同，因此单频信号输入时，其输出信号的波形没有失真，只是整个波形发生了移位，但对于方波信号，由于其中包含了各种频率的信号，因此不同频率的信号相频失真不同，因此输出波形不再是方波。
信号与系统实验报告
信号与系统实验报告
实验名称：信号与系统实验
一、实验目的：
1.了解信号与系统的基本概念
2.掌握信号的时域和频域表示方法
3.熟悉常见信号的特性及其对系统的影响
二、实验内容：
1.利用函数发生器产生不同频率的正弦信号，并通过示波器观察其时域和频域表示。
2.通过软件工具绘制不同信号的时域和频域图像。
3.利用滤波器对正弦信号进行滤波操作，并通过示波器观察滤波前后信号的变化。
三、实验结果分析：
1.通过实验仪器观察正弦信号的时域表示，可以看出信号的振幅、频率和相位信息。
2.通过实验仪器观察正弦信号的频域表示，可以看出信号的频率成分和幅度。
3.利用软件工具绘制信号的时域和频域图像，可以更直观地分析信号的特性。
4.经过滤波器处理的信号，可以通过示波器观察到滤波前后的信号波形和频谱的差异。
四、实验总结：
通过本次实验，我对信号与系统的概念有了更深入的理解，掌
握了信号的时域和频域表示方法。
通过观察实验仪器和绘制图像，我能够分析信号的特性及其对系统的影响。
此外，通过滤波器的处理，我也了解了滤波对信号的影响。
通过实验，我对信号与系统的理论知识有了更加直观的了解和应用。
南昌大学实验报告二
学生姓名： 学 号： 专业班级：实验类型：□ 验证 □ 综合 □ 设计 □ 创新 实验日期： 3.15 实验成绩： _MATLAB 基础上机训练二一、 实验项目名称：MATLAB 基础上机训练二二、 实验目的：1、熟悉掌握Matlab 的基本关系运算符，逻辑运算符，和简单的循环结构，流程控制。
2、掌握线性代数和矩阵的计算以及数据分析。
三、实验说明：第四章1、要求n ，使n ！是一个100位数字的最小值。
2、输入数据n ，判断其奇偶性。
第五章1、设⎪⎩⎪⎨⎧=+=++=++1,223,132********x x x x x x x x 利用两种方法求x 1,x 2以及x 3 。
2、设A=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡5210150110520125，试利用MATLAB 求其特征值与特征向量。
3、设B=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡1817161514131211109876543,试利用MATLAB 产生向量C 1=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡161284，C 2=[7 8 910],C 3=⎥⎦⎤⎢⎣⎡18171413,C 4=[4 8 13 18]学生姓名： 潘书敏 学 号： 6100210062 专业班级： 通信101 实验类型：□ 验证 □ 综合 □ 设计 □ 创新 实验日期： 3.15 实验成绩： _4、设D 1=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡161284，D 2=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡10987，试利用MATLAB 产生向量D 3=⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡10169128874，D 4=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡5412310987161284。
第六章1、设y 1=2x 4+4x 3+2x+1,y 2=2x 2+1,试利用MATLAB 求y 1+y 2,y 1-y 2,y 1÷y 2以及 y 1*y 2。
2、设y=2x 4+5x 3+4x+1,试利用MATLAB 求该多项式的根、dxdv以及在区间[-1,3]内100点的值并作图表示。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
实验报告课程名称：信号与系统实验指导老师：学生姓名：Du学号：6102213专业班级：通信工程133班2015年06月18 日目录实验一、利用matlab求解电路微分方程实验二、利用matlab进行傅里叶变换实验三、利用matlab简单识别男女声实验四、利用matlab画出离散系统零极点图实验五、利用matlab求出离散系统响应（1）实验六、利用matlab求出离散系统响应（2）实验七、利用matlab绘制各种信号实验八、利用matlab对信号进行频域分析和卷积实验九、软件实验期末考试实验十、周期信号频谱测试实验十一、模拟滤波器频率特性测试实验十二、连续时间系统模拟实验一：利用matlab求解电路微分方程一、实验目的：利用matlab对连续时间系统的时域进行分析，求解电路微分方程。
二、实验原理：1）先根据题目要求建立矩阵描述相应微分方程2）定义lsim函数所需系统sys为之后调用lsim做准备3）定义仿真时间，对时间进行抽样4）方法一：利用step函数求取系统单位阶跃响应：step()step()函数的用法：y=step(num,den,t)：其中num和den分别为系统传递函数描述中的分子和分母多项式系数，t为选定的仿真时间向量，一般可以由t=0:step:end等步长地产生出来。
该函数返回值y为系统在仿真时刻各个输出所组成的矩阵并绘图5）方法二：定义构造阶跃信号作为激励信号，定义t=0时u=1/2，利用lsim函数仿真阶跃响应lsim（）函数用法：lsim(sys,u,t):模拟了系统SYS对输入u持续t秒时间的响应并绘图6）方法一：利用impulse函数求取单位冲激响应：impulse（）Impulse（）函数用法Impulse（sys,t）可直接求出连续时间系统冲激响应并绘图7）方法二：定义构造冲激信号作为激励信号，定义t=0时delta=100，利用lsim函数仿真阶跃响应，从响应中减去一个delta并绘图三、实验任务：已知某系统微分方程为r’’+r’+r=e’+e分别用两种方法计算阶跃响应和冲激响应，结果对比四、实验内容：clear;clc;close all;a=[1,1,1];b=[1,1];sys=tf(b,a);t=[0:0.01:3];figure(1);subplot(2,2,1);step(sys);subplot(2,2,2);x_step=zeros(size(t));x_step(t>0)=1;x_step(t==0)=1/2;lsim(sys,x_step,t);subplot(2,2,3);[h1,t1]=impulse(sys,t);plot(t1,h1,'k');title('Impulse Response'); xlabel('Time(sec)');ylabel('Amplitude');subplot(2,2,4);x_delta=zeros(size(t));x_delta(t==0)=100;[y1,t] =lsim(sys,x_delta,t); y2=y1'-x_delta;plot(t,y2,'k');title('Impulse Response'); xlabel('Time(sec)');ylabel('Amplitude');五、实验结果：实验二、利用matlab进行傅里叶变换一、实验目的：利用matlab对指定信号进行傅里叶变换二、实验原理：1）定义信号时间并时域离散化抽样2）定义信号函数3）利用矩阵计算法进行傅里叶变换和逆变换并绘图三、实验任务：如图4.4所示锯齿波信号，分别取一个周期的抽样数据X1（t），0<=t<=1和五个周期的数据X（t），0<=t<5,计算其傅立叶变换X1（w）和X（w），比较有和不同并解释原因。
四、实验内容：clc;clear all;close all;T1=1;N1=2000;t1=linspace(0,T1-T1/N1,N1)';x1=0*t1;0<t1&t1<1;x1=sawtooth(t1*2*pi,0);OMG=16*pi;K=2000;omg=linspace(-OMG/2,OMG/2-OMG/K,K)'; F1=0*omg;fs1=0*t1;U=exp(-1j*kron(omg,t1.'));F1=T1/N1*U*x1;fs1=OMG/2/pi/K*exp(1j*kron(t1,omg.'))*F1;T2=5;N2=2000;t2=linspace(0,T2-T2/N2,N2)';x=0*t2;0<t2&t2<5;x=sawtooth(t2*2*pi,0);F2=0*omg;fs2=0*t2;U=exp(-1j*kron(omg,t2.'));F2=T2/N2*U*x;fs2=OMG/2/pi/K*exp(1j*kron(t2,omg.'))*F2;figure(1);subplot(2,2,1);plot(omg,abs(F1),'r');xlabel('Frequency');ylabel('Amplitude');title('µ¥ÖÜÆÚ¾â³Ý²¨·ùƵÌØÐÔÇúÏß'); subplot(2,2,2);plot(t1,fs1,'r');xlabel('Time');ylabel('Amplitude');title('µ¥ÖÜÆÚ¾â³Ý²¨±ä»»');subplot(2,2,3);plot(omg,abs(F2),'r');xlabel('Frequency');ylabel('Amplitude');title('ÎåÖÜÆÚ¾â³Ý²¨·ùƵÌØÐÔÇúÏß'); subplot(2,2,4);plot(t2,fs2,'r');xlabel('Time(sec)');ylabel('Amplitude');title('ÎåÖÜÆÚ¾â³Ý²¨±ä»»');五、实验结果实验三、利用matlab简单识别男女声一、实验目的：利用matlab对声音信号进行傅里叶变换，通过基频的比较区分男女声二、实验原理：调用音频文件，stream为音频数据，sample_rate为采样率取双声道中声道1，从1到音频数据总长度作为一个新数组对这个新数组进行快速傅里叶变换，得到幅度特性对频域1到中值的幅度进行归一化找到最大幅度对应频率，最大振幅的频率是基频的整数倍倍数可能是从1到frequency_max，要提取出需要的倍数最大幅度的频率除以倍数即为我们所需要的基频根据广泛调查，女声男声基频的临界值为220，进行比较，得出结论三、实验内容：close all;clear all;clc;[stream, sample_rate] = audioread('442.wav');stream = stream(1 : length(stream), 1);frequency = abs(fft(stream, sample_rate));len = floor(length(frequency) / 2);frequency = frequency(1 : len) / max(frequency);[~, frequency_max] = max(frequency);frequency_sum = zeros(frequency_max);for i = 1 : frequency_maxstep = frequency_max / i;index = ceil(step : step : frequency_max);frequency_sum(i) = sum(frequency(index)) * step;endchange = zeros(1, frequency_max);for i = 1 : length(frequency_sum) - 1change(i) = frequency_sum(i) - frequency_sum(i + 1);endchange_sum = zeros(1, length(frequency_sum));for i = 1 : length(change_sum)index = i : i : length(change_sum);change_sum(i) = sum(change(index)) * i;end[~, change_sum_max] = max(change_sum);frequency_base = ceil(frequency_max / change_sum_max); disp(frequency_base);limit=200;if frequency_base <limitdisplay '男'elsedisplay '女'end四、实验结果：1)2)实验四、利用matlab画出离散系统零极点图一、实验目的：利用matlab计算出离散系统零极点、画出离散系统零极点图二、实验任务：第八章练习题1三、实验程序：clc;clear all;close all;b=[0,-3];a=[2,-5,2];[z,p,k]=residuez(b,a)[b,a]=zp2tf(z,p,1);figure;zplane(b,a);syms z n;X=((-3)*z^(-1))/(2-5*z^(-1)+2*z^(-2));x=iztrans(X);x1=x*heaviside(n-2)三、实验结果：z =-11p =2.00000.5000k =[ ]x1 =-heaviside(n - 2)*(2^n - (1/2)^n)实验五、利用matlab求出离散系统响应一、实验目的利用matlab求出离散系统单位样值响应、零状态响应二、实验任务第八章练习题2三、实验程序clc;clear all;close all;syms n z;a=[1,3,-1];b=[1];figure;subplot(2,1,1),zplane(b,a);subplot(2,1,2),impz(b,a);x=(n+n^2)*(heaviside(n)-heaviside(n-6));X=ztrans(x);H=1/(1+3*z^(-1)-z^(-2));Y=H*X;y1=iztrans(Y);y=simplify(y1)h=iztrans(H);四、实验结果y =(12697*(-1)^n*2^(1 - n)*13^(1/2)*(13^(1/2) + 3)^(n - 1))/39 - 320*kroneckerDelta(n - 2, 0) - 93*kroneckerDelta(n - 3, 0) - 21*kroneckerDelta(n - 4, 0) - (2134*(-1)^((3*n)/2)*cos(n*(pi/2 + asinh(3/2)*i)))/3 - 3521*kroneckerDelta(n, 0) - 1065*kroneckerDelta(n - 1, 0) - (12697*(-1)^n*2^(1 - n)*13^(1/2)*(3 - 13^(1/2))^(n - 1))/39实验六、利用matlab求出离散系统响应一、实验目的：已知零极点，绘出零极点图。