MATLAB simulink 通信系统建模与仿真 第四章部分习题答案
第四章、
2、仿真一个三角波信号通过AWGN信道后的结果，分别用randn函数和awgn函数实现。
3、修改例题4.7中的Tx模块采用BPSK基带调制，重新仿真，观察仿真结果与4.7有何不同？
4、产生最大多普勒频移为120的多径瑞丽衰落信道，假设信号的抽样时间间隔为1/100000s，多径延迟为[0 6e-5 11e-5],各经增益为[0 -3 -6],所有路径的接收信号强度之和为0，画出信道的功率随时间变化曲线。
实验一 基于Matlab的控制系统模型
实验一 基于Matlab 的控制系统模型 姓名 学号 班级 一、实验目的 1) 熟悉Matlab 的使用环境，学习Matlab 软件的使用方法和编程方法。 2) 学习使用Matlab 进行各类数学变换运算的方法。 3) 学习使用Matlab 建立控制系统模型的方法。 二、实验原理 1. 香农采样定理 对一个具有有限频谱的连续信号f (t )进行连续采样，当采样频率满足ωs ≥ωmax 时，采样信号f *(t )能无失真的复现原连续信号。 (1) 作信号f (t )=5e 10t 和f *(t ) =5e 10kT 的曲线，比较采样前后的差异。 0.05 0::0.5 5*(10*) subplot(2,1,1) plot(,) grid subplot(2,1,2) stem(,) grid T t T f exp t t f t f ===- 请改变采样周期T ，观察不同的采样周期下的采样效果。
(2) 频谱曲线 50:1:50 5./(100.^2) (,)w F sqrt w plot w F grid =-=+ 若|F (j ωmax ) |=0.1|F (0)|，选择合理的采样周期T 并验加以证。 400:20:400 200 2*/05/*(1./(100.^2)) 15/*(1./(100().^2)) 25/*(1./(100().^2)) (,0,,1,,2) w ws Ts pi ws F Ts sqrt w F Ts sqrt w ws F Ts sqrt w ws plot w F w F w F grid =-===+=+-=++ 请改变采样频率，观察何时出现频谱混叠？ 2. 拉式变换和Z 变换 (1) 使用Matlab 求函数的拉氏变换 拉式变换： 反拉氏变换： ()()()()()()2 222 1exp -*123*exp -*4sin *5exp -*s 11/(1) 21/()31/4/() 51/(*(2)*(*c 3)o ) s *yms syms a w t f a t laplace f f t f t a t f s a f s ilaplace f f s a f s f w s w f s s s w t f a t w t ==+==+====++== (2) 使用Matlab 求函数的Z 变换 Z 变换： 反Z 变换：
实验七 SIMULINK仿真集成环境
实验七 SIMULINK 仿真集成环境 一、实验目的 熟悉SIMULINK 的模型窗口、熟练掌握SIMULINK 模型的创建，熟练掌握常用模块的操作及其连接。 二、实验内容 (1) SIMULINK 模型的创建和运行。 (2)一阶系统仿真 三、实验步骤 1. Simulink 模型的创建和运行 (1) 创建模型。 ① 在MATLAB 的命令窗口中输入simulink 语句，或者单击MATLAB 工具条上的SIMULINK 图标，SIMULINK 模块库浏览器。 ②在MA TLAB 菜单或库浏览器餐单中选择File|New|Model ，或者单击库浏览器的图标，即可新建一个“untitle ”的空白模型窗口。 ③打开“Sources ”模块库，选择“Sine Wave ”模块，将其拖到模型窗口，再重复一次；打开“Math Operatioins ”模块库选取“Product ”模块；打开“Sinks ”模块库选取“Scope ”模块。 (2) 设置模块参数 ① 修改模块注释。单击模块的注释处，出现虚线的编辑框，在编辑框中修改注释。 ② 双击下边“Sine Wave ”模块，弹出参数对话框，浆“Frequency ”设置为100；双击“Scope ”模块，弹出示波器窗口，然后单击示波器图标，弹出参数对话框，修改示波器的通道数“Number of axes ”为3. ③如图所示，用信号线连接模块。 (3) 启动仿真 ① 单击工具栏上的图标或者选择Simulation|Start 菜单项，启动仿真；然后双击“Scope ”模块弹出示波器窗口，可以看到波形图。 ② 修改仿真步长。在模块窗口的Simulation 菜单下选择“Configuration Parameters ”命令，把“Max step size ”设置为0.01；启动仿真，观察波形是不是比原来光滑。 ③再次修改“Max step size ”为0.001；设置仿真终止时间为10s ；启动仿真，单击示波器工具栏中的按钮，可以自动调整显示范围，可以看到波形的起点不是零点，这是因为步长改小后，数据量增大，超出了示波器的缓冲。 浆示波器的参数对话框打开，选择“Data history ”页，把“Limit data point tolast ”设置为10000；再次启动仿真，观察示波器将看到完整的波形。 2.. 一阶系统仿真 使用阶跃信号作为输入信号，经过传递函数为1 6.01 s 的一阶系统，观察其输出。 ①设置“Step ”模块的“Step time ”为0；浆仿真参数的最大步长“Max step size ”设置为0.01. 把结果数据输出到工作空间。 ②打开“Sources ”模块库，选取“Clock ”模块添加到模型窗口中。 ③代开“Sinks ”模块库，选取两个“To workspace ”模块添加到模型窗口中，两个模块分别连接输出和“Clock ”模块。
MATLAB Simulink系统建模与仿真 实验报告
MATLAB/Simulink 电力系统建模与仿真 实验报告 姓名：****** 专业：电气工程及其自动化 班级：******************* 学号：*******************
实验一无穷大功率电源供电系统三相短路仿真 1.1 无穷大功率电源供电系统仿真模型构建 运行MATLAB软件，点击Simulink模型构建，根据电路原理图，添加下列模块： （1）无穷大功率电源模块（Three-phase source） （2）三相并联RLC负荷模块（Three-Phase Parallel RLC Load） （3）三相串联RLC支路模块（Three-Phase Series RLC Branch） （4）三相双绕组变压器模块（Three-Phase Transformer (Two Windings)） （5）三相电压电流测量模块（Three-Phase V-I Measurement） （6）三相故障设置模块（Three-Phase Fault） （7）示波器模块（Scope） （8）电力系统图形用户界面（Powergui） 按电路原理图连接线路得到仿真图如下： 1.2 无穷大功率电源供电系统仿真参数设置 1.2.1 电源模块 设置三相电压110kV，相角0°，频率50Hz，接线方式为中性点接地的Y形接法，电源电阻0.00529Ω，电源电感0.000140H，参数设置如下图：
1.2.2 变压器模块 变压器模块参数采用标幺值设置，功率20MVA，频率50Hz，一次测采用Y型连接，一次测电压110kV，二次侧采用Y型连接，二次侧电压11kV，经过标幺值折算后的绕组电阻为0.0033，绕组漏感为0.052，励磁电阻为909.09，励磁电感为106.3，参数设置如下图： 1.2.3 输电线路模块 根据给定参数计算输电线路参数为：电阻8.5Ω，电感0.064L，参数设置如下图： 1.2.4 三相电压电流测量模块 此模块将在变压器低压侧测量得到的电压、电流信号转变成Simulink信号，相当于电压、电流互感器的作用，勾选“使用标签（Use a label）”以便于示波器观察波形，设置电压标签“Vabc”，电流标签“Iabc”，参数设置如下图：
控制系统的MATLAB仿真与设计课后答案
控制系统的MATLAB仿真与设计课后答案
>>z=-4*sqrt(2)*sin(t); >>plot3(x,y,z,'p'); >>title('Line in 3-D Space'); >>text(0,0,0,'origin'); >>xlabel('X'),ylable('Y'),zlable('Z');grid; 4>>theta=0:0.01:2*pi; >>rho=sin(2*theta).*cos(2*theta); >>polar(theta,rho,'k'); 5>>[x,y,z]=sphere(20); >>z1=z; >>z1(:,1:4)=NaN; >>c1=ones(size(z1)); >>surf(3*x,3*y,3*z1,c1); >>hold on >>z2=z; >>c2=2*ones(size(z2)); >>c2(:,1:4)=3*ones(size(c2(:,1:4))); >>surf(1.5*x,1.5*y,1.5*z2,c2); >>colormap([0,1,0;0.5,0,0;1,0,0]); >>grid on >>hold off 第四章 1>>for m=100:999 m1=fix(m/100); m2=rem(fix(m/10),10); m3=rem(m,10); if m==m1*m1*m1+m2*m2*m2+m3*m3*m3 disp(m) end end 2M文件：function[s,p]=fcircle(r) s=pi*r*r; p=2*pi*r; 主程序： [s,p]=fcircle(10) 3>>y=0;n=100; for i=1:n y=y+1/i/i; end >>y
基于MATLAB的变压器仿真 与分析
于MATLAB_Simulink的牵引变压器建模与仿真 基于MATLAB/Simulink的牵引变压器建模与仿真徐（西安铁路局安康供电段新陕西汉中 723000）摘要：针对多种牵引变压器接线方式，建立数学模型，基于Matlab/Simulink仿真软件，建立牵引变压器的仿真模型，并验证数学模型和仿真模型的一致性。利用所建立仿真模型对不同接线形式牵引变压器在不同条件下对公用电网产生的谐波和负序影响进行仿真试验，对研究各种类型的牵引变压器特性在我国电气化铁路的应用提供条件。关键词：牵引变压器；数学模型；仿真模型；Matlab/Simulink 中图分类号：U223.6 文献标识码：A 文章编号：1671－7597（2011）0610061－03 牵引变压器按其特性可分为平衡接线和不平衡接线。其中不平衡接线有单相接线、Vv接线和YNd11接线；平衡接线是试图实现三相两相对称变换而提出的，主要代表方式有Scott，Leblanc、Kubler、Wood-bridge、阻抗匹配接线等。本次主要总结了常用牵引变压器的特点并建立数学模型，包括每种牵引变压器的原理结构、原次边电气量关系等，基于Matlab/Simulink软件建立牵引变压器仿真模型，并对牵引变压器在不同条件下的负序、谐波特性的进行了研究. 1 牵引变压器数学模型研究 1.1 YNd11接线 YNd11变压器接线原理如下图所示，如果忽略激磁电流及其漏阻抗压降，二次侧绕组ac相与一次侧绕组A相同相，cb相与C相同相。由于变压器一次侧绕组A，B，C相与电力系统的相序一致，A相滞后C相，对应的二次侧ac也滞后cb相[2]。其中Z为牵引端口对应变压器漏抗，和β相的端口电压。 1.2 Vv接线 Vv接线牵引变压器接线原理如图2所示。为二次侧空载相即α相图2 Vv接线牵引变压器设Vv接线变压器一次侧、二次侧绕组匝数分别为可得电流输入输出关系[3]：和，电压输入输出关系如下：图1 YNd11接线牵引变压器设YNd11接线变压器一次侧、二次侧绕组匝数分别为和假设变压器原边中性点接地，可以得出一次侧三相电流。，其中为牵引端口对应变压器漏抗，为二次侧空载相即α相和β相的端口电压。 1.3 Scott接线 Scott接线变压器（又称T形接法变压器）属于能完成三相-两相变换的平衡变压器，Scott接线牵引变压器接线原理如图3所示。图3 Scott牵引变压器接线原理图 1 61 设一次侧绕组BC的匝数为次侧绕组AD的匝数为，记，二次的绕组ad、bc的匝数为，则一。可得电流输入输出关系[4]：把一次侧绕组电流用相电流替换，即为：式中，为从三相端子流进变压器的电流。输出端口电压方程为：图6 YNd11接线牵引变压器两供电臂输出电压波形从电压输出波形中可以得到α供电臂电压波形超前β供电臂电压波形120°，在对称阻性负载下，两臂电流输出波形幅值相同，相位相差120°，满足理论值。 2.2 Vv接线牵引变压器 Vv 接线牵引变压器是由两个单相牵引变压器并联而成，仿真模型如图7所示.在仿真模型中牵引变压器T1和T2的原、次边变比设置为110kV/27.5kV。对，于
自动实验一——典型环节的MATLAB仿真 报告
班级 姓名 学号 XXXXXX电子与信息工程学院实验报告册 课程名称：自动控制原理实验地点： 实验时间同组实验人： 实验题目：典型环节的MATLAB仿真 一、实验目的： 1．熟悉MATLAB桌面和命令窗口，初步了解SIMULINK功能模块的使用方法。 2．通过观察典型环节在单位阶跃信号作用下的动态特性，加深对各典型环节响应曲线的理解。 3．定性了解各参数变化对典型环节动态特性的影响。 二、实验原理及SIMULINK图形： 1．比例环节的传递函数为22 12 11 ()2100,200 Z R G s R K R K Z R =-=-=- == 其对应的模拟电路及SIMULINK图形如图1-3所示。 2．惯性环节的传递函数为 2 21 121 121 2 ()100,200,1 10.21 R Z R G s R K R K C uf Z R C s =-=-=-=== ++ 其对应的模拟电路及SIMULINK图形如图1-4所示。 3．积分环节(I)的传递函数为 uf C K R s s C R Z Z s G1 , 100 1.0 1 1 ) ( 1 1 1 1 1 2= = - = - = - = 其对应的模拟电路及SIMULINK图形如图1-5所示。 图1-5 积分环节的模拟电路及及SIMULINK图形 图1-4 惯性环节的模拟电路及SIMULINK图形
4．微分环节(D)的传递函数为 uf C K R s s C R Z Z s G 10,100)(111112==-=-=-= uf C C 01.012=<< 其对应的模拟电路及SIMULINK 图形如图1-6所示。 5．比例+微分环节（PD ）的传递函数为 )11.0()1()(111212+-=+-=-=s s C R R R Z Z s G uf C C uf C K R R 01.010,10012121=<<=== 其对应的模拟电路及SIMULINK 图形如图1-7所示。 6．比例+积分环节（PI ）的传递函数为 )11(1)(11212s R s C R Z Z s G +-=+-=-= uf C K R R 10,100121=== 其对应的模拟电路及SIMULINK 图形如图1-8所示。 三、实验设备： 计算机 Matlab 软件 四、试验内容： 按下列各典型环节的传递函数，建立相应的SIMULINK 仿真模型，观察并记录其单位阶跃响应波形。 ① 比例环节1)(1=s G 和2)(1=s G ； ② 惯性环节11)(1+= s s G 和1 5.01)(2+=s s G ③ 积分环节s s G 1)(1= ④ 微分环节s s G =)(1 ⑤ 比例+微分环节（PD ）2)(1+=s s G 和1)(2+=s s G 图1-6 微分环节的模拟电路及及SIMULINK 图形 图1-7 比例+微分环节的模拟电路及SIMULINK 图形 图1-8 比例+积分环节的模拟电路及SIMULINK 图形曲线
MATLAB仿真课后习题
第一章习题 3.请指出以下的变量名（函数名、M文件名）中，哪些是合法的？Abc 2004x lil-1 wu_2004 a&b qst.u _xyz 解：合法的变量名有：Abc wu_2004 4．指令窗操作 （1）求[12+2×（7-4）]÷32的运算结果 解：>> [12+2*(7-4)]/3^2 ans = 2 （2）输入矩阵A=[1，2，3；4，5，6；7，8，9]，观察输出。解：>> A=[1,2,3;4,5,6;7,8,9] A = 1 2 3 4 5 6 7 8 9 （3）输入以下指令，观察运算结果； clear;x=-8:0.5:8; y=x'; X=ones(size(y))*x; Y=y*ones(size(x)); R=sqrt(X.^2+Y.^2)+eps;
Z=sin(R)./R; mesh(X,Y,Z); colormap(hot) xlabel('x'),ylabel('y'),zlabel('z') 解： 7．指令行编辑 （1）依次键入以下字符并运行：y1=2*sin(0.3*pi)/(1+sqrt(5)) 解：>>y1=2*sin(0.3*pi)/(1+sqrt(5)) y1 = 0.5000 (2)通过反复按键盘的箭头键，实现指令回调和编辑，进行新的计算；y2=2*cos(0.3*pi)/(1+sqrt(5)) 解：>>y2=2*cos(0.3*pi)/(1+sqrt(5)) y2 =
0.3633 11.编写题4中（3）的M脚本文件，并运行之。解：
第二章习题 1.在指令窗中键入x=1:0.2:2和y=2:0.2:1,观察所生成的数组。 解:>> x=1:0.2:2 x = 1.0000 1.2000 1.4000 1.6000 1.8000 2.0000 >> y=2:0.2:1 y = Empty matrix: 1-by-0 2．要求在[0，2π]上产生50个等距采样数据的一维数组，试用两种不同的指令实现。 解： y1=0:2*pi/49:2*pi y2=linspace(0,2*pi,50) 3.计算e -2t sint ，其中t 为[0，2π]上生成的10个等距采样的数组。 解：>> t=linspace(0,2*pi,10); x=exp(-2*t).*sin(t) x = 0 0.1591 0.0603 0.0131 0.0013 -0.0003 -0.0002 -0.0001 -0.0000 -0.0000 4.已知A=??????4321 , B=?? ????8765,计算矩阵A 、B 乘积和点乘. 解：>> A=[1,2;3,4]; B=[5,6;7,8]; x=A*B x =
倒立摆系统的建模及Matlab仿真资料
第1 页共11 页 倒立摆系统的建模及Matlab仿真 1.系统的物理模型 考虑如图(1)所示的倒立摆系统。图中，倒立摆安装在一个小车上。这里仅考虑倒立摆在图面内运动的二维问题。 图(1)倒立摆系统 假定倒立摆系统的参数如下。 摆杆的质量：m=0.1g l=1m小车的质量：摆杆的长度：2重力加速度：g=9.8m/M=1kg s摆杆的质量在摆杆的中心。 设计一个控制系统，使得当给定任意初始条件(由干扰引起)时，最大超调量?≤10%，调节时间ts ≤4s ，通过小车的水平运动使倒立摆保持在垂直位置。 2.系统的数学模型 2.1建立倒置摆的运动方程并将其线性化。 为简化问题，在数学模型中首先假设：1)摆杆为刚体；2）忽略摆杆与支点之间的摩擦；3）忽略小车与接触面间的摩擦。 ?）,在u设小车瞬时位置为z,摆心瞬时位置为（作用下，小车及摆均产生加速远 动，sin?lz根据牛顿第二定律，在水平直线远动方向的惯性力应与u平衡，于是有 22dzd?)?sinu?M?m(zl22dtdt???2????z(M?mml?)cos?mlusin? 即：??①
绕摆轴转动的惯性力矩与重力矩平衡，因而有． 第2 页共11 页 2??d??? sin??lcosm(z?lsinmgl)??2dt?????22???????即： nis?l?ocgcosincoszs?ls??② 以上两个方程都是非线性方程，为求得解析解，需作线性化处理。由于控制的目的是保持倒立摆直?2?????且可忽略则，立，在试驾合适的外力条件下，假定θ很小，接近于零时合理的，1sincos??,项。于是有 ???M?zm?u?ml??)(③ ????g?z?l??④联立求解可得1mg?u?z????MM 1)?m(M????u??MlMl 列写系统的状态空间表达式。2.2??T xx,x,x,，选取系统变量则 xx,x,xx?,42134123xx??211mgux???x?32MM x?x?431)(M?mu?x?x? 34MlMl 即00100????z??1mg??????000?z?????d MM??Bu?Ax?xux????????00001???dt????1gm?(M)????000??????? MlMl??????Cx?0?y?xx1001代入数据计算得到：0100????000?1??????T0D,?0??1BA?,?001,C100??1000??00011?? 11 页3 页共第 3.设计控制器3.1判断系统的能控性和稳定性 1100????0011????23BBAABAB?Q?故被控对象完全可控， rank()=4,Q kk??11?0?10??011?10???22???11?。出现大于零的特征值，故被，，0 解得特征值为 0由特征方程0??11I?A?)(控对象不稳定3.2确定希望的极点， 另一对为远极点，认为系统性能主要由主导，选其中一对为主导极点和希望的极点n=4ss21极点决定，远极点只有微小影响。根据二阶系统的关系式，先确定主导极点???42??1????10.?e??t1.67?有，闭环可得；取误差带，于是取，则6.?059?0.02.?0? pns??n2????1?js??=-10.8j,远极点选择使它和原点的距离大于主导极点与原点 距离主导极点为?n,21s??15倍，取的54,33.3采用状态反馈方法使系统稳定并配置极点 ??kkkk?k；状态反馈系统的状态方程，馈状态反的控制规律为为kxu??3102?，其
实验三__SIMULINK仿真实验
实验三 SIMULINK 仿真实验 一、实验目的 1.熟悉Simulink 的操作环境并掌握绘制系统模型的方法。 2.掌握Simulink 中子系统模块的建立与封装技术。 3.对简单系统所给出的数学模型能转化为系统仿真模型并进行仿真分析。 二、实验设备及条件 计算机一台（带有MATLAB7.0软件环境）。 三、实验内容 1.建立下图5-1所示的Simulink 仿真模型并进行仿真，改变Gain 模块的增益，观察Scope 显示波形的变化。 图3-1 正弦波产生及观测模型 2．利用Simulink 仿真下列曲线，取πω2=。 t t t t t t x ωωωωωω9sin 9 17sin 715sin 513sin 31sin )(++++=。 仿真参考模型如下图3-2，Sine Wave5模块参数设置如下图3-3，请仿真其结果。
图3-2 ()x t 的仿真参考模型图 图3-3 Sine Wave5模块参数设置图 3. 已知某控制系统的传递函数如题3-4图所示。试利用SIMULINK 建模仿真，并用示波器显示该系统的阶跃响应曲线。(注：系统中e －0.4 s 环节表示的是控制中的延时环节，可用SIMULINK 的连续系统模块库中的“Transport Delay”模块表示) 图3-4 4、已知某控制系统的传递函数如题3-5图所示。 试利用SIMULINK 建模，并实现以下功能： (1) 将已建模型转化为一个名为“mysys”的子系统； (2) 将已建子系统进行适当的封装； (3) 封装完毕后双击子系统图标，在弹出的属性设置窗口中对变量进行赋值(Tm = 0.5，Tp = 1)，并在模型中加入源模块和显示模块，观察系统的阶跃响应曲线。
MATLAB的建模和仿真
课程设计说明书 题目：基于Matlab的IIR滤波器设计与仿真班级：2012 级电气五班 姓名：王璐 学号：201295014178 指导教师：张小娟 日期：2015年 1 月12日
课程设计任务书
基于MATLAB的IIR滤波器设计与仿真 前言 数字信号处理（digital signal processing，DSP）是从20世纪60年代以来，随着信息学科和计算机学科的高速发展而迅速发展起来的一门新兴学科。数字信号处理是把信号用数字或符号表示的序列，通过计算机或通用（专用）信号处理设备，用数字的数值计算方法处理（例如滤波、变换、压缩、增强、估计、识别等），以达到提取有用信息便于应用处理的目的。数字信号处理系统有精度高、灵活性高、可靠性高、容易大规模集成、时分复用、可获得高性能指标、二维与多维处理等特点。正是由于这些突出的特点，使得它在通信、语音、雷达、地震测报、声呐、遥感、生物医学、电视、仪器中得到愈来愈广泛的应用。在数字信号处理中起着重要的作用并已获得广泛应用的是数字滤波器（DF，Digital Filter），根据其单位冲激响应函数的时域特性可分为两类：无限冲激响应IIR（Infinite Impulse Response）滤波器和有限冲激响应FIR（Finite Impulse Response）滤波器。MATLAB的基本数据单位是矩阵，它的指令表达式与数学、工程中常用的形式十分相似，故用MATLAB来结算问题要比用C，FORTRAN等语言完成相同的事情简捷得多，并且MATLAB也吸收了像Maple等软件的有点，使MATLAB成为一个强大的数学软件，在新的版本中也加入了对C，FORTRAN，C++，JA V A的支持。可以直接调用，用户也可以将自己编写的实用程序导入到MATLAB函数库中方便自己以后调用。 1 数字滤波器概述 数字滤波器是对数字信号实现滤波的线性时不变系统。数字滤波实质上是一种运算过程，实现对信号的运算处理。输入数字信号（数字序列）通过特定的运算转变为输出的数字序列，因此，数字滤波器本质上是一个完成特定运算的数字计算过程，也可以理解为一台计算机。描述离散系统输出与输入关系的卷积和差分方程只是给数字信号滤波器提供运算规则，使其按照这个规则完成对输入数据的处理。时域离散系统的频域特性：Y(eωj)=X(eωj)H(eωj) 其中Y(eωj)、X(eωj)分别是数字滤波器的输出序列和输入序列的频域特性（或称为
Matlab控制系统计算机辅助设计
实验目录 实验一：Matlab环境熟悉与基本运算（设计型）实验二：Matlab语言程序设计(设计型） 实验三：控制系统模型的建立（设计型） 实验四：Simulink仿真入门（验证型） 实验五：控制系统时域仿真分析（设计型） 实验六：Simulink环境下时域仿真 实验七：控制系统根轨迹仿真分析 实验八：控制系统频域仿真分析（设计型）
1、矩阵运算（1）矩阵的乘法 A=[1 2;3 4]; B=[5 5;7 8]; y=A^2*B y = 105 115 229 251 （2）矩阵除法 A=[1 2 3;4 5 6;7 8 9]; B=[1 0 0;0 2 0;0 0 3]; y1=A\B 警告: 矩阵接近奇异值，或者缩放错误。结果可能不准确。RCOND = 1.541976e-18。y1 = 1.0e+16 * -0.4504 1.8014 -1.3511 0.9007 -3.6029 2.7022 -0.4504 1.8014 -1.3511 y2=A/B y2 = 1.0000 1.0000 1.0000 4.0000 2.5000 2.0000 7.0000 4.0000 3.0000 （3）矩阵的转置及共轭转置 A=[5+i,2-i,1;6*i,4,9-i]; y1=A.' y1 = 5.0000 + 1.0000i 0.0000 + 6.0000i 2.0000 - 1.0000i 4.0000 + 0.0000i 1.0000 + 0.0000i 9.0000 - 1.0000i y2=A' y2 = 5.0000 - 1.0000i 0.0000 - 6.0000i 2.0000 + 1.0000i 4.0000 + 0.0000i 1.0000 + 0.0000i 9.0000 + 1.0000i 实验名称：Matlab环境熟悉与基本运算（设计型）
MATLAB仿真实验全部
实验一 MATLAB 及仿真实验（控制系统的时域分析） 一、实验目的 学习利用MATLAB 进行控制系统时域分析，包括典型响应、判断系统稳定性和分析系统的动态特性； 二、预习要点 1、 系统的典型响应有哪些？ 2、 如何判断系统稳定性？ 3、 系统的动态性能指标有哪些？ 三、实验方法 （一） 四种典型响应 1、 阶跃响应： 阶跃响应常用格式： 1、)(sys step ；其中sys 可以为连续系统，也可为离散系统。 2、),(Tn sys step ；表示时间围0---Tn 。 3、),(T sys step ；表示时间围向量T 指定。 4、),(T sys step Y =；可详细了解某段时间的输入、输出情况。 2、 脉冲响应： 脉冲函数在数学上的精确定义：0 ,0)(1)(0 ?==?∞ t x f dx x f 其拉氏变换为：)()()()(1 )(s G s f s G s Y s f === 所以脉冲响应即为传函的反拉氏变换。 脉冲响应函数常用格式： ① )(sys impulse ； ② ); ,();,(T sys impulse Tn sys impulse ③ ),(T sys impulse Y = （二） 分析系统稳定性 有以下三种方法： 1、 利用pzmap 绘制连续系统的零极点图； 2、 利用tf2zp 求出系统零极点； 3、 利用roots 求分母多项式的根来确定系统的极点 （三） 系统的动态特性分析 Matlab 提供了求取连续系统的单位阶跃响应函数step 、单位脉冲响应函数impulse 、零输入响应函数initial 以及任意输入下的仿真函数lsim.
基于Matlab、Simulink 的AM通信系统仿真设计与研究
天津理工大学计算机与通信工程学院通信工程专业设计说明书 基于Matlab/Simulink 的AM通信系统仿真设计与研究 姓名杜艳玮 学号 20092177 班级 09通信-2 指导老师赵健 日期2012/12/16
目录 摘要 (3) 第一章前言 (4) 1.1专业设计任务及要求 (4) 1.2 Matlab简介 (4) 1.4 通信系统模型 (6) 第二章 AM调制原理及仿真 (7) 2.1 AM调制原理 (7) 2.1.1 AM介绍 (7) 2.1.2 AM调制原理框图 (8) 2.2 AM调制方式的Matlab仿真 (8) 2.2.1 载波信号分析 (8) 2.2.2 AM调制 (9) 2.3 AM调制方式Matlab-simulink仿真 (10) 2.3.1 仿真框图 (10) 2.3.2 仿真结果 (11) 第三章 AM解调 (13) 3.1 AM解调原理 (13) 3.2 AM解调方式Matlab仿真 (13) 3.2.1 滤波前AM解调信号波形 (13) 3.2.2 AM调制信号解调 (15) 3.3 AM解调方式的Matlab-simulink仿真 (17) 3.3.1 仿真框图 (17) 3.3.2 仿真结果 (18) 第四章结论 (19) 参考文献 (20)
摘要 学习AM调制原理，AM调制就是由调制信号去控制高频载波的幅度，使之随调制信号作线性变化的过程。在波形上，幅度已调信号的幅度随基带信号的规律而呈正比地变化。解调方法利用相干解调。解调就是实现频谱搬移，通过相乘器与载波相乘来实现。通过相干解调，通过低通滤波器得到解调信号。相干解调时，接收端必须提供一个与接受的已调载波严格同步的本地载波，它与接受的已调信号相乘后，经低通滤波器取出低频分量，得到原始的基带调制信号。通过信号的功率谱密度的公式，得到功率谱密度。利用Matlab和Matlab-Simulink仿真建立AM调制的通信系统模型，用Matlab仿真程序画出调制信号、载波、已调信号、相干解调之后信号的波形以及功率频谱密度，分析所设计系统性能。用Matlab-Simulink仿真建立基于相干解调的AM仿真模型，详细叙述模块参数的设置，分析仿真结果。 关键字：AM调制相干解调 Matlab仿真 Matlab-Simulink仿真
Simulink 比值控制系统课后习题
Simulink 比值控制课后答案 作业题目： 在例一中如系统传递函数为43 ()151s G s e s -= +，其他参数不变，试对其进行单闭环比 值控制系统仿真分析，并讨论43 ()151 s G s e s -=+分母中“15”变化10%±时控制系统的鲁棒 性。 （1）分析从动量无调节器的开环系统稳定性。 由控制理论知，开环稳定性分析是系统校正的前提。系统稳定性的分析可利用Bode 图进行，编制MATLAB Bode 图绘制程序（M-dile ）如下： clear all close all T=15;K0=3;tao=4; num=[K0];den=[T,1]; G=tf(num,den,'inputdelay',tao); margin(G) 执行该程序得系统的Bode 图如图所示，可见系统是稳定的。幅值裕量为6.77dB ，对应增益为2.2。 -40-30-20-100 10M a g n i t u d e (d B )10 10101010 -2160 -1800-1440-1080-720 -3600P h a s e (d e g ) Bode Diagram Gm = 6.77 dB (at 0.431 rad/sec) , P m = 66.3 deg (at 0.189 rad/sec) Frequency (rad/sec) （2）选择从动量控制器形式及整定其参数。 根据工程整定的论述，选择PI 形式的控制器，即() I p K G s K s =+ 。本处采用稳定边界法整定系统。先让I K =0，调整p K 使系统等幅振荡（由稳定性分析图知在p K =2.2附近时系统震荡），即使系统处于临界稳定状态。
用MATLAB处理线性系统数学模型
实验一 用MATLAB 处理线性系统数学模型 [说明] 一个控制系统主要由被控对象、测量装置、控制器和执行器四大部分构成。MATLAB 软件的应用对提高控制系统的分析、设计和应用水平起着十分重要的作用。采用MATLAB 软件仿真的关键问题之一是在MATLAB 软件平台上怎样正确表示被控对象的数学模型。 [实验目的] 1．了解MATLAB 软件的基本特点和功能； 2．掌握线性系统被控对象传递函数数学模型在MATLAB 环境下的表示方法及转换； 3．掌握多环节串联、并联、反馈连接时整体传递函数的求取方法； 4． 掌握在SIMULINK 环境下系统结构图的形成方法及整体传递函数的求取方法； 5．了解在MATLAB 环境下求取系统的输出时域表达式的方法。 [实验指导] 一、被控对象模型的建立 在线性系统理论中，一般常用的描述系统的数学模型形式有： （1）传递函数模型——有理多项式分式表达式 （2）传递函数模型——零极点增益表达式 （3）状态空间模型（系统的内部模型） 这些模型之间都有着内在的联系，可以相互进行转换。 1、传递函数模型——有理多项式分式表达式 设系统的传递函数模型为 111011 1......)()()(a s a s a s a b s b s b s b s R s C s G n n n n m m m m ++++++++= =---- 对线性定常系统，式中s 的系数均为常数，且a n 不等于零。 这时系统在MATLAB 中可以方便地由分子和分母各项系数构成的两个向量唯一地确定，这两个向量常用num 和den 表示。 num=[b m ,b m-1,…,b 1,b 0] den=[a n ,a n-1,…,a 1,a 0]
simulink仿真实验报告
电机与拖动控制实验及其MATLAB仿真： 《电机与拖动控制实验及其MATLAB仿真》是2014年11月18日清华大学出版社出版的图书，作者是曹永娟。 内容简介： 本书分上、下两篇。上篇为电机与拖动控制实验教程，针对MCL 系列电机实验教学系统进行介绍，包括变压器、同步电机、异步电机、直流电机以及直流调速系统、交流调速系统拖动控制实验内容。 目录： 上篇电机与拖动控制实验 第1章电机实验装置和基本要求 1.1MCLⅡ型电机教学实验台 1.2实验装置和挂件箱的使用 1.2.1MCLⅡ型电机实验装置交流及直流电源操作说明 1.2.2仪表的使用 1.2.3挂件箱的使用 1.2.4交直流电机的使用 1.2.5导轨、测速发电机及转速计的使用 第2章电机与拖动控制实验基本要求和安全操作规程 2.1实验基本要求 2.2实验前的准备 2.3实验的进行 2.4实验报告
2.5实验安全操作规程 第3章变压器实验 3.1单相变压器 3.1.1实验目的 3.1.2预习要点 3.1.3实验项目 3.1.4实验设备及仪器 3.1.5实验方法 3.1.6实验报告 3.2三相变压器 3.2.1实验目的 3.2.2预习要点 3.2.3实验项目 3.2.4实验设备及仪器 3.2.5实验方法 3.2.6实验报告 3.3三相变压器的连接组和不对称短路3.3.1实验目的 3.3.2预习要点 3.3.3实验项目 3.3.4实验设备及仪器 3.3.5实验方法
3.3.6实验报告 3.3.7附录 3.4三相变压器的并联运行3. 4.1实验目的 3.4.2预习要点 3.4.3实验项目 3.4.4实验设备及仪器 3.4.5实验方法 3.4.6实验报告 第4章同步电机实验 4.1三相同步发电机的运行特性4.1.1实验目的 4.1.2预习要点 4.1.3实验项目 4.1.4实验设备及仪器 4.1.5实验方法 4.1.6实验报告 4.1.7思考题 4.2三相同步发电机的并联运行4.2.1实验目的 4.2.2预习要点 4.2.3实验项目
倒立摆系统的建模及Matlab仿真
倒立摆系统的建模及Matlab 仿真 1.系统的物理模型 考虑如图(1)面内运动的二维问题。 图(1)倒立摆系统 假定倒立摆系统的参数如下。 摆杆的质量：m=0.1g 摆杆的长度：l =1m 小车的质量： M=1kg 重力加速度：g=9.8m/2s 摆杆的质量在摆杆的中心。 设计一个控制系统，使得当给定任意初始条件(由干扰引起)时，最大超调量δ ≤10%，调节时 间ts ≤4s ，通过小车的水平运动使倒立摆保持在垂直位置。 2.系统的数学模型 2.1建立倒置摆的运动方程并将其线性化。 为简化问题，在数学模型中首先假设：1)摆杆为刚体；2）忽略摆杆与支点之间的摩擦；3）忽略小车与接触面间的摩擦。 设小车瞬时位置为z,摆心瞬时位置为（θsin l z +）,在u 作用下，小车及摆均产生加速远动，根据牛顿第二定律，在水平直线远动方向的惯性力应与u 平衡，于是有 u l z dt d m dt z d M =++)sin (22 22θ 即： u ml ml z m M =-++θθθθsin cos )(2&&&&& ① 绕摆轴转动的惯性力矩与重力矩平衡，因而有
θθθsin cos )sin (22mgl l l z dt d m =??? ????+ 即： θθθθθθθsin cos sin cos cos 22g l l z =-+&&&&& ② 以上两个方程都是非线性方程，为求得解析解，需作线性化处理。由于控制的目的是保持倒立摆直 立，在试驾合适的外力条件下，假定θ很小，接近于零时合理的，则1cos ,sin ≈≈θθθ，且可忽略θ θ2&项。于是有 u ml z m M =++θ&&&& )( ③ θθg l z =+&&&& ④ 联立求解可得 u Ml Ml m M u M M mg z 1)(1 -+=+- =θθθ&&&& 2.2列写系统的状态空间表达式。 选取系统变量4321,,,x x x x ， []T x x x x x 4321,,,=则 u Ml x Ml m M x x x u M x M mg x x x 1 )(134433221-+= =+-==&&&& 即 []Cx x x y Bu Ax u Ml M x Ml g m M M mg z z dt d x ===+=?????? ? ???????-+?????????? ??? ? +- =???? ????????=000110100)(0 010 0000000 1 1θθ&&& 代入数据计算得到： [][]0,0001,1010,01100 1000010000 1 0==-=? ? ??? ? ??? ???-=D C B A T
控制系统的MATLAB仿真和设计课后答案
第二章 1>>x=[15 22 33 94 85 77 60] >>x(6) >>x([1 3 5]) >>x(4:end) >>x(find(x>70)) 2>>T=[1 -2 3 -4 2 -3] ; >>n=length(T); >>TT=T'; >>for k=n-1:-1:0 >>B(:,n-k)=TT.^k; >>end >>B >>test=vander(T) 3>>A=zeros(2,5); >>A(:)=-4:5 >>L=abs(A)>3 >>islogical(L) >>X=A(L) 4>>A=[4,15,-45,10,6;56,0,17,-45,0] >>find(A>=10&A<=20) 5>>p1=conv([1,0,2],conv([1,4],[1,1]));
>>p2=[1 0 1 1]; >>[q,r]=deconv(p1,p2); >>cq='商多项式为'; cr='余多项式为'; >>disp([cq,poly2str(q,'s')]),disp([cr,poly2str(r,'s')]) 6>>A=[11 12 13;14 15 16;17 18 19]; >>PA=poly(A) >>PPA=poly2str(PA,'s') 第三章 1>>n=(-10:10)'; >>y=abs(n); >>plot(n,y,'r.','MarkerSize',20) >>axis equal >>grid on >>xlabel('n') 2>>x=0:pi/100:2*pi; >>y=2*exp(-0.5*x).*sin(2*pi*x); >>plot(x,y),grid on; 3>>t=0:pi/50:2*pi; >>x=8*cos(t); >>y=4*sqrt(2)*sin(t); >>z=-4*sqrt(2)*sin(t); >>plot3(x,y,z,'p');
matlab控制系统传递函数模型
MATLAB及 控制系统仿真 实验 班级：智能0702
姓名：刘保卫 学号： 06074053(18) 实验四控制系统数学模型转换及MATLAB实现 一、实验目的 熟悉MATLAB 的实验环境。 掌握MATLAB 建立系统数学模型的方法。 二、实验内容 （注：实验报告只提交第2 题） 1、复习并验证相关示例。 （1）系统数学模型的建立 包括多项式模型（Transfer Function，TF），零极点增益模型(Zero-Pole，ZP)，状态空间模型 （State-space,SS）； （2）模型间的相互转换 系统多项式模型到零极点模型（tf2zp）,零极点增益模型到多项式模型（zp2tf）,状态空间模 型与多项式模型和零极点模型之间的转换（tf2ss,ss2tf,zp2ss…）； （3）模型的连接 模型串联（series）,模型并联（parallel）,反馈连接（feedback） 2、用MATLAB 做如下练习。 （1）用2 种方法建立系统的多项式模型。 程序如下： %建立系统的多项式模型（传递函数） %方法一，直接写表达式 s=tf('s') Gs1=(s+2)/(s^2+5*s+10) %方法二，由分子分母构造 num=[1 2]; den=[1 5 10]; Gs2=tf(num,den) figure pzmap(Gs1) figure pzmap(Gs1) grid on 运行结果：
易知两种方法结果一样 Transfer function: s Transfer function: s + 2 -------------- s^2 + 5 s + 10 Transfer function: s + 2 -------------- s^2 + 5 s + 10 （2）用2 种方法建立系统的零极点模型和多项式模型。 程序如下： %方法一 s=tf('s') Gs1=10*(s+1)/((s+1)*(s+5)*(s+10)) % zpk模型 ZPK=zpk(Gs1) %方法二 % tf模型 num=[10 10]; den=conv([1 1],conv([1 5],[1 10])); Gs2=tf(num,den) % zpk模型 ZPK=zpk(Gs2) figure pzmap(Gs1) figure pzmap(Gs1) grid on 运行结果： 易知两种方法结果一样 Transfer function: s Transfer function: