TU Wien:Computernumerik VU (Schranz-Kirlinger)/Übungen SS16/Beispiel 13
Zur Navigation springen
Zur Suche springen
Lösung[Bearbeiten | Quelltext bearbeiten]
Matlab Code
function Bsp3
data = csvread('datenDO.csv');
xj = data(:,1);
fj = data(:,2);
N = size(data,1);
A = [N+1 sum(xj) sum(xj.^2) ;
sum(xj) sum(xj.^2) sum(xj.^3) ;
sum(xj.^2) sum(xj.^3) sum(xj.^4)];
b = [sum(fj) sum(fj .* xj) sum(fj .* (xj.^2))]';
res = A \ b;
fprintf('%s %f %s %f %s %f %s\n', 'Polynom =', res(1), , res(2), 'x +', res(3), 'x^2');
end