信号与系统课程设计报告.docx
- 文档编号:3725340
- 上传时间:2022-11-24
- 格式:DOCX
- 页数:19
- 大小:229.05KB
信号与系统课程设计报告.docx
《信号与系统课程设计报告.docx》由会员分享,可在线阅读,更多相关《信号与系统课程设计报告.docx(19页珍藏版)》请在冰豆网上搜索。
信号与系统课程设计报告
课程设计任务书
课程名称信号与系统课程设计
题目连续LTI系统的分析
专业通信工程
班级通信1102
学号
姓名
起止日期2013.06.26
设计地点计通楼304
指导教师
职称
分数
辽宁石油化工大学计算机与通信工程学院
一、课程设计的内容和要求
用MATLAB软件设计连续LTI系统的分析。
参考题目如下:
1、已知信号f1(t)=u(t)-u(t-2),f2(t)=cos(2πt)用MATLAB绘制
f1(t)+f2(t)和f1(t)×f2(t)的波形。
2、已知某连续系统的微分方程
2y''(t)+y'(t)+8y(t)=f(t)
f(t)=e-2tu(t)
试用MATLAB绘制出该系统的冲激响应和阶跃响应的波形。
3、编程求解正弦函数f(t)=sin(t)、取样函数f(t)=Sa(t)、门函数
、双边指数函数
信号的傅立叶变换式。
4、已知某连续系统的的系统函数为
H(s)=(s2+3s+2)/(8s4+2s3+3s2+5),试用MATLAB求出该系统的零极点,画出零极点分布图
5、设有两个稳定的LTI系统,可分别由下列微分来描述:
试分别画出它们的系统频率响应的幅值和相位特性曲线。
6、信号fl(t)和f2(t)如图所示。
取t=0:
0.005:
2.5,计算信号f(t)=f1(t)+f2(t)cos(50t)的值并画出波形;
而一可实现的实际系统的H(w)为
H(w)=
用freqs函数画出H(w)的幅度和相位曲线。
用lsim函数求出信号f1(t)和、f2(t)通过系统H(w)的零状态响应y1(t)和y2(t),并根据理论知识解释所得的结果。
二、对课程设计成果的要求(包括课程设计说明书(论文)、图纸、实物样品等)
1.课程设计要求符合《课程设计说明书规范》
2.图表格式要规范。
3上机运行结果,应满足课程设计任务书的基本要求,并有明确的结论。
三、课程设计进度计划
1)知识回顾和课程设计安排,集中进行,半天
2)查阅文献,分散进行,1天
3)数字信号处理设计并实施,两个半天
4)整理文档,编写课程设计说明书,两天
系(教研室)主任审查意见
系主任(签字):
年月日
摘要
信号与系统课程设计是在学习完信号与系统以及复变函数与积分变换等课程的基础上设立的一门实践课。
通过本次设计让学生学会使用MATLAB软件,同时学会建立m文件,并且在m文件下建立MATLAB函数和语句以及对其进行调试、运行、函数图像处理保存。
进一步加深对信号的各种运算的理解,提高学生的综合分析问题和解决问题的能力,为后续课程和相关专业研究奠定基础
Abstract
Signalandsystemcourseisdesignedafterstudyingsignalsandsystems,andcomplexfunctionandintegraltransformcourseestablishedonthebasisofapracticallesson.ThroughthedesigntoletthestudentslearnhowtousetheMATLABsoftware,atthesametimelearnhowtobuildmfile,andtheMATLABfunctionandstatementinthe.mfileandcarriesonthedebuggingandrunning,thefunctionofimageprocessing.Furtherdeepentheunderstandingofthevariousoperationsofthesignal,improvingthestudents'comprehensiveproblemanalysisandproblemsolvingskills,laythefoundationforfollow-upcoursesandrelatedresearch
一、中英文摘要
二、正文
2.1设计题目
2.2试验环境
2.3实验目的
2.4课程设计原理及说明
2.5课程设计的过程及调试
2.6设计的收获及结论
2.7主要参考文献
一、设计题目
1、已知信号f1(t)=u(t)-u(t-2),f2(t)=cos(2πt)用MATLAB绘制
f1(t)+f2(t)和f1(t)×f2(t)的波形。
2、已知某连续系统的微分方程
2y''(t)+y'(t)+8y(t)=f(t)
f(t)=e-2tu(t)
试用MATLAB绘制出该系统的冲激响应和阶跃响应的波形。
3、编程求解正弦函数f(t)=sin(t)、取样函数f(t)=Sa(t)、门函数
、双边指数函数
信号的傅立叶变换式。
4、已知某连续系统的的系统函数为
H(s)=(s2+3s+2)/(8s4+2s3+3s2+5),试用MATLAB求出该系统的零极点,画出零极点分布图
5、设有两个稳定的LTI系统,可分别由下列微分来描述:
试分别画出它们的系统频率响应的幅值和相位特性曲线。
6、信号fl(t)和f2(t)如图所示。
取t=0:
0.005:
2.5,计算信号f(t)=f1(t)+f2(t)cos(50t)的值并画出波形;而一可实现的实际系统的H(w)为
H(w)=
用freqs函数画出H(w)的幅度和相位曲线。
用lsim函数求出信号f1(t)和、f2(t)通过系统H(w)的零状态响应y1(t)和y2(t),并根据理论知识解释所得的结果。
二、试验环境
计算机、MATLAB软件
三、实验目的
信号与系统课程设计是在学习完基础上设立的一门实践课。
通过本次设计让学生学会使用MATLAB软件,同时学会建立m文件,并且在m文件下建立MATLAB函数和语句以及对其进行调试、运行、函数图像处理保存。
更深层次的学习信号与系统以及复变函数与积分变换等课程,掌握信号的的各种处理及计算,同时提高学生的综合分析问题和解决问题的能力,为后续课程和相关专业研究奠定基础
四、课程设计原理及说明
1、信号的MATLAB表示
(1)用MATLAB实现函数f(t)=Sa(t),并绘制f(t)波形。
抽样函数Sa(t)在MATLAB中用sinc函数表示,定义为sinc(t)=sin(pt)/pt
其调用形式为y=sinc(t)
clearall;
t=-5:
0.1:
5;
f=sinc(t);%产生抽样函数
plot(t,f);grid;xlabel('t');ylabel('Sa(t)')
(2)正弦函数f(t)=Ksin(wt+a)
正弦信号和余弦信号在MATLAB中分别用sin和cos表示,其调用形式为
clearall;
t=-8:
0.01:
8;
k=2;w=1;a=pi/6;
f=k*sin(w*t+a);
plot(t,f);grid;xlabel('t');ylabel('f(t)')
(3)单边指数函数f(t)=
指数信号Keat在MATLAB中可用exp表示,其调用形式为yK*exp(a*t)
t=0:
0.001:
10;k=1;a=2;
f=k*exp(-a*t);
plot(t,f);grid;xlabel('t');ylabel('f(t)')
(4)单位冲激信号
clearall;
t0=0;tf=5;dt=0.01;t1=1;
t=[t0:
dt:
tf];
st=length(t);
n1=floor((t1-t0)/dt);
x1=zeros(1,st);x1(n1)=1/dt;
stairs(t,x1),gridon;axis([0,5,0,22])
(5)单位阶跃信号
利用MATLAB的可视化功能,在work目录下创建函数Heaviside(),可方便地绘出单位阶跃信号波形。
需要创建的Heaviside()函数代码如下:
functionf=Heaviside(t)
f=(t>0);
调用该函数,并运用plot命令。
%利用plot命令绘制单位阶跃信号波形
clearall;
t=-1:
0.01:
3;
f=Heaviside(t);
plot(t,f);
axis([-1,3,-0.2,1.2]);
set(gcf,'color','w');
title('单位阶跃信号')
(6)矩形脉冲信号
矩形脉冲信号在MATLAB中用rectpuls()函数表示,其调用形式为
yrectpuls(t,width)
用以产生一个幅度为1,宽度为width以t为对称的矩形波。
clearall;
t=0:
0.001:
4;T=1;
ft=rectpuls(t-2*T,2*T);%产生一个以t=t-2T为中心,宽度为2T的矩形脉冲
plot(t,ft);gridon;axis([04-0.51.5])
(7)三角波脉冲信号
三角波脉冲信号在MATLAB中用tripuls()函数表示,其调用形式为
ytripuls(t,width,skew)
用以产生一个最大幅度为1,宽度为width的三角波。
函数值的非零范围为
(-width/2,width/2)。
t=-3:
0.001:
3;%t决定函数的横坐标范围
ft=tripuls(t,4,0.5);%产生一个最大幅度为1,宽度为4,斜度为0.5的三角波
plot(t,ft);gridon;
axis([-33-0.51.5])
(7)用MATLAB实现函数x(t)e0.1tsin(2/3t),并绘制x(t)波形。
clearall;
t=0:
0.1:
30;
x=exp(-0.1*t).*sin(2/3*t);
axis([030-11]);plot(t,x);grid;xlabel('time(second)');ylabel('x(t)')
2、连续时间系统零状态响应的求解
MATLAB工具箱提供了一个用于求解零初始
条件下微分方程数值解的函数lsim(),其调用形式为ylsim(sys,f,t)
该调用格式对向量t定义的时间范围内,绘制LTI系统的时域波形,同时绘制出系统的激励信号对应的时域波形。
其中,t表示计算系统响应的抽样点向量,f是系统输入信号向量,sys是LTI系统模型。
在求解微分方程时,LTI的模型sys要借助MATLAB中的tf()函数来获得,
其调用形式为systf(b,a)
其中,b和a分别为微分方程右端和左端各项系数的向量。
例如,对于下列微分方程:
解:
(1)
ts=0;te=5;dt=0.01;
num=[1];den=[12100];
sys=tf(num,den);%调用LTI系统模型的函数
t=ts:
dt:
te;
f=10*sin(2*pi*t);
y=lsim(sys,f,t);%求零初始条件微分方程数值解
plot(t,y);
xlabel('t(sec)');ylabel('y(t)');grid
(2)
clearall;
ts=0;te=5;dt=0.01;
num=[10];den=[12100];
sys=tf(num,den);
t=ts:
dt:
te;
y=impulse(sys,t);%求解冲激响应的函数
plot(t,y);grid;xlabel('time(sec)');ylabel('h(t)')
3、求解冲激响应可用MATLAB工具箱提供的impulse()函数,求解阶跃响应可
用step()函数。
其调用形式分别为
y=step(sys,t)
y=impulse(sys,t)
其中,t表示计算系统响应的抽样点向量,sys是LTI系统模型。
4、连续时间系统的零、极点分析与MATLAB实现
利用MATLAB绘制连续系统零、极点分布图
如果连续系统的系统函数已知,则可利用MATLAB提供的函数tf、pole、zero、pzmap可以方便地求出系统函数的零、极点,并绘出其零、极点分布图。
首先,根据系统函数分子和分母多项式的系数,调用tf函数生成系统函数对象(tf对象)。
tf调用格式为systf(num,den)
其中输入参量num为系统函数分子多项式系数构成的行向量,den为系统函数分母多项式系数构成的行向量。
输出参量sys为MATLAB定义的系统函数对象。
其次,生成函数对象sys后,即可调用pole函数、zero函数、pzmap函数求出系统的零、极点,并绘制其零极点分布图。
(a)pole()函数
pole函数用于计算系统函数的极点,调用格式为p=pole(sys)
输出参量p为返回包含系统函数所有极点位置的列向量。
(b)zero()函数
zero函数用于计算系统函数的零点,调用格式为z=zero(sys)
输出参量z为返回包含系统函数所有零点位置的列向量。
(c)pzmap()函数
pzmap函数用于绘制系统函数的零、极点分布图和计算系统函数的零、极点位置。
调用格式为
(1)pzmap(sys)
调用该命令直接绘制出系统函数的零、极点分布图。
(2)[p,z]=pzmap(sys)
输出参量p,z为返回包含系统函数所有极点、零点位置的列向量。
调用该命令并不绘制系统函数的零、极点分布图。
num=[251];
den=[123];
sys=tf(num,den)
利用MATLAB求出该系统的零、极点,并画出零、极点分布图。
程序如下:
clearall;
a=[10-4];
b=[12-321];
h=tf(a,b);
p=pole(h)%计算系统函数极点位置列向量
z=zero(h)%计算系统函数零点位置列向量
pzmap(h)
5、利用MATLAB工具箱中提供的fourier()函数可直接求出其傅立叶变换。
其调用格式为:
F=fourier(f);
求f(t)=e-2|t|的傅立叶变换
解:
利用MATLAB提供的fourier()函数求解,命令如下:
symst;
f=exp(-2*abs(t));
F=fourier(f)
6、当非周期序列写成下列有理多项式的形式
可以用MATLAB中提供的freqz()函数来计算上面非周期序列的傅立叶变换值,调用格式为h=freqz(b,a,w)
上式中,b和a分别是其分子多项式和分母多项式的系数向量,即
w为抽样的频率点,h为傅立叶变换在抽样点w上的值。
例如利用MATLAB画出a=-0.9时,
程序如下:
clearall;
b=[1];a=[1-0.9];
w=linspace(0,2*pi,512);%线性均匀分0-2p的间隔,共512点
h1=freqz(b,a,w);
plot(w/pi,abs(h1),'k:
');
xlabel('\omega/\pi');
legend('\alpha=0.9','\alpha=-0.9')
五、课程设计的过程及调试
1、1f1(t)+f2(t)的程序及波形.
1、11f1(t)+f2(t)的调试程序
clearall;
t=-1:
0.01:
3;
f1=Heaviside(t)-Heaviside(t-2);f2=cos(2*pi*t);
plot(t,f1+f2);
axis([-1,2,-0.5,1.2]);
set(gcf,'color','w');
title('f1+f2')
1.12f1(t)+f2(t)调试波形
1、2f1(t)×f2(t)的程序及波形
1、21f1(t)×f2(t)的调试程序
clearall;
t=-1:
0.01:
3;
f1=Heaviside(t)-Heaviside(t-2);f2=cos(2*pi*t);
plot(t,f1.*f2);
axis([-1,3,-0.2,1.2]);
set(gcf,'color','w');
title('f1*f2')
1、22f1(t)×f2(t)调试波形
2、1函数的冲激响应程序及波形
2.11函数的冲激响应调试程序
clearall;
ts=0,te=5;dt=0.01;
num=[1];
den=[218];
sys=tf(num,den);
t=ts:
dt:
te;
y=impulse(sys,t);
plot(t,y);grid;xlabel('time(sec)');ylabel('s(t))
2.12函数的冲激响应的调试波形
2、2函数的阶跃响应及波形
2.21函数的阶跃响应的调试程序
clearall;
ts=0,te=15;dt=0.01;
num=[1];
den=[218];
sys=tf(num,den);
t=ts:
dt:
te;
y=step(sys,t);
plot(t,y);grid;xlabel('t(sec)');ylabel('y(t)');grid;
title('½×Ô¾ÏìÓ¦')
2.22函数的阶跃响应的调试波形
3、1正弦函数f(t)=sin(t)信号的傅立叶变换式
3.11f(t)=sin(t)的调试程序
symst;
f=sin(t);
F=fourier(f)
3.12f(t)=sin(t)的调试结果
F=i*pi*(-dirac(w-1)+dirac(w+1))
3、2取样函数f(t)=Sa(t)信号的傅立叶变换式
3.21取样函数f(t)=Sa(t)信号的调试程序
clearall;
t=-5:
0.1:
5;
symst;
f=sinc(t);
F=fourier(f);
grid;xlabel('t');ylabel('Sa(t)')
F=fourier(f)
3.22取样函数f(t)=Sa(t)的调试波形
F=heaviside(w+pi)-heaviside(w-pi)
3、3门函数
信号的傅立叶变换式
3.31门函数
信号的傅立叶变换式的调试程序
symst;
f=sym('Heaviside(t+1)')-sym('Heaviside(t-1)')
F=fourier(f)
3.32门函数
信号的傅立叶变换式的调试波形
F=2/w*sin(w)
3、4双边指数函数
信号的傅立叶变换式
3、41双边指数函数
信号的傅立叶变换式的程序
Symst;
f=exp(-2*abs(t));
F=fourier(f)
3、42双边指数函数
信号的傅立叶变换式的调试波形
F=4/(4+w^2)
4连续系统的的系统函数H(s)=(s2+3s+2)/(8s4+2s3+3s2+5)零极点
4、1系统函数H(s)=(s2+3s+2)/(8s4+2s3+3s2+5)零极点的调试程序
clearall
a=[132];
b=[82305];
h=tf(a,b);
p=pole(h)
z=zero(h)
pzmap(h)
4、2系统函数H(s)=(s2+3s+2)/(8s4+2s3+3s2+5)零极点的调试波形
5、1第一个稳定的LIT系统的频率响应的幅值和相位特性曲线
5.11系统的频率响应的幅值和相位的调试程序
clearall;
b=[3];
a=[13];
w=linspace(0,2*pi,512);
hl=freqz(b,a,w);
plot(w/pi,abs(hl),'k:
');
xlabel('\omega/\pi');
legend;
5.12系统的频率响应的幅值和相位的调试波形
5、2第二个稳定的LIT系统的频率响应的幅值和相位特性曲线
5.21系统的频率响应的幅值和相位的调试程序
clearall;
b=[501];
a=[143];
w=linspace(0,2*pi,512);
hl=freqz(b,a,w);
plot(w/pi,abs(hl),'k:
');
xlabel('\omega/\pi');
legend;
5.22系统的频率响应的幅值和相位的调试波形
6、1f(t)=f1(t)+f2(t)cos(50t)的值及波形
6.11f(t)=f1(t)+f2(t)cos(50t)的值的调试程序
clearall;
t=0:
0.005:
2.5;
f1=Heaviside(t)-Heaviside(t-1);
f2=tripuls(t-1,2,1);
f3=f2.*cos(50*t);
f=f1+f3;
plot(t,f);
gridon;
axis([-3,3,-0.5,2.5]);
set(gcf,'color','w')
6.12f(t)=f1(t)+f2(t)cos(50t)的值的调试波形
;
6、2函数信号f1(t)通过系统H(w)的零状态响应y1(t)
6.21函数信号f1(t)通过系统H(w)的零状态响应y1(t)的调试程序
clearall;
ts=0;te+5;dt=0.01;
num=[10000];
den=[126.131341.422613.1];
sys=tf(num,den);
t=ts:
dt:
te;
f1=Heaviside(t)-Heaviside(t-1);
y=lsim(sys,f,t);
plot(t,y);
xlabel('t(sec)');ylabel('y(t)');grid
set(gcf,'color','w');
6.22函数信号f1(t)通过系统H(w)的零状态响应y1(t)的调试波形
6、3函数信号f2(t)通过系统H(w)的零状态响应y2(t)
6、31函数信号f2(t)通过系统H(w)的零状态响应y2(t)的调试程序
clearall;
ts=0;te+5;dt=0.01;
num=[10000];
den=[126.131341.422613.1];
sys=tf(num,den);
t=ts:
dt:
te;
f=tripuls(t-1,2,1).*cos(50*t);
y=lsim(sys,f,t);
plot(t,y);
xlabel('t(sec)');ylabel('y(t)');grid
set(gcf,'color','w');
6、32函数信号f2(t)通过系统H(w)的零状态响应y2(t)的调
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 信号 系统 课程设计 报告