Budownictwo 2006 Wydział Inżynierii Lądowej Politechniki Krakowskiej    FAQ   Szukaj   Użytkownicy   Grupy   Galerie   Rejestracja   Profil   Zaloguj się, by sprawdzić wiadomości   Zaloguj  wektory wlasne laborki nr3          Forum Budownictwo 2006 Wydział Inżynierii Lądowej Politechniki Krakowskiej Strona Główna -> Metody numeryczne Zobacz poprzedni temat :: Zobacz następny temat   Autor Wiadomość rzabaa Operator Łopaty Dołączył: 20 Wrz 2007 Posty: 19 Przeczytał: 0 tematów Ostrzeżeń: 0/5 Wysłany: Wto 23:41, 08 Kwi 2008    Temat postu: wektory wlasne laborki nr3 jaki byla tresc zadania?? wszyscy mieli to samo??? Post został pochwalony 0 razy Powrót do góry Zobacz poprzedni temat :: Zobacz następny temat   Autor Wiadomość Szczur Majster Dołączył: 07 Lis 2006 Posty: 143 Przeczytał: 0 tematów Pomógł: 1 raz Ostrzeżeń: 0/5 Skąd: Żywiec Płeć: Mężczyzna Wysłany: Śro 9:47, 09 Kwi 2008    Temat postu: ma ktoś może już to napisane? Post został pochwalony 0 razy Powrót do góry 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 Wysłany: Śro 15:26, 09 Kwi 2008    Temat postu: może by ktoś napisał już coś co było? Post został pochwalony 0 razy Powrót do góry Zobacz poprzedni temat :: Zobacz następny temat   Autor Wiadomość Szczur Majster Dołączył: 07 Lis 2006 Posty: 143 Przeczytał: 0 tematów Pomógł: 1 raz Ostrzeżeń: 0/5 Skąd: Żywiec Płeć: Mężczyzna Wysłany: Śro 16:44, 09 Kwi 2008    Temat postu: trzeba było obliczyć wektor i wartości własne macierzy 3 stopnia, za pomocą metody potęgowej albo metody iteracji odwrotnej. metoda potęgowa: Cytat: function[iter,lambda,x]=met_pot(A) n=length(A); xk=[1;1;1]; lamk=0.0; c=(xk'*xk)^(1/2); uk=xk/c; epsL=1; epsR=1; iter=1; while ((epsL>1e-6)&&(epsR>1e-6))&&(iter<200) xk1=A*uk; lamk1=uk'*xk1; c=(xk1'*xk1)^(1/2); xk1=xk1/c; epsL=abs((lamk1-lamk)/lamk1); epsR=([xk1-uk]'*[xk1-uk])^(1/2); iter=iter+1; lamk=lamk1; uk=xk1; end lambda=lamk1; x=xk1; end; A = [1 2 3 ; 4 5 6 ; 7 8 9] [iteracja,lambda,x]=met_pot(a) [wektory,x] = eig(A) Tak wyglądał mój program, który spisałem z kartki, którą dostałem od kumpla, który siedział obok mnie . Nie ja to napisałem (mam nadzieje, że autor się nie obrazi), nie wiem jak to działa i nie wiem czy 3 ostatnie linijki działają. Post został pochwalony 0 razy Ostatnio zmieniony przez Szczur dnia Śro 17:10, 09 Kwi 2008, w całości zmieniany 2 razy Powrót do góry Zobacz poprzedni temat :: Zobacz następny temat   Autor Wiadomość Woodchack Kierownik Budowy Dołączył: 19 Paź 2006 Posty: 335 Przeczytał: 0 tematów Ostrzeżeń: 0/5 Skąd: Bielsko-Biala Płeć: Mężczyzna Wysłany: Śro 19:41, 09 Kwi 2008    Temat postu: Na mailu naszym gr.4 jest ten program poprawiony działający w pełni...chwila to wkleje... A=[-2 5 9;4 -13 11;1 0 -18] function [iter,lambda,x]=pys(A) n=100; xk=[1;0;0]; lamk=0; c=(xk'*xk)^(1/2); uk=xk/c; epsL=1; epsR=1; nr_iteracji=1; while(epsL>10^-6)&(epsR>10^-6) xk1=A*uk; lamk1=uk'*xk1; c=[xk1'*xk1]^(1/2); uk1=xk1/c; epsL=abs((lamk1-lamk)/lamk1); epsR=((uk1-uk)'*(uk1-uk))^(1/2); nr_iteracji=nr_iteracji+1; if (nr_iteracji>n) printf("proces rozbiezny"); return; end lamk=lamk1; uk=uk1; end iter=nr_iteracji; lambda=lamk1; x=uk1; endfunction [iter,lambda,x]=pys(A) Prosze Post został pochwalony 0 razy Ostatnio zmieniony przez Woodchack dnia Śro 19:46, 09 Kwi 2008, w całości zmieniany 1 raz Powrót do góry Wyświetl posty z ostatnich: Wszystkie Posty1 Dzień7 Dni2 Tygodnie1 Miesiąc3 Miesiące6 Miesięcy1 Rok Najpierw StarszeNajpierw Nowsze         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: Wybierz forum .....::::: Ważne!!! :::::.....----------------OgłoszeniaInne rocznikiSpecjalizacjeHydeparkMaterialy i nawierzchnie drogowe.Podstawy sytemu transportuMechanika budowliMechanika gruntow i fundamentowanieTeoria sprezystosci i plastycznosciMetody komputeroweMiernictwoJezyczki obceHydeParkJezyki obceMechanika BudowliMechanika gruntow i fundamentowanie.Fizyka budowliMetody komputeroweAutoCADMiernictwoWstęp do mostów i mosty drewnianeKonstrukcje betonoweKonstrukcje metaloweWplyw srodowiska na budowleArchiwum - Mosty .....::::: Specjalizacje :::::.....----------------Drogi KolejoweDrogi, Ulice i Autostrady - Gr. 1Drogi, Ulice i Autostrady - Gr. 2Konstrukcje Budowlane i Inżynierskie - Gr. 1Konstrukcje Budowlane i Inżynierskie - Gr. 2Konstrukcje Budowlane i Inżynierskie - Gr. 3Konstrukcje Budowlane i Inżynierskie - Gr. 4Drogi, Ulice i Autostrady Konstrukcje Budowlane i InżynierskieMosty i Budowle PodziemneTechnologia i Organizacja BudownictwaTeoria Konstrukcji Inżynierskich .....::::: Offtopic :::::.....----------------HydeparkGiełdaKosz .....::::: Grupy :::::.....----------------Grupa 1Grupa 2Grupa 3Grupa 4Grupa 5Grupa 6Grupa 7Grupa 8Grupa 9Grupa 10Grupa 11 .....::::: արխիվ :::::.....----------------Hydepark - pierwszy rok...Hydepark - drugi rok...AnalizaAlgebraBudownictwo ogólneFizykaGeologiaGeometria wykreślna/Rysunek technicznyInformatykaMateriały budowlaneMechanikaMetody statystyczneMetody numeryczneMechanika gruntówPodstawy ekonomiiTechnologia betonuWytrzymałość materiałówPozostałeTeoria sprezystosci  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 Na skróty:.....::::: Ważne!!! :::::..........::::: Offtopic :::::..........::::: արխիվ :::::.....KoszInne rocznikiOgłoszeniaAlgebraPozostałeMetody numeryczneAnalizaBudownictwo ogólneGeometria wykreślna/Rysunek technicznyWytrzymałość materiałówFizykaMechanikaMetody statystyczneInformatykaGeologiaMechanika gruntówTechnologia betonuMateriały budowlaneGiełdaPodstawy ekonomii