系统分析基础实验.docxVIP

班级 电气 姓名 学号 实验 1 2 3 地点 日期 指导教师 成绩 系统分析基础实验一信号与系统程序计算与分析 1实验目的：用M语言研究信号与系统时域与频域响应。 2实验步骤 2. 1 连续系统响应计算 已知:H(s)=1/(s2+1.4 初值 y(O)=l, f (0)=0. 5 求响应。 参考程序： Hb=[ ]; Ha=[ ]; Fb=[ ]; Fa=[ Yzsb=conv (Hb, Fb); Yzsa二conv (Ha, Fa); [rzs, pzs, kzs]二residue(Yzsb, Yzsa) y0=l; pyO=O.5; Yzib二y0*[ ]+pyO*[ ] ; Yzia = Ha ; [rzi, pzi, kzi]二residue (Yzib, Yzia) t=0:0. 01:20; yzs=impulse(Yzsb, Yzsa, t); yzi=impulse(Yzib, Yzia, t); y二yzs+yzi; plot(t, yzs, B ,t, yzi, g , t, y, s+1),输入 f (t)=sin(2. t)u(t) F(s)二 1； %H,F算子（L变换）的分了分母 % 0状态响应分子分母. %部分分式展开 % 0输入响应计算 %定义时间范围 % 0状态响应 % 0输入响应 %全响应 r ） %显示 运算结果 响应曲线 2. 2系统频率特性 (1)棊本环节频率特性 1 积分H(s)=l/s: freqs([1], [1,0]) 特点： 2 微分 H(s)=s: freqs([1,0], [1]) 特点: 3 一阶惯性 H(s)=l/(s+l): freqs([l], [1,1]) 特点： 4 二阶小阻尼 H(s)=1/(s2+0. ls+1) freqs([1], [1,0. 1,1]) 特点: ⑵4阶巴氏低通Q c二1.频率I特性 [BL, AL]=butter (4, 1, s ), freqs (BL, AL) BL AL 传 函 H(s) 频 率 特 性 ⑶ 4阶巴氏带通Q c二0. 5_, 1. 5—频率特性 [BP, AP]=butter(4, [0. 5,1.5]/ s ), freqs (BP, AP) 频 率 特 性 2.3离散系统响应让算 已知:H(z)=l/(l-0.2 z1+ 0. 5z),输入 x (n) =0. 5 nu(n) X(z) = 初值 y (-1)二 l,y(-2)二 0.5 的响应 参考程序： Hb=[ ]; Ha=[ ]； Fb=[ ]; Fa=[ Yzsb二conv(Hb, Fb);Yzsa=conv(Ha, Fa); [rzs,pzs,kzs]=residuez (Yzsb,Yzsa); rzs, abs (pzs), angle(pzs) yl=l; y2=0. 5; Yzib二y0*[ ]+py0*[ ] ; Yzia = Ha [rzi,pzi,kzi]=residuez (Yzib,Yzia); rzi, abs(pzi), angle(pzi) n二0:20; yzs=impz(Yzsb, Yzsa, n); yzi=impz (Yzib, Yzia, n); y二yzs+yzi; plot (t, yzs, J B ,t, yzi, g , t, y, r ];% H , F算子(L变换)的分子分母 % 0状态响应分子分母. %部分分式展开 % 0输入响应计算 f %部分分式展开 %定义时间范围 % 0状态响应 % 0输入响应 %全响应 ) %显示 rzs pzs yzs (n)= rzi pzi yzi (n) y (n) 响 应 曲 线 2.4离散系统频域特性 (1)基本环节频率特性 1 积分 H(z)=l/(1-z-1): freqz([l], [1,-1]) 特点： 2 微分 H(z)= (1-z) freqz([l,-l], [1])特点: 3 二点平均 H(z)= (l+zJ/2: freqz([1,1], 2) 特点: 4 一阶 H(z) = (l-0. 9)/(l-0. 9z-1) freqz([1-0. 9], [1,-0. 9]) 特点： ⑵4阶巴氏低通0 / H =0. 5_频率特性 (3) [BL, AL]=butter (4, 0. 5_), freqz (BL, AL) BL AL H(z) 频 率 特 性 ⑷4阶巴氏带通()/ Ji =0. 3_, 0. 6—频率特性 [BP, AP]=butter (4, [0. 3—, 0. 6_]), freqz (BP, AP) 频 率 特 性 *附：M语言有-淀符号运算能力 例：Lb L? ? syms t s %定义变屋符号 ? l