TU Wien:Computernumerik VU (Schranz-Kirlinger)/Übungen SS16/Beispiel 13

Aus VoWi
Zur Navigation springen Zur Suche springen

Lösung[Bearbeiten | Quelltext bearbeiten]

datenDO.csv

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

PDF Ausarbeitung