文档库 最新最全的文档下载
当前位置:文档库 › 离散时间系统的时域特性分析实验报告

离散时间系统的时域特性分析实验报告

离散时间系统的时域特性分析实验报告
离散时间系统的时域特性分析实验报告

多元时间序列建模分析

应用时间序列分析实验报告

单位根检验输出结果如下:序列x的单位根检验结果:

1967 58.8 53.4 1968 57.6 50.9 1969 59.8 47.2 1970 56.8 56.1 1971 68.5 52.4 1972 82.9 64.0 1973 116.9 103.6 1974 139.4 152.8 1975 143.0 147.4 1976 134.8 129.3 1977 139.7 132.8 1978 167.6 187.4 1979 211.7 242.9 1980 271.2 298.8 1981 367.6 367.7 1982 413.8 357.5 1983 438.3 421.8 1984 580.5 620.5 1985 808.9 1257.8 1986 1082.1 1498.3 1987 1470.0 1614.2 1988 1766.7 2055.1 1989 1956.0 2199.9 1990 2985.8 2574.3 1991 3827.1 3398.7 1992 4676.3 4443.3 1993 5284.8 5986.2 1994 10421.8 9960.1 1995 12451.8 11048.1 1996 12576.4 11557.4 1997 15160.7 11806.5 1998 15223.6 11626.1 1999 16159.8 13736.5 2000 20634.4 18638.8 2001 22024.4 20159.2 2002 26947.9 24430.3 2003 36287.9 34195.6 2004 49103.3 46435.8 2005 62648.1 54273.7 2006 77594.6 63376.9 2007 93455.6 73284.6 2008 100394.9 79526.5 run; proc gplot; plot x*t=1 y*t=2/overlay; symbol1c=black i=join v=none; symbol2c=red i=join v=none w=2l=2; run; proc arima data=example6_4; identify var=x stationarity=(adf=1); identify var=y stationarity=(adf=1); run; proc arima; identify var=y crrosscorr=x; estimate methed=ml input=x plot; forecast lead=0id=t out=out; proc aima data=out; identify varresidual stationarity=(adf=2); run;

实验四 控制系统频率特性的测试(实验报告)

实验四 控制系统频率特性的测试 一. 实验目的 认识线性定常系统的频率特性,掌握用频率特性法测试被控过程模型的原理和方法,根据开环系统的对数频率特性,确定系统组成环节的参数。 二.实验装置 (1)微型计算机。 (2)自动控制实验教学系统软件。 三.实验原理及方法 (1)基本概念 一个稳定的线性定常系统,在正弦信号的作用下,输出稳态与输入信号关系如下: 幅频特性 相频特性 (2)实验方法 设有两个正弦信号: 若以)(t x ω为横轴,以)(y t ω为纵轴,而以t ω作为参变量,则随t ω的变化,)(t x ω和 )(y t ω所确定的点的轨迹,将在 x--y 平面上描绘出一条封闭的曲线(通常是一个椭圆)。这 就是所谓“李沙育图形”。 由李沙育图形可求出Xm ,Ym ,φ,

四.实验步骤 (1)根据前面的实验步骤点击实验七、控制系统频率特性测试菜单。 (2)首先确定被测对象模型的传递函数, 预先设置好参数T1、T2、ξ、K (3)设置好各项参数后,开始仿真分析,首先做幅频测试,按所得的频率范围由低到高,及ω由小到大慢慢改变,特别是在转折频率处更应该多取几个点 五.数据处理 (一)第一种处理方法: (1)得表格如下: (2)作图如下: (二)第二种方法: 由实验模型即,由实验设置模型根据理论计算结果绘制bode图,绘制Bode图。

(三)误差分析 两图形的大体趋势一直,从而验证了理论的正确性。在拐点处有一定的差距,在某些点处也存在较大的误差。 分析: (1)在读取数据上存在较大的误差,而使得理论结果和实验结果之间存在。 (2)在数值应选取上太合适,而使得所画出的bode图形之间存在较大的差距。 (3)在实验计算相角和幅值方面本来就存在着近似,从而使得误差存在,而使得两个图形之间有差异 六.思考讨论 (1)是否可以用“李沙育”图形同时测量幅频特性和想频特性 答:可以。在实验过程中一个频率可同时记录2Xm,2Ym,2y0。 (2)讨论用“李沙育图形”测量频率特性的精度,即误差分析(说明误差的主要来源)答:用“李沙育图形”测量频率特性的精度从上面的分析处理上也可以看出是比较高的,但是在实验结果和理论的结果之间还是存在一定的差距,这些误差主要来自于从“李沙育图形”上读取数据的时候存在的误差,也可能是计算机精度方面的误差。 (3)对用频率特性测试系统数学模型方法的评测 答:用这种方法进行此次实验能够让我们更好地了解其过程,原理及方法。但本次实验的数据量很大,需要读取较多坐标,教学软件可以更智能一些,增加一些自动读取坐标的功能。 七.实验总结 通过本次实验,我加深了对线性定常系统的频率特性的认识,掌握了用频率特性法测试被控过程模型的原理和方法。使我把书本知识与实际操作联系起来,加深了对课程内容的理解。在处理数据时,需要进行一定量的计算,这要求我们要细心、耐心,作图时要注意不能用普通坐标系,而是半对数坐标系进行作图。

离散时间系统特性分析

实验五实验报告 实验名称:离散时间系统特性分析

一、实验目的: 1 。深入理解单位样值响应,离散系统的频率响应的概念; 2。 掌握通过计算机进行求得离散系统的单位样值响应,以及离散系统的频率 响应的方法。 二、实验原理: 对于离散系统的单位样值而言,在实际处理过程中,不可能选取无穷多项的取值。往往是选取有限项的取值,当然这里会产生一个截尾误差,但只要这个误差在相对小一个范围里,可以忽略不计。 另外,在一些实际的离散系统中,往往不是事先就能得到描述系统的差分方程的,而是通过得到系统的某些相应值,则此时系统的分析就需借助计算机的数值处理来进行,得到描述系统的某些特征,甚至进而得到描述系统的数学模型。 本实验首先给出描述系统的差分方程,通过迭代的方法求得系统的单位样值响应,进而求得该离散系统的频率响应。限于试验条件,虽然给出了系统方程,但处理的方法依然具有同样的实际意义。 具体的方法是: 1 在给定系统方程的条件下,选取激励信号为δ(n),系统的起始状态为零 状态,通过迭代法,求得系统的单位样值响应h(n)(n=0,…,N )。 2 利用公式 其中Ω的取值范围为0~2π 。计算系统的频率响应。 三、实验内容 1 已知系统的差分方程为 利用迭代法求得系统的单位样值响应,取N =10。 2 利用公式 其中

#include #include #define N 10 #define M 20 #define pi 3.1415926 struct pinlv{ double fu; double xiang;}; double h[N+1],x[N+1]; struct pinlv PL(double w) { double a=0, b=0,fu,xiang; int k; struct pinlv FX; for(k=0;k<=N;k++){ a=a+h[k]*cos(-k*w); b=b+h[k]*sin(-k*w);} fu=sqrt(a*a+b*b); xiang=atan(b/a); if((a<0)&(b>0)) xiang=xiang+pi; if((a<0)&(b<0)) xiang=xiang-pi; FX.fu=fu; FX.xiang=xiang; return(FX); } main() { int i,j; double w0; struct pinlv FX[M+1]; FILE *fp1,*fp2; fp1=fopen("H:\\单位样值响应.txt","w"); fp2=fopen("H:\\频率特性.txt","w"); h[-1]=0;h[-2]=0; for(i=-1;i<=N;i++) x[i]=0; x[0]=1; for(i=0;i<=N;i++) h[i]=1.3*h[i-1]-0.4*h[i-2]+x[i-1]; printf("系统的单位样值响应为\n"); fprintf(fp1,"系统的单位样值响应(从x[0]开始)为\n"); fprintf(fp1,"激励x[i] 响应y[i]\n"); for(i=0;i<=N;i++)

频率特性分析

实验三 频率特性分析 一·实验目的 1.掌握频率特性的基本概念,尤其是频率特性的几种表示方法。 2.能熟练绘制极坐标频率特性曲线(奈奎斯特曲线)和对数频率特性曲线,尤其要注意的是在非最小相位系统时曲线的绘制。 3.正确应用频率稳定判别方法,包括奈奎斯特稳定判据和对数稳定判据。 4.熟练正确计算相位裕量和幅值裕量。 5.掌握闭环频率特性的基本知识以及有关指标的近似估算方法。 二·实验内容 1增加开环传递函数零极点个数对奈奎斯特图的影响 1)改变有限极点个数n ,使n=0,1,2,3 Nyquist Diagram Real Axis I m a g i n a r y A x i s -2 -101234 -3.5-3-2.5-2-1.5-1-0.50 0.511.52n=0 n=1 n=2 n=3 2)改变原点处极点个数v ,当v=1,2,3,4, Nyquist Diagram Real Axis I m a g i n a r y A x i s -2 -1.5 -1 -0.5 00.5 1 1.5 2 -2-1.5 -1 -0.5 00.5 1 1.5 2 System: sys P hase Margin (deg): -32.9Delay Margin (sec): 4.41At frequency (rad/sec): 1.3 Closed Loop Stable? No System: sys P hase Margin (deg): -121Delay Margin (sec): 3.49At frequency (rad/sec): 1.2 Closed Loop Stable? No System: sys P hase Margin (deg): 150Delay Margin (sec): 2.28At frequency (rad/sec): 1.15Closed Loop Stable? No System: sys P hase Margin (deg): 51.8Delay Margin (sec): 0.575 At frequency (rad/sec): 1.57 Closed Loop Stable? Yes v=1 v=2 v=3 v=4

系统频率特性的测试实验报告

东南大学自动化学院课程名称:自动控制原理实验 实验名称:系统频率特性的测试 姓名:学号: 专业:实验室: 实验时间:2013年11月22日同组人员: 评定成绩:审阅教师:

一、实验目的: (1)明确测量幅频和相频特性曲线的意义; (2)掌握幅频曲线和相频特性曲线的测量方法; (3)利用幅频曲线求出系统的传递函数; 二、实验原理: 在设计控制系统时,首先要建立系统的数学模型,而建立系统的数学模型是控制系统设计的重点和难点。如果系统的各个部分都可以拆开,每个物理参数能独立得到,并能用物理公式来表达,这属机理建模方式,通常教材中用的是机理建模方式。如果系统的各个部分无法拆开或不能测量具体的物理量,不能用准确完整的物理关系式表达,真实系统往往是这样。比如“黑盒”,那只能用二端口网络纯的实验方法来建立系统的数学模型,实验建模有多种方法。此次实验采用开环频率特性测试方法,确定系统传递函数。准确的系统建模是很困难的,要用反复多次,模型还不一定建准。另外,利用系统的频率特性可用来分析和设计控制系统,用Bode 图设计控制系统就是其中一种。 幅频特性就是输出幅度随频率的变化与输入幅度之比,即)()(ωωi o U U A =。测幅频特性时, 改变正弦信号源的频率,测出输入信号的幅值或峰峰值和输输出信号的幅值或峰峰值。 测相频有两种方法: (1)双踪信号比较法:将正弦信号接系统输入端,同时用双踪示波器的Y1和Y2测量系统的输入端和输出端两个正弦波,示波器触发正确的话,可看到两个不同相位的正弦波,测出波形的周期T 和相位差Δt ,则相位差0360??=ΦT t 。这种方法直观,容易理解。就模拟示波 器而言,这种方法用于高频信号测量比较合适。 (2)李沙育图形法:将系统输入端的正弦信号接示波器的X 轴输入,将系统输出端的正弦信号接示波器的Y 轴输入,两个正弦波将合成一个椭圆。通过椭圆的切、割比值,椭圆所在的象限,椭圆轨迹的旋转方向这三个要素来决定相位差。就模拟示波器而言,这种方法用于低频信号测量比较合适。若用数字示波器或虚拟示波器,建议用双踪信号比较法。 利用幅频和相频的实验数据可以作出系统的波Bode 图和Nyquist 图。 三、预习与回答: (1)实验时,如何确定正弦信号的幅值?幅度太大会出现什么问题,幅度过小又会出现什 么问题? 答:根据实验参数,计算正弦信号幅值大致的范围,然后进行调节,具体确定调节幅值时,首先要保证输入波形不失真,同时,要保证在频率较大时输出信号衰减后人能够测量出来。如果幅度过大,波形超出线性变化区域,产生失真;如果波形过小,后续测量值过小,无法精确的测量。

实验三___离散时间系统的时域分析

实验三 离散时间系统的时域分析 1.实验目的 (1)理解离散时间信号的系统及其特性。 (2)对简单的离散时间系统进行分析,研究其时域特性。 (3)利用MATLAB对离散时间系统进行仿真,观察结果,理解其时域特性。 2.实验原理 离散时间系统,主要是用于处理离散时间信号的系统,即是将输入信号映射成的输出的某种运算,系统的框图如图所示: (1)线性系统 线性系统就是满足叠加原理的系统。如果对于一个离散系统输入信号为时,输出信号分别为,即:。 而且当该系统的输入信号为时,其中a,b为任意常数,输出为,则该系统就是一个线性离散时间系统。 (2)时不变系统 如果系统的响应与激励加于系统的时刻无关,则该系统是时不变系统。对于一个离散时间系统,若输入,产生输出为,则输入为,产生输出为,即: 若,则。 通常我们研究的是线性时不变离散系统。 3.实验内容及其步骤 (1)复习离散时间系统的主要性质,掌握其原理和意义。 (2)一个简单的非线性离散时间系统的仿真 系统方程为: x = cos(2*pi*0.05*n); x1[n] = x[n+1] x2[n] = x[n] x3[n] = x[n-1] y = x2.*x2-x1.*x3; 或者:y=x*x- x[n+1]* x[n-1] 是非线性。 参考:% Generate a sinusoidal input signal clf; n = 0:200; x = cos(2*pi*0.05*n); % Compute the output signal x1 = [x 0 0]; % x1[n] = x[n+1] x2 = [0 x 0]; % x2[n] = x[n] x3 = [0 0 x]; % x3[n] = x[n-1]

离散信号与系统时域分析

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

时间序列分析实验报告

时间序列分析实验报告 P185#1、某股票连续若干天的收盘价如表5-4 (行数据)所示。 表5-4 304 303 307 299 296 293301 293 301 295 284286 286 287 284 282278 281 278 277279 278 270 268 272 273 279 279280 275 271 277 278279 283 284 282 283279 280 280 279278 283 278 270 275 273 273 272275 273 273 272 273272 273 271 272 271273 277 274 274272 280 282 292 295 295 294 290 291 288 288 290 293 288 289 291 293 293 290 288 287 289 292 288 288 285 282 286 286 287 284 283 286 282 287 286 287 292 292 294 291 288 289 选择适当模型拟合该序列的发展,并估计下一天的收盘价。 解: (1)通过SA漱件画出上述序列的时序图如下: 程序: data example5_1; in put x@@; time=_ n_; cards ; 304 303 307 299296 293 301 293 301 295 284286286 287 284 282 278 281 278277 279 278 270 268 272 273279279 280 275 271 277 278 279283 284 282 283 279 280 280279278 283 278 270 275 273 273272 275 273 273 272 273 272273271 272 271 273 277 274 274272 280 282 292 295 295 294290291 288 288 290 293 288 289291 293 293 290 288 287 289292288 288 285 282 286 286 287284 283 286 282 287 286 287292292 294 291 288 289 proc gplot data =example5_1; plot x*time= 1; symbol1 c=black v=star i =join; run ; 上述程序所得时序图如下: 上述时序图显示,该序列具有长期趋势又含有一定的周期性,为典型的非平稳序列。又因为该序列呈现曲线形式,所以选择2阶差分。

自动控制原理控制系统的频率特性实验报告

肇庆学院 工程学院 自动控制原理实验报告 12 年级 电气一班 组员:王园园、李俊杰 实验日期 2014/6/9 姓名:李奕顺 学号:201224122130老师评定 ________________ 实验四:控制系统的频率特性 一、实验原理 1.被测系统的方块图:见图4-1 将频率特性测试仪内信号发生器产生的超低频正弦信号的频率从低到高变化, 并施加于 被测系统的输人端[r(t)],然后分别测量相应的反馈信号 [b(t)]和误差信号[e(t)]的对数幅 值和 相位。频率特性测试仪测试数据经相关运算器后在显示器中显示。 根据式(4 — 3)和式(4 — 4)分别计算出各个频率下的开环对数幅值和相位, 在半对数座标 纸上作出实验曲线:开环对数幅频曲线和相频曲线。 系统(或环节)的频率特性 幅值和相角: G (j 3)是一个复变量,可以表示成以角频率 3为参数的 G(j 3)= G(j 3)|/G(j 3) (4 — 1) 本实验应用频率特性测试仪测量系统或环节的频率特牲。 图4-1所示系统的开环频率特性为: G 1(j 3)G 2(j 3) B(j 3) 」 B(j 3) E(j 3) E(j 3) E(j 3) (4—2) 采用对数幅频特性和相频特性表示,则式( 20lgG1(j 3) G2(j 3)H(j 3)= 2 叫鵲 = 20lgB(j 3) -20lg E(j 3) (4— 3) G 1(j 3)G 2(j 3)H(j 3) 二 B(j 3)- . E(j 3) (4—4) 图4-1 被测系统方块图 4— 2 )表示 为:

根据实验开环对数幅频曲线画出开环对数幅频曲线的渐近线,再根据渐近线的斜率和转 角频确定频率特性(或传递函数)。所确定的频率特性(或传递函数)的正确性可以由测量的相频曲线来检验,对最小相位系统而言,实际测量所得的相频曲线必须与由确定的频率特牲(或传递函数)所画出的理论相频曲线在一定程度上相符,如果测量所得的相位在高频 (相 对于转角频率)时不等于-90 ° (q —p)[式中p和q分别表示传递函数分子和分母的阶次], 那么,频率特性(或传递函数)必定是一个非最小相位系统的频率特性。 2.被测系统的模拟电路图:见图4-2 图4-2被测系统 二、实验内容 (1)将U21 DAC单元的OUT端接到对象的输入端。 ⑵将测量单元的CH1 (必须拨为乘I档)接至对象的输出端。 ⑶将Ul SG单元的ST和S端断开,用排线将ST端接至U26控制信号单元中的PB0。(由于在每次测量前,应对对象进行一次回零操作,ST即为对象锁零控制端,在这里,我们用8255的PB0 口对ST进行程序控制) ⑷在PC机上分别输入角频率为1, 10,100,300,并使用“ +”、“―”键选择合适的幅值,按ENTER键后,输入的角频率开始闪烁,直至测量完毕时停止,屏幕即显示所测对象的输出及信号源,移动游标,可得到相应的幅值和相位,得到的实验波形图如图4-3到图4-10所示: 图4-3输入频率为1的波形图1

使用Multisim进行电路频率特性分析

使用Multisim进行电路频率响应分析 作者:XChuda Multisim的AC Analysis功能用于对电路中一个或多个节点的电压/电流频响特性进行分析,画出伯德图。本文基于Multisim 11.0。 1、实验电路 本例使用如图的运放电路进行试验。该放大电路采用同相输入,具有(1+100/20=)6倍的放大倍数,带300欧负载。方框部分象征信号源,以理想电压源串联电阻构成。 请不要纠结于我把120Vrms的电压源输入双15V供电的运放这样的举动是否犯二,电压源在AC Analyses中仅仅是作为一个信号入口的标识,其信号类型、幅值和频率对分析是没有贡献的,但是它的存在必不可少,否则无法得到仿真结果! 2、操作步骤 搭好上述电路后,就可以进行交流分析了。

一般设置Frequency parameters和Output两页即可,没有特殊要求的话其他选项保持默认,然后点Simulate开始仿真。切记是点Simulate,点OK的话啥都不会发生。

按照上述步骤仿真结果如下: 分析结果是一份伯德图。在上下两个图表各自区域上按右键弹出列表有若干选项,各位可自己动手试试。右键菜单中的Properties可打开属性对话框,对图表进行更为详细的设置。 3、加个电容试试 从上面伯德图分析结果看出,该电路具有高通特性,是由输入耦合电容C3造成的。现在在输入端加入一个退耦电容试试。电路如下:

在输入端加入220pF退耦电容后C1与后面的放大电路输入电阻构成低通滤波器,可滤除高频干扰。加入C1后,放大电路的输出应该具有带通特性。用AC Analysis分析加入C1后的电路频响特性: 奇怪,为什么高通不见了?一阵疑惑,我甚至动笔算了同相输入端的阻容网络复频域的特性,无论C1是否加入,从同相输入端向左看出去的阻容电路都有一个横轴为0的零点,所以幅度特性应该是从0Hz处开始上升的!对,从0Hz开始!回头看看电路加入C1前仿真的伯德图,发现竖轴范围是13dB~13.3dB! 我们尝试放大来看看。现在重新进行AC分析,将频率范围设置为0.1~10Hz,结果如下图。OK,没问题,果然是高通的,只是截止频率非常低(0.3Hz左右),刚才的仿真频率范围从1Hz开始,自然是看不到的。从中也看出,图表中数字后加小写m,是毫赫兹(mHz)的意思,而不是兆赫兹(MHz)。

频域分析实验报告

频域分析实验报告 班级: 学号: 姓名:

一、实验内容: 1利用计算机作出开环系统的波特图; 2、观察记录控制系统的开环频率特性; 3、控制系统的开环频率特性分析。 二、仿真原理: 对数频率特性图(波特图): 对数频率特性图包括了对数幅频特性图和对数相频特性图。横坐标为频率w,采用对数分度,单位为弧度/秒;纵坐标均匀分度,分别为幅值函数20lgA(w),以dB表示;相角,以度表示。MATLAB提供了函数bode()来绘制系统的波特图,其用法如下: (1)bode(num,den):可绘制出以连续时间多项式传递函数表示的系统的波特图。 (2)当带输出变量[mag,pha,w]或[mag,pha]引用函数时,可得到系统波特图相应的幅值mag、相角pha及角频率点w矢量或只是返回幅值与相角。相角以度为单位,幅值可转换为分贝单位:magdb=20×log10(mag) 二、实验验证 1、用Matlab作Bode图。要求:画出对应Bode图。 (1)G(S)=25/S2+4s+25 (7)G(S)=9(s2+0.2s+1)/s(s2+1.2s+9);

图 1 图 2 (1)G(S)=25/S2+4s+25 可以看成是一个比例环节和一个振荡环节组成,所以k=1,T1=0.04,因为v=0,所以在转折频率之前都为20lgk,因为k=1所以斜率为0,经过转折频率,分段直线斜率的变化量为-40db/dec。

(7)G(S)=9(s2+0.2s+1)/s(s2+1.2s+9); 可以看成是一个二阶微分环节和一个积分环节和一个振荡环节组成,化常数为1后,v=1,t1=1,t2=1/3,所以我们可以看到,在起始阶段是-20*vdb/dec,所以一开始斜率为-20db/dec。当经过1/3的转折频率之后分段直线的改变量为40db/dec,当经过1的转折频率之后分段直线的改变量为-40db/dec。故图像如图所示。 第二题: 典型二阶系统Gs=Wn2/s2+2ζWns+Wn2,试绘制取不同值时的Bode图。取Wn=8,ζ=0.1,0.2,0.3,,0.5,0.6; 图 3 如图所示。

离散时间系统的时域分析

第七章离散时间系统的时域分析 §7-1 概述 一、离散时间信号与离散时间系统 离散时间信号:只在某些离散的时间点上有值的 信号。 离散时间系统:处理离散时间信号的系统。 混合时间系统:既处理离散时间信号,又处理连 续时间信号的系统。 二、连续信号与离散信号 连续信号可以转换成离散信号,从而可以用离散时间系统(或数字信号处理系统)进行处理: 三、离散信号的表示方法:

1、 时间函数:f(k)<——f(kT),其中k 为序号,相当于时间。 例如:)1.0sin()(k k f = 2、 (有序)数列:将离散信号的数值按顺序排列起来。例如: f(k)={1,0.5,0.25,0.125,……,} 时间函数可以表达任意长(可能是无限长)的离散信号,可以表达单边或双边信号,但是在很多情况下难于得到;数列的方法表示比较简单,直观,但是只能表示有始、有限长度的信号。 四、典型的离散时间信号 1、 单位样值函数:? ??==其它001)(k k δ 下图表示了)(n k ?δ的波形。

这个函数与连续时间信号中的冲激函数 )(t δ相似,也有着与其相似的性质。例如: )()0()()(k f k k f δδ=, )()()()(000k k k f k k k f ?=?δδ。 2、 单位阶跃函数:? ??≥=其它001)(k k ε 这个函数与连续时间信号中的阶跃函数)(t ε相似。用它可以产生(或表示)单边信号(这里称为单边序列)。 3、 单边指数序列:)(k a k ε

比较:单边连续指数信号:)()()(t e t e t a at εε=,其 底一定大于零,不会出现负数。 (a) 0.9a = (d) 0.9a =? (b) 1a = (e) 1a =? (c) 1.1a = (f) 1.1a =?

实验六 离散时间系统的时域分析

信号与系统实验报告 实验名:离散时间信号与系统的频域分析 实验六离散时间系统的时域分析 一、实验目的 1、掌握离散时间信号与系统的频域分析方法,从频域的角度对信号与系统的特性进行分析。 2、掌握离散时间信号傅里叶变换与傅里叶逆变换的实现方法。 3、掌握离散时间傅里叶变换的特点及应用 4、掌握离散时间傅里叶变换的数值计算方法及绘制信号频谱的方法 二、预习内容 1、离散时间信号的傅里叶变换与逆变换。 2、离散时间信号频谱的物理含义。 3、离散时间系统的频率特性。 4、离散时间系统的频域分析方法。 三、实验原理 1. 离散时间系统的频率特性

2. 离散时间信号傅里叶变换的数值计算方法 3.涉及到的Matlab 函数

四、实验内容 1、离散时间系统的时域分析 1 离散时间傅里叶变换 (1)下面参考程序是如下序列在范围?4π≤ω≤ 4π的离散时间傅里叶变换 %计算离散时间傅里叶变换的频率样本 clear all; w=-4*pi:8*pi/511:4*pi; num=[2 1]; den=[1 -0.6]; h=freqz(num,den,w); subplot(2,1,1)

plot(w/pi,real(h)); grid; title(‘实部’) xlabel(‘omega/\pi’); yl abel(‘振幅’); subplot(2,1,2) plot(w/pi, imag(h)); grid; title(‘虚部’) xlabel(‘omega/\pi’); ylabel(‘振幅’); figure; subplot(2,1,1) plot(w/pi, abs(h)); grid; title(‘幅度谱’) xlabel(‘omega/\pi’); ylabel(‘振幅’); subplot(2,1,2) plot(w/pi, angle (h)); grid; title(‘相位谱’) x label(‘omega/\pi’); ylabel(‘以弧度为单位的相位’);

时间序列分析实验报告汇总.doc

《时间序列分析》课程实验报告

一、上机练习(P124) 1.拟合线性趋势 12.79 14.02 12.92 18.27 21.22 18.81 25.73 26.27 26.75 28.73 31.71 33.95 程序: data xiti1; input x@@; t=_n_; cards; 12.79 14.02 12.92 18.27 21.22 18.81 25.73 26.27 26.75 28.73 31.71 33.95 ; proc gplot data=xiti1; plot x*t; symbol c=red v=star i=join; run; proc autoreg data=xiti1; model x=t; output predicted=xhat out=out; run; proc gplot data=out; plot x*t=1 xhat*t=2/overlay; symbol2c=green v=star i=join; run; 运行结果:

分析:上图为该序列的时序图,可以看出其具有明显的线性递增趋势,故使用线性模型进行拟合:x t=a+bt+I t,t=1,2,3,…,12 分析:上图为拟合模型的参数估计值,其中a=9.7086,b=1.9829,它们的检验P值均小于 0.0001,即小于显著性水平0.05,拒绝原假设,故其参数均显著。从而所拟合模型为: x t=9.7086+1.9829t.

分析:上图中绿色的线段为线性趋势拟合线,可以看出其与原数据基本吻合。 2.拟合非线性趋势 1.85 7.48 14.29 23.02 37.42 74.27 140.72 265.81 528.23 1040.27 2064.25 4113.73 8212.21 16405.95 程序: data xiti2; input x@@; t=_n_; cards; 1.85 7.48 14.29 23.02 37.42 74.27 140.72 265.81 528.23 1040.27 2064.25 4113.73 8212.21 16405.95 ; proc gplot data=xiti2; plot x*t; symbol c=red v=star i=none; run; proc nlin method=gauss; model x=a*b**t; parameters a=0.1 b=1.1; der.a=b**t; der.b=a*t*b**(t-1); output predicted=xh out=out; run; proc gplot data=out; plot x*t=1 xh*t=2/overlay;

控制系统的频率特性分析

实验六 控制系统的频率特性分析 1.已知系统传递函数为:1 2.01)(+=s s G ,要求: (1) 使用simulink 进行仿真,改变正弦输入信号的频率,用示波器观察输 出信号,记录不同频率下输出信号与输入信号的幅值比和相位差,即 可得到系统的幅相频率特性。 F=10时 输入: 输出:

F=50时 输入:输出: (2)使用Matlab函数bode()绘制系统的对数频率特性曲线(即bode图)。 提示:a)函数bode()用来绘制系统的bode图,调用格式为: bode(sys) 其中sys为系统开环传递函数模型。 参考程序: s=tf(‘s’); %用符号表示法表示s G=1/(0.2*s+1); %定义系统开环传递函数 bode(G) %绘制系统开环对数频率特性曲线(bode图)

实验七连续系统串联校正 一.实验目的 1.加深理解串联校正装置对系统动态性能的校正作用。 2. 对给定系统进行串联校正设计,并通过matlab实验检验设计的正确性。二.实验内容 1.串联超前校正 系统设计要求见课本例题6-3,要求设计合理的超前校正环节,并完成以下内容用matlab画出系统校正前后的阶跃相应,并记录系统校正前后的超调量及调节时间 num=10; 1)figure(1) 2)hold on

3)figure(1) 4)den1=[1 1 0]; 5)Gs1=tf(num,den1); 6)G1=feedback(Gs1,1,-1); 7)Step(G1) 8) 9)k=10; 10)figure(2) 11)GO=tf([10],[1,1,0]); 12)Gc=tf([0.456,1],[1,00114]); 13)G=series(G0,Gc); 14)G1=feedback(G,1); 15)step(G1);grid

电路实验__电路频率特性的研究要点说明

东南大学电工电子实验中心 实验报告 课程名称:电路实验 第二次实验 实验名称:电路频率特性的研究 院(系):仪器科学与工程学院专业: 姓名:学号: 实验室: 实验组别: 同组人员:实验时间: 评定成绩:审阅教师:

电路频率特性的研究 一、 实验目的 1. 掌握低通、带通电路的频率特性; 2. 应用Multisim 软件测试低通、带通电路频率特性及有关参数; 3. 应用Multisim 软件中的波特仪测试电路的频率特性。 二、 实验原理 研究电路的频率特性,即是分析研究不同频率的信号作用于电路所产生的响应函数与激励函数的比值关系。通常情况下,研究具体电路的频率特性,并不需要测试构成电路所有元件上的响应与激励之间的关系,只需要研究由工作目的所决定的某个元件或支路的响应与激励之间的关系。本实验主要研究一阶RC 低通电路,二阶RLC 低通、带通电路的频率特性。 (一):网络频率特性的定义 电路在一个正弦电源激励下稳定时,各部分的响应都是同频率的正弦量,通过正弦量的相量,网络函数|()|H jw 定义为:. ().|()||()|j w Y H w H jw e X ?== 其中Y 为输出端口的响应,X 为输入端口的激励。由上式可知,网络函数是频率的函数,其中网络函数的模|()|H jw 与频率的关系称为幅频特性,网络函数的相角()w ?与频率的关系称为相频特性,后者表示了响应与激励的相位差与频率的关系。 (二):网络频率特性曲线 1. 一阶RC 低通网络 网络函数: 其模为: 辐角为: 显然,随着频率的增高,|H(j ω)|将减小, 即响应与激励的比值减小,这说明低频信 4590 (a) RC低通网络(b) 幅频特性 (c) 相频特性 ()H j ω()) RC ?ω=().0.1/1 1/1i U j c H j R j C j RC U ωωωω=== ++

spss时间序列模型

《统计软件实验报告》SPSS软件的上机实践应用 时间序列分析

数学与统计学学院 一、实验内容: 时间序列是指一个依时间顺序做成的观察资料的集合。时间序列分析过程中最常用的方法是:指数平滑、自回归、综合移动平均及季节分解。 本次实验研究就业理论中的就业人口总量问题。但人口经济的理论和实践表明,就业总量往往受到许多因素的制约,这些因素之间有着错综复杂的联系,因此,运用结构性的因果模型分析和预测就业总量往往是比较困难的。时间序列分析中的自回归求积分移动平均法(ARIMA)则是一个较好的选择。对于时间序列的短期预测来说,随机时序ARIMA是一种精度较高的模型。 我们已辽宁省历年(1969-2005)从业人员人数为数据基础建立一个就业总量的预测时间序列模型,通过spss建立模型并用此模型来预测就业总量的未来发展趋势。 二、实验目的: 1.准确理解时间序列分析的方法原理 2.学会实用SPSS建立时间序列变量 3.学会使用SPSS绘制时间序列图以反应时间序列的直观特征。

4.掌握时间序列模型的平稳化方法。 5.掌握时间序列模型的定阶方法。 6.学会使用SPSS建立时间序列模型与短期预测。 7.培养运用时间序列分析方法解决身边实际问题的能力。 三、实验分析: 总体分析: 先对数据进行必要的预处理和观察,直到它变成稳态后再用SPSS对数据进行分析。 数据的预处理阶段,将它分为三个步骤:首先,对有缺失值的数据进行修补,其次将数据资料定义为相应的时间序列,最后对时间序列数据的平稳性进行计算观察。 数据分析和建模阶段:根据时间序列的特征和分析的要求,选择恰当的模型进行数据建模和分析。 四、实验步骤: SPSS的数据准备包括数据文件的建立、时间定义和数据期间的指定。 SPSS的时间定义功能用来将数据编辑窗口中的一个或多个变量指定为时间序列变量,并给它们赋予相应的时间标志,具体操作步骤是: 1.选择菜单:Date→Define Dates,出现窗口:

一阶RC电路频率特性的研究实验报告

北京交通大学电子信息工程学院2011~2012 实验报告 实验题目:一阶RC电路频率特性的研究。 实验内容及结果: 1.低通电路的研究 实验电路: 实验数据: 低通电路数据 频率电平频率 a b 相位差100 -0.1 100 0.6 4.8 7.1808 300 -0.9 200 1.2 4.4 15.8266 475 -2 300 1.9 4.2 26.8965 560 -2.5 400 2.2 4 33.367 641 -3 500 2.4 3.6 41.8103 704 -3.5 600 2.4 3.2 51.3752 788 -4 700 2.4 3 53.1301 849 -4.5 800 2.4 2.8 58.9973 926 -5 1000 2 2.4 56.4427 1000 -5.5 2000 3.2 3.4 70.2501 1072 -6 3000 2.1 2.2 72.6586 1149 -6.5 5000 1.2 1.2 90 1240 -7 1340 -7.5 1430 -8 1520 -8.5 1600 -9 1660 -9.5 1860 -10 2400 -12 3040 -14 3780 -16 4700 -18 5000 -19

电平图: 相位差图:

2.高通电路研究 实验数据: 高通电路数据 频率电平频率 a b 相位差5000 0 5000 0.6 4.6 7.4947 1500 -0.4 4000 0.8 4.6 10.0154 1200 -0.8 3000 1 4.5 12.8396 1030 -1 2000 1.4 4.5 18.1262 899 -1.4 1000 2 4 30 740 -2 600 2.4 3.2 48.5904 663 -2.4 500 2.4 3.1 50.732 588 -3 400 2.1 2.4 61.045 532 -3.5 300 2 2.1 72.2472 481 -4 100 0.8 0.8 90 440 -4.5 400 -5 372 -5.5 344 -6 339 -6.5 325 -7 261 -8 227 -9 200 -10 165 -11 100 -16

相关文档
相关文档 最新文档