北航飞行力学大作业.
信号与系统实验5
信号与系统实验（五） 班级：通信5班姓名：刘贺洋学号：11081515 班级：通信5班姓名：章仕波学号：11081522 1.符号函数的傅里叶变换 (1)下面参考程序和运行结果是信号||2 f- t =的傅里叶变换，分析程序，判 e ) (t 断运行结果正确与否。 syms t; %时间符号 f=exp(-2*abs(t)); %符号函数 F=fourier(exp(-2*abs(t))); subplot(1,2,1); ezplot(f); subplot(1,2,2); ezplot(F); 1（1）图
(2)参考上述程序试画出信号)(32 )(3t u e t f t -=的波形及其幅频特性曲线。 1（2）源程序： syms t ; %时间符号 f=2/3*exp(-3*t)*heaviside(t); %符号函数 F=fourier(f); subplot(1,2,1); ezplot(f); subplot(1,2,2); ezplot(abs(F)); 1（2）图： 2.符号函数的傅里叶变换
(1)下面参考程序是求信号211)(ωω+=j F 的逆傅里叶变换，分析程 序，比较运行结果。 源程序2（1） syms t w; F=1/(1+w^2); f=ifourier(F,t); subplot(1,2,1); ezplot(F); subplot(1,2,2); ezplot(f); 2（1）图： （2）求信号ωωωsin 2)(=j F 的逆傅里叶变换，并用程序验证。
源程序2（2） syms t w; F=2*sin(w)/w; f=ifourier(F,t); subplot(1,2,1); ezplot(F); subplot(1,2,2); ezplot(f); 图2（2）：
北航最优化方法大作业参考
北航最优化方法大作业参考
1 流量工程问题 1.1 问题重述 定义一个有向网络G=(N,E)，其中N是节点集，E是弧集。令A是网络G的点弧关联矩阵，即N×E阶矩阵，且第l列与弧里(I,j)对应，仅第i行元素为1，第j行元素为-1，其余元素为0。再令b m=(b m1,…,b mN)T，f m=(f m1,…,f mE)T，则可将等式约束表示成： Af m=b m 本算例为一经典TE算例。算例网络有7个节点和13条弧，每条弧的容量是5个单位。此外有四个需求量均为4个单位的源一目的对，具体的源节点、目的节点信息如图所示。这里为了简单，省区了未用到的弧。此外，弧上的数字表示弧的编号。此时，c=((5,5…,5)1 )T， ×13 根据上述四个约束条件，分别求得四个情况下的最优决策变量x=((x12,x13,…,x75)1× )。 13 图 1 网络拓扑和流量需求
1.2 7节点算例求解 1.2.1 算例1（b1=[4;-4;0;0;0;0;0]T） 转化为线性规划问题： Minimize c T x1 Subject to Ax1=b1 x1>=0 利用Matlab编写对偶单纯形法程序，可求得: 最优解为x1*=[4 0 0 0 0 0 0 0 0 0 0 0 0]T 对应的最优值c T x1=20 1.2.2 算例2（b2=[4;0;-4;0;0;0;0]T） Minimize c T x2 Subject to Ax2=b2 X2>=0 利用Matlab编写对偶单纯形法程序，可求得: 最优解为x2*=[0 4 0 0 0 0 0 0 0 0 0 0 0]T 对应的最优值c T x2=20 1.2.3 算例3（b3=[0;-4;4;0;0;0;0]T） Minimize c T x3 Subject to Ax3=b3 X3>=0 利用Matlab编写对偶单纯形法程序，可求得: 最优解为x3*=[4 0 0 0 4 0 0 0 0 0 0 0 0]T 对应的最优值c T x3=40
北航电子电路设计数字部分实验报告
电子电路设计数字部分实验报告 学院： 姓名：
实验一简单组合逻辑设计 实验内容 描述一个可综合的数据比较器，比较数据a 、b的大小，若相同，则给出结果1，否则给出结果0。 实验仿真结果 实验代码 主程序 module compare(equal,a,b); input[7:0] a,b; output equal; assign equal=(a>b)1:0; endmodule 测试程序
module t; reg[7:0] a,b; reg clock,k; wire equal; initial begin a=0; b=0; clock=0; k=0; end always #50 clock = ~clock; always @ (posedge clock) begin a[0]={$random}%2; a[1]={$random}%2; a[2]={$random}%2; a[3]={$random}%2; a[4]={$random}%2; a[5]={$random}%2; a[6]={$random}%2; a[7]={$random}%2; b[0]={$random}%2; b[1]={$random}%2; b[2]={$random}%2; b[3]={$random}%2; b[4]={$random}%2;
b[5]={$random}%2; b[6]={$random}%2; b[7]={$random}%2; end initial begin #100000 $stop;end compare m(.equal(equal),.a(a),.b(b)); endmodule 实验二简单分频时序逻辑电路的设计 实验内容 用always块和@(posedge clk)或@(negedge clk)的结构表述一个1/2分频器的可综合模型，观察时序仿真结果。 实验仿真结果
2014信号与系统实验五(2)
Matlab拓展实验1：周期方波的傅立叶级数分析 一、周期方波信号的产生 1、周期方波函数：square(t,duty) （1）其中t为时间向量。 （2）duty为占空比，此参数可缺省，缺省时，duty=50。 注：周期方波的占空比为x%，是指一个周期内，高电平持续时间占整个周期时长的x%。但在square函数中，参数duty是一个0~100的数值。例如，若需产生占空比为20%的方波，则应设置duty=20。 2、功能描述：square函数产生一个周期为2π、高低电平分别为±1的周期方波，通过适当的编程可以调整为任意周期，任意幅度，任意中心值。 【例1】在时间范围[?4π,4π]产生周期的2π的周期方波。 clear;clf;clc; step=pi/200; %设置时间步长 t=-4*pi:step:4*pi; %设置时间范围 x1=square(t); %产生占空比为50%的周期方波 duty=20;x2=square(t,duty); %产生占空比为20%的周期方波 subplot(2,1,1);plot(t,x1); axis([-4*pi,4*pi,-1.5,1.5]);title('占空比为0.5的周期方波'); subplot(2,1,2);plot(t,x2); axis([-4*pi,4*pi,-1.5,1.5]);title('占空比为0.2的周期方波'); 上述程序的运行结果如下：
【思考】如将例1程序中的作图部分修改如下（修改部分以红色标出），这种处理方式称之为横坐标对pi归一化，图形结果会有何不同？这时横坐标的数值标注含义有何变化？ subplot(2,1,1);plot(t/pi,x1); axis([-4,4,-1.5,1.5]);title('占空比为0.5的周期方波'); subplot(2,1,2);plot(t/pi,x2); axis([-4,4,-1.5,1.5]);title('占空比为0.2的周期方波'); 【实践1】产生周期为1，占空比为30%，高电平为1，低电平为0的周期方波，时间范围为[?2,2]，即能观察到4个完整的周期。（时间变量的尺度变换） 二、傅立叶级数分解的计算 1、周期信号傅立叶级数展开的理论分析 周期信号的傅立叶级数有若干种形式（如三角形式、指数形式），这里以指数形式的傅立叶级数为例进行说明。 信号x(t)满足，x t=x(t+T)，其中T称为周期（或基波周期，单位s），F=1T称为频率（或基波频率，单位Hz），ω=2πT=2πF称为基波角频率，单位rad s。 （1）指数形式的傅立叶级数：
最优化方法大作业答案
1.用薄钢板制造一体积5m 3，长度不小于4m ，无上盖的货箱，要求钢板耗量最小。确定货箱的长x 1、宽x 2和高x 3。试列出问题的数学模型。 解：min 32312122x x x x x x z ++= s.t 5321=x x x 41≥x 0,,321≥x x x 2.将下面的线性规划问题表示为标准型并用单纯形法求解 max f=x 1+2x 2+x 3 s ．t ．2x 1+x 2-x 3≤2 -2x 1+x 2-5x 3≥-6 4x 1+x 2+x 3≤6 x i ≥0 i=1，2，3 解：先化标准形： Min 321x x x z -+= 224321=+-+x x x x 6525321=++-x x x x 646321=+++x x x x 列成表格：
1 2 1 610011460105122001112----- 可见此表已具备1°，2°，3°三个特点，可采用单纯形法。首先从底行中选元素-1，由2/2,6/2,6/4最小者决定选第一行第一列的元素2，标以记号，迭代一次得 1 2 1 2102310401162010021212 11-------- 再从底行中选元素-2/3，和第二列正元素1/2，迭代一次得 1 2 12 32 30 210231040116201002121211- ------ 再从底行中选元素-3，和第二列正元素2，迭代一次得 4 2 3 3 410120280114042001112--- 再迭代一次得 10 2 30 2 10 6 221023 1010213000421021013-- 选取最优解：
北航电子电路设计训练模拟分实验报告
北航电子电路设计训练模拟部分实验报告
————————————————————————————————作者：————————————————————————————————日期：
电子电路设计训练模拟部分实验 实验报告
实验一：共射放大器分析与设计 1.目的： （1）进一步了解Multisim的各项功能，熟练掌握其使用方法，为后续课程打好基础。 （2）通过使用Multisim来仿真电路，测试如图1所示的单管共射放大电路的静态工作点、电压放大倍数、输入电阻和输出电阻，并观察 静态工作点的变化对输出波形的影响。 （3）加深对放大电路工作原理的理解和参数变化对输出波形的影响。 （4）观察失真现象，了解其产生的原因。 图 1 实验一电路图 2.步骤： （1）请对该电路进行直流工作点分析,进而判断管子的工作状态。 （2）请利用软件提供的各种测量仪表测出该电路的输入电阻。 （3）请利用软件提供的各种测量仪表测出该电路的输出电阻。 （4）请利用软件提供的各种测量仪表测出该电路的幅频、相频特性曲线。 （5）请利用交流分析功能给出该电路的幅频、相频特性曲线。 （6）请分别在30Hz、1KHz、100KHz、4MHz和100MHz这5个频点利用示波器测出输入和输出的关系，并仔细观察放大倍数和相位差。 （提示：在上述实验步骤中，建议使用普通的2N2222A三极管，并请注 意信号源幅度和频率的选取，否则将得不到正确的结果。） 3.实验结果及分析： （1）根据直流工作点分析的结果，说明该电路的工作状态。 由simulate->analyses->DC operating point,可测得该电路的静态工作点为：
信号与系统-实验报告-实验五
实验五 连续信号与系统的S 域分析 学院 班级 姓名 学号 一、实验目的 1. 熟悉拉普拉斯变换的原理及性质 2. 熟悉常见信号的拉氏变换 3. 了解正/反拉氏变换的MATLAB 实现方法和利用MATLAB 绘制三维曲面图的方法 4. 了解信号的零极点分布对信号拉氏变换曲面图的影响及续信号的拉氏变换与傅氏变换的关系 二、 实验原理 拉普拉斯变换是分析连续时间信号的重要手段。对于当t ∞时信号的幅值不衰减的时间信号，即在f(t)不满足绝对可积的条件时，其傅里叶变换可能不存在，但此时可以用拉氏变换法来分析它们。连续时间信号f(t)的单边拉普拉斯变换F(s)的定义为： 拉氏反变换的定义为： 显然，上式中F(s)是复变量s 的复变函数，为了便于理解和分析F(s)随s 的变化规律，我们将F(s)写成模及相位的形式：()()()j s F s F s e ?=。其中，|F(s)|为复信号F(s)的模，而()s ?为F(s)的相位。由于复变量s=σ+jω，如果以σ为横坐标（实轴），jω为纵坐标（虚轴），这样，复变量s 就成为一个复平面，我们称之为s 平面。从三维几何空间的角度来看，|()|F s 和()s ?分别对应着复平面上的两个曲面，如果绘出它们的三维曲面图，就可以直观地分析连续信号的拉氏变换F(s)随复变量s 的变化情况，在MATLAB 语言中有专门对信号进行正反拉氏变换的函数，并且利用 MATLAB 的三维绘图功能很容易画出漂亮的三维曲面图。 ①在MATLAB 中实现拉氏变换的函数为： F=laplace( f ) 对f(t)进行拉氏变换，其结果为F(s) F=laplace (f,v) 对f(t)进行拉氏变换，其结果为F(v) F=laplace ( f,u,v) 对f(u)进行拉氏变换，其结果为F(v) ②拉氏反变换 f=ilaplace ( F ) 对F(s)进行拉氏反变换，其结果为f(t) f=ilaplace(F,u) 对F(w)进行拉氏反变换，其结果为f(u) f=ilaplace(F,v,u ) 对F(v)进行拉氏反变换，其结果为f(u) 注意： 在调用函数laplace( )及ilaplace( )之前，要用syms 命令对所有需要用到的变量（如t,u,v,w ）等进行说明，即要将这些变量说明成符号变量。对laplace( )中的f 及ilaplace( )中的F 也要用符号定义符sym 将其说明为符号表达式。具体方法参见第一部分第四章第三节。 例①：求出连续时间信号 ()sin()()f t t t ε=的拉氏变换式，并画出图形 求函数拉氏变换程序如下： syms t s %定义符号变量 ft=sym('sin(t)*Heaviside(t)'); %定义时间函数f(t)的表达式
最优化方法大作业
发动机空燃比控制器 引言：我主要从事自动化相关研究。这里介绍我曾经接触过的发动机空燃比控制器设计中的优化问题。 发动机空燃比控制器设计中的最优化问题 AFR =a f m m && (1) 空燃比由方程（1）定义，在发动机运行过程中如果控制AFR 稳定在14.7可以获 得最好的动力性能和排放性能。如果假设进入气缸的空气流量a m &可以由相关单元检测得到，则可以通过控制进入气缸的燃油流量f m &来实现空燃比的精确控制。由于实际发动机的燃油喷嘴并不是直接对气缸喷燃油，而是通过进气歧管喷燃油，这么做会在进 气歧管壁上液化形成油膜，因此不仅是喷嘴喷出的未液化部分燃油会进入气缸，油膜 蒸发部分燃油也会进入气缸，如方程（2）。这样如何更好的喷射燃油成为了一个问题。 1110101122211ττττ?? ?? -?? ??????????=+????????-????????????-???? ? ??? ?? ????????? ?f f f v X x x u x x X x y =x && (2) 其中12、,==ff fv x m x m &&=f y m &,=fi u m &这里面，表示油膜蒸发量ff m &、fv m &表示为液化部分燃油、fi m &表示喷嘴喷射的燃油，在τf 、τv 、X 都已知的情况下，由现代控制理论知识，根据系统的增广状态空间模型方程（3） 0000001 1 011011114.70ττττ????-?? ??????????=-+-??????????????? ??????????????? ?? ??=?????? f f v v a X X u +q q m y q x x x &&& (3) 其中()0 14.7?t a q = y -m &。由极点配置方法，只要设计控制器方程（4），就可以 使得y 无差的跟踪阶跃输入，那么y 也能较好的跟踪AFR *a m /&。 12-- u =K q K x (4) 这里面的12、K K 确定，可由主导极点概念降维成两个参数12C ,C ，虽然都是最终稳态无差，但是目标是使得瞬态过程中y 和阶跃输入y r 的差异尽可能的小。所以原问
信号与系统实验实验报告
信号与系统实验实验报 告 文档编制序号：[KKIDT-LLE0828-LLETD298-POI08]
实验五连续系统分析一、实验目的 深刻理解连续时间系统的系统函数在分析连续系统的时域特性、频域特性及稳定性中的重要作用及意义，掌握根据系统函数的零极点设计简单的滤波器的方法。掌握利用MATLAB分析连续系统的时域响应、频响特性和零极点的基本方法。 二、实验原理 MATLAB提供了许多可用于分析线性时不变连续系统的函数，主要包含有系统函数、系统时域响应、系统频域响应等分析函数。 三、实验内容 1.已知描述连续系统的微分方程为，输入，初始状态 ，计算该系统的响应，并与理论结果比较，列出系统响应分析的步骤。 实验代码： a=[1 10]; b=[2]; [A B C D]=tf2ss(b,a); sys=ss(A,B,C,D); t=0: :5; xt=t>0; sta=[1]; y=lsim(sys,xt,t,sta); subplot(3,1,1); plot(t,y); xlabel('t'); title('系统完全响应 y(t)'); subplot(3,1,2); plot(t,y,'-b'); hold on yt=4/5*exp(-10*t)+1/5; plot(t,yt,' : r'); legend('数值计算','理论计算'); hold off xlabel('t'); subplot(3, 1 ,3); k=y'-yt; plot(t,k); k(1) title('误差');
实验结果： 结果分析： 理论值 y(t)=0. 8*exp(-10t)*u(t)+ 程序运行出的结果与理论预期结果相差较大误差随时间增大而变小，初始值相差最大，而后两曲线基本吻合，表明该算法的系统响应在终值附近有很高的契合度，而在初值附近有较大的误差。 2.已知连续时间系统的系统函数为，求输入分别为，， 时，系统地输出，并与理论结果比较。 a=[1,3,2,0]; b=[4,1]; sys=tf(b,a); t=0: :5; x1=t>0; x2=(sin(t)).*(t>0); x3=(exp(-t)).*(t>0); y1=lsim(sys,x1,t); y2=lsim(sys,x2,t); y3=lsim(sys,x3,t); subplot(3,1,1); plot(t,y1); xlabel('t'); title('X(t)=u(t)'); subplot(3,1,2); plot(t,y2); xlabel('t'); title('X(t)=sint*u(t)'); subplot(3, 1 ,3); plot(t,y3); xlabel('t'); title('X(t)=exp(-t)u(t)'); 实验结果： 结果分析： a=[1,3,2,0]; b=[4,1]; sys=tf(b,a); t=0: :5; x1=t>0; x2=(sin(t)).*(t>0); x3=(exp(-t)).*(t>0); y1=lsim(sys,x1,t); y2=lsim(sys,x2,t); y3=lsim(sys,x3,t); subplot(3,1,1); plot(t,y1,'-b');
北航惯性导航大作业
惯性导航基础课程大作业报告（一）光纤陀螺误差建模与分析 班级：111514 姓名： 学号 2014年5月26日
一.系统误差原理图 二.系统误差的分析 （一）漂移引起的系统误差 1. εx ，εy ，εz 对东向速度误差δVx 的影响 clc;clear all; t=1:0.01:25; g=9.8; L=pi/180*39; Ws=2*pi/84.4*60; Wie=2*pi/24; R=g/(Ws)^2; e=0.1*180/pi; mcVx1=e*g*sin(L)/(Ws^2-Wie^2)*(sin(Wie*t)-Wie*sin(Ws*t)/Ws); mcVx2=e*((Ws^2-(Wie^2)*((cos(L))^2))/(Ws^2-Wie^2)*cos(Ws*t)-(Ws^2)*((sin(L))^2)*cos(Wi e*t)/(Ws^2-Wie^2)-(cos(L))^2); mcVx3=(sin(L))*(cos(L))*R*e*((Ws^2)*cos(Wie*t)/(Ws^2-Wie^2)-(Wie^2)*cos(Ws*t)/(Ws^2-Wi e^2)-1); plot(t,[mcVx1',mcVx2',mcVx3']); title('Ex,Ey,Ez 对Vx 的影响'); xlabel('时间t'); ylabel('Vx(t)'); 0,δλδL ,v v δδ
legend('Ex-mcVx1','Ey-mcVx2','Ez-mcVx3'); grid; axis square; 分析：εx，εy，εz对东向速度误差δVx均有地球自转周期的影响，εx，εy还会有舒勒周期分量的影响，其中，εy对δVx的影响较大。 2.εx，εy，εz对东向速度误差δVy的影响 clc;clear all; t=1:0.01:25; g=9.8; L=pi/180*39; Ws=2*pi/84.4*60; Wie=2*pi/24; R=g/(Ws)^2; e=0.1*180/pi; mcVy1=e*g*(cos(Wie*t)-cos(Ws*t))/(Ws^2-Wie^2); mcVy2=g*sin(L)*e/(Ws^2-Wie^2)*(sin(Wie*t)-Wie/Ws*sin(Ws*t)); mcVy3=g*cos(L)*e/(Ws^2-Wie^2)*(sin(Wie*t)-Wie/Ws*sin(Ws*t)); plot(t,[mcVy1',mcVy2',mcVy3']); title('Ex,Ey,Ez对Vy的影响'); xlabel('时间t'); ylabel('Vy(t)'); legend('Ex-mcVy1','Ey-mcVy2','Ez-mcVy3'); grid; axis square;
北航飞行器设计与应用力学系.doc
航空科学与工程学院 2016年研究生入学考试复试大纲 一、复试方式：笔试＋面试 二、复试组织： 1、笔试：由航空学院统一组织，考试科目及复试大纲另见《航空科学与工程学院2013年考研复试安排》。 2、口试：以学科专业组为单位，由3-5位硕士生导师组成面试小组（组长为教授），每位考生的面试时间为20分钟。 三、复试流程和评分标准： 1）检查并核实考生面试所必备的个人证件和材料；考生可以提供有助于证明自己背景和能力的相关材料，证件和材料完备是面试的必要条件。 2）考生用英语口述个人基本情况、兴趣等，面试小组老师就考生基本情况提问，考生用英文回答问题。 3）考生朗读一段考场指定的专业外语短文，并口头翻译成中文。 4）面试小组老师就基础理论知识提问，学生用中文回答问题。 5）面试小组老师就专业知识提问，学生用中文回答问题。 面试结束后考生退场，在3-5个工作日后见航空学院网站“招生就业”栏目的“研究生招生”，会通知出学院的拟录取名单，在7层的研究生教学橱窗也会公布。 四、考场纪律 考生准时到达指定的复试考场，遵守考场秩序，尊重考试教师。 五、各学科专业组具体复试内容及参考书： 1、飞行力学与飞行安全系2016年硕士研究生入学复试程序 方式： 由3~6位硕士生导师组成面试小组，每位考生的面试时间为20分钟。 范围： 面试范围包括英语口语能力、专业英语阅读理解能力、专业基础理论知识和专业知识。具体环节如下： 1）对考生学习背景、心理、爱好和志愿等基本情况的了解。 2）考察考生的英语阅读和口头表达能力。
3）基础理论和专业知识面试。基础理论包括自动控制原理、理论力学和材料力学。专业知识包括飞行力学、飞行安全、飞行器总体设计、空气动力学等。 参考书： 基础理论可以选用任何一本考生熟悉的《自动控制原理》、《理论力学》、《材料力学》教材。专业课可以参考《飞机飞行动力学》（熊海泉编）或《飞机飞行性能》、《飞机的稳定与控制》等方面的参考书。 面试流程和评分标准： 1）检查并核实考生面试所必备的个人证件和材料；证件和材料完备是面试的必要条件。2）考生用英语口述个人基本情况、兴趣等，面试小组老师就考生基本情况提问，考生回答问题。 3）读一段指定的专业外语，并口头翻译成中文。 4）面试小组老师就基础理论知识提问，学生回答问题。 5）面试小组老师就专业知识提问，学生回答问题。 6）问答结束后，考生退场，面试老师根据考核要求和面试情况，对考生进行评分。 7）所有考生面试结束后，面试老师根据总体情况，对所有考生进行综合评估和比较，给出面试成绩。 2、人机与环境工程/制冷及低温工程2016年硕士研究生入学复试程序 方式： 由3~5位硕士生导师组成面试小组，每位考生的面试时间为20分钟。 范围： 1）英语阅读和口头表达能力。 2）对考生心理、基本情况的了解。 3）基础理论和专业知识面试。基础理论包括：自动控制原理，理论力学，流体力学；专业知识包括工程热力学，传热学，人机工程，低温制冷。考生可以选择其中1门基础理论和1门专业课作为面试内容，或者是综合知识。 参考书： 可以选用任何一本考生熟悉的《自动控制原理》、《理论力学》、《流体力学》教材。专业课可以选用考生熟悉的《工程热力学》，《传热学》，《人机工程》，低温制冷等方面的参考书。 面试流程和评分标准： 1）检查并核实考生面试所必备的个人证件和材料；证件和材料完备是面试的必要条件. 2）考生用英语口述个人基本情况、兴趣等，面试小组老师就考生基本情况提问，考生回答问题。 3）读一段指定的专业外语，并口头翻译成中文。 4）面试小组老师就基础理论知识提问，学生回答问题。 5）面试小组老师就专业知识提问，学生回答问题。 6) 问答结束后，考生退场，面试老师根据考核要求和面试情况，对考生进行评分。
北航自动控制原理实验报告(完整版)
自动控制原理实验报告 一、实验名称：一、二阶系统的电子模拟及时域响应的动态测试 二、实验目的 1、了解一、二阶系统阶跃响应及其性能指标与系统参数之间的关系 2、学习在电子模拟机上建立典型环节系统模型的方法 3、学习阶跃响应的测试方法 三、实验内容 1、建立一阶系统的电子模型，观测并记录在不同时间常数T时的响应曲线，测定过渡过程时间T s 2、建立二阶系统电子模型，观测并记录不同阻尼比的响应曲线，并测定超调量及过渡过程时间T s 四、实验原理及实验数据 一阶系统 系统传递函数： 由电路图可得，取则K=1，T分别取：0.25, 0.5, 1 T 0.25 0.50 1.00 R2 0.25MΩ0.5M Ω1MΩ C 1μ1μ1μ T S 实测0.7930 1.5160 3.1050 T S 理论0.7473 1.4962 2.9927 阶跃响应曲线图1.1 图1.2 图1.3 误差计算与分析 （1）当T=0.25时，误差==6.12%； （2）当T=0.5时，误差==1.32%； （3）当T=1时，误差==3.58% 误差分析：由于T决定响应参数，而,在实验中R、C的取值上可能存在一定误差，另外，导线的连接上也存在一些误差以及干扰，使实验结果与理论值之间存在一定误差。但是本实验误差在较小范围内，响应曲线也反映了预期要求，所以本实验基本得到了预期结果。 实验结果说明 由本实验结果可看出，一阶系统阶跃响应是单调上升的指数曲线，特征有T确定，T越小，过度过程进行得越快，系统的快速性越好。 二阶系统 图1.1 图1.2 图1.3
系统传递函数： 令 二阶系统模拟线路 0.25 0.50 1.00 R4 210.5 C2 111 实测45.8% 16.9% 0.6% 理论44.5% 16.3% 0% T S实测13.9860 5.4895 4.8480 T S理论14.0065 5.3066 4.8243 阶跃响应曲线图2.1 图2.2 图2.3 注：T s理论根据matlab命令[os,ts,tr]=stepspecs(time,output,output(end),5)得出，否则误差较大。 误差计算及分析 1）当ξ=0.25时，超调量的相对误差= 调节时间的相对误差= 2）当ξ=0.5时，超调量的相对误差==3.7% 调节时间的相对误差==3.4% 4）当ξ=1时，超调量的绝对误差= 调节时间的相对误差==3.46% 误差分析：由于本试验中，用的参量比较多，有R1,R2,R3,R4;C1,C2;在它们的取值的实际调节中不免出现一些误差，误差再累加，导致最终结果出现了比较大的误差，另外，此实验用的导线要多一点，干扰和导线的传到误差也给实验结果造成了一定误差。但是在观察响应曲线方面，这些误差并不影响，这些曲线仍旧体现了它们本身应具有的特点，通过比较它们完全能够了解阶跃响应及其性能指标与系统参数之间的关系，不影响预期的效果。 实验结果说明 由本实验可以看出，当ωn一定时，超调量随着ξ的增加而减小，直到ξ达到某个值时没有了超调；而调节时间随ξ的增大，先减小，直到ξ达到某个值后又增大了。 经理论计算可知，当ξ=0.707时，调节时间最短，而此时的超调量也小于5%，此时的ξ为最佳阻尼比。此实验的ξ分布在0.707两侧，体现了超调量和调节时间随ξ的变化而变化的过程，达到了预期的效果。 图2.2 图2.1 图2.3
15秋北航《信号与系统》在线作业二100分答案
北航《信号与系统》在线作业二 一、单选题（共10 道试题，共30 分。） 1. 信号〔ε(t)－ε(t－2)〕的拉氏变换的收敛域为________。 A. Re[s]>0 B. Re[s]>2 C. 全S平面 D. 不存在 -----------------选择：C 2. 信号的时宽与信号的频宽之间呈________。 A. 正比关系 B. 反比关系 C. 平方关系 D. 没有关系 -----------------选择：B 3. If f1(t) ←→F1(jω)，f2(t) ←→F2(jω) Then________。 A. [a f1(t) + b f2(t) ] ←→ [a F1(jω) *b F2(jω) ] B. [a f1(t) + b f2(t) ] ←→ [a F1(jω) - b F2(jω) ] C. [a f1(t) + b f2(t) ] ←→ [a F1(jω) + b F2(jω) ] D. [a f1(t) + b f2(t) ] ←→ [a F1(jω) /b F2(jω) ] -----------------选择：D 4. 某信号的频谱是周期的离散谱，则对应的时域信号为________。 A. 连续的周期信号 B. 连续的非周期信号 C. 离散的非周期信号 D. 离散的周期信号 -----------------选择：D 5. 信号在时域拥有的总能量，________其频谱在频域内能量的总和。 A. 大于 B. 等于 C. 小于 D. 不等于 -----------------选择：B 6. 理想低通滤波器是________。 A. 因果系统 B. 物理可实现系统 C. 非因果系统 D. 响应不超前于激励发生的系统 -----------------选择：C 7. 连续周期信号的傅氏变换是________。 A. 连续的 B. 周期性的 C. 离散的 D. 与单周期的相同 -----------------选择：C
FPGA实验报告北航电气技术实验
FPGA电气技术实践 实验报告 院（系）名称宇航学院 专业名称飞行器设计与工程（航天）学生学号XXXXXXXX 学生姓名XXXXXX 指导教师XXXX 2017年11月XX日
实验一四位二进制加法计数器与一位半加器的设计实验时间：2017.11.08（周三）晚实验编号20 一、实验目的 1、熟悉QuartusII的VHDL的文本编程及图形编程流程全过程。 2、掌握简单逻辑电路的设计方法与功能仿真技巧。 3、学习并掌握VHDL语言、语法规则。 4、参照指导书实例实现四位二进制加法计数器及一位半加器的设计。 二、实验原理 .略 三、实验设备 1可编程逻辑实验箱EP3C55F484C8 一台（包含若干LED指示灯，拨码开关等）2计算机及开发软件QuartusII 一台套 四、调试步骤 1四位二进制加法计数器 (1)参照指导书实例1进行工程建立与命名。 (2)VHDL源文件编辑 由于实验箱上LED指示灯的显示性质为“高电平灭，低电平亮”，为实现预期显示效果应将原参考程序改写为减法器，且”q1<= q1+1”对应改为”q1<= q1-1”，以实现每输入一个脉冲“亮为1，灭为0”。 由于参考程序中的rst清零输入作用并未实现，所以应将程序主体部分的最外部嵌套关于rst输入是否为1的判断，且当rst为1时，给四位指示灯置数”1111”实现全灭，当rst为0时，运行原计数部分。 (3)参照指导书进行波形仿真与管脚绑定等操作，链接实验箱并生成下载文件 (4)将文件下载至实验箱运行，观察计数器工作现象，调试拨动开关查看是否清零。 可以通过改变与PIN_P20（工程中绑定为clk输入的I/O接口）相连导线的另一端所选择的实验箱频率时钟的输出口位置，改变LED灯显示变化频率。 并且对照指导书上对实验箱自带时钟频率的介绍，可以通过改变导线接口转换输入快慢，排查由于clk输入管脚损坏而可能引起的故障。
北理工信号与系统实验(5)
实验5 连续时间系统的复频域分析 一、实验目的 1．掌握拉普拉斯变换及其反变换的定义，并掌握MATLAB实现方法。 2．学习和掌握连续时间系统函数的定义及复频域分析方法。 3．掌握系统零极点的定义，加深理解系统零极点分布与系统特性的关系。 二、实验原理与方法 1.拉普拉斯变换 连续时间信号x(t)的拉普拉斯变换、反变换定义为 MATLAB中，可以采用符号数学工具箱的laplace函数和ilaplace函数进行拉氏变换和拉氏反变换。 L=laplace(F)符号表达式F的拉氏变换，F中时间变量为t，返回变量为s的结果表达式。 L=laplace(F,t)用t替换结果中的变量t。 F=ilaplace(L)以s为变量的符号表达式L的拉氏反变换，返回时间变量为t的结果表达式。 F=ilaplace(L,x)用x替换结果中的变量t。 2.连续时间系统的系统函数 连续时间系统的系统函数是系统单位冲击响应的拉氏变换 连续时间系统的系统函数还可以由输入和输出信号的拉氏变换之比得到 单位冲击响应h(t)反映了系统的固有性质，而H(s)从复频域反映了系统的固有性质。由上式描述的连续时间系统，其系统系数为s的有理函数
3.连续时间系统的零极点分析 系统的零点指使上式的分子多项式为零的点，极点指使分母多项式为零的点，零点使系统的值为零，极点使系统的值为无穷大。通常将系统函数的零极点绘在s平面上，零点用○表示，极点用×表示，得到零极点分布图。 由零点定义可知，零点和极点分别指上式分子和分母多项式的根，利用MATLAB求多项式的根可以通过roots来实现，该函数调用格式为： r=roots(c) c为多项式的系数向量，返回值r为多项式的根向量。 分别对上式的分子、分母多项式求根即可求得零极点。 此外，在MATLAB中还提供了pzmap函数来求取零极点和绘制系统函数的零极点分布图，该函数的调用格式为： Pzmap(sys) 绘出系统模型sys描述的系统的零极点分布图。 [p,z]=pzmap(sys) 这种调用方法返回极点和零点，而不绘出零极点分布图。其中sys为系统传函模型，由t命令sys=tf(b,a)实现，b、a为传递函数的分子和分母多项式的系数向量。 MATLAB还提供了两个专用函数tf2zp和zp2tf来实现系统传递函数模型和灵机电网增益模型的转换，其调用格式为 [z,p,k]=th2zp(b,a) [b,a]=th2zp(z,p,k) 其中b、a为传递函数的分子和分母多项式的系数向量，返回值z为零点列向量，p为极点列向量，k为系统函数零极点形式的增益。 三、实验内容 1.已知系统的冲激响应h(t)u(t)u(t2) =，试采用复频域的 =--，输入信号x(t)u(t) 方法求解系统的响应，编写MATLAB程序实现。 解：MATLAB程序如下： h=sym('heaviside(t)-heaviside(t-2)'); H=laplace(f) x=sym('heaviside(t)'); X=laplace(x) Y=H*X
北航数值分析大作业第二题精解
目标：使用带双步位移的QR 分解法求矩阵10*10[]ij A a =的全部特征值，并对其中的每一个实特征值求相应的特征向量。已知：sin(0.50.2)() 1.5cos( 1.2)(){i j i j ij i j i j a +≠+== (i,j=1,2, (10) 算法： 以上是程序运作的逻辑，其中具体的函数的算法，大部分都是数值分析课本上的逻辑，在这里特别写出矩阵A 的实特征值对应的一个特征向量的求法： ()[]()() []()[]()111111I 00000 i n n n B A I gause i n Q A I u Bu u λλ-?-?-=-?-?? ?-=????→=??????→= ?? ? 选主元的消元 检查知无重特征值 由于=0i A I λ- ，因此在经过选主元的高斯消元以后，i A I λ- 即B 的最后一行必然为零，左上方变 为n-1阶单位矩阵[]()()11I n n -?-，右上方变为n-1阶向量[]()11n Q ?-，然后令n u 1=-，则 ()1,2,,1j j u Q j n ==???-。
这样即求出所有A所有实特征值对应的一个特征向量。 #include 航空科学与工程学院 《飞行力学实验班》课程实验飞机典型模态特性仿真 实验报告 学生姓名:姜南 学号:11051136 专业方向:飞行器设计与工程 指导教师:王维军 (2014年 6 月29日 一、实验目的 飞机运动模态是比较抽象的概念, 是课程教学中的重点和难点。本实验针对这一问题,采用计算机动态仿真和在人-机飞行仿真实验平台上的驾驶员在环仿真实验,让学生身临其境地体会飞机响应与模态特性的关系,加深对飞机运动模态特性的理解。 二、实验内容 1.纵向摸态特性实验 计算某机在某状态下的短周期运动、长周期运动的模态参数;进行时域的非实时或实时仿真实验,操纵升降舵激发长、短周期运动模态,并由结果曲线分析比较模态参数;放宽飞机静稳定性,观察典型操纵响应曲线,并通过驾驶员在环实时仿真体验飞机的模态特性变化。 2.横航向模态特性实验 计算某机在某状态下的滚转、荷兰滚、螺旋模态参数;进行时域仿真计算,操纵副翼或方向舵,激发滚转、荷兰滚等运动模态,并由结果曲线分析比较模态参数。 三、各典型模态理论计算方法及模态参数结果表 1 纵向模态纵向小扰动运动方程 0000 1 00 0e p e p e p u w e u w q p u w q X X u u X X g Z Z w w Z Z Z q q M M M M M δδδδδ δδδθθ????????-???? ????????? ? ???????????=+??????????????????? ?????????????????? A =[ X u X ?w Z u Z w 0?g Z q 0M ?u M ?w0 M q 010] =[?0.01999980.0159027?0.0426897?0.04034850?32.2869.6279 0?0.00005547?0.001893500?0.54005010] A 的特征值方程 |λ+0.0199998?0.01590270.0426897 λ+0.0403485032.2 ?869.627900.000055470.001893500λ+0.540050 ?1λ |=0 特征根λ1,2=?0.290657205979137±1.25842158268078i λ3,4=?0.00954194402086311±0.0377636398212079i 半衰期t 1/2由公式t 1/2=? ln2λ 求得,分别为 t 1/2,1=2.38475828674173s t 1/2,3=72.6421344585972s 振荡频率ω分别为 ω1=1.25842158268078rad/s ω3=0.0377636398212079rad/s 周期T 由公式T =北航飞行力学实验班飞机典型模态特性仿真实验报告(精)