% Sample data
clear
close all
Fs=1000;
t=0:1/Fs:10;
Fsig=5
A = 500;
m=50;
d=200;

y=A*sin(2*pi*Fsig*t) + m*t + d;



% High-pass filter design
Fc   = Fsig/2;   % Hz cuttoff frequency
N    = 2;    % filter order

[b,a] = butter(N,Fc/Fs,'high');
yf = filter(b,a,y);

figure('Name','Sample File Format Data')
plot(t,y,'b',t,yf,'r')
xlabel('Time(s)');
ylabel('Encoder Counts');
legend('Raw', 'Filtered')
fileID = fopen('50_10test.txt','w');
for i=1:length(y),
    fprintf(fileID,'%f,%f\n',y(i),yf(i));
end
fclose(fileID);