CodeViz —— 一款分析C_C++源代码中函数调用关系的调用图生成工具
通过函数绘制一阶二阶传递函数伯德图
关于一阶二阶传递函数的伯德图 一阶惯性系统的通式为： 将式子两边同时除以a0得 令0 0a K b =为系统静态灵敏度； 0 1a a =τ为系统时间常数； 则有 )()()1( s KX s Y s =+τ 故有 ) 1()()()(+==s K s X s Y s H τ 以液柱式温度计为例，传递函数为 )1(1)()()(+==s s X s Y s H τ 可得频率响应函数 )1j (1)(+= τωs H )()()(001t x b t y a dt t dy a =+)()()(0001t x a b t y dt t dy a a =+
可得传递函数的幅频与相频特性 2)1(1 )()(τωωω+==j H A ωτωω?arctan )()(-=∠=j H 在MATLAB 上输入程序（此时令1=τ） num=[1]; den=[1,1]; figure sys=tf(num,den); bode(sys);grid on 可得bode 图
二阶惯性系统的通式为： 将式子两边同时除以a 0得 令0 0a K b =为系统静态灵敏度； 20n a a = ω为系统无阻尼固有频率； 1 012a a a =ξ为系统阻尼器 传递函数为 12) ()()(22++==n n s s K s X s Y s H ωξω 可得传递函数的幅频与相频特性 2222)(4)1(1 )()(2n n K j H A ωωξωωωω+-== )()()()(001222t x b t y a dt t dy a dt t y d a =++)()()()(00012202t x a b t y dt t dy a a dt t y d a a =++
5-Excel中三个查找引用函数的用法(十分有用)
Excel中三个查找引用函数的用法(十分有用) 在Excel中，我们经常会需要从某些工作表中查询有关的数据复制到另一个工作表中。比如我们需要把学生几次考试成绩从不同的工作表中汇总到一个新的工作表中，而这几个工作表中的参考人数及排列顺序是不完全相同的，并不能直接复制粘贴。此时，如果使用Excel的VLOOKUP、INDEX或者OFFSET函数就可以使这个问题变得非常简单。我们以Excel 2007为例。 图1 假定各成绩工作表如图1所示。B列为姓名，需要汇总的项目“总分”及“名次”位于H列和I列（即从B列开始的第7列和第8列）。而汇总表则如图2所示，A 列为姓名列，C、D两列分别为要汇总过来的第一次考试成绩的总分和名次。其它各次成绩依次向后排列。 图2
一、VLOOKUP函数 我们可以在“综合”工作表的C3单元格输入公式“=VLOOKUP($B3,第1次!$B$1:$I$92,7,FALSE)”，回车后就可以将第一位同学第一次考试的总分汇总过来了。 把C3单元格公式复制到D3单元格，并将公式中第三个参数“7”改成“8”，回车后，就可以得到该同学第一次考试名次。 选中C3:D3这两个单元格，向下拖动填充句柄到最后就可以得到全部同学的总分及名次了。是不是很简单呀？如图3所示。 VLOOKUP函数的用法是这样的：VLOOKUP(参数1,参数2,参数3,参数4)。“参数1”是“要查找谁？”本例中B3单元格，那就是要查找B3单元格中显示的人名。“参数2”是“在哪里查找？”本例中“第1次!$B$1:$I$92”就是告诉Excel在“第1次”工作表的B1:I92单元格区域进行查找。“参数3”是“找第几列的数据？”本例中的“7”就是指从“第1次”工作表的B列开始起，第7列的数据，即H列。本例中“参数4”即“FALSE”是指查询方式为只查询精确匹配值。 该公式先在“第1次”工作表的B!:I92单元格区域的第一列（即B1:B92单元格区域）查找B3单元格数据，找到后，返回该数据所在行从B列起第7列（H列）的数据。所以，将参数3改成“8”以后，则可以返回I列的数据。 由此可以看出，使用VLOOKUP函数时，参数1的数据必须在参数2区域的第一列中。否则是不可以查找的。
求下图所示系统的传递函数
一、求下图所示系统的传递函数)(/)(0s U s U i 。 （10分） ) 1()()(3132320+++-=CS R R R R CS R R s U s U i 一、控制系统方块图如图所示： （1）当a =0时，求系统的阻尼比ξ，无阻尼自振频率n ω和单位斜坡函数输入时的稳态误差； （2）当ξ=时，试确定系统中的a 值和单位斜坡函数输入时系统的稳态误差； 系统的开环传函为 s a s s G )82(8)(2++=闭环传函为8)82(8)()(2+++=s a s s R s Y 25.0 83.2 36.0===ss n e ωξ 4 25.0==ss e a 设某控制系统的开环传递函数为 ) 22()(2++=s s s k s G 试绘制参量k 由0变至∞时的根轨迹图，并求开环增益临界值。 （15分） 1）j p j p p --=+-==110 321 2）πππ?σ3 5,,332=-=a a （10分） 3）ω=j 2±，c k =4，开环增益临界值为K=2 设某系统的特征方程为23)(234+--+=s s s s s D ，试求该系统的特征根。 列劳斯表如下 0000220112311 2 3 4 s s s s --- （4分） 得 辅助方程为0222=+-s ，解得1,121-==s s （4分）
最后得1,243=-=s s 设某控制系统的开环传递函数为 )()(s H s G =) 10016()12.0(752+++s s s s 试绘制该系统的Bode 图，并确定剪切频率c ω的值 剪切频率为s rad c /75.0=ω 某系统的结构图和Nyquist 图如图(a)和(b)所示，图中 2)1(1)(+=s s s G 23 ) 1()(+=s s s H 试判断闭环系统稳定性，并决定闭环特征方程正实部根的个数。 (16分) 解：由系统方框图求得内环传递函数为： s s s s s s s H s G s G +++++=+23452 474)1()()(1)( (3分) 内环的特征方程：04742345=++++s s s s s (1 分) 由Routh 稳定判据： 01: 03 10 :16 :044: 171: 01234s s s s s 七、设某二阶非线性系统方框图如图所示，其中 4 , 2.0 , 2.00===K M e 及s T 1=， 试画出输入信号)(12)(t t r ?=时系统相轨迹的大致图形，设系统原处于静止状态。 （16分） 解：根据饱和非线性特性，相平面可分成三个区域，运动方程分别为
Matlab中Bode图的绘制技巧(精)
Matlab中Bode图的绘制技巧 我们经常会遇到使用Matlab画伯德图的情况，可能我们我们都知道bode这个函数是用来画bode图的，这个函数是Matlab内部提供的一个函数，我们可以很方便的用它来画伯德图，但是对于初学者来说，可能用起来就没有那么方便了。 譬如我们要画出下面这个传递函数的伯德图： 1.576e010 s^2 H(s= ------------------------------------------------------------------------------------------ s^4 + 1.775e005 s^3 + 1.579e010 s^2 + 2.804e012 s + 2.494e014 (这是一个用butter函数产生的2阶的，频率范围为[20 20K]HZ的带通滤波器。 我们可以用下面的语句： num=[1.576e010 0 0]; den=[1 1.775e005 1.579e010 2.804e012 2.494e014]; H=tf(num,den; bode(H 这样，我们就可以得到以下的伯德图： 可能我们会对这个图很不满意，第一，它的横坐标是rad/s，而我们一般希望横坐标是HZ；第二，横坐标的范围让我们看起来很不爽；第三，网格没有打开（这点当然我们可以通过在后面加上grid on解决）。 下面，我们来看看如何定制我们自己的伯德图风格： 在命令窗口中输入：bodeoptions
我们可以看到以下内容：ans = Title: [1x1 struct] XLabel: [1x1 struct] YLabel: [1x1 struct] TickLabel: [1x1 struct] Grid: 'off' XLim: {[1 10]} XLimMode: {'auto'} YLim: {[1 10]} YLimMode: {'auto'} IOGrouping: 'none' InputLabels: [1x1 struct] OutputLabels: [1x1 struct] InputVisible: {'on'} OutputVisible: {'on'} FreqUnits: 'rad/sec' FreqScale: 'log' MagUnits: 'dB' MagScale: 'linear' MagVisible: 'on' MagLowerLimMode: 'auto' MagLowerLim: 0 PhaseUnits: 'deg' PhaseVisible: 'on' PhaseWrapping: 'off'
C++函数调用过程深入分析
函数调用过程分析 刘兵QQ: 44452114 E-mail: liubing2000@https://www.wendangku.net/doc/8c8156275.html, 0. 引言 函数调用的过程实际上也就是一个中断的过程，那么C++中到底是怎样实现一个函数的调用的呢？参数入栈、函数跳转、保护现场、回复现场等又是怎样实现的呢？本文将对函数调用的过程进行深入的分析和详细解释，并在VC 6.0环境下进行演示。分析不到位或者存在错误的地方请批评指正，请与作者联系。 首先对三个常用的寄存器做一下说明，EIP是指令指针，即指向下一条即将执行的指令的地址；EBP 为基址指针，常用来指向栈底；ESP为栈指针，常用来指向栈顶。 看下面这个简单的程序并在VC 6.0中查看并分析汇编代码。 图1 1. 函数调用 g_func函数调用的汇编代码如图2： 图2 首先是三条push指令，分别将三个参数压入栈中，可以发现参数的压栈顺序是从右向左的。这时我们可以查看栈中的数据验证一下。如图3所示，从右边的实时寄存器表中我们可以看到ESP（栈顶指针）值为0x0012FEF0，然后从中间的内存表中找到内存地址0x0012FEF0处，我们可以看到内存中依次存储了0x00000001（即参数1），0x00000002（即参数2），0x00000003（即参数3），即此时栈顶存储的是三个参数值，说明压栈成功。 图3
然后可以看到call指令跳转到地址0x00401005，那么该地址处是什么呢？我们继续跟踪一下，在图4中我们看到这里又是一条跳转指令，跳转到0x00401030。我们再看一下地址0x00401030处，在图5中可以看到这才是真正的g_func函数，0x00401030是该函数的起始地址，这样就实现了到g_func函数的跳转。 图4 图5 2. 保存现场 此时我们再来查看一下栈中的数据，如图6所示，此时的ESP（栈顶）值为0x0012FEEC，在内存表中我们可以看到栈顶存放的是0x00401093，下面还是前面压栈的参数1,2,3，也就是执行了call指令后，系统默认的往栈中压入了一个数据（0x00401093），那么它究竟是什么呢？我们再看到图3，call指令后面一条指令的地址就是0x00401093，实际上就是函数调用结束后需要继续执行的指令地址，函数返回后会跳转到该地址。这也就是我们常说的函数中断前的“保护现场”。这一过程是编译器隐含完成的，实际上是将EIP（指令指针）压栈，即隐含执行了一条push eip指令，在中断函数返回时再从栈中弹出该值到EIP，程序继续往下执行。 图6 继续往下看，进入g_func函数后的第一条指令是push ebp，即将ebp入栈。因为每一个函数都有自己的栈区域，所以栈基址也是不一样的。现在进入了一个中断函数，函数执行过程中也需要ebp寄存器，而在进入函数之前的main函数的ebp值怎么办呢？为了不被覆盖，将它压入栈中保存。 下一条mov ebp, esp 将此时的栈顶地址作为该函数的栈基址，确定g_func函数的栈区域（ebp为栈底，
典型环节的Bode图
控制系统的开环频率特性 目的：掌握开环Bode 图的绘制 根据Bode 图确定最小相位系统的传递函数 重点：开环Bode 图的绘制、根据Bode 图确定最小相位系统的传递函数 1 开环伯德图手工作图的一般步骤： 1）将开环传递函数表示为时间常数表达形式，计算各个典型环节的交接频率 2）求20lgK 的值，并明确积分环节的个数ν 3）通过（1,20lgK ）绘制斜率为-20vdB/dec 低频段 4）随着频率增加，每遇到一个典型环节的交接频率，就改变一次斜率 最小相位系统定义： 递函数的零点、极点全部位于S 左半平面，同时又无纯滞后环节的系统称为最小相位系统。否则就是非最小相位系统。 对数幅频特性与相频特性之间存在确定的对应关系。对于一个最小相位系统，我们若知道了其幅频特性，它的相频特性也就唯一地确定了。也就是说：只要知道其幅频特性，就能写出此最小相位系统所对应的传递函数，而无需再画出相频特性。 非最小相位系统高频时相角迟后大，起动性能差，响应缓慢。对响应要求快的系统，不宜采用非最小相位元件。 2 典型环节的伯德图 绘制曲线在MA TLAB 中实现，利用下述的程序段： num=[b2 b1 b0]; den=[1 a2 a1 a0]; H=tf(num,den); bode(H) margin(H) hold on 2.1 比例环节 传递函数：()G s K = 频率特性：()G j K ω= 对数幅频特性：()20lg L j K ω= 对数相频特性：()0?ω= 程序段： num=[0 10]; den=[0 1]; H=tf(num,den); bode(H) margin(H) hold on 结论：放大环节的对数幅频特性是一条幅值为20lgK 分贝，且平行于横轴的直线，相频特性是一条和横轴重合的直线。 K>1时，20lgK>0dB ；K<1时，20lgK<0dB 。 2.2 惯性环节（低通滤波特性） 传递函数：1()1G s s τ= + 频率特性：()()()j G j A e ?ωωω= 对数幅频特性：2 1()20lg 1() L ωτω=+ 对数相频特性：()arctan ?ωτω=- 绘制1()10.1G s s =+的Bode 图 程序段： num=[0 1]; den=[0.1 1];H=tf(num,den); bode(H) margin(H) hold on 结论：惯性环节的对数幅频特性可以用在1ωτ= 处相交于0分贝的两条渐近直线来近似表示：当1ωτ 时，是一条0分贝的直线； 当1ωτ 时，是一条斜率为-20dB/dec 的直线。 惯性环节具有低通特性，对低频输入能精确地复现，而对高频输入要衰减，且产生相位迟后。因此，它只能复现定常或缓慢变化的信号。 2.3 积分环节 传递函数：1 ()G s s τ= 频率特性：()()()j G j A e ?ωωω= 对数幅频特性：1 ()20lg L j ωτω = 对数相频特性：()2 π?ω=- 在同一坐标中绘制1()G s s = 、1()0.1G s s = 和 1()0.01G s s = 的Bode 图 num1=[0 1];den1=[1 1];H1=tf(num1,den1); bode(H1)margin(H1)hold on
Excel中三个查找引用函数的用法(十分有用)
在Excel中，我们经常会需要从某些工作表中查询有关的数据复制到另一个工作表中。比如我们需要把学生几次考试成绩从不同的工作表中汇总到一个新的工作表中，而这几个工作表中的参考人数及排列顺序是不完全相同的，并不能直接复制粘贴。此时，如果使用Excel的VLOOKUP、INDEX或者OFFSET函数就可以使这个问题变得非常简单。我们以Excel 2007为例。 图1 假定各成绩工作表如图 1所示。B列为，需要汇总的项目“总分”及“名次”位于H列和I列（即从B列开始的第7列和第8列）。而汇总表则如图2所示，A列为列，C、D两列分别为要汇总过来的第一次考试成绩的总分和名次。其它各次成绩依次向后排列。
图2 一、 VLOOKUP函数 我们可以在“综合”工作表的C3单元格输入公式“=VLOOKUP($B3,第1次!$B$1:$I$92,7,FALSE)”，回车后就可以将第一位同学第一次考试的总分汇总过来了。 把C3单元格公式复制到D3单元格，并将公式中第三个参数“7”改成“8”，回车后，就可以得到该同学第一次考试名次。 选中C3:D3这两个单元格，向下拖动填充句柄到最后就可以得到全部同学的总分及名次了。是不是很简单呀？如图3所示。
VLOOKUP函数的用法是这样的：VLOOKUP(参数1,参数2,参数3,参数4)。“参数1”是“要查找谁？”本例中B3单元格，那就是要查找B3单元格中显示的人名。“参数2”是“在哪里查找？”本例中“第1次!$B$1:$I$92”就是告诉Excel在“第1次”工作表的B1:I92单元格区域进行查找。“参数3”是“找第几列的数据？”本例中的“7”就是指从“第1次”工作表的B列开始起，第7列的数据，即H列。本例中“参数4”即“FALSE”是指查询方式为只查询精确匹配值。 该公式先在“第1次”工作表的B!:I92单元格区域的第一列（即B1:B92单元格区域）查找B3单元格数据，找到后，返回该数据所在行从B列起第7列（H列）的数据。所以，将参数3改成“8”以后，则可以返回I列的数据。 由此可以看出，使用VLOOKUP函数时，参数1的数据必须在参数2区域的第一列中。否则是不可以查找的。 二、INDEX函数 某些情况下，VLOOKUP函数可能会无用武之地，如图4所示。“综合”工作表中，列放到了A 列，而B列要求返回该同学所在的班级。但我们看前面的工作表就知道了，“班级”列是位于“”列前面的。所以，此时我们不可能使用VLOOKUP函数来查找该同学的班级。而INDEX函数就正可以一试身手。
求下图所示系统的传递函数
一、求下图所示系统的传递函数 ) (/)(0s U s U i 。 （10分） ) 1()()(313 2320+++-=CS R R R R CS R R s U s U i 一、控制系统方块图如图所示： （1）当a =0时，求系统的阻尼比ξ，无阻尼自振频率n ω和单位斜坡函数输入时的稳态误差； （2）当ξ=0.7时，试确定系统中的a 值和单位斜坡函数输入时系统的稳态误差； 系统的开环传函为 s a s s G )82(8)(2++= 闭环传函为8)82(8 )()(2 +++=s a s s R s Y 25.0 83.2 36.0===ss n e ωξ 4 25.0==ss e a 设某控制系统的开环传递函数为 ) 22()(2 ++= s s s k s G 试绘制参量k 由0变至∞时的根轨迹图，并求开环增益临界值。 （15分） 1）j p j p p --=+-==110321 2） πππ?σ3 5 ,,332=- =a a （10分） 3）ω=j 2±，c k =4，开环增益临界值为K=2 设某系统的特征方程为23)(2 3 4 +--+=s s s s s D ，试求该系统的特征根。 列劳斯表如下 022******* 2 34 s s s s ---
得辅 助 方 程 为 222=+-s ，解得 1,121-==s s （4分） 最后得1, 243=-=s s 设某控制系统的开环传递函数为 )()(s H s G = ) 10016() 12.0(752+++s s s s 试绘制该系统的Bode 图，并确定剪切频率c ω的值 剪切频率为s rad c /75.0=ω 某系统的结构图和Nyquist 图如图(a)和(b)所示，图中 2)1(1)(+=s s s G 2 3 ) 1()(+=s s s H 试判断闭环系统稳定性，并决定闭环特征方程正实部根的个数。 (16分) 解：由系统方框图求得内环传递函数为： s s s s s s s H s G s G +++++= +23452 474)1()()(1)(
试求图示有源网络的传递函数和Bode图.docx
6-1试求图示有源网络的传递甫数和Bode 图，并说明其网络特性。 6-2已知单位反馈控制系统的开环传递函数为 G(f)二 10 5(0.25 +1) 当串联校正装置的传递函数G c ($)如下所示时: (1) G c (5)= 0.2s +1 0.05s +1 2($ +1) (10s+ 1) 1?试绘出两种校正时校正前和校正后系统Bode 图; 2.试比较两种校正方案的优缺点。 6-3已知单位反馈系统的对数幅频特性Illi 线如图屮厶)@), 串联校正装置G c (s)的对 数幅频特性如图中&9),要求: 1. 在图小画出系统校止后的对数幅频特性厶(e); 2. 写出校正后系统的开环传递函数； 3. 分析校止装置G c (5)对系统的作用。 6-4系统的结构图如图所示，试利用根轨迹法设计超前校止装置，使系统满足下列性 能指标：=0.7 , t s =1.45, K v = 。 6—5已知一单位反馈系统的开环传递函数为 习题6— 1图
试设计一?校正装置，使系统的相角裕量厂> 45° ,剪切频率0. > 50$ j 0 6-6单位反馈系统的开环传递函数为 设计一串联滞后校正装置，使系统相角裕量/ > 40° ,并保持原有的开环增益。 6-7设单位反馈系统的开环传递函数为 G(s)= --------------- ------------ 5(0.15 + 1)(0.255 + 1) 试设计--校正装置，使系统满足下列性能指标，速度误差系数K,, 相角裕量 / > 40° ,剪切频率 > 0.5s~} o 6-8单位反馈系统的开环传递函数为 若耍求校正后系统的谐振峰值=1.4,谐振频率> lor 1，试确定校正装置的形 式与参数。 6-9单位反馈系统的结构如图所示，现用速度反馈来校正系统，校正后系统具有临界 G(s) = 200 5(0.15 + 1) G() = 4 s(2s +1) G(s)= 10 5(0.255 +1)(0.055 +1) 习题6 —3图
C语言函数调用规定
在C语言中，假设我们有这样的一个函数： int function（int a，int b） 调用时只要用result = function（1，2）这样的方式就可以使用这个函数。但是，当高级语言被编译成计算机可以识别的机器码时，有一个问题就凸现出来：在CPU中，计算机没有办法知道一个函数调用需要多少个、什么样的参数，也没有硬件可以保存这些参数。也就是说，计算机不知道怎么给这个函数传递参数，传递参数的工作必须由函数调用者和函数本身来协调。为此，计算机提供了一种被称为栈的数 据结构来支持参数传递。 栈是一种先进后出的数据结构，栈有一个存储区、一个栈顶指针。栈顶指针指向堆栈中第一个可用的数据项（被称为栈顶）。用户可以在栈顶上方向栈中加入数据，这个操作被称为压栈（Push），压栈以后，栈顶自动变成新加入数据项的位置，栈顶指针也随之修改。用户也可以从堆栈中取走栈顶，称为弹出栈（pop），弹出栈后，栈顶下的一个元素变成栈顶， 栈顶指针随之修改。 函数调用时，调用者依次把参数压栈，然后调用函数，函数被调用以后，在堆栈中取得数据，并进行
计算。函数计算结束以后，或者调用者、或者函数本身修改堆栈，使堆栈恢复原装。
在参数传递中，有两个很重要的问题必须得到明 确说明： 当参数个数多于一个时，按照什么顺序把参数压 入堆栈 函数调用后，由谁来把堆栈恢复原装在高级语言中，通过函数调用约定来说明这两个问题。常见的调用约定有： stdcall cdecl fastcall thiscall naked call stdcall调用约定 stdcall很多时候被称为pascal调用约定，因为pascal是早期很常见的一种教学用计算机程序设计语言，其语法严谨，使用的函数调用约定就是stdcall.在Microsoft C++系列的C/C++编译器中，常常用PASCAL宏来声明这个调用约定，类似的宏还有 WINAPI和CALLBACK. stdcall调用约定声明的语法为（以前文的那个函 数为例）：
函数调用过程分析
1. 函数调用过程分析
1. 函数调用
我们用下面的代码来研究函数调用的过程。 例 19.1. 研究函数的调用过程
int bar(int c, int d) { int e = c + d; return e; }
int foo(int a, int b) { return bar(a, b); }
int main(void) { foo(2, 3); return 0; }
如果在编译时加上-g 选项（在第 10 章 gdb 讲过-g 选项），那么用 objdump 反汇编时可以把 C 代 码和汇编代码穿插起来显示，这样 C 代码和汇编代码的对应关系看得更清楚。反汇编的结果很长， 以下只列出我们关心的部分。
$ gcc main.c -g $ objdump -dS a.out ...
08048394
int e = c + d; 804839a: 804839d: 80483a0: 80483a2: 8b 55 0c 8b 45 08 01 d0 89 45 fc mov mov add mov 0xc(%ebp),%edx 0x8(%ebp),%eax %edx,%eax %eax,-0x4(%ebp)
return e; 80483a5: } 80483a8: 80483a9: c9 c3 leave ret 8b 45 fc mov -0x4(%ebp),%eax
080483aa
int foo(int a, int b) { 80483aa: 80483ab: 80483ad: 55 89 e5 83 ec 08 push mov sub %ebp %esp,%ebp $0x8,%esp
return bar(a, b); 80483b0: 80483b3: 80483b7: 8b 45 0c 89 44 24 04 8b 45 08 mov mov mov 0xc(%ebp),%eax %eax,0x4(%esp) 0x8(%ebp),%eax
自动控制原理基础伯德图
使用MATLAB 绘制频率特性曲线 姓名 黄勇 班级 16电气本三 学号 4702160186 一、频率特性 在定义谐波输入下，输出响应中与输入同频率的谐波分量与谐波输入的幅值之比A(ω)为幅频特性，相位之差）（ω?为相频特性，并称其指数表达形式： ()()()j G j A e ?ω ωω= 为系统的频率特性。 总结上述我们可知：频率特性由两个部分组合而成，分别是幅频特性和相频特性。稳态系统的输出信号与输入信号的相位之差我们称其为相频特性。稳态系统输出与输入的幅值之比称为幅频特性。另外频率响应对稳定系统和不稳定系统都适应，其中稳定系统的频率特性可以通过实验的方法确定。 二、频率特性的几何表示法 ? 幅相频率特性曲线 简称幅相特性曲线，或幅相特性，或极坐标图。 ? 对数频率特性曲线 又称为伯德曲线或伯德图。 ? 对数幅相曲线 又称为尼科尔斯曲线或尼科尔斯图。 三、惯性环节频率特性的绘制 惯性环节的表达式为： () 1 1G s Ts = + T 的取值分别为2、4、7，使用MATLAB 软件绘制
MATLABA的函数指令如下： 指令说明：num为分子指令；den为分母指令；此次画图调用了伯德图画法（bode指令）。绘制图如下：
T=2时。 MATLABA的函数指令如下： 绘制图如下： 同理当T=4时。MATLABA的函数指令如下：
绘制图如下： 四、振荡环节频率特性的绘制 振荡环节的传递函数为： ()22 1=21 n n G s s s ζωω+ +
在 2 01取值，本次取值分别为0.1 0.3 0.5 0.707 0.85 0.91 1。 方法一：使用伯德图画MATLAB函数程序指令如下： MATLAB图形显示如下：
函数调用
杨振平
●函数定义后，并不能自动执行，必须通过函数调用来实现函数的功能。 ●函数调用，即控制执行某个函数。 ●C++中，主函数可以调用其它子函数，而其它函数之间也可以相互调用。 ●在本节中，我们将介绍一下内容： ?函数调用的格式 ?参数的传递方式 ?为形参指定默认值 ?数组名作函数参数 ?结构体变量作函数参数
函数调用的一般格式： <函数名>（<实际参数表>）//有参调用 或<函数名>（）//无参调用 其中： ●<函数名>为要使用的函数的名字。 ●<实际参数表>是以逗号分隔的实参列表，必须放在一对圆括号中。 <实参表>与<形参表>中参数的个数、类型和次序应保持一致。 ●当调用无参函数时，函数名后的圆括号不能省略。
1.实参的几种形式 ●形参为简单类型变量，对应的实参可以是：常量，变量及表达式。 ●形参为数组，对应的实参为数组（名）。 ●形参为结构类型，对应的实参为结构类型变量。 如：调用已知三边求三角形面积的函数Area。 double Area(double,double,double); //函数声明 cout< 2.函数调用的形式 (1)函数调用作为一个独立的语句（用于无返回值的函数）调用的形式为： 函数名（实参表）；或函数名（）； 如：调用print_char函数（用户定义的无返回值函数）。 print_char(‘*’,6); //连续显示6个‘*’字符。 分析函数调用关系图(call graph)的几种方法 绘制函数调用关系图对理解大型程序大有帮助。我想大家都有过一边读源码（并在头脑中维护一个调用栈），一边在纸上画函数调用关系，然后整理成图的经历。如果运气好一点，借助调试器的单步跟踪功能和call stack窗口，能节约一些脑力。不过如果要分析的是脚本语言的代码，那多半只好老老实实用第一种方法了。如果在读代码之前，手边就有一份调用图，岂不妙哉？下面举出我知道的几种免费的分析C/C++函数调用关系的工具。 函数调用关系图（call graph）是图（graph），而且是有向图，多半还是无环图（无圈图）——如果代码中没有直接或间接的递归的话。Graphviz是专门绘制有向图和无向图的工具，所以很多call graph分析工具都以它为后端(back end)。那么前端呢？就看各家各显神通了。 调用图的分析分析大致可分为“静态”和“动态”两种，所谓静态分析是指在不运行待分析的程序的前提下进行分析，那么动态分析自然就是记录程序实际运行时的函数调用情况了。 静态分析又有两种方法，一是分析源码，二是分析编译后的目标文件。 分析源码获得的调用图的质量取决于分析工具对编程语言的理解程度，比如能不能找出正确的C++重载函数。Doxygen是源码文档化工具，也能绘制调用图，它似乎是自己分析源码获得函数调用关系的。GNU cflow也是类似的工具，不过它似乎偏重分析流程图（flowchart）。 对编程语言的理解程度最好的当然是编译器了，所以有人想出给编译器打补丁，让它在编译时顺便记录函数调用关系。CodeViz（其灵感来自Martin Devera (Devik) 的工具）就属于此类，它（1.0.9版）给GCC 3.4.1打了个补丁。另外一个工具egypt的思路更巧妙，不用大动干戈地给编译器打补丁，而是让编译器自己dump出调用关系，然后分析分析，交给Graphviz去绘图。不过也有人另起炉灶，自己写个C语言编译器（ncc），专门分析调用图，勇气可嘉。不如要是对C++语言也这么干，成本不免太高了。分析C++的调用图，还是借助编译器比较实在。 分析目标文件听起来挺高深，其实不然，反汇编的工作交给binutils的objdump 去做，只要分析一下反汇编出来的文本文件就行了。下面是Cygwin下objdump -d a.exe的部分结果： 00401050 <_main>: 401050: 55 p 习 题 6－1 试求图示有源网络的传递函数和Bode 图，并说明其网络特性。 6—2 已知单位反馈控制系统的开环传递函数为 ) 12.0(10)(+=s s s G 当串联校正装置的传递函数)(s G c 如下所示时： （1）1 05.012.0)(++=s s s G c (2))110()1(2)(++=s s s G c 1．试绘出两种校正时校正前和校正后系统Bode 图； 2 6—3 已知单位反馈系统的对数幅频特性曲线如图中)(0ωL ，串联校正装置)(s G c 的对数幅频特性如图中)(ωc L ，要求： 1．在图中画出系统校正后的对数幅频特性)(ωL ； 2 3．分析校正装置)(s G c 对系统的作用。 6—4系统的结构图如图所示，试利用根轨迹法设计超前校正装置，使系统满足下列性 能指标7.0=ζ，s t s 4.1=，12-=s K v 。 6—5 已知一单位反馈系统的开环传递函数为 ) 11.0(200)(+= s s s G 试设计一校正装置，使系统的相角裕量?≥45γ，剪切频率150-≥s c ω。 6—6 单位反馈系统的开环传递函数为 ) 12(4)(+=s s s G c 设计一串联滞后校正装置，使系统相角裕量?≥40γ，并保持原有的开环增益。 6—7 设单位反馈系统的开环传递函数为 ) 125.0)(11.0(5)(++=s s s s G 试设计一校正装置，使系统满足下列性能指标，速度误差系数15-=s K v ，相角裕量 ?≥40γ，剪切频率15.0-≥s c ω。 6—8 单位反馈系统的开环传递函数为 ) 105.0)(125.0(10)(++=s s s s G 若要求校正后系统的谐振峰值4.1=r M ，谐振频率110-≥s r ω，试确定校正装置的形 式与参数。 6—9 单位反馈系统的结构如图所示，现用速度反馈来校正系统，校正后系统具有临界 作者: Badcoffee Email: blog.oliver@https://www.wendangku.net/doc/8c8156275.html, 2004年10月 原文出处: https://www.wendangku.net/doc/8c8156275.html,/yayong 这是作者在学习X86汇编过程中的学习笔记，难免有错误和疏漏之处，欢迎指正。 1. 编译环境 OS: Axianux 1.0 Compiler: gcc 3..2.3 Linker: Solaris Link Editors 5.x Debug Tool: gdb Editor: vi 2. 最简C代码分析 为简化问题，来分析一下最简的c代码生成的汇编代码： # vi test1.c int main() { return 0; } 编译该程序，产生二进制文件： # gcc -o start start.c # file start start: ELF 32-bit LSB executable, Intel 80386, version 1 (SYSV), for GNU/Linux 2.2.5, dynamically linked (uses shared libs), not stripped start是一个ELF格式32位小端(Little Endian)的可执行文件，动态链接并且符号表没有去除。这正是Unix/Linux平台典型的可执行文件格式。 用gdb反汇编可以观察生成的汇编代码： [wqf@15h166 attack]$ gdb start GNU gdb Asianux (6.0post-0.20040223.17.1AX) Copyright 2004 Free Software Foundation, Inc. GDB is free software, covered by the GNU General Public License, and you are welcome to change it and/or distribute copies of it under certain conditions. Type "show copying" to see the conditions. There is absolutely no warranty for GDB. Type "show warranty" for details. This GDB was configured as "i386-asianux-linux-gnu"...(no debugging symbols found)https://www.wendangku.net/doc/8c8156275.html,ing host libthread_db library "/lib/tls/libthread_db.so.1". (gdb) disassemble main --->反汇编main函数 Dump of assembler code for function main: 0x08048310 一、Matlab7.0函数及运算命令 1、运算符： ＋：加，－：减，*：乘，/：除，\：左除^：幂，‘：复数的共轭转置，（）：制定运算顺序。 2、常用函数表： sin( ) 正弦（变量为弧度）Cot( ) 余切（变量为弧度）sind( ) 正弦（变量为度数）Cotd( ) 余切（变量为度数）asin( ) 反正弦（返回弧度）acot( ) 反余切（返回弧度）Asind( ) 反正弦（返回度数）acotd( ) 反余切（返回度数）cos( ) 余弦（变量为弧度）exp( ) 指数 cosd( ) 余弦（变量为度数）log( ) 对数 acos( ) 余正弦（返回弧度）log10( ) 以10为底对数acosd( ) 余正弦（返回度数）sqrt( ) 开方 tan( ) 正切（变量为弧度）realsqrt( ) 返回非负根 tand( ) 正切（变量为度数）abs( ) 取绝对值 atan( ) 反正切（返回弧度）angle( ) 返回复数的相位角atand( ) 反正切（返回度数）mod(x,y) 返回x/y的余数sum( ) 向量元素求和 3、其余函数可以用help elfun和help specfun命令获得。 4、常用常数的值： pi 3.1415926……. realmin 最小浮点数，2^-1022 i 虚数单位 realmax 最大浮点数，（2－eps）2^1022 j 虚数单位 Inf 无限值 eps 浮点相对经度＝2^-52 NaN 空值 二、常用对象操作：除了一般windows窗口的常用功能键外。 1、!dir 可以查看当前工作目录的文件。!dir& 可以在dos状态下查看。 2、who 可以查看当前工作空间变量名，whos 可以查看变量名细节。 3、功能键： 功能键快捷键说明 方向上键Ctrl+P 返回前一行输入 方向下键Ctrl+N 返回下一行输入 方向左键Ctrl+B 光标向后移一个字符 方向右键Ctrl+F 光标向前移一个字符 Ctrl+方向右键Ctrl+R 光标向右移一个字符 Ctrl+方向左键Ctrl+L 光标向左移一个字符 home Ctrl+A 光标移到行首 End Ctrl+E 光标移到行尾 栈：在函数调用时，第一个进栈的是主函数中函数调用后的下一条指令（函数调用语句的下一条可执行语句）的地址，然后是函数的各个参数，在大多数的C 编译器中，参数是由右往左入栈的，然后是函数中的局部变量。注意静态变量是不入栈的。 当本次函数调用结束后，局部变量先出栈，然后是参数，最后栈顶指针指向最开始存的地址，也就是主函数中的下一条指令，程序由该点继续运行。 当发生函数调用的时候,栈空间中存放的数据是这样的: 1、调用者函数把被调函数所需要的参数按照与被调函数的形参顺序相反的顺序压入栈中,即:从右向左依次把被调函数所需要的参数压入栈; 2、调用者函数使用call指令调用被调函数,并把call指令的下一条指令的地址当成返回地址压入栈中(这个压栈操作隐含在call指令中); 3、在被调函数中,被调函数会先保存调用者函数的栈底地址(push ebp),然后再保存调用者函数的栈顶地址,即:当前被调函数的栈底地址(mov ebp,esp); 4、在被调函数中,从ebp的位置处开始存放被调函数中的局部变量和临时变量,并且这些变量的地址按照定义时的顺序依次减小,即:这些变量的地址是按照栈的延伸方向排列的,先定义的变量先入栈,后定义的变量后入栈; 所以,发生函数调用时,入栈的顺序为: 参数N 参数N-1 参数N-2 ..... 参数3 参数2 参数1 函数返回地址 上一层调用函数的EBP/BP 局部变量1 局部变量2 .... 局部变量N 函数调用栈如下图所示: 解释: 首先,将调用者函数的EBP入栈(push ebp),然后将调用者函数的栈顶指针ESP 赋值给被调函数的EBP(作为被调函数的栈底,mov ebp,esp),此时,EBP寄存器处于一个非常重要的位置,该寄存器中存放着一个地址(原EBP入栈后的栈顶),以该地址为基准,向上(栈底方向)能获取返回地址、参数值,向下(栈顶方向)能获取函数的局部变量值,而该地址处又存放着上一层函数调用时的EBP值; 一般而言,SS: [ebp+4]处为被调函数的返回地址,SS:[EBP+8]处为传递给被调函数的第一个参数(最后一个入栈的参数,此处假设其占用4字节内存)的 值,SS:[EBP-4]处为被调函数中的第一个局部变量,SS:[EBP]处为上一层EBP值;分析函数调用关系图(call graph)的几种方法
试求图示有源网络的传递函数和Bode图
c语言函数调用详细过程
matlab基本函数调用
函数调用栈