频谱在MATLAB，终点proble

Question

我试图采取用以下行125个频段的频谱的导数微分：

dW=diff(wavelength); 
dR=diff(data); 
df=dR./dW; 

问题是下一步我想将它与比较原始光谱数字和视觉上通过绘图，但df的大小是124，但我的原始波长是125.问题是我必须删除第一个或最后一个乐队？然而一些光谱分析软件的输出并没有改变尺寸。考虑到乐队的平均水平也不起作用，它会使图表显示疯狂的行为。

Answer 1

diffis basically：

Y = [X(2)-X(1) X(3)-X(2) ... X(m)-X(m-1)] 

这意味着它比你输入一个短（你不能从一无所有减去一些东西，对吧？）。

你当然需要做什么取决于你想要做什么，但最少的“意义改变”方法（保持与采样时间有关的因果关系的类型）应该是在你的dW和dR之前添加一个任意值。

顺便说一下，如果dW有零（只要两个连续的数据值相同，就会发生），那么您的比率df=dR./dW可能有很多NaN s。