Budownictwo 2006 Wydział Inżynierii Lądowej Politechniki Krakowskiej

ewelka87

napisal ktos moze jakis algorytm??

krakoz

Dobry

function a=apr(A)
n = size(A,1);
x=A(:,1);
y=A(:,2);
for i=(1:n)
d(:,i)=x.^(i-1);
end
% for k=(1:n)
% y = A(k,1);
% end
a=d'*d\d'*y;
% a = d'*y*(d'^-1)*(d^-1)
printf("Wynik to: a=%i",a);
endfunction
function A=czyt(nazwa)
fid=fopen(nazwa,'r');
A=fscanf(fid,'%f',[2,Inf])';
endfunction
function X=mkplot(a)
n = size(a,2);
x = (0:0.1:20);
y = 0;
for k=(1:n)
y = y + a(k,1) * x.^(n-k-1);
end
% printf("Punkty wykresu:\n");
% printf("%f ",x);
% printf("\n");
% printf("%f ", y);
plot(x,y);
endfunction
A=czyt('dan12')
% a=apr(A)
% A = rand(2,2)
%A = [ 1,2,3,4,5,6,7,8,9; 10,1,10,1,10,1,10,1,10 ]
a = apr(A)
mkplot(a)

kosciel

tak to powinno końcowo wyglądać :
function A=czyt(nazwa)
fid=fopen(nazwa,'r');
A=fscanf(fid,'%f',[2,Inf])';
endfunction
function y=g(A,x)
y=0
for i=1:length(A)
y=y+A(i)*x.^(i-1)
end
end
function A=f(n,x,y)
for i=1:n+1
D(:,i)=x'.^(i-1);
end
A=D'*D\D'*y
end
A=czyt('dan3')
x=A(:,1)'
y=A(:,2)
#n - tutaj wpisz liczbę baz
n=3
A=f(n,x,y)
x1=x(1):0.1 Mad

(length(x))
for i=1:length(x1)
y1=g(A,x1)
end

plot(x,y,"@1",x1,y1,"-b;funkcja aproksymujaca;")

pause()

Maciej

a o ktorej w czwartek mozna sie udac do zioma? ;]

dyniak

z tego co wiem to tylko w środy ma konsultacje, ale słyszałem jak gadał z kościelem że chyba jutro rano koło 8.00-8.30 można do niego przyjść. zapytaj go bo ja nie jestem pewny tego co pisze ;p

saddam

Czwartek, 8.30. Ostateczny termin. Ja sie z nim tak ugadalem.

Maciej

ołki dołki