function [D] = actualdatphi(A,f,phi,sigma,time)
save('sigma.txt','sigma','-ascii')
save('time.txt','time','-ascii')
M = modtimedatphi(A,f,phi,time)
randn('state',0)
D = M + sigma*randn(time,1)
D;
save('Ddata.txt', 'D', '-ascii')