控制系统仿真实验报告
洗衣机模糊控制仿真实验报告
洗衣机模糊控制仿真实验报告一、实验目的本实验旨在通过对洗衣机运行过程的模糊控制仿真实验，帮助学生更好地了解模糊控制的基本原理和实现方法。
二、实验原理洗衣机模糊控制系统主要包括模糊控制器、模糊推理机和输出规则等三个部分。
模糊控制器是模糊系统的核心部分，其主要作用是将输入信号转化为模糊集，并将控制输出信号转化为真实输出信号。
模糊控制器的输入为洗衣机工作状态的一些参数，例如水位、温度等，输出为洗衣机运行状态的一些控制命令，例如加热、搅拌等。
模糊推理机是由一系列规则组成的系统，它负责根据输入的模糊集和一组先验规则，进行模糊推理，得到控制输出信号的模糊集，即模糊控制器的中间变量。
输出规则主要为控制输出信号的模糊集赋值，即将模糊集中各个元素映射到真实输出信号的取值范围内。
三、实验步骤1、建立洗衣机的模糊控制系统模型，包括模糊控制器、模糊推理机和输出规则等。
2、设置洗衣机的运行参数，例如水位、温度等，作为模糊控制器的输入。
3、根据洗衣机的运行状态，制定一组先验规则，作为模糊推理机的输入，并进行模糊推理。
4、根据模糊推理得到的控制输出信号的模糊集，进行输出规则的映射，得到洗衣机的真实控制命令。
5、根据洗衣机的控制命令，模拟洗衣机的工作流程。
6、对洗衣机的工作流程进行仿真实验，并记录实验结果。
四、实验结果分析经过多次实验，得到了洗衣机的模糊控制系统的优化参数，能够实现洗衣机的良好控制。
通过对实验结果的分析，可以发现，模糊控制系统可以有效地调节洗衣机的运行状态，使其在不同的工作状态下保持稳定且高效的运行。
同时，模糊控制系统也具有很强的适应性和鲁棒性，可以自适应地调节参数，应对各种不同的运行环境。
五、实验总结本实验通过模拟洗衣机的工作流程，对模糊控制系统的基本原理和实现方法进行了深入探究，能够有效地帮助学生掌握模糊控制系统的设计和应用方法。
同时，在实验过程中，也需要注意对实验数据和结论的分析和总结，以便更好地优化模糊控制系统的参数和性能，实现最佳控制效果。
电力拖动自动控制系统实验报告
电⼒拖动⾃动控制系统实验报告电⼒拖动⾃动控制系统实验实验⼀转速反馈控制直流调速系统的仿真⼀、实验⽬的1、了解MATLAB下SIMULINK软件的操作环境和使⽤⽅法。
2、对转速反馈控制直流调速系统进⾏仿真和参数的调整。
⼆、转速反馈控制直流调速系统仿真根据课本的操作步骤可得到如下的仿真框图：图 1 仿真框图1、运⾏仿真模型结果如下：图2 电枢电流随时间变化的规律图3 电机转速随时间变化的规律2、调节参数Kp=0.25 1/τ=3 系统转速的响应⽆超调但调节时间长3、调节参数Kp=0.8 1/τ=15 系统转速的响应的超调较⼤，但快速性较好实验⼩结通过本次实验初步了解了MATLAB下SIMULINK的基本功能，对仿真图的建⽴了解了相关模块的作⽤和参数设置。
并可将其⽅法推⼴到其他类型控制系统的仿真中。
实验⼆转速、电流反馈控制直流调速系统仿真⼀、实验⽬的及内容了解使⽤调节器的⼯程设计⽅法，是设计⽅法规范化，⼤⼤减少⼯作计算量，但⼯程设计是在⼀定近似条件下得到的，⽤MATLAB仿真可根据仿真结果对设计参数进⾏必要的修正和调整。
转速、电流反馈控制的直流调速系统是静、动态性能优良、应⽤最⼴泛的直流调速系统，对于需要快速正、反转运⾏的调速系统，缩短起动、制动过程的时间成为提⾼⽣产效率的关键。
为了使转速和电流两种负反馈分别起作⽤，可在系统⾥设置两个调节器，组成串级控制。
⼀、双闭环直流调速系统两个调节器的作⽤1)转速调节器的作⽤（1）使转速n跟随给定电压*mU变化，当偏差电压为零时，实现稳态⽆静差。
（2）对负载变化起抗扰作⽤。
（3）其输出限幅值决定允许的最⼤电流。
2)电流调节器的作⽤（1）在转速调节过程中，使电流跟随其给定电压*iU变化。
（2）对电⽹电压波动起及时抗扰作⽤。
（3）起动时保证获得允许的最⼤电流，使系统获得最⼤加速度起动。
（4）当电机过载甚⾄于堵转时，限制电枢电流的最⼤值，从⽽起⼤快速的安全保护作⽤。
当故障消失时，系统能够⾃动恢复正常。
控制工程基础仿真实验报告
本科生课程论文控制工程基础仿真实验报告实验一一阶系统的单位阶跃响应一、实验目的1、学会使用ATLABM编程绘制控制系统的单位阶跃响应曲线；2、掌握准确读取动态特征指标的方法；3、研究时间常数T对系统性能的影响；4、掌握一阶系统11Ts+时间响应分析的一般方法；5、通过仿真实验，直观了解各典型环节的时间响应和频率响应，巩固课程中所学的基本概念和基本原理；二、实验要求1、输入3个不同的时间常数T，观察一阶系统11Ts+的单位阶跃响应曲线的变化，绘制响应曲线图，并分析时间常数T对系统性能的影响。
2、若通过实验已测得一阶系统11Ts+的单位阶跃响应曲线，试说明如何通过该曲线确定系统的时间常数T。
三、实验内容（一）实验设备计算机；WINDOWS操作系统，并安装Matlab语言编程环境。
（二）实验原理通过对各种典型环节的仿真实验，可以直观的看到各种环节的时间响应和频率响应的图像。
通过对所得图像的分析可以得出各种参数如何影响系统的性能。
四、实验过程在Matlab平台对一阶系统11Ts+的单位阶跃响应进行仿真。
（1）输入3个不同的时间常数T，观察一阶系统单位阶跃响应曲线的变化，绘制响应曲线图，并分析时间常数T对系统性能的影响。
在Matlab中进行操作，其代码如下：1.num=1;2.den=[11];3.g=tf(num,den)4.5.g =6.7.18. -----9. s + 110.11.Continuous-time transfer function.12.13.>> step(g)14.hold on15.>> step(tf(1,[21]))16.>> step(tf(1,[41]))17.>> legend('T=1','T=2','T=4');（2）对于已测得的一阶系统的单位阶跃响应曲线，分析通过该曲线确定系统的时间常数T的方法。
PID仿真实验报告
PID仿真实验报告PID控制是一种经典的控制方法，被广泛应用于工业自动化控制系统中。
本次实验主要针对PID控制器的参数调整方法进行仿真研究。
实验目的：1.研究PID控制器的工作原理；2.了解PID参数调整的方法；3.通过仿真实验比较不同PID参数对系统控制性能的影响。
实验原理：PID控制器由比例（P）、积分（I）、微分（D）三个控制部分组成。
比例控制：输出与误差成比例，用来修正系统集成误差；积分控制：输出与误差的积分关系成比例，用来修正系统持续存在的静态误差；微分控制：输出与误差变化率成比例，用来修正系统的瞬态过程。
PID参数调整方法有很多种，常见的有经验法、Ziegler-Nichols法和优化算法等。
实验中我们使用经验法进行调整，根据系统特性来进行手动参数调整。
实验装置与步骤：实验装置：MATLAB/Simulink软件、PID控制器模型、被控对象模型。
实验步骤：1. 在Simulink中建立PID控制器模型和被控对象模型；2.设定PID控制器的初始参数；3.运行仿真模型，并记录系统的响应曲线；4.根据系统响应曲线，手动调整PID参数；5.重复第3步和第4步，直到系统的响应满足要求。
实验结果与分析：从图中可以看出，系统的响应曲线中存在较大的超调量和震荡，说明初始的PID参数对系统控制性能影响较大。
从图中可以看出，系统的响应曲线较为平稳，没有出现明显的超调和震荡。
说明手动调整后的PID参数能够使系统达到较好的控制效果。
总结与结论：通过本次实验，我们对PID控制器的参数调整方法进行了研究。
通过手动调整PID参数，我们能够改善系统的控制性能，提高系统的响应速度和稳定性。
这为工业自动化控制系统的设计和优化提供了参考。
需要注意的是，PID参数的调整是一个复杂的工作，需要结合具体的控制对象和要求进行综合考虑。
而且，不同的参数调整方法可能适用于不同的控制对象和场景。
因此，在实际应用中，需要根据具体情况选择合适的参数调整方法，并进行实验验证。
单闭环直流调速系统的设计与仿真实验报告
单闭环直流调速系统的设计与仿真实验报告摘要：本文基于基本原理和方法，设计和仿真了一个单闭环直流调速系统。
首先介绍了直流电机调速的基本原理，然后根据系统要求，设计了控制系统的结构和参数，包括PID控制器的参数调整方法。
接下来使用Matlab/Simulink软件进行系统仿真实验，对系统的性能进行评估。
最后根据仿真结果对系统进行分析和总结，并提出了可能的改进方法。
关键词：直流电机调速、单闭环控制系统、PID控制器、仿真实验一、引言直流电机广泛应用于机械传动系统中，通过调节电机的电压和电流实现电机的调速。
在实际应用中，需要确保电机能够稳定运行，并满足给定的转速要求。
因此，设计一个高性能的直流调速系统至关重要。
本文基于单闭环控制系统的原理和方法，设计和仿真了一个直流调速系统。
首先介绍了直流电机调速的基本原理，然后根据系统要求，设计了控制系统的结构和参数，并采用PID控制器进行调节。
接着使用Matlab/Simulink软件进行系统仿真实验，并对系统的性能进行评估。
最后根据仿真结果对系统进行分析和总结，并提出了可能的改进方法。
二、直流电机调速的基本原理直流电机调速是通过调节电机的电压和电流实现的。
电压变化可以改变电机的转速，而电流变化可以改变电机的转矩。
因此，通过改变电机的电压和电流可以实现电机的调速。
三、控制系统设计和参数调整根据系统的要求，设计一个单闭环控制系统，包括传感器、控制器和执行器。
传感器用于测量电机的转速，并将信息传递给控制器。
控制器根据测量的转速和给定的转速进行比较，并调节电机的电压和电流。
执行器根据控制器的输出信号来控制电机的电压和电流。
在本实验中，采用PID控制器进行调节。
PID控制器的输出信号由比例项、积分项和微分项组成，可以根据需要对各项参数进行调整。
调整PID控制器的参数可以使用试错法、频率响应法等方法。
四、系统仿真实验使用Matlab/Simulink软件进行系统仿真实验，建立直流调速系统的模型，并对系统进行性能评估。
仿真伺服程序实验报告
仿真伺服程序实验报告实验报告：仿真伺服程序一、引言伺服系统是一种用于控制和调节机械运动的自动控制系统，广泛应用于机器人、数控设备、航空航天等领域。
为了加深了解伺服系统的工作原理和参数调节方法，本实验采用仿真软件进行伺服程序的实验。
二、实验目的1. 理解伺服系统的结构和工作原理；2. 掌握伺服系统参数的调节方法；3. 利用仿真软件搭建伺服系统并进行参数调节；4. 分析伺服系统的动态性能和稳态性能。
三、实验步骤1. 搭建伺服系统仿真模型根据实验要求，采用仿真软件搭建伺服系统的数学模型。
包括伺服电机、负载、传感器等元件，并建立其间的物理联系和数学关系。
2. 参数调节方法根据实验要求，采用合适的参数调节方法对伺服系统进行参数调节。
常用的方法有试验法、根轨迹法等。
3. 仿真实验利用仿真软件对搭建好的伺服系统进行参数调节，并观察系统的动态响应和静态性能。
根据实验结果分析系统的稳定性、速度响应和位置精度等。
四、实验结果与分析通过仿真实验，在伺服系统仿真模型中分别采用试验法和根轨迹法进行参数调节。
观察到伺服系统的动态响应和静态性能都有明显的改善。
1. 试验法参数调节试验法是一种较为简单直观的参数调节方法。
通过逐渐调节参数，并观察系统的动态响应，根据实验结果进行调整。
在实验中，采用试验法进行参数调节后，伺服系统的速度响应和位置精度均有明显的提升。
2. 根轨迹法参数调节根轨迹法是一种基于系统的极点位置的参数调节方法。
通过绘制系统的根轨迹图，分析图形的形状和位置，调整参数使得系统稳定性和动态性能得到优化。
在实验中，采用根轨迹法进行参数调节后，伺服系统的速度响应和位置精度均有进一步提升。
综合分析两种参数调节方法的实验结果，可以得出伺服系统的参数调节对系统的动态响应和静态性能有着显著的影响。
合理的参数调节可以提升系统的稳定性、速度响应和位置精度，从而满足实际应用需求。
五、实验总结通过本次仿真实验，掌握了伺服系统的结构和工作原理，了解了伺服系统参数的调节方法，并利用仿真软件成功搭建了伺服系统的数学模型进行参数调节。
电力拖动自动控制系统Matlab仿真实验报告模板
电力拖动自动控制系统Matlab仿真实
验报告
电力拖动自动控制系统
---Matlab仿真实验报告
实验一二极管单相整流电路一．【实验目的】
1．经过对二极管单相整流电路的仿真，掌握由电路原理图转换成仿真电路的基本知识；
2．经过实验进一步加深理解二极管单向导通的特性。
图1-1 二极管单相整流电路仿真模型图
二．【实验步骤和内容】
1.仿真模型的建立
①打开模型编辑窗口；
②复制相关模块；
③修改模块参数；
④模块连接；
2.仿真模型的运行
①仿真过程的启动；
②仿真参数的设置；
3.观察整流输出电压、电流波形并作比较，如图1-2、1-3、1-4所示。
三．【实验总结】
由于负载为纯阻性，故输出电压与电流同相位，即波形相同，但幅值不等，如图1-4所示。
图1-2 整流电压输出波形图图1-3 整流电流输出波形图
图1-4 整形电压、电流输出波形图
实验二三相桥式半控整流电路
一．【实验目的】
1．经过对三相桥式半控整流电路的仿真，掌握由电路原理图转换成仿真电路的基本知识；
2．研究三相桥式半控整流电路整流的工作原理和全过程。
二．【实验步骤和内容】
1.仿真模型的建立：打开模型编辑窗口，复制相关模块，修改模块
参数，模块连接。
2.仿真模型的运行；仿真过程的启动，仿真参数的设置。
相应的参数设置：
（1）交流电压源参数U=100 V，f=25 Hz，三相电源相位依次延迟120°。
实验报告_控制器(3篇)
第1篇一、实验目的1. 理解控制器的基本原理和设计方法；2. 掌握常见控制器的实现方法；3. 学会使用控制器进行系统仿真和分析。
二、实验内容1. 控制器基本原理；2. 控制器设计方法；3. 控制器实现与仿真；4. 控制器性能分析。
三、实验原理控制器是自动控制系统中实现控制功能的核心部件，其作用是使被控对象的输出信号与期望信号之间保持一定的稳定关系。
控制器的设计方法主要包括比例-积分-微分（PID）控制、模糊控制、自适应控制等。
四、实验步骤1. 实验一：控制器基本原理（1）了解控制器的定义、功能及分类；（2）学习控制器的作用和基本原理；（3）分析控制器在不同控制策略下的特点。
2. 实验二：控制器设计方法（1）学习PID控制器的原理和参数整定方法；（2）学习模糊控制器的原理和设计方法；（3）分析自适应控制器的原理和特点。
3. 实验三：控制器实现与仿真（1）使用MATLAB/Simulink软件搭建控制系统模型；（2）根据控制器设计方法，实现PID控制器、模糊控制器和自适应控制器；（3）对控制器进行仿真分析，验证其性能。
4. 实验四：控制器性能分析（1）分析不同控制器在仿真过程中的性能表现；（2）比较不同控制器的优缺点，总结适用场景；（3）提出改进控制器的建议。
五、实验结果与分析1. 实验一：控制器基本原理通过学习，我们了解到控制器是自动控制系统的核心部件，其作用是实现被控对象的输出信号与期望信号之间的稳定关系。
控制器的基本原理包括比例控制、积分控制和微分控制。
2. 实验二：控制器设计方法（1）PID控制器：通过比例、积分和微分三个环节，对被控对象的误差进行控制。
在实际应用中，需要根据系统特性对PID参数进行整定，以达到最佳控制效果。
（2）模糊控制器：基于模糊逻辑，通过模糊推理实现控制。
模糊控制器具有鲁棒性强、易于实现等优点。
（3）自适应控制器：根据系统动态特性，实时调整控制器参数，以适应系统变化。
自适应控制器具有自适应性、抗干扰性强等特点。
电力拖动自动控制系统Matlab仿真实验报告
电力拖动自动控制系统---Matlab仿真实验报告实验一二极管单相整流电路一．【实验目的】1．通过对二极管单相整流电路的仿真，掌握由电路原理图转换成仿真电路的基本知识；2．通过实验进一步加深理解二极管单向导通的特性。
图1-1二极管单相整流电路仿真模型图二．【实验步骤和内容】1.仿真模型的建立1打开模型编辑窗口；2复制相关模块；3修改模块参数；4模块连接；2.仿真模型的运行1仿真过程的启动；2仿真参数的设置；3.观察整流输出电压、电流波形并作比较，如图1-2、1-3、1-4所示。
三．【实验总结】由于负载为纯阻性，故输出电压与电流同相位，即波形相同，但幅值不等，如图1-4所示。
图1-2整流电压输出波形图图1-3整流电流输出波形图图1-4整形电压、电流输出波形图实验二三相桥式半控整流电路一．【实验目的】1．通过对三相桥式半控整流电路的仿真，掌握由电路原理图转换成仿真电路的基本知识；2．研究三相桥式半控整流电路整流的工作原理和全过程。
二．【实验步骤和内容】1.仿真模型的建立：打开模型编辑窗口，复制相关模块，修改模块参数，模块连接。
2.仿真模型的运行；仿真过程的启动，仿真参数的设置。
相应的参数设置：（1）交流电压源参数U=100V，f=25Hz，三相电源相位依次延迟120°。
（2）晶闸管参数Rn=0.001Ω，Lon=0.0001H，Vf=0V，Rs=50Ω，Cs=250e-6F。
（3）负载参数R=10Ω，L=0H，C=inf。
（4）脉冲发生器的振幅为5V，周期为0.04s（即频率为25Hz），脉冲宽度为2。
图2-1三相桥式半控整流电路仿真模型图当α=0°时，设为0.0033s，0.0166s，0.0299s。
图2-2α=0°整流输出电压等波形图当α=60°时，触发信号初相位依次设为0.01s，0.0233s，0.0366s。
图2-3α=60°整流输出电压等波形图三.【实验总结】三相可控整流电路中，最基本的是三相半波可控整流电路，应用最为广泛的是三相桥式全控整流电路、双反星形可控整流电路以及十二脉波可控整流电路等，均可在三相半波的基础上进行分析。
电力拖动自动控制系统Matlab仿真实验报告
电力拖动自动控制系统---Matlab仿真实验报告实验一二极管单相整流电路一.【实验目的】1.通过对二极管单相整流电路的仿真,掌握由电路原理图转换成仿真电路的基本知识;2.通过实验进一步加深理解二极管单向导通的特性。
图1-1 二极管单相整流电路仿真模型图二．【实验步骤与内容】1.仿真模型的建立①打开模型编辑窗口;②复制相关模块;③修改模块参数;④模块连接;2.仿真模型的运行①仿真过程的启动;②仿真参数的设置;3.观察整流输出电压、电流波形并作比较,如图1-2、1-3、1-4所示。
三．【实验总结】由于负载为纯阻性,故输出电压与电流同相位,即波形相同,但幅值不等,如图1-4所示。
图1-2 整流电压输出波形图图1-3 整流电流输出波形图图1-4 整形电压、电流输出波形图实验二三相桥式半控整流电路一.【实验目的】1.通过对三相桥式半控整流电路的仿真,掌握由电路原理图转换成仿真电路的基本知识;2.研究三相桥式半控整流电路整流的工作原理与全过程。
二．【实验步骤与内容】1.仿真模型的建立:打开模型编辑窗口,复制相关模块,修改模块参数,模块连接。
2.仿真模型的运行;仿真过程的启动,仿真参数的设置。
相应的参数设置:(1)交流电压源参数U=100 V,f=25 Hz,三相电源相位依次延迟120°。
(2)晶闸管参数Rn=0、001 Ω,Lon=0、000 1 H,Vf=0 V,Rs=50 Ω,Cs=250e-6 F。
(3)负载参数R=10 Ω,L=0 H,C=inf。
(4)脉冲发生器的振幅为5 V, 周期为0、04 s ( 即频率为25 Hz), 脉冲宽度为2。
图2-1 三相桥式半控整流电路仿真模型图当α=0°时, 设为0、003 3s,0、016 6s,0、029 9 s。
图2-2 α=0°整流输出电压等波形图当α=60°时,触发信号初相位依次设为0、01s,0、0233s,0、0366s。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
控制系统仿真实验报告班级：测控1402班姓名：王玮学号：072018年01月实验一经典的连续系统仿真建模方法一实验目的:1 了解和掌握利用仿真技术对控制系统进行分析的原理和步骤。
2 掌握机理分析建模方法。
3 深入理解阶常微分方程组数值积分解法的原理和程序结构，学习用Matlab编写数值积分法仿真程序。
4 掌握和理解四阶Runge-Kutta法，加深理解仿真步长与算法稳定性的关系。
二实验内容:1. 编写四阶 Runge_Kutta 公式的计算程序，对非线性模型（3）式进行仿真。
（1）将阀位u 增大10％和减小10％，观察响应曲线的形状;（2）研究仿真步长对稳定性的影响，仿真步长取多大时RK4 算法变得不稳定（3）利用 MATLAB 中的ode45()函数进行求解，比较与（1）中的仿真结果有何区别。
2. 编写四阶 Runge_Kutta 公式的计算程序，对线性状态方程（18）式进行仿真（1）将阀位增大10％和减小10％，观察响应曲线的形状;（2）研究仿真步长对稳定性的影响，仿真步长取多大时RK4 算法变得不稳定（4）阀位增大10％和减小10％，利用MATLAB 中的ode45()函数进行求解阶跃响应，比较与（1）中的仿真结果有何区别。
三程序代码:龙格库塔:%RK4文件clccloseH=[,]';u=; h=1;TT=[];XX=[];for i=1:h:200k1=f(H,u);k2=f(H+h*k1/2,u);k3=f(H+h*k2/2,u);k4=f(H+h*k3,u);H=H+h*(k1+2*k2+2*k3+k4)/6;TT=[TT i];XX=[XX H];end;hold onplot(TT,XX(1,:),'--',TT,XX(2,:)); xlabel('time')ylabel('H')gtext('H1')gtext('H2')hold on水箱模型:function dH=f(H,u)k=;u=;Qd=;A=2;a1=;a2=;dH=zeros(2,1);dH(1)=1/A*(k*u+Qd-a1*sqrt(H(1)));dH(2)=1/A*(a1*sqrt(H(1))-a2*sqrt(H(2)));2编写四阶 Runge_Kutta 公式的计算程序，对线性状态方程（18）式进行仿真：1 阀值u对仿真结果的影响U=;h=1; U=;h=1;U=;h=1;2 步长h对仿真结果的影响:U=;h=5; U=;h=20;U=;h=39 U=;h=50由以上结果知,仿真步长越大,仿真结果越不稳定。
采用ode45算法程序如下:function dH=liu(t,H)k=;u=;Qd=;A=2;a1=;a2=;dH=zeros(2,1);dH(1)=1/A*(k*u+Qd-a1*sqrt(H(1)));dH(2)=1/A*(a1*sqrt(H(1))-a2*sqrt(H(2)));在命令窗口运行以下程序:[t,H]=ode45('liu',[1 200],[ ]);plot(t,H(:,1),['r','+'],t,H(:,2),['g','*'])u= u=u=用ode45与用龙格库塔法仿真结果基本一致。
2编写四阶 Runge_Kutta 公式的计算程序，对线性状态方程（18）式进行仿真：%RK4文件clcclearclosex=[,]';u=; h=5;TT=[];XX=[];for i=1:h:200k1=f2(x,u);k2=f2(x+h*k1/2,u);k3=f2(x+h*k2/2,u);k4=f2(x+h*k3,u);x=x+h*(k1+2*k2+2*k3+k4)/6;TT=[TT i];XX=[XX x];end;hold onplot(TT,XX(1,:),'--',TT,XX(2,:));xlabel('time')ylabel('x')gtext('x1')gtext('x2')hold on线性函数：function dx=f2(x,u)%线性Qd=;a1=;a2=;A=2;k=;dx=zeros(2,1);dx(1)=1/A*(k*u-x(1)/(2*sqrt/a1)+Qd);dx(2)=1/A*(x(1)/(2*sqrt/a1)-x(2)/(2*sqrt/a2));1 阀值u对仿真结果的影响:U=；h=1； U=；h=1；U=；h=1；2 步长h对仿真结果的影响:U=；h=5； U=；h=20；U=；h=35； U=；h=50；当步长为50时仿真结果开始不稳定，仿真步长越大,仿真结果越不稳定。
采用ode45算法程序如下:function dx=liu2(t,x)%ÏßÐÔQd=;u=;a1=;a2=;A=2;k=;dx=zeros(2,1);dx(1)=1/A*(k*u-x(1)/(2*sqrt/a1)+Qd);dx(2)=1/A*(x(1)/(2*sqrt/a1)-x(2)/(2*sqrt/a2));在命令窗口运行以下程序:>> [t,x]=ode45('liu2',[1 200],[ ]);>> plot(t,x(:,1),['r','+'],t,x(:,2),['g','*']) U=； u=；U=用ode45与用龙格库塔法仿真结果基本一致。
四思考题：1讨论仿真步长对稳定性和仿真精度的影响仿真步长越短，稳定性越高，精确度越高。
2你是怎样实现阀位增大和减小的，对于非线性模型和线性模型方法一样吗通过改变u来改变阀的开度，线性系统和非线性系统方法一样。
实验二面向结构图的仿真第一部分线性系统仿真一实验目的1．掌握理解控制系统闭环仿真技术。
2．掌握理解面向结构图的离散相似法的原理和程序结构。
3．掌握 MATLAB 中C2D 函数的用法，掌握双线性变换的原理。
二实验内容根据上面的各式，编写仿真程序，实现无扰动时给定值阶跃仿真实验1. 取K P = ，T i = 85 s T = 10s，ΔH2 S =H2set_ percent = 80, ΔQ d = 0，tend = 700，进行仿真实验，绘制响应曲线。
clcclear allA=2;ku=;H10=;alpha12 = sqrt(H10); alpha2 = sqrt(H20);R12=2*sqrt(H10)/alpha12; R2=2*sqrt(H20)/alpha2; H1SpanLo=0;H2SpanLo=0;H1SpanHi=;H2SpanHi=;Kp=;Ti=85;R12*AR12ad = 1/(A*R12);a1 = 1/(A*R12);a2 = 1/(A*R2);Kc=Kp/Ti;Kd = 1/A;K1 = ku/A;K2 = 1/(A*R12);uc(1)=0;ud(1)=0;u1(1)=0;u2(1)=0;xc(1)=0;xd(1)=0;x1(1)=0;x2(1)=0;yd(1)=0;yc(1)=0;y1(1)=0;y2(1)=0;nCounter = 70;T=10;k=1;deltaQd=0;H20_percent=(H20-H2SpanLo)/(H2SpanHi-H2SpanLo)*100; H2=80;tend = nCounter*T;for t=T:T:tendk=k+1;uc(k)= (H2 - (y2(k-1)+H20-H2SpanLo)/(H2SpanHi-H2SpanLo)*100)/100;ud(k)=deltaQd;u1(k)=yc(k-1);u2(k)=y1(k-1);xc(k) = xc(k-1) + Kc*T*uc(k-1);yc(k)=xc(k)+bc*Kc*uc(k);xd(k) = exp(-ad*T)*xd(k-1) + Kd/ad*(1-exp(-ad*T))*ud(k);yd(k)=xd(k); x1(k) = exp(-a1*T)*x1(k-1) + K1/a1*(1-exp(-a1*T))*u1(k);y1(k)=x1(k); x2(k) = exp(-a2*T)*x2(k-1) + K2/a2*(1-exp(-a2*T))*u2(k);y2(k)=x2(k); endHlevel(:,1)=(y1+H10-H1SpanLo)/(H1SpanHi-H1SpanLo)*100;Hlevel(:,2)=(y2+H20-H2SpanLo)/(H2SpanHi-H2SpanLo)*100;yc=(yc+*100;y2sp=H2*ones(size(y1'));yv=yc;textPositionH1=max(Hlevel(:,1));textPositionH2=max(Hlevel(:,2));H2Steady=Hlevel(size(Hlevel(:,1),1),1)*ones(size(y1')); xmax=max(0:T:tend);xmin=0;ymax=110;ymin=50;scrsz = get(0,'ScreenSize');gca=figure('Position',[5 10 scrsz(3)-10 scrsz(4)-90]);%gca=figure('Position',[5 10 scrsz(3)/2 scrsz(4)/])set(gca,'Color','w');plot(0:T:tend,Hlevel(:,1),'r','LineWidth',2)hold onplot(0:T:tend,Hlevel(:,2),'b','LineWidth',2)hold onplot(0:T:tend,yv,'k','LineWidth',2)hold onplot(0:T:tend,y2sp,'g','LineWidth',2)hold onplot(0:T:tend,H2Steady,'y','LineWidth',2)line([tend/2 tend/2+27],[(ymax-ymin)/2+ymin-(ymax-ymin)/10(ymax-ymin)/2+ymin-(ymax-ymin)/10],'Color','r','LineWidth',6)text(tend/2+27,(ymax-ymin)/2+ymin-(ymax-ymin)/10,' 第一个水箱的液位H1','FontSize',16)line([tend/2 tend/2+27],[(ymax-ymin)/2+ymin-(ymax-ymin)/6(ymax-ymin)/2+ymin-(ymax-ymin)/6],'Color','b','LineWidth',6)text(tend/2+27,(ymax-ymin)/2+ymin-(ymax-ymin)/6,' 第二个水箱的液位H2','FontSize',16)line([tend/2 tend/2+27],[(ymax-ymin)/2+ymin-(ymax-ymin)/(ymax-ymin)/2+ymin-(ymax-ymin)/],'Color','g','LineWidth',6)text(tend/2+27,(ymax-ymin)/2+ymin-(ymax-ymin)/,' 第二个水箱的液位给定值','FontSize',16)line([tend/2 tend/2+27],[(ymax-ymin)/2+ymin-(ymax-ymin)/(ymax-ymin)/2+ymin-(ymax-ymin)/],'Color','k','LineWidth',6)text(tend/2+27,(ymax-ymin)/2+ymin-(ymax-ymin)/,'阀位变化情况','FontSize',16)axis([xmin xmax ymin ymax]);text(tend/5,ymax+,' 实验二不考虑阀位饱和特性时的控制效果','FontSize',22) grid2. 用 MATLAB 求出从输入到输出的传递函数，并将其用c2d 函数，利用双线性变换法转换为离散模型，再用dstep()函数求离散模型的阶跃响应，阶跃幅值为3。