|
Seite 1 von 1
|
[ 10 Beiträge ] |
|
Alte Matlab Tests Programmiert
Autor |
Nachricht |
Spielli
Registriert: 10/2010 Beiträge: 22
|
Alte Matlab Tests Programmiert
Sommersemester 2008 A1: A2: A3: | | | | Code: function [ lam,Mx ] = poweriteration( A,x0,tau ) Mx=x0; lam=x0'*A*x0; test=1; i=1; while test Mx=[Mx,A*Mx(:,i)/norm(A*Mx(:,i))]; lam=[lam,Mx(:,i+1)'*A*Mx(:,i+1)]; i=i+1; if abs(lam(i)<=tau) tautest=tau; else tautest=tau*abs(lam(i)); end test=(norm(A*Mx(:,i)-lam(i)*Mx(:,i))>tau) * (abs(lam(i-1)-lam(i))>tautest); end end | | | | |
A4:
|
So 23-01-2011 16:06:47 |
|
|
Spielli
Registriert: 10/2010 Beiträge: 22
|
Alte Matlab Tests Programmiert
Wintersemester 2008/2009 A1: | | | | Code: function [L,U] = computeLU(A) n=length(A); L=eye(n); U=zeros(n); U(1,:)=A(1,:); for i=1:n-1 for j=i+1:n L(j,i)=(1/U(i,i))*(A(j,i)-L(j,1:i-1)*U(1:i-1,i)); U(i+1,j)=A(i+1,j)-L(i+1,1:i)*U(1:i,j); end end end | | | | |
A2: | | | | Code: function [ L,U ] = computeLU2( A ) n=length(A); L=eye(n); U=zeros(n); U(1,:)=A(1,:); for i=1:n-1 L(i+1:n,i)=(1/U(i,i))*(A(i+1:n,i)-L(i+1:n,1:i-1)*U(1:i-1,i)); U(i+1,i+1:n)=A(i+1,i+1:n)-L(i+1,1:i)*U(1:i,i+1:n); end end | | | | |
A3: A4:
|
So 23-01-2011 16:09:19 |
|
|
Spielli
Registriert: 10/2010 Beiträge: 22
|
Alte Matlab Tests Programmiert
Wintersemester 2009/2010 A2: A3: A4: A5: | | | | Code: function phi = richardson( f,x,h0,eps ) fx=f(x); h=h0; y(1)=(f(x+h0)-fx)/h0; nev(1)=neville(h,y,0); h=[h,h(length(h))/2]; y(2)=(f(x+h(2))-fx)/h(2); nev(2)=neville(h,y,0); i=2; while abs(nev(2)-nev(1))>eps*max(abs(nev(1)),abs(nev(2))) i=i+1; h=[h,h(length(h))/2]; y(i)=(f(x+h(i))-fx)/h(i); nev(1)=nev(2); nev(2)=neville(h',y',0); end phi=nev(2); end | | | | |
|
So 23-01-2011 16:13:56 |
|
|
Spielli
Registriert: 10/2010 Beiträge: 22
|
Alte Matlab Tests Programmiert
Sommersemester 2010 A1: A2: A3: A4:
|
So 23-01-2011 16:15:51 |
|
|
Sahin
Registriert: 10/2007 Beiträge: 25 + 7
|
Re: Alte Matlab Tests Programmiert
märz 2010 bsp 1 (8pkt) hat jemand vielleicht aufgabe 2, 3 und 4 richtig? vor allem rootpoly wäre sehr interessant. bei 3. bin ich soweit gekommen: bzw 2. aber das liefert weiß irgendwer was damit nicht stimmt?
|
Di 14-06-2011 10:34:51 |
|
|
MarK
Registriert: 10/2007 Beiträge: 29 + 498 Wohnort: Pfeilheim, Wien8
Studium: (alt) Bachelor Computerwissenschaften
|
Bitte benutzt die [ code][/code] tags wenn ihr code eingibt:
[ code]Hello World![/ code]
Hab eure Posts editiert.
_________________
|
Di 14-06-2011 11:26:10 |
|
|
Sahin
Registriert: 10/2007 Beiträge: 25 + 7
|
Re: Alte Matlab Tests Programmiert
Bs4) 24.Jun.2010 bzw. Bsp4) 4.Okt.2010
Zuletzt geändert von Sahin am Do 16-06-2011 10:38:14, insgesamt 3-mal geändert.
|
Do 16-06-2011 09:47:22 |
|
|
Sahin
Registriert: 10/2007 Beiträge: 25 + 7
|
Re: Alte Matlab Tests Programmiert
Bsp3) 24.Jun.2010 integrate
|
Do 16-06-2011 10:23:11 |
|
|
Sahin
Registriert: 10/2007 Beiträge: 25 + 7
|
Re: Alte Matlab Tests Programmiert
Bsp1) 21.Jan.2010 (mit schleifen) | | | | Code: function [L,U] = computeLU(A) n = size(A,1); L = eye(n); U = zeros(n);
for i =1:n for k = 1:n if k>=i sum = 0; for j = 1:i-1 sum = sum + L(i,j)*U(j,k); end U(i,k)=A(i,k)-sum; else sum = 0; for j = 1:i-1 sum = sum + L(i,j)*U(j,k); end L(i,k)=(1/U(k,k))*(A(i,k)-sum); end end end end | | | | |
Bsp 4) mit Vektorarithmetik schleifen eliminiert | | | | Code: function [L,U] = computeLU2(A) n = size(A,1); L = eye(n); U = zeros(n,n);
for i = 1:n for k = 1:n if k >= i U(i,k) = A(i,k)-sum([L(i,[1:i-1])*U([1:i-1],k)]); else L(i,k)=(1/U(k,k))*(A(i,k)-sum([L(i,[1:i-1])*U([1:i-1],k)])); end end end end | | | | |
|
Do 16-06-2011 11:54:47 |
|
|
Adoncia
Registriert: 10/2009 Beiträge: 25 + 7
|
Ich hätte eine Frage bzgl Nachtest WS 10/11: Bei Bsp 2 wo wir den Grad n eines Polynoms bestimmen müssen.
kann man das nicht einfacht mit n=length(p) ? ermitteln?
|
Do 16-06-2011 12:31:14 |
|
|
|
Seite 1 von 1
|
[ 10 Beiträge ] |
|
Wer ist online? |
Mitglieder in diesem Forum: 0 Mitglieder und 1 Gast |
|
Du darfst neue Themen in diesem Forum erstellen. Du darfst Antworten zu Themen in diesem Forum erstellen. Du darfst deine Beiträge in diesem Forum nicht ändern. Du darfst deine Beiträge in diesem Forum nicht löschen. Du darfst keine Dateianhänge in diesem Forum erstellen.
|
|