微机保护工作原理.pptVIP

Subplot(311) stem(n,x) title(input sequence);ylabel(x) Subplot(312) stem(n,h) title(impulse respond sequence);ylabel(h) Subplot(313) stem(n,y) title(output sequence); xlabel(n);ylabel(y) 卷积计算举例 因果系统 h(n) = 0 , ( n 0 ) 稳定系统 差分方程 ：离散系统的一种基本数学描述 系数 ak,bk为常数 , x(n),y(n) 为系统的输入、输出序列。 计算机保护中主要采用以下两种解法： 在给定输入和给定的初始条件下，用递推方法求解线性差分方程。 保护中的实时故障处理程序中的数字信号滤波与保护算法，基本上采用这种方法。 变换域方法。即利用Z变换的方法。 这种方法简便有效，是用来分析数字滤波器及保护算法的基本方法。 三、Z变换 Z变换定义 离散序列的Z变换（双边Z变换 ） 因果序列Z变换（单边Z变换 ） 收敛域、性质和计算方法，参阅相关资料 2．Z变换的应用 Z变换求解差分方程 求Z变换，解得： 两边同时取反变换，就能够解得输出序列。 系统函数（传递函数） 对于线性时不变系统，由h(n)确定输入的响应 ： y (n) = x(n) ? h(n) 取Z变换，由离散卷积定理： Y (Z) = X(Z)? H(Z) 得系统函数 H(Z) = Y (Z) / X(Z) 即：单位冲击响应的Z变换 传递函数的一般形式 因果稳定系统: 传递函数的全部极点必须在Z平面上的单位园内 ck、dk分别是Z平面上的零点和极点 传递函数形式： 零极点形式： 采样过程与离散信号的频谱 一、时域理想采样 时域采样过程: 采样脉冲信号与连续信号相乘 x s(t) = x(t) ? s(t) 实际采样脉冲信号：周期矩形脉冲 理想采样过程: 矩形脉冲的宽度趋于零 反映采样过程的本质特征 二、离散信号的频谱 x (t) 的傅立叶变换为 X(?) x s(t) 的傅立叶变换为 Xs(?) 根据频域卷积定理 理想采样信号的傅立叶变换为 因此有 由δ函数的筛选特性，即有离散采样信号的频谱： 连续信号频谱的周期延拓 注意： 以上的讨论中，并未限定连续信号为周期信号或非周期信号 因此： 周期信号的离散采样信号频谱为 周期信号频谱的周期延拓 非周期信号的离散采样信号频谱为 非周期信号频谱的周期延拓 时域离散化 ? 频域周期性 时域周期性 ? 频域离散化 离散周期信号的频域变换特征 时域离散化 ? 频域周期性 时域非周期性 ? 频域连续性 离散非周期信号的频域变换特征 离散信号的频谱分析举例 例：时域连续信号为Xa(t)=Ae-atsin(?0t), A=400, a=50?*?2,?0=a，分别以采样频率fs=1000Hz,400Hz 对Xa(t)进行等间隔采样，分析X(t)及Xs(n)的幅频特性。 说明： 严格地讲，Matlab无法计算连续函数的频谱特性。 工程上近似计算时，当采样频率足够大，则频谱混叠忽略不计，从而可对采样序列进行傅里叶变换，得到连续函数的频谱特性。 为了便于比较，需要计算幅度归一化的频谱特性。 采样序列的傅里叶变换 根据定义 用求和代替积分 设置频率序列 [ wi ] , 就可求出一系列 X(wi)=x*exp(-j*wi*n’)。 因为x是行向量，n’及exp(-j*wi*n’)是同长的列向量，两者的向量点乘就包含了逐项相乘后的连加运算，得出一个标量，即该频率点的频率响应。 不同的wi可用for循环来计算，但高效计算方法是Matlab的矩阵计算： 将w设成行向量，放在n’之后，n’*w及 exp(-j*n’*w)就成为一个矩阵，其行数与x相同，列数与w的长度相同。 用x左乘exp(-j*n’*w)以后，各列分别得出相应wi处的响应，最后得到采样序列的傅里叶变换： X(w)= x*exp(-j*n’*w) 采样序列的傅里叶变换子程序 func