用fft编程计算材料性质
 

我用一台太赫兹传输时域光谱仪（THz-TDS）测试一种材料（频率从0到3 THz），得到材料的时域测试数据（附件“waveform.csv”）。其中，A栏是时间（单位：ps）；B栏是没材料时的幅度（单位：nA）；C栏是有材料时的幅度（单位：nA）。

根据一篇学术论文（附件“Terahertz Dielectric Properties of Polymers.pdf”）Part II.2里的公式，我用Matlab编程计算材料性质（附件“Calculate_Substrate_Properties.m”）。

材料的厚度为0.5mm。Matlab计算值与理论值不同（附件“Results.pdf”）。我觉得我的程序中相位计算有误。请问有谁能帮我找出程序哪里有错？