信号处理仿真及应用 第六讲 第七讲.pptVIP

* 第六讲 电路仿真应用实例 本讲主要内容： 反复练习 * P18，一个二阶电路系统的仿真 输入为电压信号，包含电阻、电容、电感。含两个储能元件。 uc (t) iL (t) 请先复习! * 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] 请先复习! * 练习1. 试仿真得出一个幅度调制系统的输入输出波形。设输入被调制信号是一个幅度为2V，频率为1000Hz的余弦波，调制度为0.5，调制载波信号是一个幅度为5V，频率为10kHz的余弦波，所有余弦波的初相位为0. 这是一个什么系统？ * 练习2. 试仿真电容的充电过程(M文件编程)。 一个电压源通过电阻与电容充电。设t=0为初始时刻（初始时刻之前电路断开，不工作），电压源输出电压x(t)为单位阶跃函数，电容两端的电压为y(t),回路电流为并将电压源视为系统输出，电路的初始状态为y(0)。如图1给出了电路图和等价的系统模型。 ～ x(t) i(t) R C + - y(t) x(t) y(t) 图1 电容的充电电路以及等价系统 * 练习3.一个电路系统及其方框图模型如图2所示，将电压源v(t)视为输入信号f(t), 电容两端的电压u(t)视为输出信号y(t), 设系统的状态变量为电感元件上的电流以及电容元件两端的电压，即x1(t)=i(t), x2(t)=u(t)。试写出该电路的微分方程，并改写为系统状态方程，画出系统的实现方框图和信号