Skip to content
Permalink
Branch: master
Find file Copy path
Find file Copy path
Fetching contributors…
Cannot retrieve contributors at this time
33 lines (28 sloc) 827 Bytes
function [] = rueckwaertsdifferenz
% Funktion Vorwärtsdifferenz
function [dy] = D3f (func, x0, h)
dy = (func(x0) - func(x0 - h)) / h;
end
% Beispiel aus Serie 1 Aufgabe 2
% Geg: Funktion, Ableitung, x0
fn = @(x) log(x.^2);
df = @(x) (1/(x.^2))*2.*x;
x0 = 2;
% Ges: Rückwärtsdifferenz und Diskretisierungsfehler für h-Werte 10^-1
% bis 10^-17
res_d3f = [];
res_err_d3f = [];
eff = df(x0);
for i = 1:17
% Näherungswert berechnen für h = 10^-i
[dy] = D3f(fn, x0, 10^(-i));
% Diskretisierungsfehler berechnen
err = abs(dy - eff);
res_d3f(i) = dy;
res_err_d3f(i) = err;
end
fprintf("RESULTAT FÜR D3f\n");
res_d1f'
fprintf("DISKRETISIERUNGSFEHLER FÜR D3f\n");
res_err_d3f'
end
You can’t perform that action at this time.