clc;clear;
tps=621; %输入图片数
for n=1:5
k=tps-5;
P=zeros(k,1);
if n==1
x=csvread('实际位移.csv');
else
x=csvread(name);
end
name=strcat('位移滤波_',sprintf('%03d',n),'.csv');
P(1,1)=(3*x(1,1)+2*x(2,1)+x(3,1))/6;
P(2,1)=(2*x(1,1)+3*x(2,1)+2*x(3,1)+x(4,1))/8;
P(k-1,1)=(2*x(k,1)+3*x(k-1,1)+2*x(k-2,1)+x(k-3,1))/8;
P(k,1)=(3*x(k,1)+2*x(k-1,1)+x(k-2,1))/6;
for m=3:k-2;
P(m,1)=(x(m-2,1)+2*x(m-1,1)+3*x(m,1)+2*x(m+1,1)+x(m+2,1))/9;
end
csvwrite(name,P);
end
csvwrite('n次滤波后实际位移.csv',P);
t=0:0.118:0.118*(tps-6);
plot(t,P)
title('011斥力盘位移曲线')
xlabel('t/ms')
ylabel('s/mm')