clear,clc
% DATOS TOMADOS
f =[600,700,800,900,1000,1100,1200,1300,1400,1500,1600,1700]; %Valores de las frecuencias
finv=1./f; %Hallar los inversos de las frecuencias
londa=[0.577,0.493,0.432,0.382,0.3455,0.315,0.287,0.265,0.245,0.231,0.217,0.204]; %Valores de las longitudes de onda
%Gráfica de la recta
[a,b]=polyfit(finv,londa,1); %Hallar mediante una función de MATLAB la pendiente, la ordenada en el origen y el error de la mejor recta
finv2=[finv(length(finv)),finv(1)]; %Muestreo de la recta
londa2=a(1).*finv2+a(2); %Función de la recta
plot(finv,londa,'.r',finv2,londa2,'g'); %Pintar los datos cen discreto y la recta en continuo
axis square %Ejes cuadrados
title('\bfAjuste a la recta','FontSize', 24);
xlabel('\bfT [s]' ,'FontSize', 24);
ylabel('\bf\lambda [m]','FontSize', 24);
legend('\bfDatos muestreados','\bfAjuste a la recta');
gtext('\bf\lambda=1/f · v ','FontSize', 24)
%%
%Gráfica de la hipérbola
clear,clc
f =[600,700,800,900,1000,1100,1200,1300,1400,1500,1600,1700]; %Valores de las frecuencias
finv=1./f; %Hallar los inversos de las frecuencias
londa=[0.577,0.493,0.432,0.382,0.3455,0.315,0.287,0.265,0.245,0.231,0.217,0.204]; %Valores de las longitudes de onda
%Bucle para hallar el valor de v que cometa menos error en la hipérbola
v=0;
for k=1:length(f)
v=(v+f(k)*londa(k));
end
v=v/length(f);
f2=f(1):0.1:f(length(f));
londa3=v./f2;
plot(f,londa,'r.',f2,londa3)
axis square
title('\bfHipérbola','FontSize', 24);
xlabel('\bff [Hz]','FontSize', 24);
ylabel('\bf\lambda [m]','FontSize', 24);
legend('\bfDatos muestreados','\bfHipérbola');
gtext('\bf\lambda=v/f','FontSize', 24)
%%
%Hallar la velocidad del sonido
vsonido=(a(1)+v)/2;
fprintf('La velocidad del sonido hallada por la recta es ')
fprintf(num2str(vsonido))
fprintf('\n')
fprintf('La velocidad del sonido hallada mediante la hipérbola es ')
fprintf(num2str(v))
fprintf('\n')
