Interpolacion

Comandos del programa hecho para el método de Interpolación
%Métodos Numéricos
%Curva de Histéresis
%Cálculo del Magnetismo Remanente y Campo Coercitivo de un Material Ferromagnético
%DatosObtenidos mediante estudio
fprintf('B: 1 0.95 0.9 0.8 0.7 0.6 0.4 0.2 0 -0.2 -0.4 -0.6 -0.7 -0.8 -0.9 -1 -0.95 -0.9 -0.8 -0.7 -0.6 -0.4 -0.2 0.01 0.2 0.4 0.6 0.7 0.8 0.9 1')
fprintf('H: 193 80 42 2 -18-29 -40 -45 -48 -52 -58 -73 -85 -103 -135 -193 -80 -42 -2 18 29 40 45 48 52 58 73 85 103 135 193')
T= [193 80 42 2 -18 -29 -40 -45 -48 -52 -58 -73 -85 -103 -135 -193 -80 -42 -2 18 29 40 45 48 52 5873 85 103 135 193]

R= [1 0.95 0.9 0.8 0.7 0.6 0.4 0.2 0 -0.2 -0.4 -0.6 -0.7 -0.8 -0.9 -1 -0.95 -0.9 -0.8 -0.7 -0.6 -0.4 -0.2 0.01 0.2 0.4 0.6 0.7 0.8 0.9 1]

hold on
plot(T,R)
%En el análisisse toma la parte superior de la Curva, por lo tanto también de los vectores, obteniendo un nuevo par de vectores.
h= [193 80 42 2 -18 -29 -40 -45 -48 -52 -58 -73 -85 -103 -135 -193]
b= [1 0.95 0.90.8 0.7 0.6 0.4 0.2 0.01 -0.2 -0.4 -0.6 -0.7 -0.8 -0.9 -1]

plot(h,b)
grid

%Interpolación Lineal de los Vectores
n=1

%Interpolación donde “H” cambia de signo

while n<16
ifh(n)*h(n+1)<0
%M: Coordenada en H
M=0
%C: Coordenada en B interpolada
C=(((M-h(n))/(h(n+1) - (h(n))))*(b(n+1) - b(n))) + b(n)
%V: Coordenadas en H de la recta deinterpolación
V=linspace(h(n-1),h(n+1),5)
%K: Coordenadas en B de la recta de interpolación
K=linspace(b(n-1),b(n+1),5)
plot(V,K, 'R')
grid
plot(M,C, 'Y')disp(['El magnetismo remanente es:',num2str(C)]);
break
end

n=n+1
end

%Interpolación donde “B” cambia de signo

while n<16
if b(n)*b(n+1)<0%F: Coordenada en B
F=0
%D: Coordenada en H interpolada
D=(((F-b(n))/(b(n+1) - (b(n))))*(h(n+1) - h(n))) + h(n)
%L: Coordenada en H de la recta de interpolación...