信号与系统实验教程只有答案.docx
《信号与系统实验教程只有答案.docx》由会员分享，可在线阅读，更多相关《信号与系统实验教程只有答案.docx（58页珍藏版）》请在冰点文库上搜索。
信号与系统实验教程只有答案
信号与系统
实验教程（只有答案）
（实验报告）
这么玩！
实验一信号与系统的时域分析
三、实验内容及步骤
实验前，必须首先阅读本实验原理，读懂所给出的全部范例程序。
实验开始时，先在计算机上运行这些范例程序，观察所得到的信号的波形图。
并结合范例程序应该完成的工作，进一步分析程序中各个语句的作用，从而真正理解这些程序。
实验前，一定要针对下面的实验项目做好相应的实验准备工作，包括事先编写好相应的实验程序等事项。
Q1-1：
修改程序Program1_1，将dt改为0.2，再执行该程序，保存图形，看看所得图形的效果如何？
dt=0.01时的信号波形dt=0.2时的信号波形
这两幅图形有什么区别，哪一幅图形看起来与实际信号波形更像？
答：
Q1-2：
修改程序Program1_1，并以Q1_2为文件名存盘，产生实指数信号x（t）=e-0.5t。
要求在图形中加上网格线，并使用函数axis（）控制图形的时间范围在0~2秒之间。
然后执行该程序，保存所的图形。
修改Program1_1后得到的程序Q1_2如下：
信号x（t）=e-0.5t的波形图
clear,%Clearallvariables
closeall,%Closeallfigurewindows
dt=0.2;%Specifythestepoftimevariable
t=-2:
dt:
2;%Specifytheintervaloftime
x=exp（-0.5*t）;%Generatethesignal
plot（t,x）
gridon;
axis（[0201]）
title（'Sinusoidalsignalx（t）'）
xlabel（'Timet（sec）'）
Q1-3：
修改程序Program1_1，并以Q1_3为文件名存盘，使之能够仿真从键盘上任意输入的一个连续时间信号，并利用该程序仿真信号x（t）=e-2t。
修改Program1_1后得到的程序Q1_3如下：
信号x（t）=e-2t的波形图
clear,
closeall,
dt=0.2;
t=-2:
dt:
2;
x=input（'Inputx（t）:
'）;
plot（t,x）
gridon;
axis（[02-11]）
title（'Sinusoidalsignalx（t）'）
xlabel（'Timet（sec）'）
Q1-4：
将实验原理中所给的单位冲激信号和单位阶跃信号的函数文件在MATLAB文件编辑器中编写好，并分别以文件名delta和u存入work文件夹中以便于使用。
抄写函数文件delta如下：
抄写函数文件u如下：
functiony=delta（t）%Unitstepfunction
dt=0.01;functiony=u（t）
y=（u（t）-u（t-dt））/dt;y=（t>=0）;%y=1fort>0,elsey=0
Q1-5：
修改程序Program1_4，并以Q1_5为文件名存盘，利用axis（）函数，将图形窗口的横坐标范围改为-2≤n≤5，纵坐标范围改为-1.5≤x≤1.5。
修改Program1_4后得到的程序Q1_5如下：
信号的波形图
clear,
closeall,
n=-5:
5;
x=[zeros（1,4）,0.1,1.1,-1.2,0,1.3,zeros（1,2）];
stem（n,x,'.'）
gridon,
axis（[-25-1.51.5]）;
title（'Adiscrete-timesequencex[n]'）
xlabel（'Timeindexn'）
Q1-6：
仿照前面的示例程序的编写方法，编写一个MATLAB程序，以Q1_6为文件名存盘，使之能够在同一个图形窗口中的两个子图中分别绘制信号x[n]=0.5|n|和x（t）=cos（2πt）[u（t）-u（t-3）]。
要求选择的时间窗能够表现出信号的主要部分（或特征）。
编写的程序Q1_6如下：
信号x[n]=0.5|n|的波形图和信号x（t）=cos（2πt）[u（t）-u（t-3）]的波形图
clear,closeall,
t=-1:
0.01:
4;
xt=cos（2*pi*t）.*（u（t）-u（t-3））;
n=-5:
5;
xn=（0.5）.^abs（n）;
subplot（211）
plot（t,xt）
gridon,
title（'Originalsignalx（t）'）
subplot（212）
stem（n,xn,'.'）
gridon,
title（'Originalsignalx（n）'）
xlabel（'Timet（sec）'）
Q1-7：
根据示例程序的编程方法，编写一个MATLAB程序，以Q1_7为文件名存盘，由给定信号x（t）=e-0.5tu（t）求信号y（t）=x（1.5t+3），并绘制出x（t）和y（t）的图形。
编写的程序Q1_7如下：
编写产生x（t）的函数文件x.m
functiony=x（t）
y=exp（-0.5*t）.*u（t）;
clear,closeall,
t=-3:
0.01:
4;
xt=x（t）;%Generatetheoriginalsignalx（t）
yt=x（1.5*t+3）;
subplot（211）
plot（t,xt）%Plotx（t）
gridon,
title（'Originalsignalx（t）'）
subplot（212）
plot（t,yt）%Plotx（t）
gridon,
title（'Originalsignaly（t）'）
xlabel（'Timet（sec）'）
信号x（t）的波形图信号y（t）=x（1.5t+3）的波形图
Q1-8：
给定一个离散时间信号x[n]=u[n]–u[n-8]，仿照示例程序Program1_5，编写程序Q1_8，产生x[n]的左移序列x1[n]=x[n+6]和右移序列x2[n]=x[n-6]，并在同一个图形窗口的三个子图中分别绘制这三个序列的图形。
编写的程序Q1_8如下：
编写产生x（t）的函数文件xx.m
functiony=xx（n）
y=u（n）-u（n-8）;
clear,closeall,
n=-10:
15;
x=xx（n）;%Generatetheoriginalsignalx（n）
x1=xx（n+6）;%Shiftx（t）totheleftby2secondtogetx1（n+6）
x2=xx（n-6）;%Shiftx（t）totherightby2secondtogetx2（n-6）
subplot（311）
stem（n,x,'.'）%Plotx（t）
gridon,
title（'Originalsignalx（n）'）
subplot（312）
stem（n,x1,'.'）%Plotx1（t）
gridon,
title（'Leftshiftedversionofx（n）'）
subplot（313）
stem（n,x2,'.'）%Plotx2（t）
gridon,
title（'Rightshiftedversionofx（n）'）
xlabel（'Timet（sec）'）
信号波形图
Q1-9：
编写程序Q1_9，使之能够接受以键盘方式输入的定义在不同时间段的两个不同连续时间信号并完成卷积运算，分别绘制这两个信号及其卷积的结果的图形，图形按照2⨯2分割成四个子图。
编写的程序Q1_9如下：
clear;closeall;
dt=0.01;
t0=input（'Inputfirstsignalt0:
'）;t1=input（'Inputfirstfirstsignalt1:
'）;
tx=t0:
dt:
t1;
x=input（'Inputfirstsignalvariable（tx）:
'）;
t2=input（'Inputsecondsignalt0:
'）;t3=input（'Inputsecondsignalt1:
'）;
th=t2:
dt:
t3;
h=input（'Inputsecondsignalvariable（th）:
'）
y=dt*conv（x,h）;%Computetheconvolutionofx（t）andh（t）
subplot（221）
plot（tx,x）,gridon,title（'Signalx（t）'）
xlabel（'Timetsec'）
subplot（222）
plot（th,h）,gridon,title（'Signalh（t）'）
xlabel（'Timetsec'）
subplot（313）
plot（y）,gridon,title（'Theconvolutionofx（t）andh（t）'）
xlabel（'Timetsec'）信号x（t）、h（t）和x（t）*h（t）的波形图
Q1-10：
给定两个离散时间序列
x[n]=0.5n{u[n]-u[n-8]}
h[n]=u[n]-u[n-8]
编写程序Q1_10，计算它们的卷积，并分别绘制x[n]、h[n]和它们的卷积y[n]的图形。
编写的程序Q1_10如下：
n=0:
10;
x=（0.5）.^n.*（u（n）-u（n-8））;
h=u（n）-u（n-8）;
y=conv（x,h）;%Computetheconvolutionofx（t）andh（t）
subplot（221）
stem（n,x,'.'）,gridon,title（'Signalx（n）'）
subplot（222）
stem（n,h,'.'）,gridon,title（'Signalh（n）'）
subplot（212）
stem（y）,gridon,title（'Theconvolutionofx（n）andh（n）'）,
xlabel（'Timetsec'）;
信号x[n]、h[n]和y[n]的波形图
Q1-11已知一个序列为
编写MATLAB程序Q1_11，能够将x[n]以N=8为周期进行周期延拓得到一个周期为N=8的周期序列y[n]，并分别绘制x[n]和y[n]图形。
编写的程序Q1_11如下：
U4.m
functiony=u4（n）
y=n.*（u（n）-u（n-5））;
Q1——11.m
clear,closeall;
n=-16:
32
x=u4（n）;
T=8;y=0;
fork=-2:
4;
y=y+u4（n-k*T）;
end
subplot（211）
stem（n,x,'.'）;
gridon,
title（'Originalsignalx（n）'）
xlabel（'Timetsec'）
subplot（212）
stem（n,y）;
title（'periodsignalx（n）'）
xlabel（'Timetsec'）
gridon,信号x[n]的波形图信号y[n]的波形图
Q1-12仿照范例程序Program1_7，编写程序Q1_12，计算并绘制由如下微分方程表示的系统在输入信号为x（t）=（e-2t-e-3t）u（t）时的零状态响应和你手工计算得到的系统零状态响应曲线。
手工计算得到的系统零状态响应的数学表达式是：
编写的程序Q1_12如下：
用MATLAB绘制的手工计算的系统响应
clear,closeall;
num=input（'Typeintherightcoefficientvectorofdifferentialequation：
'）;
den=input（'Typeintheleftcoefficientvectorofdifferentialequation：
'）;
t=0:
0.01:
8;
x=input（'Typeintheexpressionoftheinputsignalx（t）：
'）;
y=lsim（num,den,x,t）;plot（t,y）
执行程序Q1_12得到的系统响应
Q1-13：
利用程序Q1_9，验证卷积的相关性质。
（a）验证性质：
选择信号x（t）的数学表达式为：
sin（t）
x（t）、δ（t）和x（t）*δ（t）的波形
验证所得结论是：
（b）验证性质：
选择信号x（t）的数学表达式为：
sin（t）t0=2
x（t）、δ（t-t0）和
的波形
验证所得结论是：
（c）验证性质：
选择信号x（t）的数学表达式为：
sin（t）选择的t1=2秒，t2=3秒。
执行程序Q1_9，输入信号x（t-t1）和δ（t-t2）的数学表达式，得到的信号及其卷积的波形图如下：
执行程序Q1_9，输入信号x（t-t2）和δ（t-t1）的数学表达式，得到的信号及其卷积的波形图如下：
验证所得结论是：
（d）验证性质：
选择信号x（t）（建议选择一个时限信号）的数学表达式为：
u（t）-u（t-3）
的数学表达式为：
手工绘制的
波形如下：
执行程序Q1_9，输入信号x（t）和u（t）的数学表达式，得到的信号及其卷积的波形图如下：
验证所得结论是：
（e）验证性质：
选择信号x（t）的数学表达式为：
sin（t）选择信号h（t）的数学表达式为：
sin（t）
选择的t0=：
1
执行程序Q1_9，输入信号x（t）和h（t-t0）的数学表达式，得到的信号及其卷积的波形图如下：
执行程序Q1_9，输入信号x（t-t0）和h（t）的数学表达式，得到的信号及其卷积的波形图如下：
验证所得结论是：
Q1-14：
做如下总结：
1、信号与系统分析，就是基于信号的分解，在时域中，信号主要分解成：
2、写出卷积的运算步骤，并谈谈你对卷积的一些基本性质的理解。
利用MATLAB计算卷积的函数是什么？
如何使用？
3、在时域中，描述一个连续时间LTI系统的数学模型有：
4、MATLAB是如何表示一个由微分方程描述的连续时间LTI系统的？
求解连续时间LTI系统的单位冲激响应、单位阶跃响应以及系统在某一个输入信号作用下的零状态响应的MATLAB函数有哪些？
四、实验报告要求
1、按要求完整书写你所编写的全部MATLAB程序
2、详细记录实验过程中的有关信号波形图（存于自带的U盘中），图形要有明确的标题。
全部的MATLAB图形应该用打印机打印，然后贴在本实验报告中的相应位置，禁止复印件。
3、实事求是地回答相关问题，严禁抄袭。
本实验完成时间：
年月日
实验二连续时间信号的频域分析
三、实验内容及步骤
实验前，必须首先阅读本实验原理，读懂所给出的全部范例程序。
实验开始时，先在计算机上运行这些范例程序，观察所得到的信号的波形图。
并结合范例程序应该完成的工作，进一步分析程序中各个语句的作用，从而真正理解这些程序。
实验前，一定要针对下面的实验项目做好相应的实验准备工作，包括事先编写好相应的实验程序等事项。
Q2-1编写程序Q2_1，绘制下面的信号的波形图：
其中，ω0=0.5π，要求将一个图形窗口分割成四个子图，分别绘制cos（ω0t）、cos（3ω0t）、cos（5ω0t）和x（t）的波形图，给图形加title，网格线和x坐标标签，并且程序能够接受从键盘输入的和式中的项数。
抄写程序Q2_1如下：
clear,closeall
T=2;dt=0.00001;t=-2*pi:
dt:
2*pi;w0=0.5*pi;
x1=cos（w0*t）;
x3=（-1/3）*cos（3*w0*t）;
x5=（1/5）*cos（5*w0*t）;
N=input（'TypeinthenumberoftheharmoniccomponentsN=:
'）;
y=0;
forq=1:
N;%Synthesiztheperiodicsignaly（t）fromthefiniteFourierseries
y=y+（1/q）.*sin（（q*pi）/2）.*cos（q*w0*t）;
end;
subplot（221）,
plot（t,x1）,title（'Theoriginalsignalcos（w0t）'）;
gridon;axis（[-2*pi,2*pi,-1,1]）,xlabel（'Timet'）
subplot（223）,
plot（t,x5）,title（'Theoriginalsignal（1/5）cos（5w0t）'）;
gridon;axis（[-2*pi,2*pi,-1,1]）,xlabel（'Timet'）
subplot（222）
plot（t,x3）,title（'Theoriginalsignal（-1/3）cos（3w0t）'）;
gridon;axis（[-2*pi,2*pi,-1,1]）,xlabel（'Timet'）
subplot（224）
plot（t,y）,title（'Thesynthesissignalofx（t）'）;
gridon;axis（[-10,10,-1,1]）,xlabel（'IndexN'）
执行程序Q2_1所得到的图形如下：
N=10
Q2-2给程序Program2_1增加适当的语句，并以Q2_2存盘，使之能够计算例题2-1中的周期方波信号的傅里叶级数的系数，并绘制出信号的幅度谱和相位谱的谱线图。
通过增加适当的语句修改Program2_1而成的程序Q2_2抄写如下：
clear,closeall
T=2;dt=0.00001;t=-2:
dt:
2;
x1=u（t）-u（t-1-dt）;x=0;
form=-1:
1
x=x+u（t-m*T）-u（t-1-m*T-dt）;%Periodicallyextendx1（t）toformaperiodicsignal
end
w0=2*pi/T;
N=input（'TypeinthenumberoftheharmoniccomponentsN=:
'）;
L=2*N+1;
fork=-N:
1:
N;
ak（N+1+k）=（1/T）*x1*exp（-j*k*w0*t'）*dt;
end
phi=angle（ak）;
y=0;
forq=1:
L;%Synthesiztheperiodicsignaly（t）fromthefiniteFourierseries
y=y+ak（q）*exp（j*（-（L-1）/2+q-1）*2*pi*t/T）;
end;
subplot（221）
plot（t,x）,title（'Theoriginalsignalx（t）'）,axis（[-2,2,-0.2,1.2]）,gridon;
subplot（222）
k=-N:
N;stem（ak）,title（'Theakofx（t）'）,axis（[-1,1,-0.4,0.4]）,gridon;
subplot（223）
k=-N:
N;stem（k,abs（ak）,'k.'）,title（'Theamplitude|ak|ofx（t）'）,axis（[-N,N,-0.1,0.6]）,gridon;
subplot（224）
stem（k,phi,'r.'）,title（'Thephasephi（k）ofx（t）'）,axis（[-N,N,-2,2]）,xlabel（'Indexk'）,gridon;
执行程序Q2_2得到的图形
Q2-3反复执行程序Program2_2，每次执行该程序时，输入不同的N值，并观察所合成的周期方波信号。
通过观察，你了解的吉伯斯现象的特点是：
N=5N=10
N=20N=40
1、周期信号的傅里叶级数与GIBBS现象
给定如下两个周期信号：
Q2-4分别手工计算x1（t）和x2（t）的傅里叶级数的系数。
信号x1（t）在其主周期内的数学表达式为：
t+1-11-t0计算x1（t）的傅里叶级数的系数的计算过程如下：
k=-10:
10;ak=0;
ak=1/2.*（sin（（k）*pi/2）./（（k）*pi/2））
N=-10:
10;
stem（k,ak）;
通过计算得到的x1（t）的傅里叶级数的系数的数学表达式是：
1/2.*（sin（（k）*pi/2）./（（k）*pi/2））
信号x2（t）在其主周期内的数学表达式为：
1|t|<0.2
00.2<|t|<1
计算x2（t）的傅里叶级数的系数的计算过程如下：
k=-10:
10;ak=0;
ak=sin（k*pi*0.2）./（k*pi）
N=-10:
10;
stem（k,ak）;
通过计算得到的x1（t）的傅里叶级数的系数的数学表达式是：
sin（k*pi*0.2）./（k*pi）
用MATLAB帮助你计算出你手工计算的傅里叶级数的系数ak从-10到10共21个系数。
从命令窗口上抄写x1（t）的21个系数如下：
ak=
Columns1through8
0.00000.0354-0.0000-0.04550.00000.0637-0.0000-0.1061
Columns9through16
0.00000.3183NaN0.31830.0000-0.1061-0.00000.0637
Columns17through21
0.0000-0.0455-0.00000.03540.0000
从命令窗口上抄写x2（t）的21个系数如下：
Columns1through8
-0.0000-0.0208-0.0378-0.0432-0.03120.00000.04680.1
升级会员 