23子系统创建与封装.PPT

仿真应用实例 1、电路仿真应用 2、信号的频谱分析 1、电路仿真应用 P18，一个二阶电路系统的仿真 输入为电压信号，包含电阻、电容、电感。含两个储能元件。 1）零状态响应：系统在t=-1秒时刻的初始状态为零，输入信号为阶跃信号f(t)=u(t)，输出为电容电压，试对时间t=[-1，10]秒的输出信号波形进行方针求解。 2）零输入响应：系统在t=-1秒时刻电容电压uc 为1V，当t=-1秒时刻信号输入二端闭合，试对时间t=[-1，10]秒的输出信号波形进行仿真求解。 3）全响应：系统在t=-1秒时刻电容电压为为uc =1V，输入信号为阶跃信号f(t)=u(t)，试对时间t=[-1，10]秒的输出信号波形进行仿真求解。 电路方程 KVL方程： 标准形式： 状态方程 iL (t) = x1 ， uC (t) = x2 微分方程的MATLAB函数格式 function xdot=方程函数名（t, x, flag, 附加参数） 时间变量矩阵： t 方程的状态变量矩阵：x 状态变量对时间的一阶导数矩阵：xdot 方程参数：附加参数 必需的，以求符合求解函数“ode45”的要求的标志变量： flag 微分方程组函数 function xdot=funcforex123(t,x,flag,R,L,C) xdot=zeros(2,1); xdot(1)=-R/L*x(1)-1/L*x(2)+1/L*f(t); xdot(2)=1/C*x(1); function in=f(t); In=(t0)*1; ode45的用法 [t, x]=ode45(方程函数名，仿真时间范围，状态变量初始值，算法选项，附加参数） 初始变量：电感电流，电容电压 这里 R，L，C为附加参数 1）零状态响应 输入信号：in= (t0)*1 (阶跃函数); 方程中系统状态为： [0；0]; 求解 clear; clc L=1; C=0.1; for R=[1.5 3 5] [t, x]=ode45(funcforex123,[-1,10],[0;0],[],R,L,C) figure(1); plot(t,x(:,1)); hold on; xlabel(time sec); text(0.9,017,\leftarrow i_L(t)); grid; figure(2); plot(t,x(:,2)); hold on; xlabel(time sec); text(0.5,03,\leftarrowu_C(t)); grid; end MATLAB提供的求解微分方程的指令 ode45、ode23、ode15s 适用场合不同 2）t=-1输入信号为零 修改输入信号：in=0; 系统状态改为：[0；1] 3）全响应情况 修改输入信号：in= (t0)*1 (阶跃函数); 系统状态为： [0；1] 信号的频谱分析 频谱分析原理 连续时间信号 ，其频谱分析采用连续时间傅立叶变换（CTFT）进行： 其中 离散傅立叶变换（DFT） N点离散序列 ，设其N点离散傅立叶变换（DFT）序列为 ，则 其中 利用DFT来近似计算CTFT 首先对连续信号离散化：确定采样时间间隔及截取的时间段。 若 的绝大部分能量集中在 ，且的非零值集中在时间段 则CTFT中的积分作近似求和计算： 对比DFT计算公式，显然有： 3 Simulink动态仿真 Simulink环境下的仿真运行 仿真参数对话框 点击Simullink模型窗simulation菜单下的Parameters命令，弹出仿真参数对话框如右图所示。它共有5页，用得较多的主要是Solver页和Workspace I/O页。 4 数字信号处理 交互式滤波器设计工具 滤波器设计函数 交互式界面的滤波器设计工具 fdatool 5 Smulink系统仿真实例 下图所示，为简化后的飞行控制系统，试建立此动态系统的Simulink模型并进行简单的仿真分析。系统输入input为单位阶跃，各参数如下： 简化后的飞行控制系统 具体要求如下： 1）采用自顶向下的设计思路。 2）对虚线框中的控制器采用子系统技术。 3）用同一示波器显示