[ Pobierz całość w formacie PDF ]
.3.Znaj�c warto�� w�asn� i odpowiadaj�cy jej wektor w�asny (obliczony np.przyprzyj�ciu dodatkowego warunku dotycz�cego d�ugo�ci tego wektora) mo�na równie�wyznaczy� drug� warto�� w�asn� i nast�pne warto�ci w�asne metod� redukcji,korzystaj�c z mo�liwo�ci rozk�adu dowolnego wektora wzgl�dem ró�nych mi�dzysob� wektorów w�asnych [2, 13].Rys.2.3Znajomo�� wektorów w�asnych jest szczególnie istotna dla macierzy Asymetrycznych o elementach rzeczywistych, gdy� np.dla macierzy kwadratowej,której kolumny s� wektorami w�asnymimo�na udowodni� zale�no�� [6, 7]oraz wykaza� ortogonalno�� (2.30) macierzy S.Obliczanie ekstremalnych co do modu�u warto�ci w�asnych za pomoc� wzorów(2.141), (2.143)���(2.145) oraz (2.148) odbywa si� w programie 2.4; tabulogrammodu�u Obliczenia tego programu jest nast�puj�cy:{Program 2.4}unit Obliczenia;interfaceusesWindows, Messages, SysUtils, Classes, Graphics, Controls,Forms, Dialogs, StdCtrls, Buttons;const nmax = 20;typeTabl1 = array [1.nmax,1.nmax] of Real;Tabl2 = array [1.nmax,1.2*nmax] of Real;.varForm3: TForm3;bl,det,lam,lam1,eps,r,r1,s,Spm,Spm1: Real;i,j,k,n,m,q,war: Integer;A,Am,Am1: Tabl1;nazwa: String;rep: Boolean;plik: Text;B: Tabl2;implementationuses Ustawienia, Informacje, Grafika, Podglad;label kon,omin,powt;{$R *.DFM}{procedure ElimGaussa(n,m: Integer; var A: Tabl2; var det: Real);}procedure TForm3.FormCreate(Sender: TObject);beginnazwa:='Dane1.dan';Dane.Lines.Text:='';end;procedure TForm3.BitBtn1Click(Sender: TObject);label omin,kon,powt;beginForm2.Show;AssignFile(plik,nazwa);Reset(plik);Readln(plik,n,eps);for i:=1 to n dofor j:=1 to n doRead(plik,A[i,j]);CloseFile(plik);if RadioButton1.Checked then war:=1;if RadioButton2.Checked then war:=2;if RadioButton3.Checked then war:=3;AssignFile(plik,Edit1.Text);Rewrite(plik);Writeln(plik,'PROGRAM 2.4.');if war=1 then beginWriteln(plik,'Obliczanie warto�ci w�asnej o module max');Writeln(plik,'i jej krotno�ci.');end;if war=2 then beginWriteln(plik,'Obliczanie ekstremalnych warto�ci w�asnych');Writeln(plik,'i ich krotno�ci.');end;if war=3 then beginWriteln(plik,'Obliczanie maksymalnego modu�u zespolonej');Writeln(plik,'warto�ci w�asnej.');end;Writeln(plik);Writeln(plik,'Stopie� macierzy - n = ',n:2);Writeln(plik,'Dok�adno�� oblicze� - eps = ',eps:11);Writeln(plik);Writeln(plik,'Macierz wspó�czynników:');for i:=1 to n do beginWriteln(plik,' wiersz nr ',i:2);k:=0; Write(plik,' ');for j:=1 to n do begink:=k+1;if k=5 then begink:=0; Writeln(plik);Write(plik,' ');end;Write(plik,' ',A[i,j]:13);end;Writeln(plik);end;Writeln(plik);rep:=False;powt:m:=2; Spm:=0;for i:=1 to n doSpm:=Spm+A[i,i];for i:=1 to n dofor j:=1 to n do begins:=0;for k:=1 to n dos:=s+A[i,k]*A[k,j];Am[i,j]:=s;end;Spm1:=0;for i:=1 to n doSpm1:=Spm1+Am[i,i];if not rep then lam:=Spm1/Spmelse lam:=Spm/Spm1;if war=3 then beginr:=Exp(0.5/m*Ln(0.5*Abs(Spm*Spm-Spm1)));goto omin;end;Spm:=Spm1;repeatfor i:=1 to n dofor j:=1 to n do begins:=0;for k:=1 to n dos:=s+A[i,k]*Am[k,j];Am1[i,j]:=s;end;Spm1:=0;for i:=1 to n doSpm1:=Spm1+Am1[i,i];if not rep then lam1:=Spm1/Spmelse lam1:=Spm/Spm1;bl:=Abs(lam1-lam);for i:=1 to n dofor j:=1 to n doAm[i,j]:=Am1[i,j];Spm:=Spm1; lam:=lam1;m:=m+1;until bl [ Pobierz całość w formacie PDF ]

  • zanotowane.pl
  • doc.pisz.pl
  • pdf.pisz.pl
  • listy-do-eda.opx.pl