MATLAB 信号分析的复频域处理

问题遇到的现象和发生背景

对于如图（1）所示的一个系统函数，其复频域表示为

其中a=1/20; w0=800*pi;
输入一段音频X，求这段音频经过此系统的两个相同级联后的响应？

我的解答思路和尝试过的方法

尝试（1）

a=1/20;w0=800*pi;
num=[0 2*a*w0 0];den=[w0^2 2*a*w0 1]; 
sys=tf(num,den); %求系统复频域
[x,Fs]=audioread('music.wav');
xs=laplace(x); %求输入音频的复频域
y=ilaplace(xs*sys^2);

报错：检查对函数 'laplace' 的调用中是否存在不正确的参数数据类型或缺少参数。

尝试（2）

a=1/20;w0=800*pi;
num=[0 2*a*w0 0];den=[w0^2 2*a*w0 1]; 
sys=tf(num,den); %求系统复频域
[x,Fs]=audioread('music.wav');
y=lsim(sys^2,x,1./Fs); %在采样率为Fs下的系统响应；

报错：错误使用 DynamicSystem/lsim (第 97 行)
When simulating the response to a specific input signal, the time vector must be specified for continuous-time
models.

我的问题：

（1）我的两次程序错误出现在哪里？
（2）如何正确实现该程序？
（3）lism是如何正确使用的？
（4）audioread读取结果是怎样的一个形式？

你好，我是有问必答小助手，非常抱歉，本次您提出的有问必答问题，技术专家团超时未为您做出解答


本次提问扣除的有问必答次数，已经为您补发到账户，我们后续会持续优化，扩大我们的服务范围，为您带来更好地服务。