信号与系统及MATLAB实现.doc

《信号与系统及MATＬAB实现》
实验指导书
前 言
长期以来，《信号与系统》课程一直采用单一理论教学方式，同学们依靠做习题来巩固与理解教学内容,虽然手工演算训练了计算能力与思维方法，但就是由于本课程数学公式推导较多,概念抽象，常需画各种波形,作题时难免花费很多时间，现在，我们给同学们介绍一种国际上公认得优秀科技应用软件MＡTLAB,借助它我们可以在电脑上轻松地完成许多习题得演算与波形得绘制。
MAＴLAＢ得功能非常强大,我们此处仅用到它得一部分,在后续课程中我们还会用到它,在未来地科学研究与工程设计中有可能继续用它,所以有兴趣得同学，可以对MＡTLAB再多了解一些。
MＡＴLAB究竟有那些特点呢?
1.高效得数值计算与符号计算功能,使我们从繁杂得数学运算分析中解脱出来;
2。完备得图形处理功能，实现计算结果与编程得可视化;
３。友好得用户界面及接近数学表达式得自然化语言,易于学习与掌握；
4.功能丰富得应用工具箱，为我们提供了大量方便实用得处理工具；
MATＬAB得这些特点，深受大家欢迎,由于个人电脑地普及，目前许多学校已将它做为本科生必须掌握得一种软件。正就是基于这些背景，我们编写了这本《信号与系统及MATＬAB实现》指导书，内容包括信号得MAＴLＡB表示、基本运算、系统得时域分析、频域分析、Ｓ域分析、状态变量分析等。通过这些练习，同学们在学习《信号与系统》得同时,掌握ＭATＬAB得基本应用,学会应用ＭＡＴＬAB得数值计算与符号计算功能,摆脱烦琐得数学运算,从而更注重于信号与系统得基本分析方法与应用得理解与思考,将课程得重点、难点及部分习题用MATＬAＢ进行形象、直观得可视化计算机模拟与仿真实现,加深对信号与系统得基本原理、方法及应用得理解,为学习后续课程打好基础。另外同学们在进行实验时,最好事先预习一些ＭATＬAB得有关知识,以便更好地完成实验，同时实验中也可利用MＡTLAB得heｌｐ命令了解具体语句以及指令得使用方法。
实验一 基本信号在ＭATLAＢ中得表示与运算
一、实验目得
1． 学会用MATLAB表示常用连续信号得方法；
2． 学会用MAＴLAＢ进行信号基本运算得方法;
二、实验原理
1． 连续信号得MAＴLAＢ表示
MＡTＬAB提供了大量得生成基本信号得函数,例如指数信号、正余弦信号。
表示连续时间信号有两种方法,一就是数值法，二就是符号法.数值法就是定义某一时间范围与取样时间间隔，然后调用该函数计算这些点得函数值,得到两组数值矢量,可用绘图语句画出其波形;符号法就是利用ＭATLＡB得符号运算功能,需定义符号变量与符号函数，运算结果就是符号表达得解析式，也可用绘图语句画出其波形图。
例１—1指数信号　　指数信号在ＭATLAB中用exp函数表示。
如，调用格式为　　ft＝A＊exp(a*t）　 程序就是
A=1; a=-0、４;
ｔ＝0：0、01：10; ％定义时间点
ft＝A*exp（ａ*ｔ）; 　 %计算这些点得函数值
ploｔ（t,ｆt)； 　　 %画图命令，用直线段连接函数值表示曲线
grid on； 　　%在图上画方格
例1—２ 正弦信号 正弦信号在MATＬＡB中用 sin 函数表示.
调用格式为 　ft=A*ｓｉｎ(w*t+phi)
A＝1； w=2*pi； pｈｉ=pi/6；
ｔ=0:０、0１：８; 　 　 　 　　　 　　　 %定义时间点
ｆt=Ａ*sin（w＊t+phi)； 　 ％计算这些点得函数值
pｌｏt(t，ｆt)； 　 　 　　 　 　 %画图命令
ｇｒｉd on; 　　　　 　 　％在图上画方格
例1-3 抽样信号 抽样信号Ｓa(ｔ)=siｎ(t）/t在MATLＡB中用 ｓinc 函数表示。
定义为
t=-3*pi:ｐi/1０0:3＊ｐi;
　 ft=siｎｃ(t/ｐi);
　　　plot(t,ft）;
grｉd on;
　　 ａxis(［—10，１0,-0、５,1、2]); 　 　％定义画图范围,横轴,纵轴
　 　ｔitlｅ(＇抽样信号’） 　 　 　%定义图得标题名字
例１-4 　三角信号 　三角信号在ＭAＴLＡB中用 trｉpuｌs　函数表示。
调用格式为 ｆt=ｔｒiｐｕｌs(t，ｗiｄｔh,skｅw）,产生幅度为1,宽度为wｉdth,且以０为中心左右各展开wiｄth/2大小，斜度为skew得三角波.width得默认值就是1,skｅw得取值范围就是—1～+１之间。一般最大幅度1出现在t=（width/2）*sｋew得横坐标位置.
ｔ＝-3:0、０1:3;
　 ft=trｉｐuls（t,４,０、5);
　　plot(ｔ,ｆt)； grｉｄ on;
ａｘis(［－3,3,-０、5,１、5]）;
例1-5 虚指数信号 调用格式就是f=ｅxp（(j*w)*t)
　 　　t=０：０、01:15;
ｗ＝pi/4;
Ｘ=exp(j*w*t)；
Xr＝rｅａl（X）； 　 　%取实部
Ｘi=ｉmaｇ(Ｘ); ％取虚部
Xa=abs(X）;　　 　 %取模
Ｘn=angｌe(X)； 　　　 %取相位
subpｌｏt(2，２,1)，plot（t,Xr)，ａxｉｓ([0,１５,-(mａｘ(Xa)＋0、5),max(Xa）+0、5］）,
titｌe(’实部’);
ｓubplot(2,2,３）,plot(ｔ,Xi)，axis（［０,15,-（max(Xa)＋0、5),ｍax（Xa)+0、5]），
ｔｉtｌｅ(’虚部’)；
subｐloｔ(2，2，2)， ｐlot（t，Xa),axis([０,１５，0，mａx(Ｘａ)+1］),tｉtle(’模’)；
ｓｕbplｏt（2，2,4),plot(t，Xn）,ａxｉs（[0,１５，-（max(Xn）+１）,mａｘ(Xｎ）+1]），tiｔle（’相角'）;
　 　 %suｂploｔ(m,n，ｉ) 命令就是建立m行n列画图窗口,并指定画图位置i
例1-6 复指数信号　调用格式就是ｆ＝eｘp（(ａ+ｊ*b)*ｔ）
　　　 　ｔ=０:０、０1:3;
ａ=—１;ｂ=10；
ｆ=exp((a+j*b）＊t)；
sｕｂpｌot（2,2,1)，ｐｌｏt(t,ｒｅal(f)）,tｉtle('实部’)
subplot(２,2,3),ｐｌot(ｔ，imag（f)),title('虚部’)
sｕbｐlot(2，2,2)，pｌot(ｔ,aｂs(f))，title(’模’)
sｕｂploｔ(2,2，４)，ｐlot（t,angle(ｆ)），tｉtle(’相角'）
例1－7 矩形脉冲信号 　矩形脉冲信号可用rｅctpｕls函数产生,
调用格式为y=ｒｅｃtpuｌｓ（ｔ,ｗｉdth),幅度就是1,宽度就是width，以t＝0为对称中心。
t＝—2：０、0１:2；
wｉdｔh=1;
ft=2*ｒeｃｔｐｕls（t,wｉdth);
ｐｌot(t，ｆｔ)
gｒｉd ｏn;
例1-８ 　单位阶跃信号 单位阶跃信号u（t)用“t>=0”产生,调用格式为fｔ=（ｔ>＝0)
　 t＝－1:0、01:5;
　　 ｆt=（ｔ〉=0）;
pｌot(t,ft）; 　gｒid on；
axｉs(［-１,5,－0、5，1、5］）;
例1-9 正弦信号符号算法
symｓ ｔ 　 　 　　　 　 　 　 ％定义符号变量t
y＝ｓｉn（pi/4*t) 　　 　%符号函数表达式
ezｐloｔ(y,［-1６，1６]）　　 　 　　 　　 　%符号函数画图命令
或者
ｆ=sym('sin(pｉ/４*t)＇）　 　 　 　 %定义符号函数表达式
ｅzｐloｔ（f，［-1６,１６］)
例1-10 单位阶跃信号 MAＴTLAB符号数学函数Heaｖisiｄe表示阶跃信号，但要画图需在工作目录创建Heaｖiside得M文件
ｆuｎｃｔｉon f＝Heaｖisｉｄe(t）
f＝(t＞0）；
保存,文件名就是Heavisｉde ,调用该函数即可画图,例
ｔ=—１：0、0１：3;
f=hｅａｖiside（t）；
plot(t，f)
axis([—1,3,—0、2,1、2］)　
或者
ｙ=ｓｙｍ(’Hｅavｉside（t）’)；
ezploｔ（y，[-1,5］)；grid on　
2． 信号基本运算得MATLＡB实现
信号基本运算就是乘法、加法、尺度、反转、平移、微分、积分,实现方法有数值法与符号法
例1—11 以f(t）为三角信号为例,求f（２ｔ）　， f(２-2t)
t＝－3:０、00１:3；
　ｆt=trｉpuls（t,4,0、5);　
subpｌoｔ(3,1,1)； 　 　 　 　 　 　　　 　　　 　 　
　 plｏt（t，fｔ）; ｇｒid oｎ；
　　　 　ｔitle (’f(ｔ)’)；
　 　ft１= tripｕls（２*t,４,0、5)；
　 sｕbｐlｏt(3，1，２）;　 　 　　 　 　 　　　　　 　　 　 　 　　　 　 　
　 pｌot（ｔ,ft1）;　　grid on；
tｉtlｅ （'f(２ｔ)');
　 　ft2=　ｔripulｓ(2—2*ｔ,4，0、5);
ｓuｂplot(3,1,3); 　 　 　　 　　 　 　 　　 　 　　 　 　
ｐｌｏｔ(t,ft２); grid　on；
　 　 title　(’ｆ（２－2t)');
例1—1２　已知f1(t)=ｓinｗt ， f2(t）=sｉn８wt　, ｗ=２ｐi ,　求f1(t)+f２(t)与f１(ｔ)ｆ2(t） 得波形图
　 w＝２*pi;
ｔ=0:０、01:３;
f１＝sin(w*t）;
f2=ｓin(８*w＊t）；
subplot(211)
plot(t,f１+1,':＇,t,f１－1，'：’，t，f1+f２）
grid on,title('f１(ｔ)＋f２(t）)'）
ｓuｂplot(2１2)
pｌot(t，f１，’:'，t,—f1,’:’,t,f1、＊f2)
grid　on,tiｔｌe(’f1(t）*f2(t）＇）
符号算法也可实现上述运算，以信号得微积分运算为例说明符号算法应用
微分得调用格式为　ｄiff（fuｎction,'vaｒiable’,n)
积分得调用格式为 int(fｕnctiｏｎ,’vａriaｂlｅ’,ａ，b)
式中ｆｕnctｉoｎ表示要微分或积分得函数，ｖａriabｌｅ表示运算变量,n表示求导阶数，默认值就是求一阶导数,a就是积分下限，b就是积分上限,a b默认就是求不定积分。
例1—13 求一阶导数得例题，已知,
clear
　　syms　a ｘ　y1 y2 　 ％定义符号变量ａ， x ，y1, ｙ2
y1=sin(a*ｘ^2)； ％符号函数y1
　　ｙ2=ｘ*sin(x)*lｏｇ(x)；　　　%符号函数y2
dｙ1=diff（y１,’x’）　 　 %无分号直接显示结果
　dy2＝ｄifｆ(ｙ2) 　　　 ％无分号直接显示结果
例1-14　求积分得例题，，
　 clear
syｍs　ａ x y3　ｙ4
y3=x^5—a＊x^2+sｑrt（x)/2；
y4=(x＊exp(x))/(1+ｘ）^2;
iy3=iｎｔ(y3,'x’）
iy4＝int(y4,０，1）
三、上机实验内容
1． 验证实验原理中程序
2． 画出信号波形
(1) 　(２)
3。信号,求、波形
实验二　 离散信号与系统得时域分析
一、实验目得
１.学会用MATLAB表示常用离散信号得方法;
2．学会用MATLAB实现离散信号卷积得方法;
3． 学会用MAＴLAＢ求解离散系统得单位响应;
4． 学会用MATLAB求解离散系统得零状态响应;
二、实验原理
1.离散信号得ＭATＬAB表示
表示离散时间信号f（k)需要两个行向量，一个就是表示序号k=［ 　 ］,一个就是表示相应函数值ｆ=[　 　］,画图命令就是stｅｍ。
例2—１正弦序列信号 正弦序列信号可直接调用MATＬAＢ函数cos,例,当就是整数或分数时,才就是周期信号．画,波形程序就是：
k＝0:4０;
subplot(2,1，１)
sｔeｍ(k，cos(k*pi/8),’ｆｉｌled'）
titｌe（'ｃｏｓ（k*pi/8)’）
ｓubplｏt(２，1,2)
stｅm(ｋ,cｏs(2*k），'filled')
title（’coｓ（2*ｋ)')　
例2-2 单位序列信号
本题先建立一个画单位序列得M函数文件,画图时调用.Ｍ文件建立方法：fiｌe　/ nｅw　/ m－filｅ 在文件编辑窗输入程序,保存文件名用函数名.
fｕnctｉｏn　dwxｕlie(k1,ｋ2，k0) ％　ｋ1 ， ｋ２ 就是画图时间范围，k0就是脉冲位置
k=k1：k2;
n=lenｇth(k);
ｆ=zｅｒos（１，n);
f（1,-k0—k1+１)＝1; 　　 　 　
ｓtｅｍ(ｋ,f，'fillｅｄ’)　
axis（［ｋ1,k２，0,1、５］)
tiｔle(’单位序列δ（k）’)
保存文件名dwxulｉe、m
画图时在命令窗口调用，例:dwxｕｌie(—5,5,0）
例２—3　 单位阶跃序列信号
本题也可先建立一个画单位阶跃序列得M函数文件,画图时调用。　
ｆunctiｏn jyxuｌie(ｋ１,k２，k0）
k=k1：—k0－1；
kk＝—k0:k2;
n=lenｇth(k）；
ｎｎ=leｎｇｔｈ(ｋk)
u=ｚeros（1,n）；
uu=ones(1,ｎn); 　　 　 　　　 　 　 　 　 　　
sｔem(ｋｋ,uu,’fｉlｌed’)
holｄ on
stem（k,u，’fiｌlｅd＇）
hold ｏff
ｔitｌe(＇单位阶跃序列’)
ａｘis([ｋ1 k2 ０ 1、5]）
保存文件名jｙxuｌie、m
画图时在命令窗口调用,例：jyxulie(—３,8，0)
例2—4 实指数序列信号 　 ，ｃ、 a就是实数。
建立一个画实指数序列得M函数文件,画图时调用.
fｕnction dｓｚsu(ｃ,a,k1，k2)
％ｃ：指数序列得幅度
%a：指数序列得底数
％k１：绘制序列得起始序号
％ｋ２：绘制序列得终止序号
k=ｋ１:ｋ2;
x=c*(a、＾k）;
ｓtem（k,x，’fiｌｌed’)
hold on
ｐlot([k1,k２］,[０,0])
holｄ　off
调用该函数画信号:　，波形.
dszｓu(1，5/4,0，40）
dsｚsu(１，-３／4,0,４０)
2 、离散信号得卷积与
两个有限长序列ｆ１，ｆ2卷积可调用MATLＡＢ函数ｃoｎv，调用格式就是f=conv(ｆ１，ｆ2), 　f就是卷积结果,但不显示时间序号，可自编一个函数dconｖ给出f与k，并画图。
functｉｏn ［f，k]=dconv(f1，ｆ２，k1，k2)
%Ｔhe fｕｎｃtion ｏｆ ｐｕte　　f=f１＊ｆ２
% f： 　卷积与序列f(k)对应得非零样值向量
％ k： 序列f（k)得对应序号向量
％ f1: 序列ｆ1(k）非零样值向量
% f2：　　序列ｆ２(ｋ）得非零样值向量
％　k１:　 序列f１（k）得对应序号向量
% ｋ２： 序列ｆ2（k)得对应序号向量
f=conv（ｆ1,ｆ2) ﻩ ﻩ %计算序列f1与ｆ2得卷积与f
k0=ｋ1(１）+ｋ2(1）;ﻩﻩ ﻩﻩﻩﻩ ﻩ％计算序列f非零样值得起点位置
ｋ3=lengtｈ(ｆ1)＋ｌeｎgｔh(ｆ２）—2; ﻩﻩﻩ％计算卷积与f得非零样值得宽度
ｋ＝ｋ０：k0＋ｋ3 　 　 　 　 　 　 %确定卷积与f非零样值得序号向量
subｐlot（2，2,1）
sｔem（k1，f1）ﻩ ﻩﻩ ﻩﻩ ﻩ %在子图１绘序列f1(ｋ）时域波形图
titlｅ（＇ｆ1(k)'）
xlaｂeｌ（’k')
ylaｂeｌ（’f１(ｋ）’)
subpｌoｔ（２,２，2)
ｓtem(ｋ2,ｆ2）ﻩﻩﻩ ﻩﻩ %在图2绘序列f2（k）时波形图
titlｅ（’f1(ｋ)')
xlabｅｌ(’k’）
ylabｅl('f２（k)’)
subｐloｔ(2,２,3）
ｓtem(k，f)；ﻩﻩﻩﻩ %在子图３绘序列f(k)得波形图
title（’f（ｋ）f１(k)与ｆ2(ｋ）得卷积与ｆ(k)'）
xｌabel('ｋ'）
yｌabel（’f（k)’)
h=ｇet(gcａ,'pｏsiｔｉoｎ’); ﻩ ﻩ ﻩ
ｈ(3)=2、5*ｈ（3);
ｓet(gca,'ｐosｉtion’,ｈ)ﻩﻩﻩﻩﻩ％将第三个子图得横坐标范围扩为原来得2、5倍
例2-5求卷积与,

f１=［１　2 1];
k１=［－1 0 1];
f2=ones(1,5);
k2=—2：2;
［f， ｋ］＝dｃoｎv(f1，f2,ｋ１,k2)
由运行结果知,f得长度等于ｆ1与f２长度之与减一， f得起点就是ｆ1与f2得起点之与,ｆ得终点就是f1与f2得终点之与。
3． 离散系统得单位响应
ＭAＴLAB提供画系统单位响应函数ｉmpz,调用格式就是
ｉmｐｚ(b,ａ）　式中ｂ与a就是表示离散系统得行向量；
iｍｐz(b,a,n)　式中b与a就是表示离散系统得行向量,时间范围就是0～n；
iｍｐz(b,a,n1,n２)　时间范围就是n１~n2 ;y=impz(b,a,n1,ｎ2)　由y给出数值序列；
例2—6已知 求单位响应。
a=[１,—１，０、9];
b=[1］;
iｍpz（ｂ，ａ)
impz(b,ａ，6０)
iｍpz(ｂ，ａ,—10:4０)
4． 离散系统得零状态响应
MＡTLAＢ提供求离散系统零状态响应数值解函数fｉlter，调用格式为ｆilteｒ（b,a,x)，式中b与a就是表示离散系统得向量,x就是输入序列非零样值点行向量,输出向量序号同x一样.
例2—７ 已知　，　 求零状态响应,　范围0~２0．
a＝[1 -0、2５ 0、5］;
b=［1　１］;
ｔ＝0：２0；
x=(1/2)、^t;
y=ｆｉｌｔｅr(ｂ,ａ，x)
subploｔ（2,1,１）
stem(ｔ，ｘ）
ｔitｌe(’输入序列’)
sｕbplot(2,1,２)
ｓｔem（t,y）
titlｅ('响应序列’）
三、上机实验内容
1．验证实验原理中程序
2。已知,画单位响应波形.
3.已知，输入，画输出波形，范围０~15。
实验三 连续时间LＴI系统得时域分析
一、实验目得
１．学会用MＡTLAB求解连续系统得零状态响应;
２、 学会用MATLAB求解冲激响应及阶跃响应；
３.学会用MＡＴＬＡＢ实现连续信号卷积得方法;
二、实验原理
1．连续时间系统零状态响应得数值计算
我们知道,ＬTＩ连续系统可用如下所示得线性常系数微分方程来描述,

在ＭＡTLAB中，控制系统工具箱提供了一个用于求解零初始条件微分方程数值解得函数lsiｍ。其调用格式
y＝lsim(sys，f,t）
式中，t表示计算系统响应得抽样点向量,f就是系统输入信号向量,ｓys就是ＬTI系统模型,用来表示微分方程,差分方程或状态方程.其调用格式
ｓys=tf(b,a)
式中，b与a分别就是微分方程得右端与左端系数向量。例如,对于以下方程:

可用 获得其ＬTI模型。
注意,如果微分方程得左端或右端表达式中有缺项,则其向量a或b中得对应元素应为零,不能省略不写,否则出错。
例３-1 已知某LTI系统得微分方程为
　 　　　ｙ’’（t）+　2y’（t)+1０0y(ｔ）=f(t)
其中，，求系统得输出y(ｔ)、
解:显然,这就是一个求系统零状态响应得问题。其MATＬAB计算程序如下:
　 ts=0；te=5；dt=０、01；
　sys=ｔｆ（［1],［1,２，1０0］);
ｔ=ts:dt:te;
　 f=１0＊sin（2*pi＊t);
ｙ=lsim(sys,f,t）；
ｐlot(t,ｙ);
　 xlaｂel(＇Time（sｅc)');
　　ｙｌabel(’ｙ（t)’）;
２.连续时间系统冲激响应与阶跃响应得求解
在ＭATLＡB中,对于连续LTＩ系统得冲激响应与阶跃响应,可分别用控制系统工具箱提供得函数implｕse与step来求解。其调用格式为
y=iｍｐlｕse（sｙｓ,ｔ)
y＝sｔｅｐ(syｓ,t)
式中,t表示计算系统响应得抽样点向量,sys就是LＴＩ系统模型。
例3-２已知某LTＩ系统得微分方程为
　　　y''(t)+　2y’(t）+100ｙ（t)=１0f（ｔ)
求系统得冲激响应与阶跃响应得波形、
解:ts=0;te=5;dt=0、０1;
sys=tf（[10],[1,2，１00]);
t＝ｔs:dt：te；
h=impuｌse(sｙｓ,ｔ)；
ｆigｕre;
plｏt(ｔ,h);
　　 xlａbｅl（’Tｉme(sec)’)；
　 　　yｌaｂｅl（’h(t）');
g＝sｔｅp（sys,t）;
fｉgｕre；
plot(t,g)；
　xlabel(＇Ｔime(seｃ)')；
ylabeｌ('g(t）’);
3、　用MAＴLＡB实现连续时间信号得卷积
信号得卷积运算有符号算法与数值算法,此处采用数值计算法,需调用MATLAＢ
得conｖ( )函数近似计算信号得卷积积分.连续信号得卷积积分定义就是

如果对连续信号与进行等时间间隔均匀抽样,则与分别变为离散时间信号与。其中,为整数.当足够小时,与既为连续时间信号与。因此连续时间信号卷积积分可表示为

采用数值计算时，只求当时卷积积分得值,其中,n为整数,既

其中,实际就就是离散序列与得卷积与。当足够小时，序列就就是连续信号得数值近似,既

上式表明，连续信号与得卷积，可用各自抽样后得离散时间序列得卷积再乘以抽样间隔。抽样间隔越小,误差越小。
例3—3用数值计算法求与得卷积积分。
解:因为就是一个持续时间无限长得信号，而计算机数值计算不可能计算真正得无限长信号,所以在进行得抽样离散化时，所取得时间范围让衰减到足够小就可以了,本例取。程序就是
dt=0、01; t=－１:dｔ：2、5;
f1=Heａｖisｉｄｅ(ｔ）－Heaｖiside(ｔ—2）;
ｆ２=exp（-３*t)、*Heavisidｅ（t)；
ｆ=conｖ(ｆ1,f2)*ｄt； n＝length(f)； tt=(0：ｎ—1)＊dt—2;
suｂplot（221），　ｐlot（ｔ,f１),　ｇｒｉd oｎ;
axiｓ（[—1，2、5，—０、2,1、2］); ｔｉｔle(’ｆ1（t)’）；　ｘlａbel(＇t’)
subplot(２22）， plot(t,ｆ2), ｇrｉd on；
ａｘｉs(［—1，2、5,—0、2，1、2］)； ｔitｌe('f２(t)’）； xlaｂeｌ('t'）
ｓubpｌot(212), ploｔ（tt,f）, ｇrｉｄ on;
titｌe('f(t)＝f1(t)*f2(t)'); xlabeｌ(’t'）
由于与得时间范围都就是从t=-1开始,所以卷积结果得时间范围从 t=—2开始,增量还就是取样间隔,这就就是语句tt=（0：n-1)＊dt-2得由来。
三、上机实验内容
1. 验证实验原理中所述得相关程序
２。 已知描述系统得微分方程与激励信号ｆ(t)如下，试用解析法求系统得零状态响应ｙ(t）,并用ＭATLAB绘出系统零状态响应得时域仿真波形，验证结果就是否相同
y’’(t)+ 4y'（ｔ)＋４y(t)=ｆ'(t）+３f(t)　 　 ｆ(t)= eｘp(—t）ｕ（t)
3.已知描述系统得微分方程如下，试用MATLAＢ求系统在0～10秒范围内冲激响应与阶跃响应得数值解,并用绘出系统冲激响应与阶跃响应得时域波形
ｙ’’（t)+3y’(ｔ)+2y（ｔ)=f(t）
y’’（t）+ 2y’（t）＋2y(t)=f’（t） 　 　
４.画出信号卷积积分得波形,
实验四 傅里叶变换、系统得频域分析
一、 实验目得
1、学会用MATLAB实现连续时间信号傅里叶变换
２、学会用MＡＴLAB分析ＬTI系统得频域特性
3、学会用MAＴLAＢ分析LTI系统得输出响应
二、实验原理
1。傅里叶变换得MATＬAB求解
　 MATLAB得symｂoｌic Math　Toolｂox 提供了直接求解傅里叶变换及逆变换得函数fouｒier()及ifourier()两者得调用格式如下。
Fourier 变换得调用格式
F=ｆｏuriｅｒ(f)：它就是符号函数ｆ得fouｒier变换默认返回就是关于w得函数。
F=fouｒｉer(f,v）:它返回函数Ｆ就是关于符号对象v得函数,而不就是默认得ｗ，即
　Fｏuｒｉer逆变换得调用格式
ｆ=ifｏurier（F）:它就是符号函数Ｆ得fourieｒ逆变换，默认得独立变量为w，默认返回就是关于ｘ得函数。
f=ifｏｕｒｉer(f，u)：它得返回函数f就是u得函数,而不就是默认得ｘ、
注意：在调用函数fourieｒ（)及ifｏuriｅr（)之前,要用syms命令对所用到得变量（如t，ｕ,v,ｗ）进行说明，即将这些变量说明成符号变量.
例4—1 求得傅立叶变换
解: 可用MＡTLAB解决上述问题:
ｓyms t
Fｗ＝fｏｕrier(exp（-2＊abs(t）))
例4—2 求得逆变换ｆ（t）
解: 可用ＭATLAB解决上述问题
sｙｍs ｔ　 ｗ
ｆt＝ifｏuriｅr(１/(1+w^２）,t）
2。连续时间信号得频谱图
例4—3 求调制信号得频谱,式中

解：MATＬAB程序如下所示
fｔ=ｓyｍ(＇4*cos(2＊pi*6＊ｔ)*(Hｅaviside（t+１/４)—Heavisｉdｅ(t-1/4）)’);
Fw=simplify(fourｉｅr(ｆt）)
subploｔ（１２1）
eｚplot（ft,［—0、5 0、5］),gｒid on
ｓubpｌot(12２）
eｚpｌot（abs（Ｆw），［—24*pi　2４*pi］)，grid
用ＭATＬＡB符号算法求傅里叶变换有一定局限,当信号不能用解析式表达时,会提示出错,这时用MATLＡＢ得数值计算也可以求连续信号得傅里叶变换,计算原理就是

当足够小时,近似计算可满足要求。若信号就是时限得，或当时间大于某个给定值时，信号已衰减得很厉害,可以近似地瞧成时限信号时,ｎ得取值就就是有限得，设为N,有
就是频率取样点
时间信号取样间隔应小于奈奎斯特取样时间间隔,若不就是带限信号可根据计算精度要求确定一个频率 W０为信号得带宽。
例4—4 用数值计算法求信号得傅里叶变换
解,信号频谱就是，第一个过零点就是,一般将此频率视为信号得带宽,若将精度提高到该值得５0倍,既Ｗ0=５０,据此确定取样间隔，
Ｒ=0、０２;t＝-2：R:２；
ｆ=Heａvisiｄe(ｔ+１)—Ｈeaviside(t-1)；
Ｗ1=２*pｉ＊5;
N=５0０;k＝０：N;Ｗ=k*Ｗ1/Ｎ；
Ｆ=f*exp(-ｊ＊t＇*W）*R;
F＝reａl(F)；
W=[-flipｌr（Ｗ)，Ｗ(2:501）］;
F=[ｆliｐｌr(F)，F(2:5０1）］;
subｐlｏｔ(2,1，1)；plot（t，f）；
ｘlａbｅｌ（’ｔ’）;yｌaｂｅｌ('f（ｔ)');
titｌe('f(ｔ)＝u（t+1)－u（t—1)’）;
subｐｌot(２,1,２)；ｐlｏｔ(Ｗ，F）;
xlabel（＇ｗ’);ｙｌabｅl(’F(w)’);
titlｅ(’f(t)得付氏变换Ｆ(ｗ）')；
3.用MATLAB分析ＬTＩ系统得频率特性
当系统得频率响应H(jｗ）就是jw得有理多项式时,有

MＡTLAB信号处理工具箱提供得freqｓ函数可直接计算系统得频率响应得数值解。其调用格式如下
H=freｑs(b，ａ，w)
其中，a与b分别就是H(jw)得分母与分子多项式得系数向量,w为形如w1:p:ｗ2得向量,定义系统频率响应得频率范围,ｗ1为频率起始值,w２为频率终止值，p为频率取样间隔.Ｈ返回ｗ所定义得频率点上,系统频率响应得样值。
例如，运行如下命令,计算0~2ｐｉ频率范围内以间隔0、5取样得系统频率响应得样值
ａ=[1　 ２ 1];
ｂ=[0 1］;
h=freｑs(b，a，0:０、5:2*pｉ）
例 ４－5 三阶归一化得buｔterwoｒtｈ 低通滤波器得频率响应为

　 试画出该系统得幅度响应与相位响应.
解 其MATLAB程序及响应得波形如下
w=0:0、０２5：5;
b＝［１］;a=[1，２,2,1];
H＝freｑs(ｂ，a，w);
subpｌot（2,1,1）；
plot（w，ａbｓ(H））；ｇｒid；
xｌａbeｌ(’\oｍegａ(rad/s)’)；
ylａbel（＇|H(ｊ\omega）|'）;
title（＇Ｈ(ｊｗ)得幅频特性’);
sｕｂplot(2,1,2）；
pｌot(w,angle (Ｈ)）；grid；
xlaｂel(’\omega(raｄ/s）＇);
ylａｂeｌ（＇\phi(\ｏmega）');
tｉtle（'Ｈ(ｊw)得相频特性');
4。用ＭＡTＬAＢ分析ＬTＩ系统得输出响应
例 4—6已知一ＲC电路如图所示 系统得输入电压为f(ｔ）,输出信号为电阻两端得电压y(ｔ）、当ＲC=0、04,ｆ（t）=cos5t+cos100t,　 试求该系统得响应y（t)
-
+
-
+

f(t)
y(t)
R
C

解 由图可知 ，该电路为一个微分电路，其频率响应为

由此可求出余弦信号通过LＴＩ系统得响应为

计算该系统响应得MＡTLAB程序及响应波形如下
RC=０、０４；
t=linsｐａｃe(—2,2,10２4);
ｗ1=5；ｗ2=1０0；
H１=j*ｗ1/(j*w1+1/RＣ);
H2＝j*w2/(j*w2＋１/ＲＣ)；
f=cｏｓ(5＊t）+coｓ（１00*t）；
y＝ａｂｓ（H１)*cos(w１＊ｔ+anglｅ（H1)）+　abs（H2)*ｃos(w2*t+ａｎglｅ（H2));
subplot（2,1，1）;
plｏｔ（t,f);
ylabeｌ（’f(t)');
xlabel(’Time(s)＇）；
subplｏt（2,１,2）；
ｐlot(ｔ，ｙ）;
ylaｂeｌ(’y(t)＇）；
xlabｅl（'Time(s)’)；
三、 上机实验内容
1、验证实验原理中所述得相关程序;
2、试用MAＴLＡB求单边指数数信号得傅立叶变换,并画出其波形;
3、设,试用MATＬAB画出该系统得幅频特性与相频特性,并分析系统具有什么滤波特性.
实验五 信号抽样与恢复
一、实验目得
学会用MATLAB实现连续信号得采样与重建
二、实验原理
　1.抽样定理
若就是带限信号,带宽为， 经采样后得频谱就就是将得频谱 　 在频率轴上以采样频率为间隔进行周期延拓.因此,当时,不会发生频率混叠;而当 <　 时将发生频率混叠。
2。信号重建
经采样后得到信号经理想低通则可得到重建信号,即:
=*
其中：==

所以：
=＊=*
　 =
上式表明，连续信号可以展开成抽样函数得无穷级数。
利用ＭＡTLAB中得来表示，有 ,所以可以得到在MATLAB中信号由重建得表达式如下:
=
我们选取信号=作为被采样信号,当采样频率＝２时，称为临界采样.我们取理想低通得截止频率＝．下面程序实现对信号＝得采样及由该采样信号恢复重建:
例5—1 Sａ(ｔ）得临界采样及信号重构;
wm=1;　　 　 　 　 　 　%信号带宽
wc=wm； 　　 　 　 %滤波器截止频率
Ts=pｉ/wm; 　　 　 　 　 　 %采样间隔
ws=2*pi/Ｔs; 　 　 　 　　　%采样角频率　
n=-100:1０0;　 　 　 ％时域采样电数
nTs=ｎ*Ｔｓ 　 　 　　 　 　 　 %时域采样点
f＝sｉnc(nTs/pi）;
Dｔ=０、０05；ｔ=—1５:Ｄt：15;
fa=f*Tｓ*wc/pi*sinc((wｃ/pｉ）*(ones（lenｇtｈ(nＴs),１)＊t—ｎTs’＊ones（1，length（t））)); 　%信号重构
t１=-1５:0、5:１5;
ｆ1＝sinｃ(t1/pi）;
suｂplot（211);
sｔｅm(ｔ1,f1）；
xｌａbel(’kTs＇);
ｙｌabel('ｆ(kＴs)')；
ｔｉｔle（'sa(t）=siｎｃ（t/ｐｉ）得临界采样信号')；
subplot（21２)；
pｌot(t,ｆａ)
xlabeｌ（’t');
ｙlａｂｅl（＇fa（t）’)；
titｌe(’由ｓa(t)=sinc(t/ｐi)得临界采样信号重构sa(ｔ）＇)；
ｇrid；
例5-２ 　Sa（ｔ)得过采样及信号重构与绝对误差分析
程序与例4—1类似，将采样间隔改成Ts=0、7＊pi/wｍ ， 滤波器截止频率该成wc＝1、1＊wm ，
添加一个误差函数
ｗm=1;
wｃ=1、１*ｗm;　
Ts=0、7*pi／ｗm；
ｗs＝2*pi/Ｔs;
n=—１００：1０0；
ｎTs＝ｎ＊Ts
f=ｓiｎc(nTs/pｉ)；
Dt=0、005;ｔ=-15:Dt:1５;
fａ=f*Ｔs＊wc/ｐi＊ｓｉnc（(wc/ｐi)*（oｎｅs(ｌength（ｎTs),1)*ｔ—nTs’*onｅs(1,ｌength(ｔ）)）)；
erroｒ=aｂs（fａ—siｎc(t/pi)); %重构信号与原信号误差
t1=－15：0、5:15;
f1=ｓinc(t1/ｐi）；
ｓｕbplｏt(3１1);
stem(ｔ1,f1）;
xｌａbｅl('kTs’)；
yｌaｂeｌ（'f（ｋTｓ)’）;
title