电气系统仿真课件7 系统仿真(含电路)幻灯片课件.pptVIP

第7章 系统仿真 常微分方程的数值求解 微分方程模型;7.1 微分方程的求解 在现在数学研究和工程实践中，很多数学模型都是用微分方程确定的，很多基本方程本身就是一个微分方程，因此求微分方程非常重要，但是大部分的微分方程目前难以求得其解析解，因此人们只有利用计算机强大的计算功能来求其数值解。MATLAB主要使用龙格-库塔法求解微分方程。 在控制系统仿真中，常用的求微分方程数值解的函数是ode23和ode45。;1. ode23 在MATLAB中，函数ode23采用2-3阶龙格-库塔法求解微分方程。 [t,y]=ode23(odefun,tspan,y0) [t,y]=ode23(odefun,tspan,y0,options) odefun：定义微分方程的形式y’=f(t,y)； tspan=[t0,tfinal]：表示微分方程的积分限从t0（始 值）到tfinal（终值），该积分限也 可以是一些离散的点； y0:初始状态列向量； options：积分参数，包括‘RelTol’（相对误差）和 ‘AbsTol’（绝对误差），可省略。;例1：使用ode23函数求解常微分方程y’=-y+x2+4x+1，x=1时，y=1。 解：首先创建函数fun1.m function f=fun1(x,y) f=-y+x^2+4*x+1; 在命令窗口中输入 [x,y]=ode23(fun1,[1,4],1); dy=-y+x.^2+4*x+1; plot(x,y,x,dy); legend(y,dy);2. ode45 在MATLAB中，函数ode45采用普通4-5阶龙格-库塔法求解微分方程。其使用方法与ode23函数的使用方法基本相同。 ode45函数是大部分场合的首选算法，ode23函数主要适用于精度较低的情形。 例2：解经典非线性方程，范得波（Van der Pol）微分方程 （w=2）。;7.2 微分方程模型 7.2.1 方法描述 在MATLAB中，使用ode函数建立微分方程模型。 微分方程模型是数学模型的一种主要形式。当采用一阶微分方程的数值积分法进行数值计算时，应把高阶微分方程变换成n个一阶微分方程形式。对于微分方程而言，除了少数可以得到解析解外，大多数只能采用数值解法。 ;（1）分析：电容电压和电流的关系;（2）建立导数函数 function dy=cap(t,y) Vi=20; R=5; C=70e-6; dy=(Vi-y)/(R*C); （3）使用ode函数进行仿真，仿真时间0～0.006s，Vc初始值为0V。 [t,y]=ode45(cap,[0,0.006],0); plot(t,y) axis([0 0.006 0 25]) title(Vc-Time) xlabel(Time/sec) ylabel(Vc/V); 当电压源为直流电压源时，加载在电容上的电压随时间呈抛物线增大，稳态值为电源电压。电容电压在t=0时取得最小值，最小值为0；电容电压在t=0.0023s时达到最大值，为20V。;例4：用简单的LC谐振电路组成滤波器电路，其电路方程是二阶微分方程。观察该RLC电路中 的时域变化情况。假设电源为直流电压源Vi=20V，电阻R=5Ω，电容C=70μF，电感L=70mH。;（1）分析：根据电路分析，由基尔霍夫定律可以得出微分方程;（2）创建导数函数 function dy=RLC(t,y) Vi=20; R=5; C=70e-06; L=70e-03; dy=zeros(2,1); dy(1)=y(2)/C; dy(2)=(Vi-y(1)-R*y(2))/L; （3）使用ode函数进行仿真，仿真时间0～0.12s，Vc初始值为0V，I初始值为0A。;[t,y]=ode45(RLC,[0 0.12],[0;0]); figure(1) subplot(2,1,1); plot(t,y(:,1)); title(Vc-Time) xlabel(Time/sec) ylabel(Vc/V) sub