Comprimento do arco da curva (MATLAB)

 


% Definir a função a ser integrada

fun = @(x) sqrt(1 + (9*x^2)/(4*x^(3/2)));


% Definir os limites de integração

a = 1;

b = 4;


% Número de pontos de amostra

n = 100;


% Calcular o espaçamento do passo

h = (b-a)/(n-1);


% Calcular o valor da integral usando o método de Simpson

L = fun(a) + fun(b);


for i = 2:2:n-1

    L = L + 4*fun(a+i*h);

end


for i = 3:2:n-2

    L = L + 2*fun(a+i*h);

end


L = L*h/3;


% Imprimir o resultado

fprintf("O comprimento do arco é: %.4f unidades de comprimento.\n", L);


SOLUÇÃO


O comprimento do arco é: 6.2927 unidades de comprimento.

(Exemplo 02) Representação para o comprimento do arco em Matlab.

ALGORITMO


% Passo 1: Definir a função que representa a curva

f = @(x) [sqrt(x.^3); -sqrt(x.^3)];


% Passo 2: Plotar a curva no intervalo desejado

x = linspace(1, 4, 1000);

y = f(x);

plot(x, y, 'LineWidth', 2);

hold on;

plot([1, 4], [1, 8], 'ro', 'MarkerSize', 10, 'LineWidth', 2);

hold off;

grid on;

xlabel('x');

ylabel('y');

title('Curva y^2 = x^3 entre os pontos (1,1) e (4,8)');


% Passo 3: Calcular o comprimento do arco

a = 1;

b = 4;

dx = 0.001;

x = a:dx:b;

y = f(x);

dy_dx = gradient(y, dx);

L = trapz(sqrt(1 + dy_dx.^2), dx);


% Passo 4: Plotar o comprimento do arco na curva

t = linspace(1, 4, 100);

xt = interp1(x, t, 'linear');

yt = interp1(y, t, 'linear');

Lx = xt(1:end-1);

Ly = yt(1:end-1);

Larco = [0 cumsum(sqrt(1 + diff(yt).^2))];

plot(x, y, 'LineWidth', 2);

hold on;

plot(Lx, Ly, 'r', 'LineWidth', 2);

hold off;

grid on;

xlabel('x');

ylabel('y');

title(sprintf('Curva y^2 = x^3 entre os pontos (1,1) e (4,8)\nComprimento do arco: %.4f', L));



REPRESENTAÇÃO GRÁFICA






Postar um comentário

Postagem Anterior Próxima Postagem