傅里叶分析和滤波

目前问题：

由状态平均方程转传递函数时候，syms带符号的不能参与ss2tf的运算，因为没有实值。
所以很困难。

傅里叶分析和滤波

[b,a] = ss2tf(A,B,C,D) 将方程组的状态空间表示形式转换为等同的传递函数。ss2tf 返回连续时间方程组的拉普拉斯变换传递函数和离散时间方程组的 Z 变换传递函数。

变换和滤波器是用于处理和分析离散数据的工具，常用在信号处理应用和计算数学中。当数据表示为时间或空间的函数时，傅里叶变换会将数据分解为频率分量。fft 函数使用快速傅里叶变换算法，相对于其他直接实现，这种方式能够减少计算成本。有关傅里叶分析的更多详细介绍，请参阅傅里叶变换。在使用传递函数修改输入数据的幅值或相位时，conv 和 filter 函数也是很有用的工具。

[b,a] = ss2tf(A,B,C,D,ni) 返回当具有多个输入的方程组的第 ni 个输入受单位脉冲影响时所生成的传递函数。

Fs = 5;               %频率5Hz
dt = 1/Fs;            %采样间隔
N = 50;                %生成 50 个时间样本
t = dt*(0:N-1);
%% 振荡器可以通过状态空间方程描述。
A = [cos(dt) sin(dt);-sin(dt) cos(dt)];
B = [1-cos(dt);sin(dt)];
C = [-1 0];
D = 1;
%% 系统使用正方向的单位脉冲进行刺激。使用该状态空间模型计算系统从全零的初始状态开始的时间演进。
u = [1 zeros(1,N-1)]    

x = [0;