[ 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