
上面是我的题目,我按照filtic函数的使用方法进行了编程
以下是我的代码
using TyPlot # 用于绘图
using TySignalProcessing
include("D:/DSP_Lab_projects/myfuns.jl")
b = [0.55,0.5,-1];#y[n]前系数
a = [1,-0.5,-0.45];#x[n]前系数
x0 = 3;#过去输入序列,注意顺序最近的在前面
y0 = 1;#过去输出序列,注意顺序最近的在前面
n = 0:20; # 时域索引(n从0到20,覆盖足够多的输出点)
x = 0.7.^n;
xic = filtic(b,a,y0,x0);
y, = filter1(b,a,x,xic);
但是还是报了一下的错误
ERROR: BoundsError: attempt to access 2-element Vector{Float64} at index [1, 2]
[2] filtic(b::Vector{Float64}, a::Vector{Float64}, ypast::Int64, xpast::Int64)
[2] filtic(b::Vector{Float64}, a::Vector{Float64}, ypast::Int64, xpast::Int64)
@ TySignalProcessing.Internal C:\Users\Public\TongYuan.julia\packages\TySignalProcessing\0UfKE\src\DigitalAndAnalogFilters\DigitalFiltering\filtic.jl:55
[3] top-level scope
@ d:\DSP_Lab_projects\实验4 LTI系统时域与Z域分析\l1 求系统的输出序列.jl:15

