這裏有兩個不同的東西 - 你怎麼發現的急劇轉變，你怎麼過濾。

讓我們依次採取這些。

尖銳過渡的特點是大麴率 - 我們可以通過取輸入曲線的diff來輕鬆檢測到這一點。使用第二個參數= 2需要兩次diff，並導致「類似」二階導數 - 但它被偏移一。所以當我們發現diff(sharkfin,2)很大的點時，我們需要抵消1來得到角點。

接下來，平滑本身。有很多技巧 - 我用一個盒子函數展示一個簡單的卷積。這樣做兩次給出了一個平滑的輸入版本。通過選擇「遠離不連續性」的原始文本和「接近不連續性」的過濾版本，我們得到了你所要求的。如果你願意，你可以「混合」點 - 使用加權版本的過濾和未過濾，取決於你與角點的距離。我沒有表現出明確，但它應該很容易看到如何擴大我已經寫的代碼：

% generate a "shark fin" function:

fin = exp(-linspace(0,4,60));

shark = [fin (1-fin+fin(end))];

shark = repmat(shark, [1 3]);

D2 = diff(shark, 2);

roundMe = find(abs(D2)>0.1*max(D2))+1; % offset by 1 because second derivative

figure;

subplot(3,1,1)

plot(shark); title 'shark plot'

hold on;

plot(roundMe, shark(roundMe),'r*')

legend('input','corners found')

% take N points on either side of the sharp corners:

N = 3;

% boxplot filtered version of the curve

boxFilt = ones(1, 2*N+1)/(2*N+1);

smoothShark1 = convn(shark, boxFilt, 'same'); % box plot

% second filter - smoother

smoothShark2 = convn(smoothShark1, boxFilt, 'same');

% plot the filtered results:

subplot(3,1,2)

plot(shark)

hold on

plot(smoothShark1);

hold on

plot(smoothShark2);

xlim([114 126])

ylim([0.8,1.1])

legend('original','box','box x2')

title 'smoothed everywhere'

% Now apply filtering only to points near the discontinuity

smoothMe = zeros(size(shark));

smoothMe(roundMe)=1;

smoothMe = convn(smoothMe, boxFilt, 'same');

smoothMe(smoothMe>0)=1; % this finds N points on either side of the corner

subplot(3,1,3)

plot(shark)

finalPlot=shark;

hold on

smoothIndx = find(smoothMe);

finalPlot(smoothIndx)=smoothShark2(smoothIndx);

plot(finalPlot,'g')

plot(smoothIndx, finalPlot(smoothIndx), 'r*')

xlim([114 126])

ylim([0.8,1.1])

legend('original','smoothed','changed')

title 'smoothed only near discontinuity'

輸出：