本代码基于MATLAB实现了在分析黑体辐射的光谱图时将通过迈克尔逊干涉仪得到的干涉图经过傅里叶变换转换为光谱图的功能.zip

# MATLAB-DFT 本代码基于MATLAB实现了在分析黑体辐射的光谱图时将通过迈克尔逊干涉仪得到的干涉图转换为光谱图的功能 This code is based on MATLAB which implements the function of the transformation from interferogram to spectrogram during the analysis of radiation of blank-body. 众所周知，迈克尔逊干涉仪可以将一束光分为两部分进行干涉产生干涉图样，通过移动动镜来改变两束光之间的光程差，通过获得不同光程差时的光强绘制出干涉图。通过在计算机中将干涉图（光强随位置即时间的变化图像）进行傅里叶变换转化为光谱图（光强随频率即波长的变化图像）然而，传统的傅里叶变换要求积分区间为无穷，即要求光程差从负无穷变为无穷，同时要求连续采样，即采样的间隔要求为无穷小。这显然是在计算机上无法实现的。但我们知道，对于干涉，在光程差较小时有较大光强，光程差较大时光程差可以忽略不计；同时，当采样间隔足够小时，连续积分近似于离散求和。所以我们退而求其次，选择有限区间，同时间隔采样，将积分变为求和，得到近似于连续傅里叶变换的离散结果，使得这种操作在计算机上可以进行。对于庞大数据量，这种方式显示出了其快速、高效的特点。 在code.m部分展示了在MATLAB上的代码，其中第一部分为预处理，第二部分为傅里叶变换，第三部分为图像显示。开头的注释为当前的一些问题，结尾的注释为添加高斯噪声后的结果。 方波.CSV和参考.CSV中为表格形式的数据，第一列为采样时间（仅供参考），第二列为光强（相对值） 方波.fig和参考.fig中为代码运行后的结果图像，第一个为原信号干涉图，第二个为光谱图。 注：本代码原型为本人学长的一次大物实验作业，我在参考此代码基础上进行修改，将其应用到这个问题上。