几种数字仿真的物理意义与代码实现

任何的系统都包括输入与输出，小到一片最小的代码片，大道航天工业系统，都是由最基本的输入输出与中间环节构成。如何处理中间的环节就显得异常重要。现如今的系统大多为数字仿真系统，也就是大多使用计算机来处理输入与反馈信号。也信号又分为连续信号与离散信号，今天我们来谈一谈利用计算机仿真来处理信号的几种典型的方法及实现的MATLAB代码。
（1）数值积分法：
我们知道，从数学的角度上考虑，y(k+1)=y(k)+∫f(t,y(t))dt这样得来的确实是最优解，如果一直这样计算下去，得到的永远都是最精确地值。但是我们的计算机却无法做到这样的计算，因为两时刻之间的斜率为变化值。于是我们利用简化的方式去代替积分，只要采样时刻取得足够小，我们就可以理所当然的人为t（k）时刻与t（k+1）时刻之间是一条直线，于是我们就认为两者之间的斜率就是不变的值。这样看来就是好办多了只要曲积分为t(k)时刻的斜率，就解决了。下面为实例实现代码：
已知：y’ = -y; y(0) = 1;

function y = euler(y0,h,f)
y = y0 + h*f;
%代码如下 ：
h =