我最近在用matlab编写霍布金森杆数据处理软件，导入数据、自动对齐波形、截波、手动对齐波形的功能都实现了，但在在应该最简单的应力应变曲线计算这步出了问题，找了好久都不知道问题出在哪。希望各位前辈有知道的帮解答一下
下面是具体情况
程序运行界面如图

       根据现实的三波分离的对齐的效果还可以，也有前辈根据一样的数据一样的对齐波形的情况算出来应力应变曲线。
       在matlab中分步运行查看工作区变量的内容后，计算出来的应力值比较理想，没有什么不合适的，但是应变率就比较奇怪。还有就是根据某一视频中方法，在图中情况下，红色的透射波在最后一段比入射波还大，这一段应该是不符合一维应力波假设的，但是视频中也算出了可靠性比较高的应力应变曲线。

下边是二波法计算应力应变的公式

下边是我计算应力应变的代码

        变量jieguo为4列n行（具体行数不定，由采样频率影响）的矩阵，第一列存储应力信息，第二列存储应变信息，第三列存储应变率信息，第四列存储时间信息（为三波分离之后图像横坐标的值）。new_time_stain变量中存储的是三波分离对齐后的数据，1、2列是入射波时间和入射波数据；3、4列为反射波；5、6列为透射波其他的固定数据保证是的对应公式里的数据，没有问题。

        主要是求应变的公式是积分计算，不知道我所理解中的积分上限时间t是否正确。请了解的前辈为我看看哪里出问题了。