信号系统实验报告
- 格式:doc
- 大小:661.50 KB
- 文档页数:12
电气学科大类级《信号与控制综合实验》课程实验报告(基本实验一:信号与系统基本实验)姓名学号专业班号同组者1 学号专业班号同组者2 学号专业班号指导教师日期实验成绩评阅人综合实验和实验报告要求信号与控制综合实验,是集多门技术基础课程以及其它延伸课程理论于一体的综合性实验课程,需要综合多门学科理论知识和实验方法来体现,因此,实验目的不是简单的课程理论验证和练习,而是综合应用、研究开发、设计创新。
应采用尽可能好的设计,使所设计的电路和系统达到要实现的功能,步骤和方案自行拟定,实现对设计思路的实验验证。
完成多个实验项目的,应将实验内容整理综合后写成一份总报告,以利于锻炼整理归纳和总结能力,在总报告中以第二级标题形式依次写下所完成的实验项目、内容及实验设计过程。
实验报告按“题目、目录、正文(分所完成的各实验项目)、结论、心得与自我评价、参考文献”6个部分撰写;正文主要包括以下几个内容:任务和目标、总体方案设计(原理分析与方案设计特点,选择依据和确定)、方案实现和具体设计(过程)、实验设计与实验结果、结果分析与讨论。
(格式方面请注意:每个图应该有图号和图名,位于图的下方,同一图号的分图应在同一页,不要跨页;每个表应该有表号和表名,位于表的上方,表号表名与表(数据)也应在同一页,不要跨页;建议各部分题目采用四号黑体、设计报告内容文字采用小四号宋体)注:报告中涉及实验指导书或教材内容,只需注明引用位置,不必在报告中再加以阐述。
不得不加引用标记地抄袭任何资料。
每一基本实验部分按计划学时100分成绩计算(100%),需要完成60分的实验项目;实验报告、设计部分和创新研究内容另外计分(分别为10%、20%和10%)。
再按照学时比例与本课程其它部分实验综合成为总实验成绩。
每一部分实验均为:基本实验:0~60分,考核基本理论的掌握和基本操作技能、实验室道德规范;实验报告:0~10分,考核思考和总结表述能力;完成设计性实验:0~20分,评价设计能力;完成创新性实验:0~10分,鼓励创新。
实验四 时域抽样与频域抽样一、实验目的加深理解连续时间信号的离散化过程中的数学概念和物理概念,掌握时域抽样定理的基本内容。
掌握由抽样序列重建原连续信号的基本原理与实现方法,理解其工程概念。
加深理解频谱离散化过程中的数学概念和物理概念,掌握频域抽样定理的基本内容。
二、 实验原理时域抽样定理给出了连续信号抽样过程中信号不失真的约束条件:对于基带信号,信号抽样频率sam f 大于等于2倍的信号最高频率m f ,即m sam f f 2≥。
时域抽样是把连续信号x (t )变成适于数字系统处理的离散信号x [k ] ;信号重建是将离散信号x [k ]转换为连续时间信号x (t )。
非周期离散信号的频谱是连续的周期谱。
计算机在分析离散信号的频谱时,必须将其连续频谱离散化。
频域抽样定理给出了连续频谱抽样过程中信号不失真的约束条件。
三.实验内容1. 为了观察连续信号时域抽样时抽样频率对抽样过程的影响,在[0,0.1]区间上以50Hz 的抽样频率对下列3个信号分别进行抽样,试画出抽样后序列的波形,并分析产生不同波形的原因,提出改进措施。
)102cos()(1t t x ⨯=π答: 函数代码为: t0 = 0:0.001:0.1;x0 =cos(2*pi*10*t0);plot(t0,x0,'r')hold onFs =50;t=0:1/Fs:0.1;x=cos(2*pi*10*t); stem(t,x); hold offtitle('连续信号及其抽样信号')函数图像为:)502cos()(2t t x ⨯=π同理,函数图像为:)0102cos()(3t t x ⨯=π同理,函数图像为:由以上的三图可知,第一个图的离散序列,基本可以显示出原来信号,可以通过低通滤波恢复,因为信号的频率为20HZ,而采样频率为50>2*20,故可以恢复,但是第二个和第三个信号的评论分别为50和100HZ,因此理论上是不能够恢复的,需要增大采样频率,解决的方案为,第二个信号的采样频率改为400HZ,而第三个的采样频率改为1000HZ,这样可以很好的采样,如下图所示:2. 产生幅度调制信号)200cos()2cos()(t t t x ππ=,推导其频率特性,确定抽样频率,并绘制波形。
《信号与系统》课程实验报告《信号与系统》课程实验报告一图1-1 向量表示法仿真图形2.符号运算表示法若一个连续时间信号可用一个符号表达式来表示,则可用ezplot命令来画出该信号的时域波形。
上例可用下面的命令来实现(在命令窗口中输入,每行结束按回车键)。
t=-10:0.5:10;f=sym('sin((pi/4)*t)');ezplot(f,[-16,16]);仿真图形如下:图1-2 符号运算表示法仿真图形三、实验内容利用MATLAB实现信号的时域表示。
三、实验步骤该仿真提供了7种典型连续时间信号。
用鼠标点击图0-3目录界面中的“仿真一”按钮,进入图1-3。
图1-3 “信号的时域表示”仿真界面图1-3所示的是“信号的时域表示”仿真界面。
界面的主体分为两部分:1) 两个轴组成的坐标平面(横轴是时间,纵轴是信号值);2) 界面右侧的控制框。
控制框里主要有波形选择按钮和“返回目录”按钮,点击各波形选择按钮可选择波形,点击“返回目录”按钮可直接回到目录界面。
图1-4 峰值为8V,频率为0.5Hz,相位为180°的正弦信号图1-4所示的是正弦波的参数设置及显示界面。
在这个界面内提供了三个滑动条,改变滑块的位置,滑块上方实时显示滑块位置代表的数值,对应正弦波的三个参数:幅度、频率、相位;坐标平面内实时地显示随参数变化后的波形。
在七种信号中,除抽样函数信号外,对其它六种波形均提供了参数设置。
矩形波信号、指数函数信号、斜坡信号、阶跃信号、锯齿波信号和抽样函数信号的波形分别如图1-5~图1-10所示。
图1-5 峰值为8V,频率为1Hz,占空比为50%的矩形波信号图1-6 衰减指数为2的指数函数信号图1-7 斜率=1的斜坡信号图1-8 幅度为5V,滞后时间为5秒的阶跃信号图1-9 峰值为8V,频率为0.5Hz的锯齿波信号图1-10 抽样函数信号仿真途中,通过对滑动块的控制修改信号的幅度、频率、相位,观察波形的变化。
信号与系统实验报告一、实验目的(1) 理解周期信号的傅里叶分解,掌握傅里叶系数的计算方法;(2)深刻理解和掌握非周期信号的傅里叶变换及其计算方法;(3) 熟悉傅里叶变换的性质,并能应用其性质实现信号的幅度调制;(4) 理解连续时间系统的频域分析原理和方法,掌握连续系统的频率响应求解方法,并画出相应的幅频、相频响应曲线。
二、实验原理、原理图及电路图(1) 周期信号的傅里叶分解设有连续时间周期信号()f t ,它的周期为T ,角频率22fT,且满足狄里赫利条件,则该周期信号可以展开成傅里叶级数,即可表示为一系列不同频率的正弦或复指数信号之和。
傅里叶级数有三角形式和指数形式两种。
1)三角形式的傅里叶级数:01212011()cos()cos(2)sin()sin(2)2cos()sin()2n n n n a f t a t a t b t b t a a n t b n t 式中系数n a ,n b 称为傅里叶系数,可由下式求得:222222()cos(),()sin()T T T T nna f t n t dtb f t n t dtTT2)指数形式的傅里叶级数:()jn tn nf t F e式中系数n F 称为傅里叶复系数,可由下式求得:221()T jn tT nF f t edtT周期信号的傅里叶分解用Matlab进行计算时,本质上是对信号进行数值积分运算。
Matlab中进行数值积分运算的函数有quad函数和int函数。
其中int函数主要用于符号运算,而quad函数(包括quad8,quadl)可以直接对信号进行积分运算。
因此利用Matlab进行周期信号的傅里叶分解可以直接对信号进行运算,也可以采用符号运算方法。
quadl函数(quad系)的调用形式为:y=quadl(‘func’,a,b)或y=quadl(@myfun,a,b)。
其中func是一个字符串,表示被积函数的.m文件名(函数名);a、b分别表示定积分的下限和上限。
1、求门函数g4(t)的傅里叶变换,并画出其频率特性曲线。
程序如下:w=linspace(-3*pi,3*pi,512);N=length(w);F=zeros(1,N);for k=1:NF(k)=quadl('sf1',-3,3,[],[],w(k));endfigure(1);plot(w,real(F));xlabel('\omega');ylabel('F(j\omega)');figure(2);plot(w,real(F)-sinc(w/2/pi).^2);xlabel('\omega');title('计算误差');运行结果如图一、图二所示:图一g4(t)的傅里叶变换图二g4(t)傅里叶变换的计算误差2、已知频率特性函数为:H(jw)=2(jw)3+jw+4(jw)4+3(jw)3+2(jw)2+5jw+2求其幅频特性和相频特性。
程序如下:w=linspace(0,5,200);b=[2 0 1 4];a=[1 3 2 5 2];H=freqs(b,a,w);subplot(2,1,1);plot(w,abs(H));set(gca,'xtick',[0 1 2 3 4 5]);set(gca,'ytick',[0 0.4 0.707 1]);xlabel('\omega');ylabel('|H(j\omega)|');subplot(2,1,2);plot(w,angle(H));set(gca,'xtick',[0 1 2 3 4 5]);xlabel('\omega');ylabel('\phi(\omega)');运行结果如图三所示:图三幅频特性和相频特性3、设H s=s(s−p1)(s−p2)设①p1=-2,p2=-30;②p1=-2,p2=3(1)针对极点参数①②,画出系统零、极点分布图,判断该系统稳定性。
《信号与系统》实验报告湖南工业大学电气与信息工程学院实验一用同时分析法观测50Hz非正弦周期信号的分解与合成一、实验目的1、用同时分析法观测50Hz非正弦周期信号的频谱,并与傅立叶级数各项的频率与系数作比较。
2、观测基波和其谐波的合成。
二、实验设备1、信号与系统实验箱:TKSS -A型或TKSS -B 型TKSS -C 型;2、双踪示波器三、实验原理1、 一个非正弦周期函数可以用一系列频率成整数倍的正弦函数来表示,其中与非正弦具有相同频率的成分称为基波或一次谐波,其他成分则根据其频率为基波频率的2、3、4、…、n 等倍数分别称为二次、三次、四次、…、n 次谐波,其幅度将随着谐波次数的增加而减小,直至无穷小。
2、 不同频率的谐波可以合成一个非正弦周期波,反过来,一个非正弦周期波也可以分解为无限个不同频率的谐波成分,3、 一个非正弦周期函数可以用傅立叶级数来表示,级数各项系数之间的关系可用一个频谱来表示,不同的非正弦周期函数具有不同的频谱图,各种不同波形及其傅氏级数表达式见表2-1,方波频谱图如图2-1表示Um1351/91/51/71/3790ωωωωωω图1-1 方波频谱图表2-1 各种不同波形的傅立叶级数表达式UmtTU 2τ方波Um0TU 2τ正弦整流全波UmTU 2τ三角波Um0T2τ正弦整流半波t tUm0tT U 2τ矩形波U1、方波 ())7s i n 715s i n 513s i n 31(s i n 4 ++++=t t t t u t u mωωωωπ 2、三角波())5s i n 2513sin 91(sin 82++-=t t t u t u mωωωπ3、半波())4c o s 1512cos 31sin 421(2 +--+=t t t u t u m ωωωππ 4、全波 ())6c o s 3514cos 1512cos 3121(4 +---=t t t u t u m ωωωπ5、 矩形波())3cos 3sin 312cos 2sin 21cos (sin 2 ++++=t T t T t T U T U t u m m ωτπωτπωτππτ实验装置的结构如图1-2所示DC20f f f f f f 3456图1-2信号分解于合成实验装置结构框图图中LPF 为低通滤波器,可分解出非正弦周期函数的直流分量。
福建xx大学xx学院信息工程类实验报告课程名称:信号与系统姓名:XXX系:电子信息工程专业:电子信息工程年级:XXX级学号:XXX指导教师:XXX职称:讲师XXX年12 月17日实验项目列表XXX 信息工程类实验报告系: 电子信息工程 专业: 电子信息工程 年级: 07级 姓名: XXX 学号: XXX 实验课程: 信号与系统实验室号:_信号与系统实验室 实验设备号: 01 实验时间: 11.19 指导教师签字: 成绩:实验一 函数信号发生器1、 实验目的1)了解单片多功能集成电路函数信号发生器的功能及特点。
2)熟悉信号与系统实验箱信号产生和测试的方法。
2、 实验仪器1)信号与系统实验箱一台。
2)20MHz 双踪示波器一台。
3、 实验原理ICL8038是单片机集成函数信号发生器,其内部框图如图1.1所示。
它由恒流源1I 和2I 、电压比较器A 和B 、触发器、缓冲器和三角波变正弦波电路等组-V EE图1.1 ICL8038原理方框图外接电容C 由两个恒流源充电和放电,电压比较器A 、B 的阀值分别为电源电压(指EE cc U U +)的2/3和1/3。
恒流源1I 和2I 的大小可通过外接电阻调节,但必须12I I >。
当触发器的输出为低电平时,恒流源2I 断开,恒流源1I 给C 充电,它的两端电压UC 随时间线性上升,当UC 达到电源电压的2/3时,电压比较器A 的输出电压发生跳变,使触发器输出由低电平变为高电平,恒流源C 接通,由于12I I >(设122I I =),恒流源2I 将电流21I 加到C 上反充电,相当于C 由一个净电流I 放电,C 两端的电压UC 又转为直线下降。
当它下降到电源电压的1/3时,电压比较器B 的输出电压发生跳变,使触发器的输出由高电平跳变为原来的低电平,恒流源2I 断开,1I 再给C 充电,…如此周而复始,产生振荡。
若调整电路,使122I I =,则触发器输出为方波,经反相缓冲器由管脚⑨输出方波信号。
《信号与系统》上机实验实验一连续时间信号的表示及可视化一.实验目的熟练掌握连续时间信号的表示及可视化处理。
二.实验源程序δf(t))=)(tf=@(t)dirac(t) %定义函数ezplot(f,[-5:5]); %利用ezplot( )命令绘制连续图形xlabel('(t)'); %横坐标ylabel('(f)'); %纵坐标f(t)= ε(t)(f=Heaviside(n))f=@(t)heaviside(t) %定义函数ezplot(f,[-5:5]); %利用ezplot( )命令绘制连续图形xlabel('(t)'); %横坐标ylabel('(f)'); %纵坐标f(t)=at e(分别取a>0及a<0)a=1时f=@(t)exp(t) %定义函数ezplot(f,[-5:5]); %利用ezplot( )命令绘制连续图形xlabel('(t)'); %横坐标ylabel('(f)'); %纵坐标a=-1时f=@(t)exp(-t) %定义函数ezplot(f,[-5:5]); %利用ezplot( )命令绘制连续图形xlabel('(t)'); %横坐标ylabel('(f)'); %纵坐标f(t)=R(t)t=-5:0.01:5; %设定时间变量t的范围及步长y=rectpuls(t,2); %用rectpuls(t a)命令表示门函数,默认以零点为中心,宽度为aplot(t,y); %用plot函数绘制连续函数grid on; %显示网格命令title('门函数'); %用title函数设置图形的名称axis([-5 5 -0.5 1.5]);f(t)=Sa(wt)w=5时,f=Sa(5*t)f=@(t)Sinc(5*t) %定义函数ezplot(f,[-5:5]); %利用ezplot( )命令绘制连续图形xlabel('(t)'); %横坐标ylabel('(f)'); %纵坐标axis([-5 5 -1.2 1.2])w=8时,f=Sa(8*t)f=@(t)sinc(8*t) %定义函数ezplot(f,[-4:4]); %利用eaplot( )命令绘制连续图形xlabel('(t)'); %横坐标ylabel('(f)'); %纵坐标f(t)=Sin(2πft)(分别画出不同周期个数的波形)f(t)=Sin(t)f=@(t)sin(t) %定义函数ezplot(f,[-15:15]); %利用eaplot( )命令绘制连续图形xlabel('(t)'); %横坐标ylabel('(f)'); %纵坐标axis([-15 15 -1.2 1.2])三.程序运行结果(1)(2)(3)-5-4-3-2-1012345-1-0.8-0.6-0.4-0.200.20.40.60.81(t)dirac(t)(f )-5-4-3-2-101234500.20.40.60.81(t)heav iside(t)(f )(4)-5-4-3-2-1012345010********607080(t)exp(t)(f )-5-4-3-2-1012345010********607080(t)exp(-t)(f )(5)-5-4-3-2-1012345-1-0.8-0.6-0.4-0.200.20.40.60.81(t)Sinc(5 t)(f )(6)-4-3-2-101234 -1-0.8-0.6-0.4-0.20.20.40.60.81(t)sinc(8 t)(f)-15-10-5051015 -1-0.8-0.6-0.4-0.20.20.40.60.81(t)sin(t)(f)实验二离散时间信号的表示及可视化一.实验目的学会对离散时间信号进行标识和可视化处理。
信号与系统实验报告目录1. 内容概要 (2)1.1 研究背景 (3)1.2 研究目的 (4)1.3 研究意义 (4)2. 实验原理 (5)2.1 信号与系统基本概念 (7)2.2 信号的分类与表示 (8)2.3 系统的分类与表示 (9)2.4 信号与系统的运算法则 (11)3. 实验内容及步骤 (12)3.1 实验一 (13)3.1.1 实验目的 (14)3.1.2 实验仪器和设备 (15)3.1.4 实验数据记录与分析 (16)3.2 实验二 (16)3.2.1 实验目的 (17)3.2.2 实验仪器和设备 (18)3.2.3 实验步骤 (19)3.2.4 实验数据记录与分析 (19)3.3 实验三 (20)3.3.1 实验目的 (21)3.3.2 实验仪器和设备 (22)3.3.3 实验步骤 (23)3.3.4 实验数据记录与分析 (24)3.4 实验四 (26)3.4.1 实验目的 (27)3.4.2 实验仪器和设备 (27)3.4.4 实验数据记录与分析 (29)4. 结果与讨论 (29)4.1 实验结果汇总 (31)4.2 结果分析与讨论 (32)4.3 结果与理论知识的对比与验证 (33)1. 内容概要本实验报告旨在总结和回顾在信号与系统课程中所进行的实验内容,通过实践操作加深对理论知识的理解和应用能力。
实验涵盖了信号分析、信号处理方法以及系统响应等多个方面。
实验一:信号的基本特性与运算。
学生掌握了信号的表示方法,包括连续时间信号和离散时间信号,以及信号的基本运算规则,如加法、减法、乘法和除法。
实验二:信号的时间域分析。
在本实验中,学生学习了信号的波形变换、信号的卷积以及信号的频谱分析等基本概念和方法,利用MATLAB工具进行了实际的信号处理。
实验三:系统的时域分析。
学生了解了线性时不变系统的动态响应特性,包括零状态响应、阶跃响应以及脉冲响应,并学会了利用MATLAB进行系统响应的计算和分析。
信号与系统实验报告一、信号的时域基本运算1.连续时间信号的时域基本运算两实验之一实验分析:输出信号值就等于两输入信号相加(乘)。
由于b=2,故平移量为2时,实际是右移1,符合平移性质。
两实验之二心得体会:时域中的基本运算具有连续性,当输入信号为连续时,输出信号也为连续。
平移,伸缩变化都会导致输出结果相对应的平移伸缩。
2.离散时间信号的时域基本运算两实验之一实验分析:输出信号的值是对应输入信号在每个n值所对应的运算值,当进行拉伸变化后,n值数量不会变,但范围会拉伸所输入的拉伸系数。
两实验之二心得体会:离散时间信号可以看做对连续时间信号的采样,而得到的输出信号值,也可以看成是连续信号所得之后的采样值。
二、连续信号卷积与系统的时域分析1.连续信号卷积积分两实验之一实验分析:当两相互卷积函数为冲激函数时,所卷积得到的也是一个冲激函数,且该函数的冲激t值为函数x,函数y冲激t值之和。
两实验之二心得体会:连续卷积函数每个t值所对应的卷积和可以看成其中一个在k值取得的函数与另外一个函数相乘得到的一个分量函数,并一直移动k值直至最后,最后累和出来的最终函数便是所得到的卷积函数。
3.RC电路时域积分两实验之一实验分析:全响应结果正好等于零状态响应与零输入响应之和。
两实验之二心得体会:具体学习了零状态,零输入,全响应过程的状态及变化,与之前所学的电路知识联系在一起了。
三、离散信号卷积与系统的时域分析1.离散信号卷积求和两实验之一实验分析:输出结果的n值是输入结果的k号与另一个n-k的累和两实验之二心得体会:直观地观察到卷积和的产生,可以看成连续卷积的采样形式,从这个方面去想,更能深入地理解卷积以及采样的知识。
2.离散差分方程求解两实验之一实验分析:其零状态响应序列为0 0 4 5 7.5,零输入响应序列为2 4 5 5.5 5.75,全状态响应序列为2 4 9 10.5 13.25,即全状态=零输入+零状态。
两实验之二心得体会:求差分方程时,可以根据全状态响应是由零输入输入以及零状态相加所得,分开来求,同时也加深了自己对差分方程的求解问题的理解。
信号实验报告信息科学与工程学院电子信息工程姓名:学号:实验一一、实验内容:1、连续信号的MATLAB表示:MATLAB提供了大量的产生基本信号的函数。
最常用的指数信号、正弦信号是MATLAB的内部函数,即不安装任何工具箱就可调用的函数。
2、编程实现题2-2已知信号x(t)如图m2-2图所示,分别用MATLAB表示信号x(t)、x(t)cos(50t),并画出波形,程序如下:function yt=x2_2(t)yt=1*(t>=-1&t<0)+4*t.*(t>=0&t<=0.5)+(-4)*(t-1).*(t>0.5&t<=1);%example2_2t=-2:0.001:2;subplot(2,1,1)plot(t,x2_2(t))title('x(t)')subplot(2,1,2)A=1;w0=50;t=-2:0.001:2;phi=0;xt=[A*cos(w0*t+phi)].*[1*(t>=-1&t<0)+4*t.*(t>=0&t<=0.5)+(-4)*(t-1).*(t>0.5&t<=1) ];plot(t,xt)title('x(t)cos(50t)')运行结果如图所示:2-3(1)编写表示题M2-3图所示的信号波形x(t)的MATLAB函数(2)试画出x(t)、x(0.5t)、和x(2-0.5t)的波形图。
程序如下:function yt=x2_3(t)yt=1*t.*(t>=0&t<2)+2*(t>=2&t<3)+(-1)*(t>=3&t<5);%example2_3t=-4:0.01:10;subplot(3,1,1)plot(t,x2_3(t))title('x(t)')axis([-4,10,-2,3])t=-4:0.01:10;subplot(3,1,2)plot(t,x2_3(0.5*t))title('x(0.5t)')axis([-4,10,-2,3])t=-6:0.01:8;subplot(3,1,3)plot(t,x2_3(2-0.5*t))title('x(2-0.5t)')axis([-6,8,-2,3]运行结果如有右图:2-4 画出题M2-4图所示信号的奇分量和偶分量。
程序如下:2-7已知连续信号x1(t)=cos(6∏t)、x2(t)=cos(14∏t)、x3(t)=cos(26∏t),以抽样函数fs=10Hz对这三个信号进行抽样得离散序列x1[k]、x2[k]、x3[k].试在同一图上画出连续信号和其他对应的离散序列。
程序如下:%program2-7A=1;w1=6*pi;phi=0;t1=0:0.0001:1;xt1=A*cos(w1*t1+phi);subplot(2,3,1)plot(t1,xt1)title('x1(t)')k1=0:20;xk1=cos(3*pi/5*k1);subplot(2,3,4)stem(k1,xk1)title('x1[k]')A=1;w2=14*pi;phi=0;t2=0:0.0001:1;xt2=A*cos(w2*t2+phi);subplot(2,3,2)plot(t2,xt2)title('x2(t)')k2=0:20;xk2=cos(7*pi/5*k2);subplot(2,3,5)stem(k2,xk2)title('x2[k]')A=1;w3=26*pi;phi=0;t3=0:0.0001:1;xt3=A*cos(w3*t3+phi);subplot(2,3,3)plot(t3,xt3)title('x3(t)')k3=0:20;xk3=cos(13*pi/5*k3); subplot(2,3,6)stem(k3,xk3)title('x3[k]')运行结果如右上图:2-8分别用square函数和Swtooth画出题M2-8图所示的离散周期序列。
程序如下:%program2-81k=0:1:40;A=1;T=20;w0=2*pi/T;xk=A*sawtooth(w0*k,0.5); stem(k,xk)运行结果如右图所示:%program2-82k=0:1:40;A=-1;T=10;w0=2*pi/T;xk=A*square(w0*k,50); stem(k,xk)运行结果如右图:实验二3-1一个连续时间LTI系统满足的微分方程为:y″(t)+3y′(t)+2y(t) =2x′(t)+x(t)%program3_1ts=0;te=10;dt=0.001; sys=tf([2 1],[1 3 2]); t=ts:dt:te;x=1*exp(-3*t)*1;y=lsim(sys,x,t);plot(t,y);xlabel('Time(sec)') ylabel('y(t)')3-2程序如下:%program3-2ts=0;te=5;dt=0.01;sys=tf([6],[1 3]);t=ts:dt:te;y=impulse(sys,t);plot(t,y);xlabel('Time(sec)') ylabel('h(t)')运行结果如有图所示:%program3-2ts=0;te=5;dt=0.01;sys=tf([6],[1 3]);t=ts:dt:te;y=step(sys,t);plot(t,y);xlabel('Time(sec)')ylabel('h(t)')运行结果如有图所示:3-4程序如下:%program3_4x=[0.85,0.53,0.21,0.67,0.84,0.12]; y=[0.68,0.37,0.83,0.52,0.71];z=conv(x,y);n=length(z);stem(-3:n-4,z);运行结果如有图所示:3-7某离散时间LTI系统满足差分方程Y[k]+0.42y[k-1]-0.19y[k-2]=0.31x[k]+0.68x[k-1],若x(t)=0.6ku(k)试求系统的零状态响应Yzs(k),并画出前20点图。
程序如下%program3_7k=0:100;a=[1,0.42,-0.19];b=[0.31,0.68];A=1;a=0.6;xk=A*a.^k;stem(k,xk)y=filter(b,a,xk);axis([0,20,0,1])xlabel('k');ylabel('Yzs[k]')运行结果如右图所示:3-8某离散时间LTI系统满足差分方程Y[k]+0.7y[k-1]-0.45y[k-2]-0.6y[k-3]=0.8x[k]-0.44x[k-1]+0.36x[k-2]+0.02x[k],试利用impz函数求其单位脉冲响应,,并画出前30点图。
程序如下:%program3_8k=0:100;a=[1 0.7 0.45 0.6];b=[0.8 -0.44 0.36 0.02];h=impz(b,a,k);stem(k,h)axis([0,30,-1,0.8])运行结果如右图所示:实验三M4-1 试求题M4-1图所示周期矩形信号和周期三角波信号的频谱,并画出频谱图,取A=1,T=2。
得如下程序4_1_1程序如下k=-4:0.01:4;A=1;T=2;w0=2*pi/T;xk=A*square(w0*k,50); k1=-4:4;N=length(k);k1=(0:N)*0.5*T;x=A*square(w0*k1,50);X=fft(x);subplot(3,1,1)stem(k1,real(X))xlabel('k');title('X[k]的实部');axis([0,40,-20,60])subplot(3,1,2)stem(k1,imag(X));xlabel('k');title('X[k]的虚部');axis([0,40,-20,60])xr=ifft(X);subplot(3,1,3)stem(k1,xr)axis([0,40,-2,2])xlabel('k');title('重建的X[k]');4-1-2 程序如下:k1=-4:0.001:4;N=length(k);k1=(0:N)*0.5*T;x=0.5*sawtooth(w0*k1,0.5)+0.5;X=fft(x);subplot(3,1,1)stem(k1,real(X))xlabel('k');title('X[k]的实部');axis([0,40,-20,60])subplot(3,1,2)stem(k1,imag(X));xlabel('k');title('X[k]的虚部');axis([0,40,-20,60])xr=ifft(X);subplot(3,1,3)stem(k1,xr)axis([0,40,-2,2])xlabel('k');title('重建的X[k]');运行结果如有图所示:M4-4 试利用MATLAB计算下列连续信号的频谱X(jw),画出频谱图4-4-1程序如下:t u=t u-t-t xt ut u=t x)[)((2)];2cos()()1(-)(()()2();2-function yt=u(t)yt=0*(t<0)+1*(t>=0);function yt=xx1(t,w);yt=(u(t)-u(t-2)).*exp(-j*w*t);w=linspace(-6*pi,6*pi,512);N=length(w);X=zeros(1,N);for k=1:NX(k)=quadl('xx1',0,2,[],[],w(k));endplot(w,real(X));xlabel('\omega');ylabel('X(j\omega)');运行结果如有图所示:4-4-2程序如下:function yt=u(t)yt=0*(t<0)+1*(t>=0);function yt=xx2(t,w);yt=cos(2*t).*(u(t)-u(t-2)).*exp(-j*w*t);w=linspace(-6*pi,6*pi,512);N=length(w);X=zeros(1,N);for k=1:NX(k)=quadl('xx2',0,2,[],[],w(k));endplot(w,real(X));xlabel('\omega ');ylabel('X(j\omega)');运行结果如下图所示:M4-2 (1)若以90.0||2||1220>=∑+=PC C Nn n 定义信号的有效带宽,试确定题M4-1图所示信号的有效带宽0w N ,取A=1,T=2;(2)画出有效带宽内有限项谐波合成的近似波形,并对结果加以讨论和比较;(3)增加谐波的项数,观察其合成的近似波形,并对结果加以讨论和比较;得如下程序::function y=x1(t,w);y=cos(2*t).*(t>=0&t<=2).*exp(-j*w*t);w=linspace(-6*pi,6*pi,1024);N=length(w);X=zeros(1,N);for k=1:NX(k)=quadl('x1',0,2,[],[],w(k)); end figure(1)plot(w,real(X));xlabel('\omega ');ylabel('X(j\omega )'运行结果如右图所示:M4-3 求题M4-3图所示周期矩形信号的幅度频谱,并画出频谱图。