我最近在用matlab編寫霍布金森桿數據處理軟件，導入數據、自動對齊波形、截波、手動對齊波形的功能都實現了，但在在應該最簡單的應力應變曲線計算這步出了問題，找了好久都不知道問題出在哪。希望各位前輩有知道的幫解答一下
下面是具體情況
程序運行界面如圖

       根據現實的三波分離的對齊的效果還可以，也有前輩根據一樣的數據一樣的對齊波形的情況算出來應力應變曲線。
       在matlab中分步運行查看工作區變量的內容后，計算出來的應力值比較理想，沒有什么不合適的，但是應變率就比較奇怪。還有就是根據某一視頻中方法，在圖中情況下，紅色的透射波在最后一段比入射波還大，這一段應該是不符合一維應力波假設的，但是視頻中也算出了可靠性比較高的應力應變曲線。

下邊是二波法計算應力應變的公式

下邊是我計算應力應變的代碼

        變量jieguo為4列n行（具體行數不定，由采樣頻率影響）的矩陣，第一列存儲應力信息，第二列存儲應變信息，第三列存儲應變率信息，第四列存儲時間信息（為三波分離之后圖像橫坐標的值）。new_time_stain變量中存儲的是三波分離對齊后的數據，1、2列是入射波時間和入射波數據；3、4列為反射波；5、6列為透射波其他的固定數據保證是的對應公式里的數據，沒有問題。

        主要是求應變的公式是積分計算，不知道我所理解中的積分上限時間t是否正確。請了解的前輩為我看看哪里出問題了。