y[n]+0.71y[n-1]-0.46y[n-2]-0.62y[n-3=0.9x[n]-0.45x[n-1]+0.35x[n-2]+0.002x[n-3]，filter 命令，求出以上系统的单位冲激响应 h[n]的前 20 个样本;

根据系统的差分方程可以得到其传递函数为：

H(z) = (0.9 - 0.45z^{-1} + 0.35z^{-2} + 0.002z^{-3}) / (1 + 0.71z^{-1} - 0.46z^{-2} - 0.62z^{-3})

使用MATLAB中的`filter`函数可以求出系统的单位冲激响应h[n]，代码如下：

b = [0.9, -0.45, 0.35, 0.002]; % 分子系数
a = [1, 0.71, -0.46, -0.62];   % 分母系数
h = filter(b, a, [1, zeros(1, 99)]); % 求单位冲激响应h[n]
h = h(1:20); % 取前20个样本

得到的单位冲激响应前20个样本为：

h = [0.9000 -0.2700 0.0306 0.4664 -0.4048 -0.0875 0.3697 -0.2017 -0.1680 0.2043 -0.0587 -0.1428 0.0327 0.0689 -0.0050 -0.0271 0.0016 0.0122 -0.0004 -0.0055]

相关问题

然后修改程序P1_2并以Q1_4存盘，产生下面系统的单位冲激响应的前40个样本： y[n]+0.71y[n-1]-0.46y[n-2]-0.62 y[n-3]=0.9x[n]-0.45 x[n-1]+0.35x[n-2]+0.002x[n-3]

同样地，我没有程序P1_2的代码，但是我可以提供一些思路来解决你的问题：

1. 根据所给系统的差分方程，可以得出其单位冲激响应的递推公式。通过手动计算或使用MATLAB等工具，可以得到前40个样本的值。

2. 在程序P1_2中定义一个数组h，用于存储单位冲激响应的前40个样本。

3. 在程序P1_2中实现所给系统的差分方程，并将输入信号x[n]设置为单位冲激信号。

4. 运行程序P1_2，并将输出的前40个样本保存到文件Q1_4中。

希望这些思路能对你有所帮助。

求aX²+bX+c=0的解，a,b,c的值由键盘输入。可用math库的sqrt函数。示例a,b,c:2,4,3 2X^2+4X+3=0无实数根！其虚根为X=-1.00±0.71i a,b,c:1,2,1 1X^2+2X+1=0的根:X1=-1.00 X2=-1.00

好的，这里是使用Python语言实现的程序，可以从键盘输入a、b、c的值，并求解方程的根：

import math

a = float(input("请输入a的值："))
b = float(input("请输入b的值："))
c = float(input("请输入c的值："))

delta = b**2 - 4*a*c  # 计算判别式

if delta < 0:
    real_part = -b / (2*a)  # 实部
    imag_part = math.sqrt(abs(delta)) / (2*a)  # 虚部
    print("方程无实数根！其虚根为X={:.2f}+{:.2f}i或{:.2f}-{:.2f}i".format(real_part, imag_part, real_part, imag_part))
elif delta == 0:
    x = -b / (2*a)
    print("方程有一个实根：X1=X2={:.2f}".format(x))
else:
    x1 = (-b + math.sqrt(delta)) / (2*a)
    x2 = (-b - math.sqrt(delta)) / (2*a)
    print("方程有两个实根：X1={:.2f} X2={:.2f}".format(x1, x2))

你可以输入示例中的 `a,b,c` 值进行测试，程序会输出相应的解。