abc Matlab - примјери и објашњења
5. СИСТЕМИ ЛИНЕАРНИХ АЛГЕБАРСКИХ ЈЕДНАЧИНА

5.2. Системи линеарних алгебарских једначина - Примјери

ПРИМЈЕР 1: Креирати m-фајл Cramer.m за рјешавање система линеарних алгебарских једначина користећи Крамерово правило:

%Rјesavanje sistema AX=B Cramerovim pravilom
%
function X=Cramer(A,B)
[m,n]=size(A);
if m~=n,
    error('Matrica nije kvadratna'),
end
if det(A)==0,
    error('Matrica je singularna'),
end
for j=1:n,
  C=A; C(:,j)=B; X(j)=det(C)/det(A);
end
X=X';

ПРИМЈЕР 2: Ријешити систем једначина матричном методом и користећи креирани фајл Cramer.m. Упоредити добијена рјешења.

>> M=[2,-3,1;-3,5,-2;1,-2,3];
N=[11;-19;14];
>> X1=inv(M)*N
X1 =
    1.0000
   -2.0000
    3.0000
>> X2=Cramer(M,N)
X2 =
     1
    -2
     3
 

ПРИМЈЕР 3: Ријешити систем једначина матричном методом:

>> syms p q
A=[-p q 1;1 -q -p;-1 p*q 1]
A =
[  -p,   q,   1]
[   1,  -q,  -p]
[  -1, p*q,   1]
>> B=[1;1;q]
B =
 1
 1
 q
>> X=inv(A)*B
X =
 -(p+1)/(p^2+p-2)-1/(p^2+p-2)+1/(p^2+p-2)*q
                            (p+1)/(p^2+p-2)
 -1/(p^2+p-2)-(p+1)/(p^2+p-2)-1/(p^2+p-2)*q
    

ПРИМЈЕР 4: Ријешити систем једначина користећи фајл Cramer.m

>> Z=[-2 3 1;1 3 -2;1 -6 1]
Z =
    -2     3     1
     1     3    -2
     1    -6     1
>> Z1=[1;1;3]
Z1 =
     1
     1
     3
>> Cramer(Z,Z1)
??? Error using ==> Cramer
Matrica je singularna

Системи линеарних алгебарских једначина    <    Index    >    Вјежбе