%Ai es el valor de amplitud de cada armonico
%phii es el valor de defasaje en radianes de dicha componente
A1=0;
phi1=2*pi*0/360;
A3=0;
phi3=2*pi*0/360;
A5=0;
phi5=2*pi*0/360;
A7=0;
phi7=2*pi*0/360;
A9=0;
phi9=2*pi*0/360;
A11=0;
phi11=2*pi*0/360;
A13=0;
phi13=2*pi*0/360;
A15=0;
phi15=2*pi*0/360;
% vector tiempo con pasos de 0.2 ms
t=0:0.0002:0.1;
v=sqrt(A1^2+A3^2+A5^2+A7^2+A9^2+A11^2+A13^2+A15^2);
% creo una funcion sin ruido
fsr=A1*sin(2*pi*50*t+phi1)+A3*sin(2*pi*3*50*t+phi3)+A5*sin(2*pi*5*50*t+phi5)+A7*sin(2*pi*7*50*t+phi7)+A9*sin(2*pi*9*50*t+phi9)+A11*sin(2*pi*11*50*t+phi11)+A13*sin(2*pi*13*50*t+phi13)+A15*sin(2*pi*15*50*t+phi15);
% agrego ruido a la se;al
for i=1:501
f(i)=fsr(i)+0.1*v*randn();
end