clear
clc

%打開數據文件----------------------
[k,s]=fopen('12.txt');
a=fscanf(k,'%g');
a1=a;
fclose all;


%去奇點（異常的噪點）------------------------
b=size(a);
for i=2:b-1;
 if a(i)>=a(i+1)+10&a(i)>=a(i-1)+10;
     
    a(i)=(a(i-1)+a(i+1))./2;
     
 end
 if a(i)+8<=a(i+1)&a(i)+8<=a(i-1);
     
     a(i)=(a(i-1)+a(i+1))./2;
 end
end
 a2=a;
 
 %平滑（去毛刺）----------------------------
b=size(a);
for l=1:4;
for i=2:b-1;
if a(i)>=a(i+1)&a(i)>=a(i-1);
     c=(a(i)+a(i-1)+a(i+1))./3;
     a(i)=c;
     a(i+1)=c;
     a(i-1)=c;
 end
 if a(i)<=a(i+1)&a(i)<=a(i-1);
     c=(a(i)+a(i-1)+a(i+1))./3;
     a(i)=c;
     a(i+1)=c;
     a(i-1)=c;
 end
 
 
 
 
end
end

%畫出處理前后的對比曲線圖-----------------------
subplot(311);
plot(a1);
subplot(312);
plot(a2);
subplot(313);
plot(a);