Forum Budownictwo 2006 Wydział Inżynierii Lądowej Politechniki Krakowskiej Strona Główna Budownictwo 2006 Wydział Inżynierii Lądowej Politechniki Krakowskiej

 
 FAQFAQ   SzukajSzukaj   UżytkownicyUżytkownicy   GrupyGrupy   GalerieGalerie   RejestracjaRejestracja 
 ProfilProfil   Zaloguj się, by sprawdzić wiadomościZaloguj się, by sprawdzić wiadomości   ZalogujZaloguj 

Aproksymacja

 
Napisz nowy temat   Odpowiedz do tematu    Forum Budownictwo 2006 Wydział Inżynierii Lądowej Politechniki Krakowskiej Strona Główna -> Metody numeryczne
Zobacz poprzedni temat :: Zobacz następny temat  
Autor Wiadomość
ewelka87
Kierownik Budowy



Dołączył: 05 Lis 2006
Posty: 305
Przeczytał: 0 tematów

Pomógł: 3 razy
Ostrzeżeń: 0/5

PostWysłany: Nie 23:35, 11 Lis 2007    Temat postu: Aproksymacja

napisal ktos moze jakis algorytm??

Post został pochwalony 0 razy
Powrót do góry
Zobacz profil autora
Zobacz poprzedni temat :: Zobacz następny temat  
Autor Wiadomość
krakoz
Kierownik Budowy



Dołączył: 05 Paź 2006
Posty: 418
Przeczytał: 0 tematów

Pomógł: 10 razy
Ostrzeżeń: 0/5

PostWysłany: Śro 0:38, 21 Lis 2007    Temat postu:

x=[0.279814 0.921852 1.453676 1.892098 2.390504 2.844429 3.514704 3.765632 4.520695 5.098113]
y=[-0.088762; 0.386523; 1.350046; 2.333342; 3.627730; 5.220384; 8.725832; 10.462116; 17.756209; 27.297747]
%%%można wstawić funkcję która sama czyta z innego pliku X iY%%

n=3;
for i=1:n+1
D(:,i)=x'.^(i-1);
end
D
A=D'*D\D'*y

%wielomian ma wyglądac tak:
%z=A1+A2*x+A3*x^2....

z=0;
for i=1:n+1
z=z+A(i,Smile*x.^(i-1);
end

plot(x,z,x,y);
pause();



Dostałem za to pozytywną ocenę. Trzeba tam jedynie pododawać kolory i kszytałty wykresu w gnuplocie


Post został pochwalony 0 razy

Ostatnio zmieniony przez krakoz dnia Śro 14:13, 21 Lis 2007, w całości zmieniany 1 raz
Powrót do góry
Zobacz profil autora
Zobacz poprzedni temat :: Zobacz następny temat  
Autor Wiadomość
Dobry
Prezes



Dołączył: 16 Paź 2006
Posty: 1273
Przeczytał: 0 tematów

Pomógł: 13 razy
Ostrzeżeń: 0/5
Skąd: Festung Krakau

PostWysłany: Śro 0:40, 21 Lis 2007    Temat postu:

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)


Post został pochwalony 0 razy
Powrót do góry
Zobacz profil autora
Zobacz poprzedni temat :: Zobacz następny temat  
Autor Wiadomość
kosciel
Kierownik Budowy



Dołączył: 11 Gru 2006
Posty: 352
Przeczytał: 0 tematów

Pomógł: 18 razy
Ostrzeżeń: 0/5
Skąd: Straconka

PostWysłany: Śro 14:21, 21 Lis 2007    Temat postu:

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.1Mad(length(x))
for i=1:length(x1)
y1=g(A,x1)
end

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

pause()


Post został pochwalony 1 raz
Powrót do góry
Zobacz profil autora
Zobacz poprzedni temat :: Zobacz następny temat  
Autor Wiadomość
Maciej
Kierownik Budowy



Dołączył: 20 Paź 2006
Posty: 485
Przeczytał: 0 tematów

Pomógł: 6 razy
Ostrzeżeń: 1/5
Skąd: Grupa 4

PostWysłany: Śro 18:08, 21 Lis 2007    Temat postu:

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

Post został pochwalony 0 razy
Powrót do góry
Zobacz profil autora
Zobacz poprzedni temat :: Zobacz następny temat  
Autor Wiadomość
dyniak
Kierownik Budowy



Dołączył: 04 Paź 2007
Posty: 407
Przeczytał: 0 tematów

Pomógł: 11 razy
Ostrzeżeń: 0/5

PostWysłany: Śro 20:35, 21 Lis 2007    Temat postu:

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

Post został pochwalony 0 razy
Powrót do góry
Zobacz profil autora
Zobacz poprzedni temat :: Zobacz następny temat  
Autor Wiadomość
saddam
Kierownik Budowy



Dołączył: 26 Lut 2007
Posty: 341
Przeczytał: 0 tematów

Pomógł: 11 razy
Ostrzeżeń: 0/5
Skąd: że znowu

PostWysłany: Śro 23:09, 21 Lis 2007    Temat postu:

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

Post został pochwalony 0 razy
Powrót do góry
Zobacz profil autora
Zobacz poprzedni temat :: Zobacz następny temat  
Autor Wiadomość
Maciej
Kierownik Budowy



Dołączył: 20 Paź 2006
Posty: 485
Przeczytał: 0 tematów

Pomógł: 6 razy
Ostrzeżeń: 1/5
Skąd: Grupa 4

PostWysłany: Śro 23:17, 21 Lis 2007    Temat postu:

ołki dołki

Post został pochwalony 0 razy
Powrót do góry
Zobacz profil autora
Wyświetl posty z ostatnich:   
Napisz nowy temat   Odpowiedz do tematu    Forum Budownictwo 2006 Wydział Inżynierii Lądowej Politechniki Krakowskiej Strona Główna -> Metody numeryczne Wszystkie czasy w strefie EET (Europa)
Strona 1 z 1

 
Skocz do:  
Nie możesz pisać nowych tematów
Nie możesz odpowiadać w tematach
Nie możesz zmieniać swoich postów
Nie możesz usuwać swoich postów
Nie możesz głosować w ankietach

fora.pl - załóż własne forum dyskusyjne za darmo
Powered by phpBB © 2001, 2005 phpBB Group
Regulamin