029

Tworzenie skryptów i współpraca z plikami danych 29

10xi + 2x₂ + 3x;, = 1,

4xi + 50x₂ + 6x3 = 2,

7xi + 8x₂ + 90x3 = 3.

W skrypcie wykorzystano zapis macierzowy układu równań w postaci

Ax = b,

gdzie:

A - macierz kwadratowa o wymiarze n razy n, b - wektor kolumnowy o n elementach.

Przykładowy skrypt czytający i zapisujący zmienne na dysku w kodzie ASCII

% axb

% Skrypt jest zapamiętany w pliku axb.m

% Uruchomienie skryptu polega na wprowadzeniu polecenia % axb

% Skrypt axb.m, macierze A, b występują w tym samym katalogu. Opis skryptu jest wyswietalny na ekranie

Skrypt rozwiązuje układ liniowych rownan A*x=b');

A - macierz kwadratowa, b - wektor kolumnowy');

%

disp ( disp ( disp ( disp ( disp ( disp ( disp ( disp ( disp ( disp ( disp ( disp ( disp ( disp ( disp ( disp ( disp ( disp ( load

x - wektor kolumnowy rozwiązań');

Rozwiązanie jest wyznaczane za pomocą polecenia');

x = A\b');

Przykład ¹) ;

10*xl +    2 *x2 +    3 *x3 = 1') ;

4*xl + 50*x2 +    6*x3 = 2');

7 *xl +    8*x2    + 90*x4 = 3');

Macierz A zapamiętano w kodzie ASCII w pliku A.dat' 10    2    3');

6 ' ) ;

90’);

zapamietano w kodzie ASCII w pliku b.dat')

50

4

7    8

Wektor b 1' ) ;

2 ' ) ;

3 ' ) ; Macierz A.dat; ' ) ;

A jest wczytana jako plik ASCII A.dat');

b

x

j est j est

wczytany

zapisany

j ako jako

plik ASCII plik ASCII

b.dat' x.dat¹

)

disp(' Wektor disp(' Wektor % Polecenia skryptu

load A.dat; % wczytany plik zostanie zapamiętany jako disp(' Wczytana macierz współczynników układu rownan A disp(A); % macierz A zostanie wyświetlona na ekranie load b.dat; % wczytany plik zostanie zapamiętany jako wektor b disp(' Wczytany wektor prawej strony układu rownan b'); disp(b);% wektor b zostanie wyświetlony na ekranie x = A\b;

macierz A

) ;