尝试movmin信封()。然后试着Saivtky-Golay过滤器如果你想顺利一些。
%初始化步骤:
clc;%清除命令窗口。
关闭所有;%关闭所有数据(除了imtool。)
clearvars;
工作空间;%确定工作区面板显示。
格式长g;
格式紧凑的;
字形大小= 15;
s =负载(“SI_and_lambda.mat”)
如果= s.SI;
λ= s.lambda;
minval = 0;
maxval = 950;
次要情节(1,1);
情节(λ,是的,“- - -”)
包含(“波长”,“字形大小”字形大小)
ylabel (的信号强度,“字形大小”字形大小)
ylim ([minval maxval])
标题(原始信号的,“字形大小”字形大小)
%全屏放大图。
集(gcf,“单位”,“归一化”,“Outerposition”,(0、0.05、1、0.95));
%过滤器和一个移动的最小窗口。
windowWidth = 27
windowWidth filteredSignal = movmin (SI);
次要情节(3、1、2);
情节(λ,filteredSignal“- - -”,“线宽”,2)
包含(“波长”,“字形大小”字形大小)
ylabel (的信号强度,“字形大小”字形大小)
ylim ([minval maxval])
标题(过滤信号移动最小的,“字形大小”字形大小)
网格在;
% Savitzky-Golay过滤器过滤。
windowWidthSG = 55
filteredSignal = sgolayfilt (filteredSignal 2 windowWidthSG);
次要情节(3,1,3);
情节(λ,filteredSignal“- - -”,“线宽”,2)
包含(“波长”,“字形大小”字形大小)
ylabel (的信号强度,“字形大小”字形大小)
ylim ([minval maxval])
标题(与移动二次Savitzky-Golay过滤器过滤信号的,“字形大小”字形大小)
网格在;
%全屏放大图。
集(gcf,“单位”,“归一化”,“Outerposition”,(0、0.05、1、0.95));
