Scribd
Upload
40 views16 pages

TD Matlab

The document contains examples of operations on matrices in MATLAB/Octave including finding eigenvalues, inverses, rank, null space, and LU decompositions of matrices. Various matrix multiplications and extractions of elements, rows and columns are also demonstrated.

Uploaded by

jalilpadawou55
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
40 views16 pages

TD Matlab

The document contains examples of operations on matrices in MATLAB/Octave including finding eigenvalues, inverses, rank, null space, and LU decompositions of matrices. Various matrix multiplications and extractions of elements, rows and columns are also demonstrated.

Uploaded by

jalilpadawou55
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 16

>> EXERCICE3

A=[3 -1 0;0 1 0;0 0 2]

A=

3 -1 0

0 1 0

0 0 2

>> B=[3 5 1;7 -2 4;-6 3 2]

B=

3 5 1

7 -2 4

-6 3 2

>> C=[5 6 5;6 10 9;5 9 10]

C=

5 6 5

6 10 9

5 9 10

>> a=eig(A)

a=

>> b=eig(B)

b=

-7.3538 + 0.0000i

5.1769 + 2.0832i

5.1769 - 2.0832i

>> c=eig(C)

c=

0.5660

1.9659

22.4681
>> [V,L]=eig(A)

V=

1.0000 0.4472 0

0 0.8944 0

0 0 1.0000

L=

3 0 0

0 1 0

0 0 2

>> [V,L]=eig(A)

V=

1.0000 0.4472 0

0 0.8944 0

0 0 1.0000

L=

3 0 0

0 1 0

0 0 2

>> A+C

ans =

8 5 5

6 11 9

5 9 12

>> A*C

ans =

9 8 6

6 10 9

10 18 20
>> A.*C

ans =

15 -6 0

0 10 0

0 0 20

>> v=[1;2;3]

v=

>> R=A+v

Error using +

Matrix dimensions must agree.

>> B=[3 5 1 1;7 -2 4 2;-6 3 2 3]

B=

3 5 1 1

7 -2 4 2

-6 3 2 3

>> k=null(B)

k=

0.1191

-0.0959

-0.6349

0.7573

>>EXERCICE4

>> A=[2 1 -3;1 -2 1];

>> b=[1 2]';

>> rank(A)

ans =

2
>> B=[2 1 -3 1;1 -2 1 2]

B=

2 1 -3 1

1 -2 1 2

>> rank(B)

ans =

>> c=[1 1 1]';

>> y=A*C

Undefined function or variable 'C'.

Did you mean:

>> y=A*c

y=

>> X=inv(A)*b

Error using inv

Matrix must be square.

>>EXO 6

>> A=[1 2 3;2 3 1;3 1 2]

A=

1 2 3

2 3 1

3 1 2

>> A([2 3],[1 3])

ans =

2 1

3 2

>> A([2 3],1:2)


ans =

2 3

3 1

>> A([2 3],:)

ans =

2 3 1

3 1 2

>> A([2 3],end)

ans =

>> A(:)

ans =

>> A(5)

ans =

>> reshape(A(:),size(A))

ans =
1 2 3

2 3 1

3 1 2

>>>> % EXO7

>> A=[1 2 3; 2 3 1;3 1 2]

A=

1 2 3

2 3 1

3 1 2

>> C=[A,zeros(3,2);zeros(2,3),eye(2)]

C=

1 2 3 0 0

2 3 1 0 0

3 1 2 0 0

0 0 0 1 0

0 0 0 0 1

>>>> TP1

1)

>> A=magic(5)

A=

17 24 1 8 15

23 5 7 14 16

4 6 13 20 22

10 12 19 21 3

11 18 25 2 9

>> sum_lignes=sum(A,2)

sum_lignes =
65

65

65

65

65

>> sum_colonnes=sum(A,1)

sum_colonnes =

65 65 65 65 65

>> sum_diagonale1=sum(diag(A))

sum_diagonale1 =

65

>> sum_diagonale2=sum(diag(fliplr(A)))

sum_diagonale2 =

65

>> A(3,:)

ans =

4 6 13 20 22

>> A(1,2)

ans =

24

>>2)

>> p=[3 -1 7]

p=

3 -1 7

>> roots(p)

ans =

0.1667 + 1.5184i

0.1667 - 1.5184i
>>

>> %plus a propose des matrices et tableaux

>> A=[10 7 8 7;7 5 6 7;8 6 10 9;7 5 9 10]

A=

10 7 8 7

7 5 6 5

8 6 10 9

7 5 9 10

>> b1=[32 23 33 31]'

>> b1 =

32

23

33

31

>> b2=[32.1 22.9 33.1 30.9]'

b2 =

32.1

22.9

33.1

30.9

>> x=inv(A)*b1

x=

1.000

1.000

1.000

1.000

>> x2=inv(A)*b2

x2 =

4.9048

-5.4762

2.7190
-0.0524

>> y=eig(A)

y=

30.6789

3.3601

0.3157

0.6453

>> >>TP4

>> A=rand(7)

A=

0.8147 0.5469 0.8003 0.0357 0.6555 0.8235 0.7655

0.9058 0.9575 0.1419 0.8491 0.1712 0.6948 0.7952

0.1270 0.9649 0.4218 0.9340 0.7060 0.3171 0.1869

0.9134 0.1576 0.9157 0.6787 0.0318 0.9502 0.4898

0.6324 0.9706 0.7922 0.7577 0.2769 0.0344 0.4456

0.0975 0.9572 0.9595 0.7431 0.0462 0.4387 0.6463

0.2785 0.4854 0.6557 0.3922 0.0971 0.3816 0.7094

>> inv(A)

ans =

0.2150 0.3124 -0.5303 0.2803 0.9697 -0.5554 -0.7391

0.7551 0.6244 -0.6416 -0.6979 0.1878 1.9664 -2.7736

0.1938 -0.8510 -0.1891 0.3615 0.4972 0.4637 -0.1897

-1.2070 -0.0659 0.9475 0.6014 -0.0246 -1.1940 1.8149

0.2754 -0.5342 1.2076 -0.1488 -0.0479 -1.5358 1.5155

0.5698 0.4188 0.0176 0.3729 -1.2926 1.4092 -1.8184

-0.4570 0.1211 0.1233 -0.4794 -0.2533 -1.4436 3.5402

>> eye(A)

Error using eye

N-dimensional arrays are not supported.

>> inv(A)*A
ans =

1.0000 0.0000 0.0000 0.0000 -0.0000 0.0000 -0.0000

0.0000 1.0000 0.0000 -0.0000 0.0000 0.0000 0.0000

-0.0000 -0.0000 1.0000 0.0000 0.0000 0.0000 -0.0000

-0.0000 -0.0000 -0.0000 1.0000 -0.0000 -0.0000 -0.0000

0.0000 -0.0000 -0.0000 -0.0000 1.0000 -0.0000 0.0000

0.0000 0.0000 0.0000 0.0000 0.0000 1.0000 0.0000

-0.0000 0.0000 -0.0000 0.0000 -0.0000 -0.0000 1.0000

>> [L,U,P]=lu(A,1)

Error using lu

Permutation flag must be 'matrix' or 'vector'.

>> [L,U,P]=lu(A)

L=

1.0000 0 0 0 0 0 0

0.1390 1.0000 0 0 0 0 0

0.9917 0.8497 1.0000 0 0 0 0

0.8920 0.4309 0.1411 1.0000 0 0 0

0.1068 0.9972 -0.5589 0.5457 1.0000 0 0

0.6923 0.9136 0.1090 0.4916 0.4693 1.0000 0

0.3049 0.4638 -0.2361 0.3868 0.4414 0.0716 1.0000

U=

0.9134 0.1576 0.9157 0.6787 0.0318 0.9502 0.4898

0 0.9430 0.2944 0.8396 0.7016 0.1850 0.1188

0 0 -1.0164 -0.5374 -0.4565 -0.4047 0.2086

0 0 0 -0.8557 0.3892 -0.0467 0.2480

0 0 0 0 -1.1244 -0.0479 0.4568

0 0 0 0 0 -0.7029 -0.3610

0 0 0 0 0 0 0.2825
P=

0 0 0 1 0 0 0

0 0 1 0 0 0 0

0 1 0 0 0 0 0

1 0 0 0 0 0 0

0 0 0 0 0 1 0

0 0 0 0 1 0 0

0 0 0 0 0 0 1

>> det(A)=prod(diag(U))*(-1)^P

Subscript indices must either be real positive integers or logicals.

det(A)=prod(diag(U))*(1)^P

Subscript indices must either be real positive integers or logicals.

>> [piv,L,U]=lu(A)

piv =

1.0000 0 0 0 0 0 0

0.1390 1.0000 0 0 0 0 0

0.9917 0.8497 1.0000 0 0 0 0

0.8920 0.4309 0.1411 1.0000 0 0 0

0.1068 0.9972 -0.5589 0.5457 1.0000 0 0

0.6923 0.9136 0.1090 0.4916 0.4693 1.0000 0

0.3049 0.4638 -0.2361 0.3868 0.4414 0.0716 1.0000

L=

0.9134 0.1576 0.9157 0.6787 0.0318 0.9502 0.4898

0 0.9430 0.2944 0.8396 0.7016 0.1850 0.1188

0 0 -1.0164 -0.5374 -0.4565 -0.4047 0.2086

0 0 0 -0.8557 0.3892 -0.0467 0.2480

0 0 0 0 -1.1244 -0.0479 0.4568

0 0 0 0 0 -0.7029 -0.3610

0 0 0 0 0 0 0.2825
U=

0 0 0 1 0 0 0

0 0 1 0 0 0 0

0 1 0 0 0 0 0

1 0 0 0 0 0 0

0 0 0 0 0 1 0

0 0 0 0 1 0 0

0 0 0 0 0 0 1
b=[1 2 3 4 5 6 7]';
y=zeros(size(b));
x=zeros(size(b));
%Foward substitution
for i=1:size(A,1)
y(i)=b(piv(i));
for j=1:i-1
y(i)=y(i)-L(i,j)*y(j);
end
end

TP2

1)
x=linspace(-2,2);
y=exp(x);
plot(x,y);
grid on
title('courbe exp(x)')
xlabel('x');
ylabel('y');
2)
R=3;
w=linspace(0,3*R);
x=R*(w-sin(w));
y=R*(1-cos(w));
plot(x,y);
hold on
grid on
xlabel('x')
ylabel('y')
title('trace de la courbe en cycloide');
3)
x=linspace(-10,10);
y=linspace(-10,10);
r=sqrt(x.*x+y.*y);
z=10.*cos(r)./r+2;
plot3(x,y,z)
grid('on')
hold('on')
xlabel('X')
ylabel('Y')
zlabel('Z')
title('trace de la courbe 3D')
3)

x=linspace(-0.1,0.1);

v=x.*sin(1./x);

>> plot(x,v,'r');

>> grid on

>> hold on

>> plot(x,-x,'b.')

plot(x,x,'b.')

>> xlabel('axe des x')

>> ylabel('axe des y')

>> gtext('y=x')

>> gtext('y=-x')

>> gtext('y=x*sin(1/x)')
>>

You might also like