内蒙工业大学数字信号处理实验一
数字信号处理实验报告
实验一 信号、系统及系统响应一、实验目的1、熟悉理想采样的性质，了解信号采样前后的频谱变化，加深对时域采样定理的理解。
2、熟悉离散信号和系统的时域特性。
3、熟悉线性卷积的计算编程方法：利用卷积的方法，观察、分析系统响应的时域特性。
4、掌握序列傅里叶变换的计算机实现方法，利用序列的傅里叶变换对离散信号、系统及其系统响应进行频域分析。
二、 实验原理1.理想采样序列：对信号x a (t)=A e −αt sin(Ω0t )u(t)进行理想采样，可以得到一个理想的采样信号序列x a (t)=A e −αt sin(Ω0nT ),0≤n ≤50，其中A 为幅度因子，α是衰减因子，Ω0是频率，T 是采样周期。
2.对一个连续时间信号x a (t)进行理想采样可以表示为该信号与一个周期冲激脉冲的乘积，即x ̂a (t)= x a (t)M(t),其中x ̂a (t)是连续信号x a (t)的理想采样；M(t)是周期冲激M(t)=∑δ+∞−∞(t-nT)=1T ∑e jm Ωs t +∞−∞,其中T 为采样周期，Ωs =2π/T 是采样角频率。
信号理想采样的傅里叶变换为X ̂a (j Ω)=1T ∑X a +∞−∞[j(Ω−k Ωs )],由此式可知：信号理想采样后的频谱是原信号频谱的周期延拓，其延拓周期为Ωs =2π/T 。
根据时域采样定理，如果原信号是带限信号，且采样频率高于原信号最高频率分量的2倍，则采样以后不会发生频率混叠现象。
三、简明步骤产生理想采样信号序列x a (n),使A=444.128，α=50√2π，Ω0=50√2π。
（1） 首先选用采样频率为1000HZ ，T=1/1000，观察所得理想采样信号的幅频特性，在折叠频率以内和给定的理想幅频特性无明显差异，并做记录；（2） 改变采样频率为300HZ ，T=1/300，观察所得到的频谱特性曲线的变化，并做记录；（3） 进一步减小采样频率为200HZ ，T=1/200，观察频谱混淆现象是否明显存在，说明原因，并记录这时候的幅频特性曲线。
数字信号处理实验
数字信号处理实验实验一信号、系统及系统响应1、实验目的认真复习采样理论、离散信号与系统、线性卷积、序列的z 变换及性质等有关内容；掌握离散时间序列的产生与基本运算，理解离散时间系统的时域特性与差分方程的求解方法，掌握离散信号的绘图方法；熟悉序列的z 变换及性质，理解理想采样前后信号频谱的变化。
2、实验内容a. 产生长度为500 的在[0,1]之间均匀分布的随机序列，产生长度为500 的均值为0 单位方差的高斯分布序列。
b. 线性时不变系统单位脉冲响应为h(n)=(0.9)nu(n)，当系统输入为x(n)=R10(n)时，求系统的零状态响应，并绘制波形图。
c. 描述系统的差分方程为：y(n)-y(n-1)+0.9y(n-2)=x(n)，其中x(n)为激励，y(n)为响应。
计算并绘制n=20,30,40,50,60,70,80,90,100 时的系统单位脉冲响应h(n)；计算并绘制n=20,30,40,50,60,70,80,90,100 时的系统单位阶跃响应s(n)；由h(n)表征的这个系统是稳定系统吗？d. 序列x(n)=(0.8)nu(n)，求DTFT[x(n)]，并画出它幅度、相位，实部、虚部的波形图。
观察它是否具有周期性？e. 线性时不变系统的差分方程为y(n)=0.7y(n-1)+x(n)，求系统的频率响应H(ejω)，如果系统输入为x(n)=cos(0.05πn)u(n)，求系统的稳态响应并绘图。
f. 设连续时间信号x(t)=e-1000|t|，计算并绘制它的傅立叶变换；如果用采样频率为每秒5000 样本对x(t)进行采样得到x1(n)，计算并绘制X1(ejω)，用x1(n)重建连续信号x(t)，并对结果进行讨论；如果用采样频率为每秒1000 样本对x(t)进行采样得到x2(n)，计算并绘制X2(ejω)，用x2(n)重建连续信号x(t)，并对结果进行讨论。
加深对采样定理的理解。
g. 设X1(z)=z+2+3z-1，X2(z)=2z2+4z+3+5z-1，用卷积方法计算X1(z)X2(z)。
数字信号处理实验报告_五个实验
实验一 信号、系统及系统响应一、 实验目的1、熟悉连续信号经理想采样前后的频谱变化关系，加深对时域采样定理的理解；2、熟悉时域离散系统的时域特性；3、利用卷积方法观察分析系统的时域特性；4、掌握序列傅立叶变换的计算机实现方法，利用序列的傅立叶变换对连续信号、离散信号及系统响应进行频域分析。
二、 实验原理及方法采样是连续信号数字处理的第一个关键环节。
对采样过程的研究不仅可以了解采样前后信号时域和频域特性发生变化以及信号信息不丢失的条件，而且可以加深对傅立叶变换、Z 变换和序列傅立叶变换之间关系式的理解。
对一个连续信号)(t x a 进行理想采样的过程可用下式表示：)()()(^t p t t xx aa=其中)(^t x a 为)(t x a 的理想采样，p(t)为周期脉冲，即∑∞-∞=-=m nT t t p )()(δ)(^t x a的傅立叶变换为)]([1)(^s m a m j X T j a XΩ-Ω=Ω∑∞-∞=上式表明^)(Ωj Xa为)(Ωj Xa的周期延拓。
其延拓周期为采样角频率（T /2π=Ω）。
只有满足采样定理时，才不会发生频率混叠失真。
在实验时可以用序列的傅立叶变换来计算^)(Ωj X a 。
公式如下：Tw jw ae X j X Ω==Ω|)()(^离散信号和系统在时域均可用序列来表示。
为了在实验中观察分析各种序列的频域特性，通常对)(jw e X 在[0，2π]上进行M 点采样来观察分析。
对长度为N 的有限长序列x(n)，有：n jw N n jw k ke m x eX--=∑=)()(1其中，k Mk πω2=，k=0,1,……M-1 时域离散线性非移变系统的输入/输出关系为 ∑∞-∞=-==m m n h m x n h n x n y )()()(*)()(上述卷积运算也可在频域实现)()()(ωωωj j j e H e X eY =三、 实验程序s=yesinput(Please Select The Step Of Experiment:\n 一.(1时域采样序列分析 s=str2num(s); close all;Xb=impseq(0,0,1); Ha=stepseq(1,1,10);Hb=impseq(0,0,3)+2.5*impseq(1,0,3)+2.2*impseq(2,0,3)+impseq(3,0,3); i=0;while(s);%时域采样序列分析 if(s==1) l=1; k=0;while(1)if(k==0)A=yesinput('please input the Amplitude:\n',...444.128,[100,1000]); a=yesinput('please input the Attenuation Coefficient:\n',...222.144,[100,600]); w=yesinput('please input the Angle Frequence(rad/s):\n',...222.144,[100,600]); end k=k+1;fs=yesinput('please input the sample frequence:\n',...1000,[100,1200]); Xa=FF(A,a,w,fs); i=i+1;string+['fs=',num2str(fs)]; figure(i)DFT(Xa,50,string); 1=yesinput 1=str2num(1); end%系统和响应分析else if(s==2)kk=str2num(kk);while(kk)if(kk==1)m=conv(Xb,Hb);N=5;i=i+1;figure(i)string=('hb(n)');Hs=DFT(Hb,4,string);i=i+1;figure(i)string('xb(n)');DFT(Xb,2,string);string=('y(n)=xb(n)*hb(n)');else if (kk==2)m=conv(Ha,Ha);N=19;string=('y(n)=ha(n)*(ha(n)');else if (kk==3)Xc=stepseq(1,1,5);m=conv(Xc,Ha);N=14;string=('y(n)=xc(n)*ha(n)');endendendi=i+1;figure(i)DFT(m,N,string);kk=yesinputkk=str2num(kk);end卷积定理的验证else if(s==3)A=1;a=0.5;w=2,0734;fs=1;Xal=FF(A,a,w,fs);i=i+1;figure(i)string=('The xal(n)(A=1,a=0.4,T=1)'); [Xa,w]DFT(Xal,50,string);i=i+1;figure(i)string =('hb(n)');Hs=DFT(Hb,4,string);Ys=Xs.*Hs;y=conv(Xal,Hb);N=53;i=i+1;figure(i)string=('y(n)=xa(n)*hb(n)');[yy,w]=DFT(y,N,string);i=i+1;figure(i)subplot(2,2,1)plot(w/pi,abs(yy));axis([-2 2 0 2]);xlabel('w/pi');ylabel('|Ys(jw)|');title(FT[x(n)*h(n)]');subplot(2,2,3)plot(w/pi,abs(Ys));axis([-2 2 0 2]);xlabel('w/pi');ylabel('|Ys(jw)|');title('FT[xs(n)].FT[h(n)]');endendend子函数:离散傅立叶变换及X(n),FT[x(n)]的绘图函数function[c,l]=DFT(x,N,str)n=0:N-1;k=-200:200;w=(pi/100)*k;l=w;c=x*Xc=stepseq(1,1,5);子函数：产生信号function c=FF(A,a,w,fs)n=o:50-1;c=A*exp((-a)*n/fs).*sin(w*n/fs).*stepseq(0,0,49); 子函数：产生脉冲信号function [x,n]=impseq(n0,n1,n2)n=[n1:n2];x=[(n-n0)==0];子函数：产生矩形框信号function [x,n]=stepseq(n0,n1,n2) n=[n1:n2];x=[(n-n0>=0)];四、 实验内容及步骤1、认真复习采样理论，离散信号与系统，线性卷积，序列的傅立叶变换及性质等有关内容，阅读本实验原理与方法。
数字信号处理实验(1-7)原始实验内容文档(含代码)
实验要求1.每个实验进行之前须充分预习准备，实验完成后一周内提交实验报告；2.填写实验报告时，分为实验题目、实验目的、实验内容、实验结果、实验小结五项；3.实验报告要求：实验题目、实验目的、实验内容、实验结果四项都可打印；但每次实验的实验内容中的重要代码（或关键函数）后面要用手工解释其作用。
实验小结必须手写！（针对以前同学书写实验报告时候抄写代码太费时间的现象，本期实验报告进行以上改革）。
实验一信号、系统及系统响应实验目的：1. 掌握使用MATLAB进行函数、子程序、文件编辑等基本操作；2. 编写一些数字信号处理中常用序列的3. 掌握函数调用的方法。
实验内容：1.在数字信号处理的基本理论和MATLAB信号处理工具箱函数的基础上，可以自己编写一些子程序以便调用。
（1）单位抽样序列δ(n-n0)的生成函数impseq.m（2）单位阶跃序列u(n-n0)的生成函数stepseq.m（3）两个信号相加的生成函数sigadd.m（4）两个信号相乘的生成函数sigmult.m（5）序列移位y(n)=x(n-n0)的生成函数sigshift.m（6）序列翻褶y(n)=x(-n)生成函数sigfold.m（7）奇偶综合函数evenodd.m（8）求卷积和2.产生系列序列，并绘出离散图。
（1） x1(n)=3δ(n-2)-δ(n+4) -5≤n≤5（2） x3(n)=cos(0.04πn)+0.2w(n) 0≤n≤50其中：w(n)是均值为0，方差为1 的白噪声序列。
3.设线性移不变系统的抽样响应h(n)=(0.9)^n*u(n),输入序列x(n)=u(n)-u(n-10),求系统的输出y(n).实验二 系统响应及系统稳定性1.实验目的（1）掌握 求系统响应的方法。
（2）掌握时域离散系统的时域特性。
（3）分析、观察及检验系统的稳定性。
2.实验原理与方法在时域中，描写系统特性的方法是差分方程和单位脉冲响应，在频域可以用系统函数描述系统特性。
数字信号处理实验报告
数字信号处理实验报告引言数字信号处理（Digital Signal Processing，DSP）是一门研究数字信号的获取、分析、处理和控制的学科。
在现代科技发展中，数字信号处理在通信、图像处理、音频处理等领域起着重要的作用。
本次实验旨在通过实际操作，深入了解数字信号处理的基本原理和实践技巧。
实验一：离散时间信号的生成与显示在实验开始之前，我们首先需要了解信号的生成与显示方法。
通过数字信号处理器（Digital Signal Processor，DSP）可以轻松生成和显示各种类型的离散时间信号。
实验设置如下：1. 设置采样频率为8kHz。
2. 生成一个正弦信号：频率为1kHz，振幅为1。
3. 生成一个方波信号：频率为1kHz，振幅为1。
4. 将生成的信号通过DAC（Digital-to-Analog Converter）输出到示波器上进行显示。
实验结果如下图所示：（插入示波器显示的正弦信号和方波信号的图片）实验分析：通过示波器的显示结果可以看出，正弦信号在时域上呈现周期性的波形，而方波信号则具有稳定的上下跳变。
这体现了正弦信号和方波信号在时域上的不同特征。
实验二：信号的采样和重构在数字信号处理中，信号的采样是将连续时间信号转化为离散时间信号的过程，信号的重构则是将离散时间信号还原为连续时间信号的过程。
在实际应用中，信号的采样和重构对信号处理的准确性至关重要。
实验设置如下：1. 生成一个正弦信号：频率为1kHz，振幅为1。
2. 设置采样频率为8kHz。
3. 对正弦信号进行采样，得到离散时间信号。
4. 对离散时间信号进行重构，得到连续时间信号。
5. 将重构的信号通过DAC输出到示波器上进行显示。
实验结果如下图所示：（插入示波器显示的连续时间信号和重构信号的图片）实验分析：通过示波器的显示结果可以看出，重构的信号与原信号非常接近，并且能够还原出原信号的形状和特征。
这说明信号的采样和重构方法对于信号处理的准确性有着重要影响。
《数字信号处理》上机实验指导书
《数字信号处理》上机实验指导书实验1 离散时间信号的产生1．实验目的数字信号处理系统中的信号都是以离散时间形态存在，所以对离散时间信号的研究是数字信号处理的基本所在。
而要研究离散时间信号，首先需要产生出各种离散时间信号。
MATLAB 是一套功能强大的工程计算及数据处理软件，广泛应用于工业，电子，医疗和建筑等众多领域。
使用MATLAB软件可以很方便地产生各种常见的离散时间信号，而且它还具有强大的绘图功能，便于用户直观地输出处理结果。
通过本实验，学生将学习如何用MATLAB产生一些常见的离散时间信号，并通过MATLAB中的绘图工具对产生的信号进行观察，加深对常用离散信号的理解。
2．实验要求本实验要求学生运用MATLAB编程产生一些基本的离散时间信号，并通过MATLAB的几种绘图指令画出这些图形，以加深对相关教学内容的理解，同时也通过这些简单的函数练习了MATLAB的使用。
3．实验原理（1）常见的离散时间信号1）单位抽样序列，或称为离散时间冲激，单位冲激：?(n)???1?0n?0 n?0如果?(n)在时间轴上延迟了k个单位，得到?(n?k)即：?1n?k ?(n?k)??0n?0?2）单位阶跃序列n?0?1 u(n)?n?0?0在MATLAB中可以利用ones( )函数实现。
x?ones(1,N);3）正弦序列x(n)?Acos(?0n??)这里，A,?0,和?都是实数，它们分别称为本正弦信号x(n)的振幅，角频率和初始相位。
f0??02?为频率。
x(n)?ej?n4）复正弦序列5）实指数序列x(n)?A?n（2）MATLAB编程介绍MATLAB是一套功能强大，但使用方便的工程计算及数据处理软件。
其编程风格很简洁，没有太多的语法限制，所以使用起来非常方便，尤其对初学者来说，可以避免去阅读大量的指令系统，以便很快上手编程。
值得注意得就是，MATLAB中把所有参与处理的数据都视为矩阵，并且其函数众多，希望同学注意查看帮助，经过一段时间的训练就会慢慢熟练使用本软件了。
课程大作业——数字信号处理实验报告
实验一 信号、系统及系统响应一．实验目的1．熟悉理想采样的性质，了解信号采用前后的频谱变化，加深对采样定理的理解。
2．熟悉离散信号和系统的时域特性。
3．熟悉线性卷积的计算编程方法：利用卷积的方法，观察、分析系统响应的时域特性。
4．掌握序列傅氏变换的计算机实现方法，利用序列的傅氏变换对离散信号、系统及系统响应进行频域分析。
二．实验原理1．连续时间信号的采样采样是从连续时间信号到离散时间信号的过渡桥梁，对采样过程的研究不仅可以了解采样前后信号时域和频域特性发生的变化以及信号内容不丢失的条件，而且有助于加深对拉氏变换、傅氏变换、z 变换和序列傅氏变换之间关系的理解。
对一个连续时间信号进行理想采样的过程可以表示为该信号和个周期冲激脉冲的乘积，即)()()(ˆt M t x t xa a = （1－1） 其中)(ˆt xa 是连续信号)(t x a 的理想采样，)(t M 是周期冲激脉冲 ∑+∞-∞=-=n nT t t M )()(δ （1－2）它也可以用傅立叶级数表示为：∑+∞-∞=Ω=n tjm s e T t M 1)( （1－3）其中T 为采样周期，T s /2π=Ω是采样角频率。
设)(s X a 是连续时间信号)(t x a 的双边拉氏变换，即有：⎰+∞∞--=dt e t xs X st aa )()( （1－4）此时理想采样信号)(ˆt xa 的拉氏变换为 ∑⎰+∞-∞=+∞∞--Ω-===m s a sta a jm s X T dt e t x s X )(1)(ˆ)(ˆ （1－5）作为拉氏变换的一种特例，信号理想采样的傅立叶变换[]∑+∞-∞=Ω-Ω=Ωm s a a m j X T j X )(1)(ˆ （1－6）由式（1－5）和式（1－6）可知，信号理想采样后的频谱是原信号频谱的周期延拓，其延拓周期等于采样频率。
根据Shannon 采样定理，如果原信号是带限信号，且采样频率高于原信号最高频率分量的2倍，则采样以后不会发生频率混淆现象。
数字信号处理实验一
数字信号处理实验一实验目的：掌握利用Matlab产生各种离散时间信号，实现信号的相加、相乘及卷积运算实验函数：参考课本77-19页，注意式(2.11.1)的表达与各matlab子函数间的关系。
1、stem(x,y) % 绘制以x为横轴，y为纵轴的离散序列图形2、[h ,t] = impz(b, a) % 求解数字系统的冲激响应h，取样点数为缺省值[h, t] = impz(b, a, n) % 求解数字系统的冲激响应h，取样点数为nimpz(b, a) % 在当前窗口用stem(t, h)函数出图3、[h ,t] = dstep(b, a) % 求解数字系统的阶跃响应h，取样点数为缺省值[h, t] = dstep (b, a, n) % 求解数字系统的阶跃响应h，取样点数为ndstep (b, a) % 在当前窗口用stairs(t, h)函数出图4、y = filter(b,a,x) % 在已知系统差分方程或转移函数的情况下求系统输出实验原理：一、常用的时域离散信号及其程序1、产生单位抽样函数δ(n)n1 = -5;n2 = 5;n0 = 0;n = n1:n2;x = [n==n0]; % x在n=n0时为1，其余为0stem(n,x,'filled'); %filled:序列圆心处用实心圆表示axis([n1,n2,0,1.1*max(x)])title('单位抽样序列')xlabel('time(n)')ylabel('Amplitude:x(n)')2、产生单位阶跃序列u(n)n1 = -2;n2 = 8;n0 = 0;n = n1:n2;x = [n>=n0]; % x在n>=n0时为1，其余为0stem(n,x,'filled');axis([n1,n2,0,1.1*max(x)])title('单位阶跃序列')xlabel('time(n)')ylabel('Amplitude:x(n)')3、复指数序列复指数序列的表示式为()(),00,0j n e n x n n σω+⎧≥⎪=⎨<⎪⎩，当0ω=时，()x n 为实指数序列；当0σ=时，()x n 为虚指数序列，即()()cos sin j n e n j n ωωω=+，即其实部为余弦序列，虚部为正弦序列。
数字信号处理实验一
实验一 离散时间信号分析一、实验目的1．熟悉MATLAB 应用环境，常用窗口的功能和使用方法。
2．掌握各种常用的序列，理解其数学表达式和波形表示。
3．掌握在计算机中生成及绘制数字信号波形的方法。
4．掌握序列的相加、相乘、移位、反褶、卷积等基本运算及计算机实现。
5．通过编程，上机调试程序，进一步增强使用计算机解决问题的能力。
二、实验原理1．序列的基本概念离散时间信号是指在离散时刻才有定义的信号，简称离散信号，或者序列。
离散时间信号在数学上可用时间序列)}({n x 来表示，其中)(n x 代表序列的第n 个数字，n 代表时间的序列，n 的取值范围为∞<<∞-n 的整数，n 取其它值)(n x 没有意义。
离散时间信号可以是由模拟信号通过采样得到，例如对模拟信号)(t a x 进行等间隔采样，采样间隔为T ，得到一个有序的数字序列)}({nT x a 就是离散时间信号，简称序列。
2．常用序列常用序列有：单位脉冲序列（单位抽样）)(n δ、单位阶跃序列)(n u 、矩形序列)(n R N 、实指数序列、复指数序列、正弦型序列等。
3．序列的基本运算序列的运算包括移位、反褶、和、积、点乘、累加、差分运算、卷积等。
4．序列的卷积运算)()()()()(n h n x m n h m x n y m *=-=∑∞-∞=上式的运算关系称为卷积运算，式中*代表两个序列卷积运算。
两个序列的卷积是一个序列与另一个序列反褶后逐次移位乘积之和，故称为离散卷积，也称两序列的线性卷积。
其计算的过程包括以下4个步骤（1）反褶：先将)(n x 和)(n h 的变量n 换成m ，变成)(m x 和)(m h ，再将)(m h 以纵轴为对称轴反褶成)(m h -。
（2）移位：将)(m h -移位n ，得)(m n h -。
当n 为正数时，右移n 位；当n 为负数时，左移n 位。
（3）相乘：将)(m n h -和)(m x 的对应点的值相乘。
数字信号处理实验1--5含代码
数字信号处理实验1--5含代码实验一离散时间信号的时域分析 1. 在MATLAB中利用逻辑关系式n,,0来实现序列，显示范围。
(产生如下，，,n,nn,n,n012图所示的单位脉冲信号的函数为impseq(n0,n1,n2)，程序如示例所示),3,n,10并利用impseq函数实现序列:; ，，，，，，yn,2,n,3，,n,6，，xn1nnnn120源代码:impseq.mfunction y=impseq(n0,n1,n2)n=[n1:n2]y=[(n-n0)==0]exp01-1.mfunction impseq(n0,n1,n2)n=-3:1:10y=2*impseq(3,-3,10)+impseq(6,-3,10);stem(n,y)n,,0，，2. 在MATLAB中利用逻辑关系式来实现序列，显示范围。
(自己编写un,nn,n,n012产生单位阶跃信号的函数，函数命名为stepseq(n0,n1,n2)) 并利用编写的stepseq函数实现序列: ，，，，，，yn,un，2，un,2,5,n,10源代码:stepseq.mfunction y=stepseq(n0,n1,n2)n=n1:1:n2y=[(n-n0)>=0]exp01-2.mfunction stepseq(n0,n1,n2)n=-5:1:20y=stepseq(-2,-5,20)+stepseq(2,-5,20)stem(n,y)3. 在MATLAB中利用数组运算符“.^”来实现一个实指数序列。
如: n ，，，，xn,0.30,n,15源代码:n=0:1:15;x=0.3.^nstem(n,x)4. 在MATLAB中调用函数sin或cos产生正余弦序列，如:π,, ，，，，xn,3sin0.4πn，，5cos0.3πn0,n,20,,5,,源代码:n=0:1:20x=11*sin(0.3*pi*n+pi/5)+5*cos(0.3*pi*n)stem(n,x)思考题:1.在MATLAB环境下产生单位脉冲序列和单位阶跃序列各有几种方法,如何使用,2.在MATLAB环境下进行序列的相乘运算时应注意什么问题,实验二离散时间系统的时域分析1. 在MATLAB中利用内部函数conv来计算两个有限长序列的卷积。
数字信号处理实验一
实验一：熟悉MATLAB环境一、思考题：1、比较实验内容第2题中的第4和第5小题的结果，试说明对于周期性信号应当如何采样才能保证周期扩展后与原信号保持一致？答：对于周期性信号，在进行采样时，其采样周期必须满足采样定理，即采样频率应该大于信号最高频率的两倍，这样才能避免迭混，以便采样后仍能准确的恢复原信号。
2、对于有限长序列，如何用MATLAB计算其DTFT？答：用函数freqz可以计算序列在给定的离散频率点上的DTFT，该变换序列是以形如式的有理函数来描述的。
这个函数的表达形式有H=freqz(num,den,w)、[H，w]=freqz(num,den,k)。
函数freqz返回的频率响应值为向量H。
在H=freqz(num,den,w)中，0到π之间指定的频率集由向量w给出。
freqz函数的自变量k就是频率点的总数。
3、对于由两个子系统级联或并联的系统，如何用MATLAB计算他们的幅频响应与相频响应？答：系统的级联或并联实现涉及到了因式分解。
在MATLAB中，我们可以用函数roots来实现多项式的因式分解。
例如，函数r=roots(h)会返回多项式向量h 的根向量。
向量h是以的升幂表示的多项式的系数。
通过计算所得的根向量，可以求出二次因式的系数。
更简单的方法是用从以给定的传输函数H(z)直接求出二阶因式的函数zp2sos。
函数sos=zp2sos(z,p,k)产生以零—极点形式确定的等效传输函数H(z)的每个二阶部分系数的矩阵sos。
二、实验内容：第一到四题源程序：第二题图：第五题：%函数命名：function [x1]=stepshift(n0,n1,n2) n=[n1:n2];x1=[(n-n0)>=0];%源程序：>> [x1]=stepshift(2,1,10);>> n=1:10;>> stem(n,x1);第六题：>> b=[1,sqrt(2),1];>> a=[1,-0.67,0.9];>> [h,w]=freqz(b,a);>> am=20*log10(abs(h)); %求幅频特性>> subplot(2,1,1);plot(w,am);>> xlabel('w');ylabel('am');>> ph=angle(h); %求相频特性>> subplot(2,1,2);plot(w,ph);>> xlabel('w');ylabel('ph');第七题：源程序：>> a=[8 -2 -1 2 3];>> b=[2 3 -1 -3];>> c=conv(a,b); %求a、b的线性卷积>> m=length(c)-1;>> n=0:1:m;>> stem(n,c);>> xlabel('n');ylabel('幅度');第八题：源程序：>> n=50;>> a=[1 -2];>> b=[1 0.1 -0.06];>> x=[1 zeros(1,n-1)];>> k=0:1:n-1;>> y=filter(a,b,x);>> stem(k,y);>> xlabel('n');ylabel('幅度');三、实验总结：通过本次实验，熟悉并掌握MATLAB的主要命令操作，比如序列的简单运算、矩阵的输入和计算等，能熟练编写绘图程序，在计算卷积和绘制幅频响应和相频响应的过程中，充分地巩固了数字信号处理学的理论知识，总之，收获颇多。
数字信号处理实验报告_完整版
实验1 利用DFT 分析信号频谱一、实验目的1.加深对DFT 原理的理解。
2.应用DFT 分析信号的频谱。
3.深刻理解利用DFT 分析信号频谱的原理，分析实现过程中出现的现象及解决方法。
二、实验设备与环境 计算机、MATLAB 软件环境 三、实验基础理论1.DFT 与DTFT 的关系有限长序列 的离散时间傅里叶变换 在频率区间 的N 个等间隔分布的点 上的N 个取样值可以由下式表示：212/0()|()()01N jkn j Nk N k X e x n eX k k N πωωπ--====≤≤-∑由上式可知，序列 的N 点DFT ,实际上就是 序列的DTFT 在N 个等间隔频率点 上样本 。
2.利用DFT 求DTFT方法1：由恢复出的方法如下：由图2.1所示流程可知：101()()()N j j nkn j nN n n k X e x n eX k W e N ωωω∞∞----=-∞=-∞=⎡⎤==⎢⎥⎣⎦∑∑∑ 由上式可以得到：IDFTDTFT( )12()()()Nj k kX e X k Nωπφω==-∑ 其中为内插函数12sin(/2)()sin(/2)N j N x eN ωωφω--= 方法2：实际在MATLAB 计算中，上述插值运算不见得是最好的办法。
由于DFT 是DTFT 的取样值，其相邻两个频率样本点的间距为2π/N ，所以如果我们增加数据的长度N ，使得到的DFT 谱线就更加精细，其包络就越接近DTFT 的结果，这样就可以利用DFT 计算DTFT 。
如果没有更多的数据，可以通过补零来增加数据长度。
3.利用DFT 分析连续信号的频谱采用计算机分析连续时间信号的频谱，第一步就是把连续信号离散化，这里需要进行两个操作：一是采样，二是截断。
对于连续时间非周期信号,按采样间隔T 进行采样，阶段长度M ，那么：1()()()M j tj nT a a a n X j x t edt T x nT e ∞--Ω-Ω=-∞Ω==∑⎰对进行N 点频域采样，得到2120()|()()M jkn Na a M kn NTX j T x nT eTX k ππ--Ω==Ω==∑因此，可以将利用DFT 分析连续非周期信号频谱的步骤归纳如下： （1）确定时域采样间隔T ，得到离散序列（2）确定截取长度M ，得到M 点离散序列,这里为窗函数。
数字信号处理实验文档
实验一：离散信号的MATLAB实现一、实验目的：1、掌握离散时间信号的一般表示方法。
2、熟悉连续信号经理想采样后的频谱变化关系，加深对时域采样定理的理解。
3、掌握离散信号序列的操作。
二、实验内容：M1-1 已知g1(t)=cos(6*pi*t), g1(t)=co 14*pi*t), g1(t)=cos(26*pi*t),以抽样频率fsam=10Hz对上述三个信号进行抽样。
在同一张图上画出g1(t),g2(t)和g3(t)及其抽样点，对所得结果进行讨论。
解：代码如下：100:100)*1/100;g1t=cos(6*pi*t);g2t=cos(14*pi*t);g3t=cos(26*pi*t);subplot(3,1,1);plot(g1t);subplot(3,1,2);plot(g2t);subplot(3,1,3);plot(g3t);绘出的图形如图1_1所示:图1_1采样频率为fsam=10Hz，采样时间为0.1s，而f1=3Hz，f2=7Hz，f3=13Hz，使得三个信号的采样图形相似，这样不能很好还原原来的信号图像。
所以对信号的采样频率应足够大，应满足fsam>=2fm.M1-2利用MATLAB的filter函数，求出下列系统的单位脉冲响应，并判断系统是否稳定。
讨论题所获得的结果。
代码1：k=1:300;x=zeros(1,300);x(1)=1;b1=[1];a1=[1,-1.845,0.850586];y1=filter(b1,a1,x);subplot(1,2,1);plot(k,y1);xlabel('k');ylabel('幅度y1');b2=[1];a2=[1,-1.85,0.85];y2=filter(b2,a2,x);subplot(1,2,2);plot(k,y2);xlabel('k');ylabel('幅度y2');图1_2_1代码2：x=zeros(1,500);x(1)=1;b1=[1];a1=[1,-1.845,0.850586];y1=filter(b1,a1,x);plot(k,y1);b2=[1];a2=[1,-1.85,0.85];y2=filter(b2,a2,x);plot(k,y1,k,y2,':');xlabel('k');ylabel('幅度');legend('y1''y2');图1_2_2结论：H1（z）的两个极点都在单位圆内，所以系统稳定，从图中可以看出响应曲线升高后有回落，系统最终趋向于0；H2（z）的一个极点在单位圆内，另一个在单位圆上，所以系统最终临界稳定，从图中可以看出响应曲线上升后没有回落，系统最终趋向于6.7左右。
数字信号处理实验一 实验报告
数字信号处理实验一1.完成本文档内容的自学阅读和其中各例题后子问题；Q1.1运行程序P1.1，以产生单位样本序列u[n]并显示它。
答: clf;n=-10:20;u=[zeros(1,10) 1 zeros(1,20)];stem(n,u);xlabel('时间序号n');ylabel('振幅');title('单位样本序列');axis([-10 20 0 1.2])Q1.2命令clf，axis，title，xlabel和ylabel的作用是什么？答：clf清除图对象，axis 控制轴刻度和风格的高层指令，title 设置图名，xlabel和ylabel设置横纵坐标轴名称。
Q1.3修改程序P1.1以产生带有延时11个单位样本的延迟单位样本序列ud[n]。
运行修改的程序并显示产生的序列。
答：clf;n=0:30;ud=[zeros(1,11) 1 zeros(1,19)];stem(n,ud);xlabel('时间序号n');ylabel('振幅');title('单位样本序列');axis([0 30 0 1.2])Q1.4修改程序P1.1以产生单位步长序列s[n].运行修改后程序并显示产生的序列。
答：clf;n = 0:30;u = [1.*n];stem(n,u);title('Unit Sample Sequence');axis([0 30 0 30])Q1.5修改程序P1.1，以产生带有超前7个样本的延时单位阶跃序列sd[n]。
运行修改后的程序并显示产生的序列。
答：clf;n = -15:30;s=[zeros(1,8) ones(1,38)];stem(n,s);xlabel('Time index n');ylabel('Amplitude'); title('Unit Sample Sequence');axis([-15 30 0 1.2]);Q1.6 运行程序P1.2，以产生复数值的指数序列。
数字信号处理实验报告
实验一：频谱分析与采样定理 subplot(3,1,1),stem(t,x2);title('指数信号'); subplot(3,1,2),stem(f1,y2);title('指数信号频谱'); subplot(3,1,3),plot(f2,y21);title('指数信号频谱'); %%%%%%%%%%%%%%%%%%%%%%%%% x3=x1.*x2; y3=T*abs(fft(x3)); y31=fftshift(y3); figure(3), subplot(3,1,1),stem(t,x3);title('两信号相乘'); subplot(3,1,2),stem(f1,y3);title('两信号相乘频谱'); subplot(3,1,3),plot(f2,y31);title('两信号相乘频谱'); 实验结果： T=1/10000，������������ =10000，L=0.10
1/ 5
实验二：卷积定理 Y2=fft(y2); Z2=X2.*Y2; z2=ifft(Z2); figure(3), subplot(321),stem(x2);title('x2'); subplot(322),stem(real(X2));title('X2'); subplot(323),stem(y2);title('y2'); subplot(324),stem(real(Y2));title('Y2'); subplot(325),stem(z2);title('z2'); subplot(326),stem(real(Z2));title('Z2'); N=6; x3=[x zeros(1,N-length(x))]; y3=[y zeros(1,N-length(y))]; X3=fft(x3); Y3=fft(y3); Z3=X3.*Y3; z3=ifft(Z3); figure(4), subplot(321),stem(x3);title('x3'); subplot(322),stem(real(X3));title('X3'); subplot(323),stem(y3);title('y3'); subplot(324),stem(real(Y3));title('Y3'); subplot(325),stem(z3);title('z3'); subplot(326),stem(real(Z3));title('Z3'); N=8; x4=[x zeros(1,N-length(x))]; y4=[y zeros(1,N-length(y))]; X4=fft(x4); Y4=fft(y4); Z4=X4.*Y4; z4=ifft(Z4); figure(5), subplot(321),stem(x4);title('x4'); subplot(322),stem(real(X4));title('X4'); subplot(323),stem(y4);title('y4'); subplot(324),stem(real(Y4));title('Y4'); subplot(325),stem(z4);title('z4'); subplot(326),stem(real(Z4));title('Z4'); %N=6 时
数字信号处理实验报告
实验报告课程名称：数字信号处理授课班级：学号：姓名：指导老师：实验一离散时间信号及系统的时域分析实验类别：基础性实验1实验目的：(1）了解MA TLAB 程序设计语言的基本特点，熟悉MA TLAB软件运行环境。
(2）熟悉MA TLAB中产生信号和绘制信号的基本命令，学会用MA TLAB在时域中产生一些基本的离散时间信号，并对这些信号进行一些基本的运算。
(3）通过MA TLAB仿真一些简单的离散时间系统，并研究它们的时域特性。
(4）通过MA TLAB进行卷积运算，利用卷积方法观察分析系统的时域特性。
2. 实验报告要求●简述实验原理及目的。
●结合实验中所得给定典型序列幅频特性曲线，与理论结果比较，并分析说明误差产生的原因以及用FFT作谱分析时有关参数的选择方法。
●记录调试运行情况及所遇问题的解决方法。
3.实验内容：思考题：9.2.1 运行程序P9.2.1，哪个参数控制该序列的增长或衰减：哪个参数控制该序列的振幅？若需产生实指数序列，应对程序作何修改？9.2.2运行程序P9.2.1，该序列的频率是多少？怎样改变它？哪个参数控制该序列的相位？哪个参数可以控制该序列的振幅？该序列的周期是多少？9.2.3 运行程序P9.2.3，对加权输入得到的y(n)与在相同权系数下输出y1(n)和y2(n)相加得到的yt(n)进行比较，这两个序列是否相等？该系统是线性系统吗？9.2.4 假定另一个系统为y(n)=x(n)x(n-1)修改程序，计算这个系统的输出序列y1(n),y2(n)和y(n)。
比较有y(n)和yt(n)。
这两个序列是否相等？该系统是线性系统吗？(提高部分)9.2.5运行程序P9.2.4，并比较输出序列y(n)和yd(n-10)。
这两个序列之间有什么关系？该系统是时不变系统吗？9.2.6 考虑另一个系统：修改程序,以仿真上面的系统并确定该系统是否为时不变系统。
（选做）n = 0:40; D = 10;a = 3.0;b = -2;x = a*cos(2*pi*0.1*n) + b*cos(2*pi*0.4*n);[x1,n1]=sigmult(n,n,x,n)[x2,n2]=sigshift(x,n,1)[y,ny1]= sigadd(x1,n1,x2,n2)[y1,ny11]= sigshift(y,ny1,D)[sx,sn]= sigshift(x,n,D)[sx1,sn1]=sigmult(n,n, sx,sn)[sx2,sn2]=sigshift(sx,sn,1)[y2,ny2]= sigadd(sx1,sn1,sx2,sn2)D= sigadd(y1,ny11,y2,ny22)六、实验心得体会：实验时间批阅老师实验成绩实验二 FFT 实现数字滤波实验类别：提高性实验 1.实验目的(1) 通过这一实验，加深理解FFT 在实现数字滤波（或快速卷积）中的重要作用，更好的利用FFT 进行数字信号处理。
数字信号处理实验报告(实验1-4)
实验一 MATLAB 仿真软件的基本操作命令和使用方法实验容1、帮助命令使用 help 命令，查找 sqrt （开方）函数的使用方法；2、MATLAB 命令窗口（1）在MATLAB 命令窗口直接输入命令行计算31)5.0sin(21+=πy 的值；（2）求多项式 p(x) = x3 + 2x+ 4的根；3、矩阵运算（1）矩阵的乘法已知 A=[1 2;3 4]， B=[5 5;7 8]，求 A^2*B（2）矩阵的行列式已知 A=[1 2 3;4 5 6;7 8 9]，求A（3）矩阵的转置及共轭转置已知A=[1 2 3;4 5 6;7 8 9]，求A'已知 B=[5+i,2-i,1;6*i,4,9-i], 求 B.' , B'（4）特征值、特征向量、特征多项式已知 A=[1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4] ，求矩阵 A的特征值、特征向量、特征多项式；（5）使用冒号选出指定元素已知： A=[1 2 3;4 5 6;7 8 9]；求 A 中第 3 列前 2 个元素；A 中所有列第 2，3 行的元素；4、Matlab 基本编程方法（1）编写命令文件：计算 1+2+…+n<2000 时的最大 n 值；（2）编写函数文件：分别用 for 和 while 循环结构编写程序，求 2 的 0 到 15 次幂的和。
5、MATLAB基本绘图命令（1）绘制余弦曲线 y=cos(t)，t∈[0，2π]（2）在同一坐标系中绘制余弦曲线 y=cos(t-0.25)和正弦曲线 y=sin(t-0.5)， t∈[0，2π]（3）绘制[0，4π]区间上的 x1=10sint 曲线，并要求：（a）线形为点划线、颜色为红色、数据点标记为加号；（b）坐标轴控制：显示围、刻度线、比例、网络线（c）标注控制：坐标轴名称、标题、相应文本；>> clear;t=0:pi/10:4*pi;y=10*sin(t);plot(t,y);plot(t,y,'-+r');grid>> xlabel('X'),ylabel('Y');>> title('Plot:y=10*sin(t)');>> text(14,10,'完整图形');实验二常见离散信号的MATLAB产生和图形显示实验容与步骤１. 写出延迟了np个单位的单位脉冲函数impseq，单位阶跃函数stepseq, n=ns:nf function [x,n]=impseq[np,ns,nf];function [x,n]=stepseq[np,ns,nf];2. 产生一个单位样本序列x1(n)，起点为ns= -10, 终点为nf=20, 在n0=0时有一单位脉冲并显示它。
数字信号处理实验1
实验1 离散时间信号的时域分析一、实验目的（1）了解MATLAB语言的主要特点及作用；（2）熟悉MATLAB主界面，初步掌握MATLAB命令窗和编辑窗的操作方法；（3）了解常用时域离散信号及其特点；（4）掌握MATLAB产生常用时域离散信号的方法；（5）掌握MATLAB中时域离散信号的基本运算方法；（6）学习简单的数组赋值、数组运算、绘图的程序编写。
二、知识点提示本章节的主要知识点是利用MATLAB产生数字信号处理的几种常用典型序列及数字序列的基本运算；重点是单位脉冲、单位阶跃、正（余）弦信号的产生；难点是MATLAB关系运算符“==、>=”的使用。
三、实验原理1．时域离散信号的概念在MATLAB中，时域的离散信号可以通过编写程序直接生成，也可以通过对连续信号等间隔抽样获得。
离散序列的时域运算主要为信号的相加和相乘，信号的时域变换包括移位、反转及尺度变换。
2．用MATLAB生成离散信号需注意的问题（1）有关数组与下标MATLAB中处理的数组，其下标默认从1开始递增，例如x=[9 8 7]，表示x(1)=9；x(2)=8；x(3)=7。
要表示一个下标不从1开始的数组，一般需要采用两个矢量，如：n=[-3:1:2]，x=[9 8 7 6 5 4]，则有x(-3)=9；x(-2)=8；x(2)=4。
（2）信号的图形绘制从本质上来讲，MATLAB及其任何计算机语言处理的信号都是离散信号。
当我们把信号的样点值取的足够密，作图时采用特殊的命令，就可以把信号近似看成连续信号。
在MATLAB中，离散信号与连续信号有时在程序编写上是一致的，只是在作图时选择不同的绘图函数而已。
连续信号作图用plot函数，绘制线形图；离散信号作图使用stem函数，绘制脉冲图。
3．常用时域离散信号常用时域离散信号有单位脉冲序列、单位阶跃序列、实指数序列、复指数序列、正弦序列、矩形序列以及随机序列等。
下面以单位脉冲序列的产生为例做详细说明。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
实验一离散时间信号与系统的分析实验二：离散系统的时域分析一．实验原理：二：实验程序：1》单位冲击响应1. 用filter函数：>> a1=[1,0.75,0.125];b1=[1,-1];n=0:20;x1=[1 zeros(1,20)];y1filter=filter(b1,a1,x1);stem(n,y1filter);title('y1filter');xlabel('x');ylabel('y');>>y1filterx2. 用conv函数>> a1=[1,0.75,0.125];b1=[1,-1];x1=[1 zeros(1,10)];[h]=impz(b1,a1,10);y1conv=conv(h,x1);n=0:19;stem(n,y1conv,'filled')>>3. 用impz函数>> a1=[1,0.75,0.125]; b1=[1,-1];impz(b1,a1,21);n (samples)A m p l i t u d eImpulse Response1》单位阶跃响应1. 用filter 函数 a1=[1,0.75,0.125]; b1=[1,-1]; n=0:20;x2=ones(1,21);y1filter=filter(b1,a1,x2); stem(n,y1filter); title('y1filter_step'); xlabel('x'); ylabel('y');2. 用conv 函数>> a1=[1,0.75,0.125]; b1=[1,-1]; x2=ones(1,21);[h,m]=impz(b1,a1,20); y1=conv(h,x2);y1conv=y1(1:21); n1=0:20; stem(n1,y1conv,'filled'); title('y1conv'); xlabel('n'); ylabel('y1[n]'); >>ny 1[n ]3. 用impz 函数 >> a=[1,0.75,0.125]; b=[1 -1]; impz(b,a) >>n (samples)A m p l i t u d e2》单位冲击响应1. 用filter函数：>> a1=[1];b1=[0 0.25 0.25 0.25 0.25];n=0:20;x1=[1 zeros(1,20)];y1filter=filter(b1,a1,x1);stem(n,y1filter);title('y1filter');xlabel('x');ylabel('y');>>yx2. 用conv函数>> a1=[1];b1=[0 0.25 0.25 0.25 0.25];x1=[1 zeros(1,10)];[h]=impz(b1,a1,10);y1conv=conv(h,x1);n=0:19;stem(n,y1conv,'filled')>>3. 用impz 函数 >> a1=[1];b1=[0 0.25 0.25 0.25 0.25]; impz(b1,a1,21);n (samples)A m p l i t u d e2>>阶跃响应：1.用filter 得到的阶跃响应： a2=[1];b2=[0.25 0.25 0.25 0.25]; n=0:20;x2=ones(1,21); y2=filter(b2,a2,x2); stem(n,y2);title('用filter 得到的系统阶跃响应’); xlabel('x2'); ylabel('y2');2.用conv得到的阶跃响应：a2=[1];b2=[0 0.25 0.25 0.25 0.25];x2=ones(1,21);[h2]=impz(b2,a2,20);y2=conv(h2,x2);y3=y2(1:21);stem(y3);girdtitle('用conv得到的系统阶跃响应'); xlabel('x3');ylabel('y3');3.用impz函数的阶跃响应：a2=[1];b2=[ 0.25 0.25 0.25 0.25];y3=impz(b2,a2,25);stem(y3);girdtitle('用impz得到的系统阶跃响应'); xlabel('x3');ylabel('y3');实验三一．实验原理：二．实验程序及其结果：两个线性卷积及其结果：>> x1=[1 2 3 4 5];x2=[1 2 1 2];L=length(x1)+length(x2)-1;XE=fft(x1,L);HE=fft(x2,L);y=ifft(XE.*HE);n=0:L-1;stem(n,y);xlable('n');ylable('Amplitude');title('xian xin juan ji');5、6、7、8点圆周卷积程序及其结果：function yc=circonv2(x1,x2,N)if length(x1)>Nerror('N must not be less than length of x1' ) endif length(x2)>Nerror('N must not be less than length of x2' ) endx1=[x1,zeros(1,N-length(x1))];x2=[x2,zeros(1,N-length(x2))];n=[0:1:N-1];x2=x2(mod(-n,N)+1);H=zeros(N,N);for n=1:1:NH (n,:)=cirshiftd(x2,n-1,N);endyc=x1*H;function y=cirshiftd(x,m,N)if length(x)>Nerror('length of x must be less than N '); endx=[x,zeros(1,N -length(x))];n=[0:1:N-1];y=x(mod (n-m,N)+1);﹪m文件x1=[1 2 3 4 5];x2=[1 2 1 2 ];ycn1=circonv2(x1,x2,5);ny1=[0:1:length(ycn1)-1];subplot(2,2,1);stem(ny1,ycn1);title('5点长的循环卷积');ycn2=circonv2(x1,x2,6);ny2=[0:1:length(ycn2)-1];subplot(2,2,2);stem(ny2,ycn2);title('6点长的循环卷积');ycn3=circonv2(x1,x2,7);ny3=[0:1:length(ycn3)-1];subplot(2,2,3);stem(ny3,ycn3);title('7点长的循环卷积');ycn4=circonv2(x1,x2,8);ny4=[0:1:length(ycn4)-1];subplot(2,2,4);stem(ny4,ycn4);title('8点长的循环卷积');线性卷积通用程序：x1=input(‘输入序列x1=’)x2=input(‘输入序列x2=’)L=length(x1)+length(x2)-1; XE=fft(x1,L);HE=fft(x2,L);y=ifft(XE.*HE);disp(‘线性卷积‘)；disp(y);n=0:L-1;stem(n,y)xlable(n,y);ylable(‘Amplitude’);title(‘线性卷积‘)圆周卷积通用程序：X1=input(‘输入序列x1=’)X2=input(‘输入序列x2=’)L=length（x1）+length(x2)-1) y=circonv(x1,x2);stem(y);xlable(‘n’);ylable(‘y’);实验五一．实验原理：1. 零、极点图，幅度、频率响应和相位响应num=[0,0528 0.0797 0.1295 0.1295 0.797 0.0582]; den=[1 -1.81072.4947 -1.8801 0.9537 -0.2336]; [z,p,k]=tf2zp(num,den);m=abs(p);disp('零点');disp(z);disp('极点');disp(z);disp('增益系数');disp(k);sos=zp2sos(z,p,k);disp('二阶节点');disp(real(sos));subplot(2,1,1);zplane(num,den);k=256;w=0:pi/k:pi);h=frepz(num,den,w);subplot(2,2,3);plot(w/pi,abs(h));gridtitle('幅度谱')xlable('\omega/\pi');ylale('幅值');subplot(2,2,4);plot(w/pi,angle(h));gridtitle('相位谱');xlable('\omega/\pi');ylable('弧度');用tf2zp实现算法和分解为二阶：。