Use FFT to Calculate Polymer Properties
3 次查看(过去 30 天)
显示 更早的评论
I have a polymer characterized by a transmission terahertz time-domain spectroscopy (THz-TDS) in the frequency range extending from 0 to 3.0 THz. And I have measured time-domain data in a .csv file (the attachment “waveform.csv”), with column A (Delay, unit: ps), column B (Reference Amplitude, unit: nA) and column C (Sample Amplitube, unit: nA).
I tried to use Matlab to calculate polymer properties (the attachment “Calculate_Substrate_Properties.m”), according to Part II.2 in a published article (Y.-S. Jin, G.-J. Kim, and S.-G. Jeon, “Terahertz dielectric properties of polymers,” J. Korean Phys. Soc., vol. 49, no. 2, pp. 513–517, 2006) (the attachment “Terahertz Dielectric Properties of Polymers.pdf”).
The polymer thickness is 0.5mm. Calculated results are not as same as expectations in the attachment (the attachment “Results.pdf”). I think it is due to wrong phase calculations. Could anyone help and tell me what problem there is in my program?
2 个评论
Rik
2020-11-18
Before you attempt to edit away your question: I made a backup here, so anyone will be able to restore your edit.
What do you think makes your question not appropriate?
Star Strider
2020-11-18
Rik — Thank you!
It turns out that I also saved the ‘Ryan waveform.CSV’ file and am attaching it to this Comment.
回答(2 个)
Star Strider
2016-9-1
See the R2015a documentation for fft for some guidance. It discusses appropriate ways to define the amplitude and frequency of the Fourier-transformed data.
This aspect of physical chemistry is not an area of my expertise, so I’m not going to attempt to understand what you’re doing, even after making an effort to understand equations (1)-(6). However some of the magnitude differences between your data and your expected data could be explained by your not normalising your Fourier transformed data by the length of the data.
See if this helps:
data_len = size(waveform.data,1); % Row Length Of Data Matrix
reference_spectra = fft(waveform.data(:, 2))/data_len;
Do the same normalisation by ‘data_len’ for your other fft calls.
Beyond that, I’m lost. I’ll leave you to the physical chemistry. (I was an undergraduate Chemistry major, but that was back in the phlogiston era.)
5 个评论
Star Strider
2016-9-5
My pleasure.
The data extraction step and code is to help anyone who has the necessary background to understand what you’re doing without the extra explanation that I require.
另请参阅
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!