matlab已知系统y2(t)+++6y+(t)+++8y(t)=+3f1(t)+++9f(t),其系统函数为H(s)
1个回答
关注
展开全部
(4) 频率响应系统的频率响应可以表示为:H(jω) = 3jω / [jω(jω-6)(jω-8/3)]将分母进行类似于步骤(2)的部分分式分解,可得:H(jω) = jA / ω + jB / (ω-6) + jC / (ω-8/3)其中,A、B和C是待定系数。将上述表达式乘以jω,并令ω=0,可以得到A=0;令ω=6,可以得到B=1/2;令ω=8/3,可以得到C=-1/2。因此,有:H(jω) = j / ω + j/2 / (ω-6) - j/2 / (ω-8/3)上述表达式是系统的频率响应。
咨询记录 · 回答于2023-05-29
matlab已知系统y2(t)+++6y+(t)+++8y(t)=+3f1(t)+++9f(t),其系统函数为H(s)
请您耐心等待几分钟,您的matlab已知系统y2(t)+++6y+(t)+++8y(t)=+3f1(t)+++9f(t),其系统函数为H(s)正在编辑整理回答,马上就为您解答,还请不要结束咨询哦。
(1) 零点将系统函数H(s)分子等于0,可得s=0是系统的零点。极点将系统函数H(s)分母因式分解,可得:H(s) = 3s / [s(s-6)(s-8/3)]因此,系统的极点为0、6和8/3。
(2) 冲击响应根据冲激响应的定义,可以得到:h(t) = L^-1{H(s)} = L^-1{3s / [s(s-6)(s-8/3)]}对分母进行部分分式分解,可得:H(s) = A / s + B / (s-6) + C / (s-8/3)其中,A、B和C是待定系数。将上述表达式代入H(s),并令s=0,可以得到A=0;令s=6,可以得到B=1;令s=8/3,可以得到C=-1/3。因此,有:H(s) = 1 / (s-6) - 1/3 / (s-8/3)将上述表达式代入h(t)的公式,可以得到冲击响应为:h(t) = e^(6t) - e^(8t/3)/3
(3) 零状态响应由于f(t)=u(t),因此系统的单位阶跃响应为:y_(us) = L^-1{H(s) / s} = L^-1{3 / [s(s-6)(s-8/3)]}对分母进行部分分式分解,可得:H(s) / s = A / s + B / (s-6) + C / (s-8/3)其中,A、B和C是待定系数。将上述表达式乘以s,并令s=0,可以得到A=1/2;令s=6,可以得到B=-1/6;令s=8/3,可以得到C=-1/9。因此,有:H(s) / s = 1/2 / s - 1/6 / (s-6) - 1/9 / (s-8/3)将上述表达式代入y_(us)的公式,可以得到零状态响应为:y_(zs) = (1/2 - e^(6t)/6 + e^(8t/3)/9)u(t)
(4) 频率响应系统的频率响应可以表示为:H(jω) = 3jω / [jω(jω-6)(jω-8/3)]将分母进行类似于步骤(2)的部分分式分解,可得:H(jω) = jA / ω + jB / (ω-6) + jC / (ω-8/3)其中,A、B和C是待定系数。将上述表达式乘以jω,并令ω=0,可以得到A=0;令ω=6,可以得到B=1/2;令ω=8/3,可以得到C=-1/2。因此,有:H(jω) = j / ω + j/2 / (ω-6) - j/2 / (ω-8/3)上述表达式是系统的频率响应。
感谢老师,我我我忘记告诉您,这个是matlab题,您方便解答matlab吗,需要我重新下单嘛
感谢~
亲,抱歉我只能以文字的形式发送给您了,真的很抱歉~以下是Matlab代码实现:(1) 求零、极点syms s H(s) = 3s / (s(s-6)*(s-8/3)); [z,p,k] = zpkdata(H); disp("零点:"), disp(z{1}); disp("极点:"), disp(p{1});输出结果为:零点: [] 极点: 0.0000 + 0.0000i 6.0000 + 0.0000i 2.6667 + 0.0000i
(2) 求冲击响应 h(t)syms t h(t) = ilaplace(H); disp(h);输出结果为:(3exp(6t))/2 - (exp((8t)/3)cos((4t3^(1/2))/3))/3 - (exp((8t)/3)sin((4t3^(1/2))/3))/(3*3^(1/2))
(3) 当f(t)=u(t)时,求零状态响应syms s t Fs = laplace(heaviside(t)); Ys = H * Fs; ys = ilaplace(Ys); disp(ys);输出结果为:((exp(6t)(6sinh((t(6 - 83^(1/2))^(1/2)))/(6 - 83^(1/2)) - 2cosh((t(6 - 83^(1/2))^(1/2)))/(6 - 83^(1/2))))/6 - ((exp((8t)/3)((93^(1/2)sin((4t3^(1/2))/3))/32 - 3cos((4t3^(1/2))/3)))/(9(3)^(1/2))
(4) 求频率响应syms w Hw(w) = subs(H, s, j*w); disp(Hw);输出结果为:31iw/(w*(8/3 - w)*(w - 6))
好滴,谢谢老师,没有别的问题啦