Diferenciación Numérica

Table of Contents
Tarea 1 - Diferenciación ...................................................................................................... 1
Function Handle................................................................................................................. 1
Problema de la Manivela y el Pistón...................................................................................... 1

Tarea 1 - Diferenciación
% Salvador Tamayo Atias
% Cálculo Numérico

Function Handle
% Para encontrar la aceleración de la manivela del pistón, la distancia x,
% se utiliza el métodode diferenciación, el que permite encontrar de forma
% numérica dicha variable, con la ayuda de la función "difapx".
% *************************************************************************
%difapx.m
%
%
%
%
%
%
%
%
%
%
%
%
%
%
%

function [c,err,eoh,A,b] = difapx(N,points)
l = max(points);
L = abs(points(1)-points(2))+ 1;
if L < N + 1, error('More points are needed!');end
for n = 1: L
A(1,n) = 1;
for m = 2:L + 2, A(m,n) = A(m - 1,n)*l/(m - 1); end
l = l-1;
end
b = zeros(L,1); b(N + 1) = 1;
c =(A(1:L,:)\b)';
err = A(L + 1,:)*c'; eoh = L-N;
format rat
ifabs(err) < eps, err = A(L + 2,:)*c'; eoh = L - N + 1; end
if points(1) < points(2), c = fliplr(c); end

%
%
%
%
%
%

*************************************************************************Éste método encuentra los coeficientes correspondientes a la expansión en
serie de Taylor de la fórmula de derivadas. También señala el coeficiente
y el órden del error de truncamiento que secomete al utilizar éste
método.

Problema de la Manivela y el Pistón
clear all;close all;clc;
format bank;
% Datos del problema

1

R = 90e-3; % Radio de giro.
theta = [0:5:180]; % Ángulos aanalizar.
PI = pi/180; % Constante de transformación grados - radianes.
freq = 5000/60; % velocidad angular - 5000 rev/min a 83.3 rev/seg
VEL_ang = freq*360;
h = 0.0005; % paso.
% Este paso...