连续时间信号和系统时域分析报告及MATLAB实现

连续时间信号和系统时域分析报告及MATLAB实现
连续时间信号和系统时域分析报告及MATLAB实现

MATLAB课程设计任务书

:王** 学号:2010******

题目:

连续时间信号和系统时域分析及MATLAB实现

初始条件:

MATLAB 7.5.0 ,Windows XP系统

实验任务:

一、用MATLAB实现常用连续时间信号的时域波形(通过改变参数,分析其时域特性)。

1、单位阶跃信号,

2、单位冲激信号,

3、正弦信号,

4、实指数信号,

5、虚指数信号,

6、复指数信号。

二、用MATLAB实现信号的时域运算

1、相加,

2、相乘,

3、数乘,

4、微分,

5、积分

三、用MATLAB实现信号的时域变换(参数变化,分析波形变化)

1、反转,

2、使移(超时,延时),

3、展缩,

4、倒相,

5、综合变化

四、用MATLAB实现信号简单的时域分解

1、信号的交直流分解,

2、信号的奇偶分解

五、用MATLAB实现连续时间系统的卷积积分的仿真波形

给出几个典型例子,对每个例子,要求画出对应波形。

六、用MATLAB实现连续时间系统的冲激响应、阶跃响应的仿真波形。

给出几个典型例子,四种调用格式。

七、利用MATLAB实现连续时间系统对正弦信号、实指数信号的零状态响应的仿真波形。

目录

1 MATLAB简介 (1)

1.1 MATLAB设计目的 (1)

1.2 MATLAB语言特点 (1)

2常用连续时间信号的时域波形 (1)

2.1单位阶跃信号 (1)

2.2单位冲激信号 (2)

2.3正弦信号 (3)

2.4实指数信号 (4)

2.5虚指数信号 (5)

2.6复指数信号 (6)

3 连续时间信号的时域运算 (7)

3.1相加 (7)

3.2相乘 (8)

3.3数乘 (9)

3.4微分 (10)

3.5积分 (11)

4.1反转 (12)

4.2时移 (13)

4.3展缩 (14)

4.4倒相 (15)

4.5综合变化 (16)

5连续时间信号简单的时域分解 (17)

5.1信号的交直流分解 (17)

5.2信号的奇偶分解 (19)

6连续时间系统的卷积积分的仿真波形 (20)

7连续时间系统的冲激响应、阶跃响应的仿真波形 (23)

7.1 IMPULSE()函数 (24)

7.2 STEP()函数 (27)

8连续时间系统对正弦信号、实指数信号的零状态响应的仿真波形 (31)

8.1 正弦信号的零状态响应 (31)

8.2 实指数信号的零状态响应 (32)

9小结 (34)

1 MATLAB简介

1.1 MATLAB设计目的

深入研究连续时间信号和系统时域分析的理论知识。利用MATLAB强大的图形处理功能、符号运算功能以及数值计算功能,实现连续时间信号和系统时域分析的仿真波形

1.2 MATLAB语言特点

MATLAB和Mathematica、Maple并称为三大数学软件。它在数学类科技应用软件中在数值计算方面首屈一指。MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的程序等,主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。MATLAB的最重要特征使他拥有解决特定应用问题的程序组,也就是TOOLBOX(工具箱),如信号处理工具箱,控制系统工具箱、神经网络工具箱、模糊逻辑工具箱、通信工具箱和数据采集工具箱等许多专用工具箱,对大多数用户来说,要想灵活、高效地运用这些工具箱,通常都需要学习相应的专业知识。

2常用连续时间信号的时域波形

连续信号又称为模拟信号,其信号存在于整个时间围,包括单位阶跃信号,单位冲激信号,正弦信号,实指数信号,虚指数信号,复指数信号。

2.1单位阶跃信号

单位阶跃信号的定义如下:

单位阶跃信号的MATLAB程序:

>> t=-0.5:0.01:5;

>> t0=1.0;

>> q=stepfun(t,t0);

>> plot(t,q);

>> axis equal

信号与系统matlab实验及答案

产生离散衰减正弦序列()π0.8sin 4n x n n ?? = ??? , 010n ≤≤,并画出其波形图。 n=0:10; x=sin(pi/4*n).*0.8.^n; stem(n,x);xlabel( 'n' );ylabel( 'x(n)' ); 用MATLAB 生成信号()0sinc at t -, a 和0t 都是实数,410t -<<,画波形图。观察并分析a 和0t 的变化对波形的影响。 t=linspace(-4,7); a=1;

t0=2; y=sinc(a*t-t0); plot(t,y); t=linspace(-4,7); a=2; t0=2; y=sinc(a*t-t0); plot(t,y);

t=linspace(-4,7); a=1; t0=2; y=sinc(a*t-t0); plot(t,y);

三组对比可得a 越大最大值越小,t0越大图像对称轴越往右移 某频率为f 的正弦波可表示为()()cos 2πa x t ft =,对其进行等间隔抽样,得到的离散样值序列可表示为()()a t nT x n x t ==,其中T 称为抽样间隔,代表相邻样值间的时间间隔,1 s f T = 表示抽样频率,即单位时间内抽取样值的个数。抽样频率取40 Hz s f =,信号频率f 分别取5Hz, 10Hz, 20Hz 和30Hz 。请在同一张图中同时画出连续信号()a x t t 和序列()x n nT 的波形图,并观察和对比分析样值序列的变化。可能用到的函数为plot, stem, hold on 。 fs = 40; t = 0 : 1/fs : 1 ; % ?μ?ê·?±e?a5Hz,10Hz,20Hz,30Hz f1=5; xa = cos(2*pi*f1*t) ; subplot(1, 2, 1) ;

北京交通大学信号与系统时域分析

北京交通大学信号与系统时域分析

【研讨题目2】 信号与系统时域分析专题研讨 【目的】 1.研究用离散方法近似计算连续信号的卷积积分; 2.通过分析近似计算卷积积分过程中出现的问题,锻炼学生分析问题和解决问题的能力; 【知识点】 信号时域分析,卷积积分,卷积和 【研讨题目】连续信号卷积积分的数值近似计算 两个连续信号的卷积积分定义为 τττd )()()(-= ?∞ ∞ -t h x t y 为了能用数值方法进行计算,需对连续信号进行抽样。记x [k ]=x (k ), h [k ]=h (k ),为 进行数值计算所选定的抽样间隔,可以证明连续信号卷积积分可近似的表示为 (Δ)Δ([][])y k x k h k ≈?* (1) 由式(1)可知,可以利用Matlab 提供的conv 函数近似计算连续信号的卷积积分。 一、(*)理论分析 为了对近似计算的结果进行分析,用解析的方法计算下列卷积积分,推出卷积积分的解析表达式; (1) 时限信号卷积积分 x 1(t )=u (t )-u (t -1),y 1(t )=x 1(t )*x 1(t ); 卷积结果为:y1(t)= x 1(t )*x 1(t )=r(t)-2*r(t-1)+r(t-2) (2) 分段常数信号卷积积分 x 2(t )= x 1(t )+2 x 1(t -1)+ x 1(t -2),h 2(t )= x 1(t )- x 1(t -1), y 2(t )=x 2(t )*h 2(t ); 卷积结果为:y2(t)= x 2(t )*h 2(t ) =y1(t)+y1(t-1)-y1(t-2)-y1(t-3) =r(t)-r(t-1)-2*r(t-2)+2*r(t-3)+r(t-4)-r(t-5) (3) 非时限信号卷积积分 x 3(t )=u (t ),h 3(t )=e -t u (t ), y 3(t )=x 3(t )*h 3(t ) 卷积结果为:y3= x 3(t )*h 3(t ) =[1-exp(-t)]*u(t) 二、(*)时限信号卷积积分的近似计算 取不同的△值,用Matlab 函数conv 近似计算卷积积分y 1(t )并画出其波形,讨论的取值对计算结果的影响。

信号与系统——MATLAB基本实验

《信号与系统MATLAB实践》第一次上机作业 实验一、熟悉MATLAB基本操作 三、基本序列运算 1.数组的加减乘除和乘方运算 A=[1 2 3]; B=[4 5 6]; C=A+B; D=A-B; E=A.*B; F=A./B; G=A.^B; subplot(2,4,1);stem(A) subplot(2,4,2);stem(B) subplot(2,4,3);stem(C) subplot(2,4,4);stem(D) subplot(2,4,5);stem(E) subplot(2,4,6);stem(F) subplot(2,4,7);stem(G) 2.绘制函数波形 (1)t=0:0.001:10

x=3-exp(-t); plot(t,x) ylabel('f(t)'); xlabel('t'); title('(1)'); (2)t=0:0.001:10 x=5*exp(-t)+3*exp(-2*t); plot(t,x) ylabel('f(t)'); xlabel('t'); title('(2)');

(3)t=0:0.001:3 x=exp(-t).*sin(2*pi*t); plot(t,x) ylabel('f(t)'); xlabel('t'); title('(3)'); (4)t=0:0.001:3 x=sin(3*t)./(3*t);

plot(t,x) ylabel('f(t)'); xlabel('t'); title('(4)'); (5)k=1:1:6 x=(-2).^(-k); stem(k) xlabel('k'); ylabel('f(k)'); title('(5)');

信号与系统MATLAB实验报告

《信号与系统》MATLAB实验报告 院系:专业: 年级:班号: 姓名:学号: 实验时间: 实验地点:

实验一 连续时间信号的表示及可视化 实验题目: )()(t t f δ=;)()(t t f ε=;at e t f =)((分别取00<>a a 及); )()(t R t f =;)()(t Sa t f ω=;)2()(ft Sin t f π=(分别画出不同周期个数 的波形)。 解题分析: 以上各类连续函数,先运用t = t1: p:t2的命令定义时间范围向量,然后调用对应的函数,建立f 与t 的关系,最后调用plot ()函数绘制图像,并用axis ()函数限制其坐标范围。 实验程序: (1) )()(t t f δ= t=-1:0.01:3 %设定时间变量t 的范围及步长 f=dirac(t) %调用冲激函数dirac () plot(t,f) %用plot 函数绘制连续函数 axis([-1,3,-0.5,1.5]) %用axis 函数规定横纵坐标的范围 (2) )()(t t f ε= t=-1:0.01:3 %设定时间变量t 的范围及步长 f=heaviside(t) %调用阶跃函数heaviside () plot(t,f) %用plot 函数绘制连续函数 title('f(t)=heaviside(t)') %用title 函数设置图形的名称 axis([-1,3,-0.5,1.5]) %用axis 函数规定横纵坐标的范围 (3) at e t f =)( a=1时: t=-5:0.01:5 %设定时间变量t 的范围及步长 f=exp(t) %调用指数函数exp ()

离散信号与系统时域分析

目录 第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) 绘制结果图形并对仿真结果进行详细的分析。

信号与系统 matlab答案

M2-3 (1) function yt=x(t) yt=(t).*(t>=0&t<=2)+2*(t>=2&t<=3)-1*(t>=3&t<=5); (2)function yt=x (t) yt=(t).*(t>=0&t<=2)+2*(t>=2&t<=3)-1*(t>=3&t<=5); t=0:0.001:6; subplot(3,1,1) plot(t,x2_3(t)) title('x(t)') axis([0,6,-2,3]) subplot(3,1,2) plot(t,x2_3(0.5*t)) title('x(0.5t)') axis([0,11,-2,3]) subplot(3,1,3) plot(t,x2_3(2-0.5*t)) title('x(2-0.5t)') axis([-6,5,-2,3]) 图像为:

M2-5 (3) function y=un(k) y=(k>=0) untiled3.m k=[-2:10] xk=10*(0.5).^k.*un(k); stem(k,xk) title('x[k]') axis([-3,12,0,11])

M2-5 (6) k=[-10:10] xk=5*(0.8).^k.*cos((0.9)*pi*k) stem(k,xk) title('x[k]') grid on M2-7 A=1; t=-5:0.001:5; w0=6*pi; xt=A*cos(w0*t); plot(t,xt) hold on A=1; k=-5:5; w0=6*pi; xk=A*cos(w0*0.1*k); stem(k,xk) axis([-5.5,5.5,-1.2,1.2]) title('x1=cos(6*pi*t)&x1[k]')

信号与系统实验报告实验一 信号与系统的时域分析

实验一信号与系统的时域分析 一、实验目的 1、熟悉与掌握常用的用于信号与系统时域仿真分析的MA TLAB函数; 2、掌握连续时间与离散时间信号的MA TLAB产生,掌握用周期延拓的方法将一个非周期信号进行周期信号延拓形成一个周期信号的MA TLAB编程; 3、牢固掌握系统的单位冲激响应的概念,掌握LTI系统的卷积表达式及其物理意义,掌握卷积的计算方法、卷积的基本性质; 4、掌握利用MA TLAB计算卷积的编程方法,并利用所编写的MA TLAB程序验证卷积的常用基本性质; 掌握MA TLAB描述LTI系统的常用方法及有关函数,并学会利用MATLAB求解LTI系统响应,绘制相应曲线。 基本要求:掌握用MA TLAB描述连续时间信号与离散时间信号的方法,能够编写MATLAB程序,实现各种信号的时域变换与运算,并且以图形的方式再现各种信号的波形。掌握线性时不变连续系统的时域数学模型用MA TLAB描述的方法,掌握卷积运算、线性常系数微分方程的求解编程。 二、实验原理 信号(Signal)一般都就是随某一个或某几个独立变量的变化而变化的,例如,温度、压力、 声音,还有股票市场的日收盘指数等,这些信号都就是随时间的变化而变化的,还有一些信号,例如在研究地球结构时,地下某处的密度就就是随着海拔高度的变化而变化的。一幅图片中的每一个象素点的位置取决于两个坐标轴,即横轴与纵轴,因此,图像信号具有两个或两个以上的独立变量。 在《信号与系统》课程中,我们只关注这种只有一个独立变量(Independent variable)的信号,并且把这个独立变量统称为时间变量(Time variable),不管这个独立变量就是否就是时间变量。 在自然界中,大多数信号的时间变量都就是连续变化的,因此这种信号被称为连续时间信号(Continuous-Time Signals)或模拟信号(Analog Signals),例如前面提到的温度、压力与声音 信号就就是连续时间信号的例子。但就是,还有一些信号的独立时间变量就是离散变化的,这种信号称为离散时间信号。前面提到的股票市场的日收盘指数,由于相邻两个交易日的日收盘指数相隔24小时,这意味着日收盘指数的时间变量就是不连续的,因此日收盘指数就是离散时间信号。 而系统则用于对信号进行运算或处理,或者从信号中提取有用的信息,或者滤出信号中某些无用的成分,如滤波,从而产生人们所希望的新的信号。系统通常就是由若干部件或单元组成的一个整体(Entity)。系统可分为很多不同的类型,例如,根据系统所处理的信号的不同,系统可分为连续时间系统(Continuous-time system)与离散时间系统(Discrete-time system),根

matlab信号与系统代码整理

连续时间系统 (1) 离散时间系统 (2) 拉普拉斯变换 (4) Z变换 (5) 傅里叶 (7) 连续时间系统 %%%%%%%%%%向量法%%%%%%%%%%%%%%%% t1=-2:0.01:5; f1=4*sin(2*pi*t1-pi/4); figure(1) subplot(2,2,1),plot(t1,f1),grid on %%%%%%%%%符号运算法%%%%%%%%%%%% syms t f1=sym('4*sin(2*pi*t-pi/4)'); figure(2) subplot(2,2,1),ezplot(f1,[-2 5])跟plot相比,ezplot不用指定t,自动生成。axis([-5,5,-0.1,1])控制坐标轴的范围xx,yy; 求一个函数的各种响应 Y’’(t)+4y’(t)+2y(t)=f”(t)+3f(t) %P187 第一题 %(2) clear all; a1=[1 4 2]; b1=[1 0 3]; [A1,B1,C1,D1]=tf2ss(b1,a1); t1=0:0.01:10; x1=exp(-t1).*Heaviside(t1); rc1=[2 1];(起始条件) figure(1) subplot(3,1,1),initial(A1,B1,C1,D1,rc1,t1);title('零输入响应') subplot(3,1,2),lsim(A1,B1,C1,D1,x1,t1);title('零状态响应') subplot(3,1,3),lsim(A1,B1,C1,D1,x1,t1,rc1);title('全响应') Y=lsim(A1,B1,C1,D1,x1,t1,rc1);title('全响应')则是输出数值解 subplot(2,1,1),impulse(b1,a1,t1:t:t2可加),grid on,title('冲激响应') subplot(2,1,2),step(b1,a1,t1:t:t2可加),grid on,title('阶跃响应') 卷积 %第九题 P189 clear all; %(1) t1=-1:0.01:3;

(完整word版)信号与系统matlab实验

习题三 绘制典型信号及其频谱图 1.更改参数,调试程序,绘制单边指数信号的波形图和频谱图。观察参数a对信号波形 及其频谱的影响。 程序代码: close all; E=1;a=1; t=0:0.01:4; w=-30:0.01:30; f=E*exp(-a*t); F=1./(a+j*w); plot(t,f);xlabel('t');ylabel('f(t)'); figure; plot(w,abs(F));xlabel('\omega');ylabel('|F(\omega)|';

E=1,a=1,波形图频谱图更改参数E=2,a=1;

更改参数a,对信号波形及其频谱的影响。(保持E=2)上图为a=1图像 a=2时

a=4时 随着a的增大,f(t)曲线变得越来越陡,更快的逼近0,而对于频谱图,随着a增大,图像渐渐向两边张开,峰值减小,陡度减小,图像整体变得更加平缓。 2.矩形脉冲信号 程序代码: close all; E=1;tao=1; t=-4:0.1:4; w=-30:0.1:30;

f=E*(t>-tao/2&tao/2)+0*(t<=-tao/2&t>=tao/2); F=(2*E./w).*sin(w*tao/2); plot(t,f);xlabel('t');ylabel('f(t)'); figure; plot(w,abs(F));xlabel('\omega');ylabel('|F(\omega)|') ; figure; plot(w,20*log10(abs(F))); xlabel('\omega');ylabel('|F(\omega)| in dB'); figure; plot(w,angle(F));xlabel('\omega');ylabel('\phi(\omega )');

信号与系统MATLAB实验

《信号与系统及MATLAB实现》实验指导书

前言 长期以来,《信号与系统》课程一直采用单一理论教学方式,同学们依靠做习题来巩固和理解教学内容,虽然手工演算训练了计算能力和思维方法,但是由于本课程数学公式推导较多,概念抽象,常需画各种波形,作题时难免花费很多时间,现在,我们给同学们介绍一种国际上公认的优秀科技应用软件MATLAB,借助它我们可以在电脑上轻松地完成许多习题的演算和波形的绘制。 MATLAB的功能非常强大,我们此处仅用到它的一部分,在后续课程中我们还会用到它,在未来地科学研究和工程设计中有可能继续用它,所以有兴趣的同学,可以对MATLAB 再多了解一些。 MATLAB究竟有那些特点呢? 1.高效的数值计算和符号计算功能,使我们从繁杂的数学运算分析中解脱出来; 2.完备的图形处理功能,实现计算结果和编程的可视化; 3.友好的用户界面及接近数学表达式的自然化语言,易于学习和掌握; 4.功能丰富的应用工具箱,为我们提供了大量方便实用的处理工具; MATLAB的这些特点,深受大家欢迎,由于个人电脑地普及,目前许多学校已将它做为本科生必须掌握的一种软件。正是基于这些背景,我们编写了这本《信号与系统及MATLAB实现》指导书,内容包括信号的MATLAB表示、基本运算、系统的时域分析、频域分析、S域分析、状态变量分析等。通过这些练习,同学们在学习《信号与系统》的同时,掌握MATLAB的基本应用,学会应用MATLAB的数值计算和符号计算功能,摆脱烦琐的数学运算,从而更注重于信号与系统的基本分析方法和应用的理解与思考,将课程的重点、

难点及部分习题用MATLAB进行形象、直观的可视化计算机模拟与仿真实现,加深对信号与系统的基本原理、方法及应用的理解,为学习后续课程打好基础。另外同学们在进行实验时,最好事先预习一些MATLAB的有关知识,以便更好地完成实验,同时实验中也可利用MATLAB的help命令了解具体语句以及指令的使用方法。 实验一基本信号在MATLAB中的表示和运算 一、实验目的 1.学会用MATLAB表示常用连续信号的方法; 2.学会用MATLAB进行信号基本运算的方法; 二、实验原理 1.连续信号的MATLAB表示 MATLAB提供了大量的生成基本信号的函数,例如指数信号、正余弦信号。 表示连续时间信号有两种方法,一是数值法,二是符号法。数值法是定义某一时间范围和取样时间间隔,然后调用该函数计算这些点的函数值,得到两组数值矢量,可用绘图语句画出其波形;符号法是利用MATLAB的符号运算功能,需定义符号变量和符号函数,运算结果是符号表达的解析式,也可用绘图语句画出其波形图。 例1-1指数信号指数信号在MATLAB中用exp函数表示。 如at )(,调用格式为ft=A*exp(a*t) 程序是 f t Ae

信号与系统MATLAB实验报告

实验报告 实验课程:信号与系统—Matlab综合实验学生姓名: 学号: 专业班级: 2012年5月20日

基本编程与simulink仿真实验 1—1编写函数(function)∑=m n k n 1并调用地址求和∑∑∑===++100 11-8015012 n n n n n n 。实验程序: Function sum=qiuhe(m,k)Sum=0For i=1:m Sum=sum+i^k End 实验结果; qiuhe(50,2)+qiuhe(80,1)+qiuhe(100,-1) ans=4.6170e+004。 1-2试利用两种方式求解微分方程响应 (1)用simulink对下列微分方程进行系统仿真并得到输出波形。(2)编程求解(转移函数tf)利用plot函数画图,比较simulink图和plot图。)()(4)(6)(5)(d 22t e t e d d t r t r d d t r d t t t +=++在e(t)分别取u(t)、S(t)和sin(20пt)时的情况! 试验过程 (1)

(2) a=[1,5,6]; b=[4,1]; sys=tf(b,a); t=[0:0.1:10]; step(sys)

连续时间系统的时域分析3-1、已知某系统的微分方程:)()()()()(d 2t e t e d t r t r d t r t t t +=++分别用两种方法计算其冲激响应和阶跃响应,对比理论结果进行验证。 实验程序: a=[1,1,1];b=[1,1];sys=tf(b,a);t=[0:0.01:10];figure;subplot(2,2,1);step(sys);subplot(2,2,2);x_step=zeros(size(t));x_step(t>0)=1;x_step(t==0)=1/2;lsim(sys,x_step,t);subplot(2,2,3);impulse(sys,t);title('Impulse Response');xlabel('Time(sec)');ylabel('Amplitude');subplot(2,2,4);x_delta=zeros(size(t));x_delta(t==0)=100;[y1,t]=lsim(sys,x_delta,t);y2=y1;plot(t,y2);title('Impulse Response');

信号与系统作业

实验一常用连续时间信号的实现 一、实验目的 (1)了解连续时间信号的特点; (2)掌握连续时间信号表示的向量法和符号法; (3)熟悉MA TLABPlot函数等的应用。 二、涉及的MATLAB函数 1.plot函数 功能:在X轴和Y轴方向都按线性比例绘制成二维图形。 2.ezplot函数 功能:绘制符号函数在一定范围内的二维图形,简易绘制函数曲线。 3.Sym函数 功能:定义信号为符号变量。 4.subplot函数 功能:产生多个绘图区间。 三、实验内容与方法 1.正弦交流信号f(t)=sin(ωt+φ) (1)符号推理法生成正弦交流信号。 MATLAB程序:. t=-0:0.001:1; f=sym('sin(2*pi*t)'); ezplot(f,[0,1]); xlabel('时间(t)'); ylabei('幅值(f)'); title(‘正弦交流信号'); 用符号法生成的正弦交流信号如图所示:

(2)数值法生成正弦交流信号。 MATLAB程序:. t=-0:0.001:1; y=sin(2*pi*t); plot(t,y,'k'); xlabel('时间(t)'); ylabei('幅值(f)'); title('正弦交流信号'); 用数值法生成的正弦交流信号如图所示: 2.单边衰减指数信号. MATLAB程序: t1=-1;t2=10;dt=0.1; t=t1:dt:t2; A1=1; %斜率 a1=0.5; %斜率 n=A1*exp(-a1*t); plot(t,n); axis([t1,t2,0,1]); xlabel('时间(t)'); ylabel('幅值(f)'); title('单边衰减指数信号'); 用数值法生成的单边衰减指数信号如图所示:

基于MATLAB的信号与系统仿真及应用

本科毕业(论文) 题 目 (中、英文 ) in The Signal System 分类 号 学号 密级 公开 学校代码 1107044431 TN911.6 基于MATLAB 的信号系统仿真及应用 The Application of MATLAB in The Signal System 工科 作者姓名 指导教师 学科门类 专业名称 电气工程及其自动化 提交论文日期 成绩评定 二零一五年五月

摘要 当前的科学信息技术正在日新月异的高速发展,而通过应用数字信号处理的方法,已成为一个非常重要的技术手段被广泛应用在通信、音频和图像、遥感,视频等领域。为了更好地了解信号与系统的基本理论和掌握其方法,从而更好地理解和掌握数字信号处理的理论知识,因此在实验过程中我们就需要通过MATLAB 计算机辅助设计平台。 本论文主要探究MATALB在信号与系统中的连续信号和离散信号中的应用,主要从连续和离散两方面入手,进一步掌握信号系统中的相关知识。同时引进计算机软件—MATLAB,对信号系统二阶系统的时域和频域分析,通过它在计算机上对程序进行仿真,阐述信号与系统理论应用与实际相联系。以此激发学习兴趣,变被动接受为主动探知,从而提升学习效果,培养主动思维,学以致用的思维习惯,也可以让人们进一步了解MATLAB软件 关键词:采样定理;MATLAB;信号与系统;抽样定理

Abstract Current, the rapid development of science and information technology are changing and through the application of digital signal processing method, has become a very important technology is widely used in communication, audio and video, remote sensing, video, etc. In order to better understand the basic theory of signal and system, and grasp the method, to better understand and master the theoretical knowledge of digital signal processing, so we need in the process of experiment by MATLAB computer aided design platform. This thesis mainly explores MATALB in signal and system, the application of discrete and continuous signals, mainly from the two aspects of the continuous and discrete, further to master relevant knowledge of signal system. Introduction of computer software - MATAB at the same time, the signal system of second order system time domain and frequency domain analysis, through its d on program on computer simulation, signal and system theory associated with the actual application. To stimulate interest in learning, change passive accept to active detection, so as to improve learning effect, active thinking, to practice habits of thinking, also can let people learn more about MATLAB software. Key words:Sampling theorem; MATLAB; Signals and systems; The sampling theorem

实验项目五:表示信号与系统的MATLAB函数、工具箱

电子科技大学 实 验 报 告 学生姓名: 学号: 指导老师: 日期:2016年 12月25 日

一、实验室名称: 科研楼a306 二、实验项目名称: 实验项目五:表示信号与系统的MATLAB 函数、工具箱 三、实验原理: 利用MATLAB 强大的数值处理工具来实现信号的分析和处理,首先就是要学会应用MATLAB 函数来构成信号。常见的基本信号可以简要归纳如下: 1、单位抽样序列 ???=01 )(n δ 00≠=n n 在MATLAB 中可以利用zeros()函数实现。 ; 1)1();,1(==x N zeros x 如果)(n δ在时间轴上延迟了k 个单位,得到)(k n -δ即: ???=-01)(k n δ 0≠=n k n 2、单位阶跃序列 ???0 1)(n u 00<≥n n 在MATLAB 中可以利用ones()函数实现。 );,1(N ones x = 3、正弦序列 )/2sin()(?π+=Fs fn A n x 采用MATLAB 实现 )/***2sin(*1:0fai Fs n f pi A x N n +=-= 4、复正弦序列

n j e n x ?=)( 采用MATLAB 实现 )**exp(1 :0n w j x N n =-= 5、指数序列 n a n x =)( 采用MATLAB 实现 n a x N n .^1 :0=-= 四、实验目的: 目的:1、加深对常用离散信号的理解; 2、熟悉表示信号的基本MATLAB 函数。 任务:基本MATLAB 函数产生离散信号;基本信号之间的简单运算;判断信 号周期。 五、实验内容: MATLAB 仿真 实验步骤: 1、编制程序产生上述5种信号(长度可输入确定),并绘出其图形。 2、在310≤≤n 内画出下面每一个信号: 1223[]sin()cos() 44[]cos ()4 []sin()cos()48n n x n n x n n n x n πππππ=== 六、实验器材: 计算机、matlab 软件、C++软件等。 七、实验数据及结果分析: 实验1: 单位抽样序列

信号与系统课后matlab作业.

(1) t=-2:0.001:4; T=2; xt=rectpuls(t-1,T); plot(t,xt) axis([-2,4,-0.5,1.5]) 图象为: (2) t=sym('t'); y=Heaviside(t); ezplot(y,[-1,1]); grid on axis([-1 1 -0.1 1.1]) 图象为:

A=10;a=-1;B=5;b=-2; t=0:0.001:10; xt=A*exp(a*t)-B*exp(b*t); plot(t,xt) 图象为: (4) t=sym('t'); y=t*Heaviside(t); ezplot(y,[-1,3]); grid on axis([-1 3 -0.1 3.1]) 图象为:

A=2;w0=10*pi;phi=pi/6; t=0:0.001:0.5; xt=abs(A*sin(w0*t+phi)); plot(t,xt) 图象为: (6) A=1;w0=1;B=1;w1=2*pi; t=0:0.001:20; xt=A*cos(w0*t)+B*sin(w1*t); plot(t,xt) 图象为:

A=4;a=-0.5;w0=2*pi; t=0:0.001:10; xt=A*exp(a*t).*cos(w0*t); plot(t,xt) 图象为: (8) w0=30; t=-15:0.001:15; xt=cos(w0*t).*sinc(t/pi); plot(t,xt) axis([-15,15,-1.1,1.1]) 图象为:

(1)function yt=x2_3(t) yt=(t).*(t>=0&t<=2)+2*(t>=2&t<=3)-1*(t>=3&t<=5); (2)function yt=x2_3(t) yt=(t).*(t>=0&t<=2)+2*(t>=2&t<=3)-1*(t>=3&t<=5); t=0:0.001:6; subplot(3,1,1) plot(t,x2_3(t)) title('x(t)') axis([0,6,-2,3]) subplot(3,1,2) plot(t,x2_3(0.5*t)) title('x(0.5t)') axis([0,11,-2,3]) subplot(3,1,3) plot(t,x2_3(2-0.5*t)) title('x(2-0.5t)') axis([-6,5,-2,3]) 图像为:

信号与系统实验(MATLAB版) (1)

《信号与系统MATLAB实现》实验指导书 电气信息工程学院 2014年2月

长期以来,《信号与系统》课程一直采用单一理论教学方式,同学们依靠做习题来巩固和理解教学内容,虽然手工演算训练了计算能力和思维方法,但是由于本课程数学公式推导较多,概念抽象,常需画各种波形,作题时难免花费很多时间,现在,我们给同学们介绍一种国际上公认的优秀科技应用软件MA TLAB,借助它我们可以在电脑上轻松地完成许多习题的演算和波形的绘制。 MATLAB的功能非常强大,我们此处仅用到它的一部分,在后续课程中我们还会用到它,在未来地科学研究和工程设计中有可能继续用它,所以有兴趣的同学,可以对MATLAB 再多了解一些。 MATLAB究竟有那些特点呢? 1.高效的数值计算和符号计算功能,使我们从繁杂的数学运算分析中解脱出来; 2.完备的图形处理功能,实现计算结果和编程的可视化; 3.友好的用户界面及接近数学表达式的自然化语言,易于学习和掌握; 4.功能丰富的应用工具箱,为我们提供了大量方便实用的处理工具; MATLAB的这些特点,深受大家欢迎,由于个人电脑地普及,目前许多学校已将它做为本科生必须掌握的一种软件。正是基于这些背景,我们编写了这本《信号与系统及MATLAB实现》指导书,内容包括信号的MA TLAB表示、基本运算、系统的时域分析、频域分析、S域分析、状态变量分析等。通过这些练习,同学们在学习《信号与系统》的同时,掌握MATLAB的基本应用,学会应用MATLAB的数值计算和符号计算功能,摆脱烦琐的数学运算,从而更注重于信号与系统的基本分析方法和应用的理解与思考,将课程的重点、难点及部分习题用MATLAB进行形象、直观的可视化计算机模拟与仿真实现,加深对信号与系统的基本原理、方法及应用的理解,为学习后续课程打好基础。另外同学们在进行实验时,最好事先预习一些MATLAB的有关知识,以便更好地完成实验,同时实验中也可利用MATLAB的help命令了解具体语句以及指令的使用方法。

信号与系统MATLAB仿真题目

考核人数______ 考核班次_______________ 任课教员_________ 出题教员签名________ 任课教研室主任签名_______日期_______ 队别__________ 教学班次___________ 学号___________ 姓名____________ …………………………密………………………………封………………………………线……………………………………… 通信系统仿真题目 1.学习电路时已知LC 谐振电路具有选择频率的作用,当输入正弦信号频率与LC 电路的谐 振频率一致时,将产生较强的输出响应,而当输入信号频率适当偏离时,输出响应相对值很弱,几乎为零(相当于窄带通滤波器)。利用这一原理可以从非正弦周期信号中选择所需的正弦频率成分。题图所示RLC 并联电路和电流1()i t 都是理想模型。已知电路的谐振频率为 0100f kHz = =,100R k =Ω谐振电路品质因素Q 足够高(可滤除邻近频率成分) 1()i t 为周期矩形波,幅度为1 mA 当1()i t 的参数(,)T τ为下列情况时,粗略地画出输出电压 2()t υ的波形,并注明幅度值。 (1)510s T s τμμ== (2)1020s T s τμμ== (3)1530s T s τμμ== 2.设()x n 为一限长序列,当0n <和n N ≥时,()0x n =,且N 等于偶数。已知[()]DFT x n = ()X k ,试用()X k 表示以下各序列的DFT 。 (1)1()(1)x n x N n =-- (2)2()(1)()n x n x n =- (3) 3() (01)()()(21)0()x n n N x n x n N N n N n ≤≤-?? =-≤≤-??? 为其他值 (4) 4()()(01) ()2 2 () N N x n x n n x n n ?≠+≤≤ -? =???为其他值 (5) 5()(01)()0 (21)0() x n n N x n N n N n ≤≤-?? =≤≤-??? 为其他值 (6) 6() ()20()n x n x n n ??? ? ?=????? 为偶数为奇数 (DFT 有限长度取2N ,k 取偶数。) (7) 7()(2)x n x n =(DFT 有限长度取 2 N )。 3.已知三角脉冲1() f t 的傅里叶变换为21()24E F Sa τωτω??= ??? 试利用有关定理求210()cos()2f t f t t τω?? =- ??? 的傅里叶变换2()F ω。1()f t 、2()f t 的波形如下图所示。 4.求下图所示半波余弦信号的傅里叶级数。若E=10V ,f=10kHz ,大致画出幅度谱。 5.求下图所示()F ω的傅里叶逆变换()f t 。

信号与系统 信号的时域分析 及Matlab实现

信号与系统 课程实验报告 实验四 实验名称 信号的时域分析 及Matlab 实现 系 别 教师姓名 实验地点 5309 实验日期 2011-06-20 学生姓名 学号 一、实验内容 1.预习实验原理。 2.对实验内容编写程序(M 文件),上机运行。 3.记录并整理实验数据。 二、实验目的 1.掌握用Matlab 分析系统时间响应的方法。 2.掌握用Matlab 分析连续、离散系统的冲激响应的方法。 3.理解系统零、极点分布与系统稳定性关系。 三、涉及实验的相关情况介绍(包含使用软件或实验设备等情况) 计算机一台(安装MATLAB6.5版本或以上版本) 四、实验试做记录(含程序、数据记录及分析) 1.设) 2)(1()(p s p s s s H --= 设①p1=-2,p2=-30; ②p1=-2,p2=3 (1)针对极点参数①②,画出系统零、极点分布图,判断该系统稳定性。 程序: clear num=[1]; den=[1,32,60]; [r,poles,k]=residue(num,den); p=roots(den);

z=roots(num); subplot(2,2,1);plot(real(p),imag(p),'*');hold on; plot(real(z),imag(z),'o');grid on T=0:0.1:10; y1=impulse(num,den,T); subplot(2,2,2);plot(T,y1);grid;title('脉冲响应') 结果: 位于s左半平面,稳定 (2)针对极点参数①②,绘出系统的脉冲响应曲线,并观察t→∞时,脉冲响应变化趋势。 程序: clear num=[1]; den=[1,-1,-6]; [r,poles,k]=residue(num,den); p=roots(den); z=roots(num); subplot(2,2,1);plot(real(p),imag(p),'*');hold on; plot(real(z),imag(z),'o');grid on T=0:0.1:10; y1=impulse(num,den,T); subplot(2,2,2);plot(T,y1);grid;title('脉冲响应') 结果:

实验一 连续时间信号与系统的时域分析及MATLAB实现

实验一 连续时间信号与系统的时域分析及MATLAB 实现 实验目的: 1. 掌握MATLAB 相关函数的调用格式及作用; 2. 会用MATLAB 生成和实现连续时间信号,并对连续系统进行分析。 实验内容: 1. 先在MATLAB 帮助菜单下查看以下函数的调用格式及作用; Plot,ezplot,sym,subplot,stepfun,diff,int,conv,impulse,step,lsim,roots. 2.绘出下列信号的时域波形。 (1)()sin t f t t =;(2)()5sin 4f t t π??= ??? ; (3)()52t f t e -=;(4)()42j t f t e π=(实部、虚部、相角及模) 3.完成下列信号相应的运算和变换。 (1)已知信号()()()11222f t t u t u t ? ?=+?+--?? ?? ???,绘制时域波形并用MA TLAB 求()()()()2,,2,f t f t f t f t +--及其波形。 (2)已知信号()()()()144f t t u t u t =-+?--????和()2sin(2)f t t π=,用MATLAB 计算()()12f t f t +;()1d f t dt ,()1t f d ττ-∞?并绘制波形。; (3)计算信号()()()12f t u t u t =--和()()()24f t u t u t =--的卷积。 4.已知描述某连续连续系统的微分方程为: ()()()()()''''282y t y t y t f t f t ++=+ 试用MATLAB 绘出该系统的冲激响应和阶跃响应,并求出系统在()()2t f t e u t -=的激励下的零状态响应。

相关文档
最新文档