问题遇到的现象和发生背景
对于如图(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读取结果是怎样的一个形式?