Ahoj, mám zadané hodnoty x a y a mám aproximovat fci f(x) = c1x^M + c2x^M−1 + ... + cM+1, M=1 pomocí metody nejmenších čtverců a určit součet čtverců odchylek. Něco jsem se snažila vypracovat, ale to znázornění 3D a vykreslení vrstevnice asi není dobře, mohl by mi někdo poradit, co s tím?
clear all;close all;clc;
%Aproximace polynomem
x=[0.1 0.3 0.4 0.6 0.9]';
y=[1.1 2.9 2.9 4.3 3.5]';
c=polyfit(x,y,1);
xx=[0:0.05:1.5]';f=polyval(c,xx);
figure(1);
plot(xx,f)
hold on;
%stem(x,y,'r');
hold off;
a1=2:0.1:4;
a2=2.5:0.1:4.5;
for i=1:21;
for j=1:21
S(i,j)=sum((a1(j)*x+a2(i)-y).^2)
end
end
subplot(2,2,2);plot(xx,f);hold on;stem(x,y)
hold off;
subplot(1,2,1);meshc(a2,a1,S);grid on;
subplot(2,2,4); contour(a2,a1,S);