离散时间系统的时域特性分析(数字信号处理实验1答案)
————————————————————————————————作者：————————————————————————————————日期：
2
二程序框图clf;
n=0:299;
x=cos(20*pi*n/256)+cos(200*pi*n/256);
num=[0.5 0.27 0.77];
den=[1 0 0];
y=filter(num,den,x);
subplot(2,1,1);
plot(n,x);
xlabel('时间信号n');
ylabel('信号幅度');
title('输入信号');
subplot(2,1,2);
plot(y);
xlabel('时间信号n');
ylabel('信号幅度');
title('输出信号');
n=0:299;
x=cos(20*pi*n/256)+cos(200*pi*n/256); num=[0.45 0.5 0.45];
den=[1 -0.53 0.46];
y=filter(num,den,x);
subplot(2,1,1);
plot(n,x);
xlabel('时间信号n');
ylabel('信号幅度');
title('输入信号');
subplot(2,1,2);
plot(y);
xlabel('时间信号n');
ylabel('信号幅度');
title('输出信号');
num=[0.5 0.27 0.77]; den=[1 0 0];
y=impz(num,den,n); stem(y);
xlabel('时间信号n'); ylabel('信号幅度'); title('冲激响应'); grid;
num=[0.45 0.5 0.45]; den=[1 -0.53 0.46]; y=impz(num,den,n); stem(y);
xlabel('时间信号n'); ylabel('信号幅度'); title('冲激响应'); grid;
a=2;
b=3;
x1=cos(20*pi*n/256); x2=cos(200*pi*n/256); x=a*x1+b*x2;
num=[0.45 0.5 0.45]; den=[1 -0.53 0.46];
y1=filter(num,den,x1); y2=filter(num,den,x2); y=filter(num,den,x);
yt=a*y1+b*y2; subplot(2,1,1);
stem(n,y);
ylabel('信号幅度'); subplot(2,1,2);
stem(n,yt);
ylabel('信号幅度');
D=10;
a=1;
b=1;
x=a*cos(20*pi*n/256)+b*cos(200*pi*n/256); xd=[zeros(1,D) x];
num=[0.45 0.5 0.45];
den=[1 -0.53 0.46];
ic=[0 0];
y=filter(num,den,x,ic);
yd=filter(num,den,xd,ic);
N=length(y);
d=y-yd(1+D:N+D);
subplot(3,1,1);
stem(n,y);
ylabel('信号幅度');
title('输出y[n]');
grid;
subplot(3,1,2);
stem(n,yd(1:length(y)));
ylabel('信号幅度');
subplot(3,1,3);
stem(n,d);
xlable('时间序号n');
ylable('信号幅度');
title('差值信号');
grid;
目录 第1章设计任务及要求 (1) 1.1课程设计内容 (1) 1.2课程设计要求 (1) 第2章设计原理 (2) 2.1离散信号与系统的时域分析设计 (2) 2.1.1描写系统特性的方法介绍 (2) 2.1.2系统的时域特性 (2) 第3章设计实现 (3) 3.1实验内容与方法 (3) 3.1.1实验内容 (3) 第4章设计结果及分析 (3) 4.1程序设计结果及分析 (4) 总结 (7) 参考文献： (7) 附录： (8)
第1章 设计任务及要求 1.1课程设计内容 编制Matlab 程序，完成以下功能，产生系统输入信号；根据系统差分方程求解单位脉冲响应序列；根据输入信号求解输出响应；用实验方法检查系统是否稳定；绘制相关信号的波形。具体要求如下： (1) 给定一个低通滤波器的差分方程为 ()0.05()0.05(1)0.9(1)y n x n x n y n =+-+- 输入信号分别为182()=()()()x n R n x n u n =， ① 分别求出系统响应，并画出其波形。 ② 求出系统的单位脉冲响应，画出其波形。 (2) 给定系统的单位脉冲响应为1102()=()()() 2.5(1) 2.5(2)(3)h n R n h n n n n n δδδδ=+-+-+-，用线性卷积法求18()=()x n R n 分别对系统h1(n)和h2(n)的输出响应，并画出波形。 (3) 给定一谐振器的差分方程为() 1.8237(1)-0.9801(2)()(2)o o y n y n y n b x n b x n =--++-令b0=1/100.49,谐振器的谐振频率为0.4rad 。 1) 用实验方法检查系统是否稳定。输入信号为u(n)时,画出系统输出波形。 2) 给定输入信号为()=sin(0.014)sin(0.4)x n n n +求出系统的输出响应，并画出其波形。 1.2课程设计要求 1. 要求独立完成设计任务。 2. 课程设计说明书封面格式要求见《天津城市建设学院课程设计教学工作规范》附表1 3. 课程设计的说明书要求简洁、通顺，计算正确，图纸表达内容完整、清楚、规范。 4. 简述离散系统时域分析方法和通过实验判断系统稳定性的方法；完成以上设计实验并对结果进行分析和解释；打印程序清单和要求画出的信号波形；写出本次课程设计的收获和体会。 5. 课设说明书要求： 1) 说明题目的设计原理和思路、采用方法及设计流程。 2) 详细介绍运用的理论知识和主要的Matlab 程序。 3) 绘制结果图形并对仿真结果进行详细的分析。
信号与系统实验报告 实验名：离散时间信号与系统的频域分析 实验六离散时间系统的时域分析 一、实验目的 1、掌握离散时间信号与系统的频域分析方法，从频域的角度对信号与系统的特性进行分析。 2、掌握离散时间信号傅里叶变换与傅里叶逆变换的实现方法。 3、掌握离散时间傅里叶变换的特点及应用 4、掌握离散时间傅里叶变换的数值计算方法及绘制信号频谱的方法 二、预习内容 1、离散时间信号的傅里叶变换与逆变换。 2、离散时间信号频谱的物理含义。 3、离散时间系统的频率特性。 4、离散时间系统的频域分析方法。 三、实验原理 1. 离散时间系统的频率特性
2. 离散时间信号傅里叶变换的数值计算方法 3.涉及到的Matlab 函数
四、实验内容 1、离散时间系统的时域分析 1 离散时间傅里叶变换 (1)下面参考程序是如下序列在范围?4π≤ω≤ 4π的离散时间傅里叶变换 %计算离散时间傅里叶变换的频率样本 clear all; w=-4*pi:8*pi/511:4*pi; num=[2 1]; den=[1 -0.6]; h=freqz(num,den,w); subplot(2,1,1)
plot(w/pi,real(h)); grid; title(‘实部’) xlabel(‘omega/\pi’); yl abel(‘振幅’); subplot(2,1,2) plot(w/pi, imag(h)); grid; title(‘虚部’) xlabel(‘omega/\pi’); ylabel(‘振幅’); figure; subplot(2,1,1) plot(w/pi, abs(h)); grid; title(‘幅度谱’) xlabel(‘omega/\pi’); ylabel(‘振幅’); subplot(2,1,2) plot(w/pi, angle (h)); grid; title(‘相位谱’) x label(‘omega/\pi’); ylabel(‘以弧度为单位的相位’);
实验2 离散系统的时域分析 一、实验目的 1、熟悉并掌握离散系统的差分方程表示法； 2、加深对冲激响应和卷积分析方法的理解。 二、实验原理 在时域中，离散时间系统对输入信号或者延迟信号进行运算处理，生成具有所需特性的输出信号，具体框图如下: 其输入、输出关系可用以下差分方程描述： 输入信号分解为冲激信号, 记系统单位冲激响应，则系统响应为如下的卷积计算式： 当时，h[n]是有限长度的（），称系统为FIR系统；反之，称系统为IIR系统。 三、实验内容
1、用MATLAB 求系统响应 1） 卷积的实现 线性移不变系统可由它的单位脉冲响应来表征。若已知了单位脉冲响应和系统激励就 可通过卷积运算来求取系统响应，即)(*)()(n h n x n y 程序： x=input(‘Type in the input sequence=’); %输入x h=input(‘Type in the impulse response sequence=’); %输入h y=conv(x,h); % 对x ，h 进行卷积 N=length(y)-1; %求出N 的值 n=0:1:N; %n 从0开始，间隔为1的取值取到N 为止 disp(‘output sequence=’); disp(y); %输出y stem(n,y); %画出n 为横轴，y 为纵轴的离散图 xlabel(‘Time index n ’); ylable(‘Amplitude ’); % 规定x 轴y 轴的标签 输入为： x=[-2 0 1 -1 3] h=[1 2 0 -1] 图形: 2） 单位脉冲响应的求取 线性时不变因果系统可用MA TLAB 的函数filter 来仿真 y=filter(b,a,x); 其中，x 和y 是长度相等的两个矢量。矢量x 表示激励，矢量a ，b 表示系统函数形式 滤波器的分子和分母系数，得到的响应为矢量y 。例如计算以下系统的单位脉冲响应 y(n)+0.7y(n-1)-0.45y(y-2)-0.6y(y-3)=0.8x(n)-0.44x(n-1)+0.36x(n-2)+0.02x(n-3) 程序: N=input(‘Desired impuse response length=’); b=input(‘Type in the vector b=’); a=input(‘Type in the vector a=’); x=[1 zeros(1,N-1)]; y=filter(b,a,x);
第七章离散时间系统的时域分析 §7-1 概述 一、离散时间信号与离散时间系统 离散时间信号：只在某些离散的时间点上有值的 信号。 离散时间系统：处理离散时间信号的系统。 混合时间系统：既处理离散时间信号，又处理连 续时间信号的系统。 二、连续信号与离散信号 连续信号可以转换成离散信号，从而可以用离散时间系统（或数字信号处理系统）进行处理： 三、离散信号的表示方法：
1、 时间函数：f(k)<——f(kT)，其中k 为序号，相当于时间。 例如：)1.0sin()(k k f = 2、 （有序）数列：将离散信号的数值按顺序排列起来。例如： f(k)={1,0.5,0.25,0.125,……,} 时间函数可以表达任意长（可能是无限长）的离散信号，可以表达单边或双边信号，但是在很多情况下难于得到；数列的方法表示比较简单，直观，但是只能表示有始、有限长度的信号。 四、典型的离散时间信号 1、 单位样值函数：? ??==其它001)(k k δ 下图表示了)(n k ?δ的波形。
这个函数与连续时间信号中的冲激函数 )(t δ相似，也有着与其相似的性质。例如： )()0()()(k f k k f δδ=， )()()()(000k k k f k k k f ?=?δδ。 2、 单位阶跃函数：? ??≥=其它001)(k k ε 这个函数与连续时间信号中的阶跃函数)(t ε相似。用它可以产生（或表示）单边信号（这里称为单边序列）。 3、 单边指数序列：)(k a k ε
比较：单边连续指数信号：)()()(t e t e t a at εε=，其 底一定大于零，不会出现负数。 (a) 0.9a = (d) 0.9a =? (b) 1a = (e) 1a =? (c) 1.1a = (f) 1.1a =?
. ... 实验二：离散LSI系统的时域分析 一、实验内容 1.知描述某离散LSI系统的差分方程为2y(n)-3y(n-1)+y(n-2)=x(n-1)，分别用impz 和dstep函数、filtic和filter函数两种方法求解系统的单位序列响应和单位阶跃响应。 用impz和dstep函数求解系统的单位序列响应和单位阶跃响应如下 a=[1,-3/2,1/2]; b=[0,1/2,0]; N=32; n=0:N-1; hn=impz(b,a,n); gn=dstep(b,a,n); subplot(1,2,1);stem(n,hn,'k'); title('系统的单位序列响应'); ylabel('h(n)');xlabel('n'); axis([0,N,1.1*min(hn),1.1*max(hn)]); subplot(1,2,2);stem(n,gn,'k'); title('系统的单位阶跃响应'); ylabel('g(n)');xlabel('n'); axis([0,N,1.1*min(gn),1.1*max(gn)]); 课程名称数字信号 实验成绩 指导教师实验报告.
... 010203000.10.20.0.0.0.0.0.0.1系统的单位序列响应h(n) n01020300112230系统的单位阶跃响应g(n)n 用函数filtic和filter求解离散系统的单位序列响应和单位阶跃
解：x01=0;y01=0; a=[1,-3/2,1/2]; b=[1/2,0,0]; N=32;n=0:N-1; xi=filtic(b,a,0); x1=[n==0]; hn=filter(b,a,x1,xi); x2=[n>=0]; gn=filter(b,a,x2,xi); subplot(1,2,1);stem(n,hn,'k'); title('系统的单位序列响应'); ylabel('h(n)');xlabel('n'); axis([0,N,1.1*min(hn),1.1*max(hn)]); . ... subplot(1,2,2);stem(n,gn,'k'); title('系统的单位阶跃响应'); ylabel('g(n)');xlabel('n'); axis([0,N,1.1*min(gn),1.1*max(gn)]); 01020300.550.60.650.70.750.80.850.90.951
实验二离散时间系统的时域分析实验
数字信号处理——实验二 武汉工程大学电气信息学院通信工程 红烧大白兔 一、实验目的 1、在时域中仿真离散时间系统，进而理解离散 时间系统对输入信号或延时信号进行简单运算处理，生成具有所需特性的输出信号的方法。 2、仿真并理解线性与非线性、时变与时不变等 离散时间系统。 3、掌握线性时不变系统的冲激响应的计算并 用计算机仿真实现。 4、仿真并理解线性时不变系统的级联、验证线 性时不变系统的稳定特性。 二、实验设备 计算机，MATLAB语言环境 三、实验基础理论 1、系统的线性 线性性质表现为系统满足线性叠加原理：若某一输入是由N个信号的加权和组成的，则输出就是由系统对这N个信号中每一个的响应的相应加权和组成的。设x1（n）和 x2（n）
分别作为系统的输入序列，其输出分别用y1(n)和y2(n)表示，即 Y1(n)=T[x1(n)], y2(n)=T[x2(n)] 若满足T[a1x1(n)+a2x2(n)]=a1y1(n)+a2y2(n) 则该系统服从线性叠加原理，或者称为该系统为线性系统。 2、系统的时不变特性 若系统的变换关系不随时间变化而变化，或者说系统的输出随输入的移位而相应移位但形状不变，则称该系统为时不变系统。对于时不变系统， 若y(n)=T[x(n)] 则T[x(n-m)]=y(n-m) 3、系统的因果性 系统的因果性既系统的可实现性。如果系统n 时刻的输出取决于n时刻及n时刻以前的输入，而和以后的输入无关，则该系统是可实现的，是因果系统。系统具有因果性的充分必要条件是h(n)=0,n<0 4、系统的稳定性 稳定系统是指有界输入产生有界输出（BIBO)的系统。如果对于输入序列x(n),存在一个不变的正有限值M，对于所有n值满足|x(n)|
实验一 常见离散信号的MATLAB 产生和图形显示 一、 实验目的 加深对常见离散信号的理解 二、实验原理 1、单位抽样序列的产生 ,10,0{=≠=n n n ）（δ 在MATLAB 中可以用zeros()函数实现 x=[1,zeros(1,N-1)]; 或x=zeros(1，N)； x(1)=1; 2、单位阶跃序列的产生 0,10,0{u ≥<=n n n ）（ 在MATLAB 中可以用ones()函数实现 x=one(1,N); 3、正弦序列的产生 在MATLAB 中实现方法如下： N=0:N-1 X=A*sin(2*pi*f*n/fs+fai) 4、复正弦序列的产生jwn e A n x *)(= 在MATLAB 中实现方法如下： n) *w *exp(j *A 1 :0=-=x N n 5、实指数序列的产生n a A n x *)(= 在MATLAB 中实现方法如下： n a A x N n .^*1 :0=-= 三、实验内容及步骤 编制程序产生以下信号，并绘出其图形。 1）产生64点的单位抽样序列)(n δ
N=64 x=[1,zeros(1,N-1)] stem(x) 2）产生64点并移位20位的单位抽样序列)20(-n δ N=64 x=[0,zeros(1,N-1)] x(20)=1 stem(x) 3）任意序列)5(7.0)4(9.2)3(6.5)2(8.1)1(4.3)(0.8)(-+-+-+-+-+=n n n n n n n f δδδδδδ b=[1]; a=[8,3.4,1.8,5.6,2.9,0.7]; xh=[1,zeros(1,20)]; h=filter(b,a,xh) figure(1); n=0:20; stem(n,h,) legend('冲激')
第3章 离散时间LTI 系统的时域分析 3.1 实验目的 ● 学会运用MATLAB 求解离散时间系统的零状态响应； ● 学会运用MATLAB 求解离散时间系统的单位取样响应； ● 学会运用MATLAB 求解离散时间系统的卷积和。 3.2 实验原理及实例分析 3.2.1 离散时间系统的响应 离散时间LTI 系统可用线性常系数差分方程来描述，即 ∑∑==-=-M j j N i i j n x b i n y a 0 )()( （3-1） 其中，i a （0=i ，1，…，N ）和j b （0=j ，1，…，M ）为实常数。 MATLAB 中函数filter 可对式（13-1）的差分方程在指定时间范围内的输入序列所产生的响应进行求解。函数filter 的语句格式为 y=filter(b,a,x) 其中，x 为输入的离散序列；y 为输出的离散序列；y 的长度与x 的长度一样；b 与a 分别为差分方程右端与左端的系数向量。 【实例3-1】 已知某LTI 系统的差分方程为 )1(2)()2(2)1(4)(3-+=-+--n x n x n y n y n y 试用MATLAB 命令绘出当激励信号为)()2/1()(n u n x n =时，该系统的零状态响应。 解：MATLAB 源程序为 >>a=[3 -4 2]; >>b=[1 2]; >>n=0:30; >>x=(1/2).^n; >>y=filter(b,a,x); >>stem(n,y,'fill'),grid on >>xlabel('n'),title('系统响应y(n)') 程序运行结果如图3-1所示。
3.2.2 离散时间系统的单位取样响应 系统的单位取样响应定义为系统在)(n 激励下系统的零状态响应，用)(n h 表示。MATLAB 求解单位取样响应可利用函数filter ，并将激励设为前面所定义的impDT 函数。例如，求解实例13-1中系统的单位取样响应时，MATLAB 源程序为 >>a=[3 -4 2]; >>b=[1 2]; >>n=0:30; >>x=impDT(n); >>h=filter(b,a,x); >>stem(n,h,'fill'),grid on >>xlabel('n'),title('系统单位取样响应h(n)') 程序运行结果如图3-2所示。 图3-1 实例3-1系统的零状态响应 图3-2 实例13-1的系统单位取样响应
实验2 离散系统的时域分析
实验2 离散系统的时域分析 一、实验目的： 加深对离散系统的差分方程、单位抽样响应和卷积分析方法的理解。 二、实验原理： 离散系统 其输入、输出关系可用以下差分方程描述： ∑∑==-= -M m m N k n m n x b k n y a )()( 输入信号分解为冲激信号， ∑∞ -∞ =-= m m n m x n x )()()(δ 系统单位抽样序列h （n ）， 则系统响应为如下的卷积计算式： ∑∞ -∞ =-= *=m m n h m x n h n x n y )()()()()( 当00≠a N k a k ,...2,1,0==时，h(n)是有限长度的（n ：[0，M]），称系统为FIR 系统；反之，称系统为IIR 系统。 在MATLAB 中，可以用函数y=filter(b,a,x)实现差分方程的仿真，也可以用函数 y=conv(x,h)计算卷积，利用函数h=impz(b,a,n)计算画出脉冲响应。 （在实验报告中对这三种函数的使用方法及参数含义做出说明，这一部分手写） 三 、实验容 1.编制程序求解下列两个系统的单位抽样响应，并绘出其图形。 （1）)1()()2(125.0)1(75.0)(--=-+-+n x n x n y n y n y 程序： N=21; b=[1 -1]; a=[1 0.75 0.125]; x=[1 zeros(1,N-1)]; n=0:1:N-1; y=filter(b,a,x);
stem(n,y); xlabel('n'); ylabel('y(n)'); 图像： （2）)]4()3()2()1([25.0)(-+-+-+-=n x n x n x n x n y 程序： N=21; b=[0 0.25 0.25 0.25 0.25]; a=[1]; x=[1 zeros(1,N-1)]; n=0:1:N-1; y=filter(b,a,x); stem(n,y); xlabel('n'); ylabel('y(n)'); 图像：
实验报告 课程名称数字信号 实验成绩______________ 指导教师 实验二：离散LSI系统的时域分析 一、实验内容 1.知描述某离散LSI系统的差分方程为2y(n)-3y(n-1)+y(n-2)=x(n-1)，分别用impz和dstep函数、filtic和filter函数两种方法求解系统的单位序列响应和单位阶跃响应。 用impz和dstep函数求解系统的单位序列响应和单位阶跃响应如下 a=[1,-3/2,1/2]; b=[0,1/2,0]; N=32; n=0:N-1; hn=imp z(b,a ,n); gn=dste p(b,a, n); sub plot(1,2,1);stem( n,hn ,'k'); title('系统的单位序列响应'); ylabel('h( n)');xlabel(' n'); axis([0,N,1.1*mi n(h n),1.1*max(h n)]); sub plot(1,2,2);stem( n,g n,'k'); title('系统的单位阶跃响应'); ylabel('g( n)');xlabel(' n'); axis([0,N,1.1*mi n(g n),1.1*max(g n)]);
1 sub plot(1,2,2);stem( n,g n,'k'); 用函数filtic 和filter 求解离散系统的单位序列响应和单位阶跃 解：x01=0;y01=0; a=[1,-3/2,1/2]; b=[1/2,0,0]; N=32 ;n=0:N-1; xi=filtic(b,a,0); x1=[ n==0]; hn=filter(b,a,x1,xi); x2=[n >=0]; gn=filter(b,a,x2,xi); sub plot(1,2,1);stem( n,h n,'k'); title （'系统的单位序列响应'）; ylabel('h( n)');xlabel(' n'); axis([0,N,1.1*mi n(h n),1.1*max(h n)]); title （'系统的单位阶跃响应'）; ylabel('g( n)');xlabel(' n'); 0.9 0.8 0.7 n 0.6 0.5 0.4 0.3 0.2 0.1 系统的单位序列响应 系统的单位阶跃响应 30 25 20 15 10 10 20 n 30 0rLliiJ 0 10 20 n 30
第4章 时域离散系统的网络结构 习题 1. 已知系统用下面差分方程描述: )1(31)()2(81)1(43)(?+??n x n x n y n y n y ＋－＝ 试分别画出系统的直接型、 级联型和并联型结构。 式中x (n )和y (n )分别表示系统的输入和输出信号。 2． 设数字滤波器的差分方程为 )2(41)1(31)1()()(?+?+?+=n y n y n x n x n y 试画出系统的直接型结构。 3. 设系统的差分方程为 y (n )=(a +b )y (n －1)－aby (n －2)+x (n －2)+(a +b )x (n －1)+ab 式中, |a |<1， |b |<1, x (n )和y (n )分别表示系统的输入和输出信号, 试画出系统的直接型和级联型结构。 4. 设系统的系统函数为 )81.09.01)(5.01()414.11)(1(4)(211211??????++?+?+=z z z z z z z H 试画出各种可能的级联型结构， 并指出哪一种最好。 5． 题 5图中画出了四个系统， 试用各子系统的单位脉冲响应分别表示各总系统的单位脉冲响应， 并求其总系统函数。 题 5图
6． 题6图中画出了10种不同的流图， 试分别写出它们的系统函数及差分方程。 题6图 7. 假设滤波器的单位脉冲响应为 h (n )=a n u (n ) 0＜a ＜1 求出滤波器的系统函数， 并画出它的直接型结构。 8. 已知系统的单位脉冲响应为 h (n )=δ(n )+2δ(n －1)+0.3δ(n －2)+2.5δ(n －3)+0.5δ(n －5) 试写出系统的系统函数， 并画出它的直接型结构。 9. 已知FIR 滤波器的系统函数为 )9.01.29.01(101)(4321????++++=z z z z z H 试画出该滤波器的直接型结构和线性相位结构。 10． 已知FIR 滤波器的单位脉冲响应为:
数字信号处理 实验三 离散时间系统的时域分析 学院：信息与通信学院 专业：电子信息工程 学号：0900220418 姓名：梁芝铭
1.实验目的 （1）理解离散时间信号的系统及其特性。 （2）对简单的离散时间系统进行分析，研究其时域特性。 （3）利用MATLAB 对离散时间系统进行仿真，观察结果，理解其时域特性。 2.实验原理 离散时间系统，主要是用于处理离散时间信号的系统，即是将输入信号映射成的输出的某种运算，系统的框图如图所示： （1）线性系统 当该系统的输入信号为12()()ax n bx n +时，其中a,b 为任意常数，输出为 121212[()()][()][()]()()T a x n b x n a T x n b T x n a y n b y n +=+=+ （2）时不变系统 若()[()]y n T x n =，则[()]()T x n k y n k -=-。 3.实验内容及其步骤 （1）复习离散时间系统的主要性质，掌握其原理和意义。 （2）一个简单的非线性离散时间系统的仿真 在MATLAB 中输入：n = 0:200; x = cos(2*pi*0.05*n); x1 = [x 0 0]; x2 = [0 x 0]; x3 = [0 0 x]; y = x2.*x2-x1.*x3; y = y(2:202); subplot(2,1,1); plot(n, x); xlabel('Time index n'); ylabel('Amplitude');title('Input Signal');subplot(2,1,2);plot(n,y);xlabel('Time index n'); ylabel('Amplitude');title('Output signal'); 结果如下：