Faculté Sultan Moulay Slimane
Faculté Polydisciplinaire-Beni Mellal
LP ENERGIES RENOUVELABLES ER
COMPTE-RENDU
Réalisé par : Encadré par :
Mohammed amine el ansary [Link] GAGA
Mohammed Kartit
ANNEE UNIVERSITAIRE 2021/2022
1
Table de matières :
Introduction……………………………………………………………3
TP N°1…………………………………………………………………………04
TP N°2…………………………………………………………………………18
TP N°3…………………………………………………………………………24
Méthode numérique……………………………………………28
Méthode de dichotomie…………………………………28
Méthode de newton……………………………………………29
Intégrale…………………………………………………………………31
Série de fourrier……………………………………………33
Conclusion………………………………………………………………35
2
Partout dans le monde, des millions d'ingénieurs et de scientifiques
utilisent MATLAB pour analyser et concevoir les systèmes et produits de demain.
MATLAB est présent dans des systèmes automobiles de sécurité active, des véhicules
spatiaux, des appareils de surveillance médicale, des réseaux électriques intelligents et
des réseaux mobiles LTE. Il est utilisé dans les domaines de l'apprentissage
automatique, le traitement du signal, la vision par ordinateur, les communications, la
finance computationnelle, la conception de contrôleurs, la robotique et bien plus .
La plate-forme MATLAB est optimisée pour résoudre les problèmes
scientifiques et techniques. Le langage MATLAB, basé sur les matrices, est le moyen le
plus naturel au monde pour exprimer les mathématiques computationnelles. Les
graphiques intégrés permettent de visualiser facilement les données afin d'en dégager
des informations. Grâce à la vaste bibliothèque de boîtes à outils prédéfinies, vous
pouvez commencer directement par les algorithmes essentiels à votre domaine.
L'environnent bureau encourage l'expérimentation, l'exploration et la découverte. Les
outils et les fonctionnalités MATLAB sont tous testés rigoureusement. Ils sont conçus
pour fonctionner conjointement.
➢ Intérêts :
Programmation infiniment plus rapide pour le calcul et pour l’affichage.
Une librairie très riche.
Possibilité d’inclure un programme en C/C++.
Langage interprété : Pas de compilation donc pas d’attente pour compiler.
Possibilité d’exécuter du code en dehors du programme.
Code facile à comprendre et très lisible.
Une aide très bien faite.
➢ Inconvénients :
Vitesse de calcul moins rapide qu’en C/C++.
Payant (équivalent gratuit scilab).
Application auto-exécutable peu pratique.
3
Cette première séance du TP a pour
but de nous familiariser avec
l’usage de MATLAB et de prendre les
bases de programmation sous ce
logiciel. On a commencé par écrire
un programme pour faire des calcule
sur les matrices et les vecteurs.
4
➢ Exercice 1 :
%script 1 premier script
a=1 ;
b=2 ;
c=a+b;
>> premierscript
a=1
b=2
c=3
Interprétation des résultats :
Dans ce premier programme en fait la somme de deux nombres.
➢ Exercice 2 :
clear all
clc
A=[1 2 -1 1;-1 1 0 3]
B=find(A>0)
find(A(B))
A=
1 2 -1 1
-1 1 0 3
B=
1
3
4
7
8
ans =
5
1
2
3
4
5
Interprétation des résultats :
Dans ce exercice en utiliser la fonction find pour trouver les nombres qui
supérieur de 0.
➢ Exercice 3:
% Question 1:
clear all
clc
u=[1;-1;2]'
v=[10;-1;3]'
w=[5;-1;4]'
A=3*u
B=norm(u,2)
C=2*u-v+5*w
D=norm(2*u-v+5*w,1)
E=norm(w-4*v,inf)
% question 2 :
u= 1 -1 2
v= 10 -1 3
w= 5 -1 4
A= 3 -3 6
B= 2.4495
C= 17 -6 21
D= 44
6
Interprétation des résultats :
Ce exercice fait les opération sur les vecteur. La fonction norm
➢ Exercice 4 :
clc
clear all;
u=11-7i;
v=-1+3i;
A=abs(u)
B=abs(v)
C=u*conj(v)
D=v*conj(u)
E=real(u^3+v^2)
F=imag(u^3+v^2)
A = 13.0384
B = 3.1623
C = -32.0000 -26.0000i
D = -32.0000 +26.0000i
E = -294
F = -2204
Interprétation des résultats :
Ce exercice appliquer des opération sur les nombre complexe. Les fonctions img et
real retournent la partie imaginaire et la partie réelle et les fonctions abs et
retournent le module.
➢ Exercice 5 :
clear all
clc
% question 1 :
A=[1 2 3;-5 3 1;-10 0 3]
B=[1 -2 5;6 1 -1;1 0 -2]
C=[10 -5;3 1]
A*B
B*A
A*B'
% question 2 :
D=eye(3)-B*B'
7
% question 3 :
S=det(A)
F=det(B)
J=det(C)
H=det(D)
E=A*A'
det(E)
% question 4 :
inv(A)
inv(B)
inv(C)
inv(D)
F=A*A'
inv(F)
% question 5 :
eig(F)
% question 6 :
eig(A)
[V]=eig(A)
A =
1 2 3
-5 3 1
-10 0 3
B =
1 -2 5
6 1 -1
1 0 -2
C =
10 -5
3 1
ans =
16 0 -3
14 13 -30
-7 20 -56
ans =
-39 -4 16
11 15 16
21 2 -3
ans =
12 5 -5
-6 -28 -7
5 -63 -16
D =
8
-29 1 9
1 -37 -8
9 -8 -4
S = 109
F = -29
J = 25
H = 421.0000
E =
14 4 -1
4 35 53
-1 53 109
ans = 1.1881e+04
ans =
0.0826 -0.0550 -0.0642
0.0459 0.3028 -0.1468
0.2752 -0.1835 0.1193
ans =
0.0690 0.1379 0.1034
-0.3793 0.2414 -1.0690
0.0345 0.0690 -0.4483
ans =
0.0400 0.2000
-0.1200 0.4000
ans =
0.1995 -0.1615 0.7720
-0.1615 0.0831 -0.5297
0.7720 -0.5297 2.5463
F =
14 4 -1
4 35 53
-1 53 109
ans =
0.0847 -0.0412 0.0208
-0.0412 0.1284 -0.0628
0.0208 -0.0628 0.0399
ans =
5.4755
15.8795
136.6450
ans =
2.2548 + 6.2197i
2.2548 - 6.2197i
2.4903 + 0.0000i
V =
2.2548 + 6.2197i
9
2.2548 - 6.2197i
2.4903 + 0.0000i
Interprétation des résultats :
Dans ce exercice en faire des opérations sur les matrices .1ere on déclarer
les matrices en suite en faire le fonction det et inv pour trever le
determin et la invers de la matrice .
➢ Exercice 6 :
clear all
clc
A=[1 -1 7;-4 2 11;8 0 3];
B=[3 -2 -1;7 8 6;5 1 3];
a=3*A
b=A.*B
c=A*B
d=A./B
e=cos(A)
f=cos(B)
g=exp(B)
a =
3 -3 21
-12 6 33
24 0 9
b =
3 2 -7
-28 16 66
40 0 9
c =
31 -3 14
57 35 49
39 -13 1
d =
0.3333 0.5000 -7.0000
-0.5714 0.2500 1.8333
1.6000 0 1.0000
e =
10
0.5403 0.5403 0.7539
-0.6536 -0.4161 0.0044
-0.1455 1.0000 -0.9900
f =
-0.9900 -0.4161 0.5403
0.7539 -0.1455 0.9602
0.2837 0.5403 -0.9900
g =
1.0e+03 *
0.0201 0.0001 0.0004
1.0966 2.9810 0.4034
0.1484 0.0027 0.0201
➢ Exercice 8 :
clear all
clc
A=[1 -1 7;-4 2 11;8 0 3];
B=[3 -2 -1;7 8 6;5 1 3];
v=[1 2 4 7]
A(v)=A(v)+0.01
B(v)=abs(B(v))
v = 1 2 4 7
A =
1.0100 -0.9900 7.0100
-3.9900 2.0000 11.0000
8.0000 0 3.0000
B =
3 2 1
7 8 6
5 1 3
➢ Exercice 9 :
clear all
clc
% Question 1:
V1=linspace(-5,5,11)
% Question 2:
v2=linspace(-500,500,1001)
% Question 3:
v3=linspace(0,pi,10)
11
V1 = -5 -4 -3 -2 -1 0 1 2 3 4 5
v2 =
Columns 1 through 23
-500 -499 -498 -497 -496 -495 -494 -493 -492 -491 -490 -489 -488 -487 -486
-485 -484 -483 -482 -481 -480 -479 -478
Columns 24 through 46
-477 -476 -475 -474 -473 -472 -471 -470 -469 -468 -467 -466 -465 -464 -463
-462 -461 -460 -459 -458 -457 -456 -455
Columns 47 through 69
-454 -453 -452 -451 -450 -449 -448 -447 -446 -445 -444 -443 -442 -441 -440
-439 -438 -437 -436 -435 -434 -433 -432
Columns 70 through 92
-431 -430 -429 -428 -427 -426 -425 -424 -423 -422 -421 -420 -419 -418 -417
-416 -415 -414 -413 -412 -411 -410 -409
Columns 93 through 115
-408 -407 -406 -405 -404 -403 -402 -401 -400 -399 -398 -397 -396 -395 -394
-393 -392 -391 -390 -389 -388 -387 -386
Columns 116 through 138
-385 -384 -383 -382 -381 -380 -379 -378 -377 -376 -375 -374 -373 -372 -371
-370 -369 -368 -367 -366 -365 -364 -363
Columns 139 through 161
-362 -361 -360 -359 -358 -357 -356 -355 -354 -353 -352 -351 -350 -349 -348
-347 -346 -345 -344 -343 -342 -341 -340
Columns 162 through 184
-339 -338 -337 -336 -335 -334 -333 -332 -331 -330 -329 -328 -327 -326 -325
-324 -323 -322 -321 -320 -319 -318 -317
Columns 185 through 207
-316 -315 -314 -313 -312 -311 -310 -309 -308 -307 -306 -305 -304 -303 -302
-301 -300 -299 -298 -297 -296 -295 -294
Columns 208 through 230
12
-293 -292 -291 -290 -289 -288 -287 -286 -285 -284 -283 -282 -281 -280 -279
-278 -277 -276 -275 -274 -273 -272 -271
Columns 231 through 253
-270 -269 -268 -267 -266 -265 -264 -263 -262 -261 -260 -259 -258 -257 -256
-255 -254 -253 -252 -251 -250 -249 -248
Columns 254 through 276
-247 -246 -245 -244 -243 -242 -241 -240 -239 -238 -237 -236 -235 -234 -233
-232 -231 -230 -229 -228 -227 -226 -225
Columns 277 through 299
-224 -223 -222 -221 -220 -219 -218 -217 -216 -215 -214 -213 -212 -211 -210
-209 -208 -207 -206 -205 -204 -203 -202
Columns 300 through 322
-201 -200 -199 -198 -197 -196 -195 -194 -193 -192 -191 -190 -189 -188 -187
-186 -185 -184 -183 -182 -181 -180 -179
Columns 323 through 345
-178 -177 -176 -175 -174 -173 -172 -171 -170 -169 -168 -167 -166 -165 -164
-163 -162 -161 -160 -159 -158 -157 -156
Columns 346 through 368
-155 -154 -153 -152 -151 -150 -149 -148 -147 -146 -145 -144 -143 -142 -141
-140 -139 -138 -137 -136 -135 -134 -133
Columns 369 through 391
-132 -131 -130 -129 -128 -127 -126 -125 -124 -123 -122 -121 -120 -119 -118
-117 -116 -115 -114 -113 -112 -111 -110
Columns 392 through 414
-109 -108 -107 -106 -105 -104 -103 -102 -101 -100 -99 -98 -97 -96 -95 -94
-93 -92 -91 -90 -89 -88 -87
Columns 415 through 437
-86 -85 -84 -83 -82 -81 -80 -79 -78 -77 -76 -75 -74 -73 -72 -71 -70
-69 -68 -67 -66 -65 -64
Columns 438 through 460
13
-63 -62 -61 -60 -59 -58 -57 -56 -55 -54 -53 -52 -51 -50 -49 -48 -47
-46 -45 -44 -43 -42 -41
Columns 461 through 483
-40 -39 -38 -37 -36 -35 -34 -33 -32 -31 -30 -29 -28 -27 -26 -25 -24
-23 -22 -21 -20 -19 -18
Columns 484 through 506
-17 -16 -15 -14 -13 -12 -11 -10 -9 -8 -7 -6 -5 -4 -3 -2 -1 0 1
2 3 4 5
Columns 507 through 529
6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23
24 25 26 27 28
Columns 530 through 552
29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46
47 48 49 50 51
Columns 553 through 575
52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69
70 71 72 73 74
Columns 576 through 598
75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92
93 94 95 96 97
Columns 599 through 621
98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113
114 115 116 117 118 119 120
Columns 622 through 644
121 122 123 124 125 126 127 128 129 130 131 132 133 134 135
136 137 138 139 140 141 142 143
Columns 645 through 667
144 145 146 147 148 149 150 151 152 153 154 155 156 157 158
159 160 161 162 163 164 165 166
Columns 668 through 690
14
167 168 169 170 171 172 173 174 175 176 177 178 179 180 181
182 183 184 185 186 187 188 189
Columns 691 through 713
190 191 192 193 194 195 196 197 198 199 200 201 202 203 204
205 206 207 208 209 210 211 212
Columns 714 through 736
213 214 215 216 217 218 219 220 221 222 223 224 225 226 227
228 229 230 231 232 233 234 235
Columns 737 through 759
236 237 238 239 240 241 242 243 244 245 246 247 248 249 250
251 252 253 254 255 256 257 258
Columns 760 through 782
259 260 261 262 263 264 265 266 267 268 269 270 271 272 273
274 275 276 277 278 279 280 281
Columns 783 through 805
282 283 284 285 286 287 288 289 290 291 292 293 294 295 296
297 298 299 300 301 302 303 304
Columns 806 through 828
305 306 307 308 309 310 311 312 313 314 315 316 317 318 319
320 321 322 323 324 325 326 327
Columns 829 through 851
328 329 330 331 332 333 334 335 336 337 338 339 340 341 342 343
344 345 346 347 348 349 350
Columns 852 through 874
351 352 353 354 355 356 357 358 359 360 361 362 363 364 365
366 367 368 369 370 371 372 373
Columns 875 through 897
374 375 376 377 378 379 380 381 382 383 384 385 386 387 388
389 390 391 392 393 394 395 396
Columns 898 through 920
15
397 398 399 400 401 402 403 404 405 406 407 408 409 410 411
412 413 414 415 416 417 418 419
Columns 921 through 943
420 421 422 423 424 425 426 427 428 429 430 431 432 433 434
435 436 437 438 439 440 441 442
Columns 944 through 966
443 444 445 446 447 448 449 450 451 452 453 454 455 456 457
458 459 460 461 462 463 464 465
Columns 967 through 989
466 467 468 469 470 471 472 473 474 475 476 477 478 479 480
481 482 483 484 485 486 487 488
Columns 990 through 1 001
489 490 491 492 493 494 495 496 497 498 499 500
v3 =
0 0.3491 0.6981 1.0472 1.3963 1.7453 2.0944 2.4435 2.7925
3.1416
Interprétation des résultats :
La fonction linspace permet de définir un vecteur de longeur
x .
➢ Exercice 10:
clear all
clc
% Question 1:
n=0:10;
v=2.^n;
disp(v)
% Question 2:
R=input('Entrer un réel');
n=0;
while (2^n)<=R
n=n+1
end
16
disp('le plus petit entier')
disp(n)
% Question 3:
l=input('entrer le nomber de linge :');
c=input('entrer le nomber de colonne :');
A=ones(l,c);
for i=1;l
for j=1;c
A(i,j)=1/(i+j-1);
end
end
disp(A)
Columns 1 through 10
1 2 4 8 16 32 64 128 256 512
1024
Entrer un réel10
le plus petit entier
4
entrer le nomber de linge :4
entrer le nomber de colonne :4
1 1 1 1
1 1 1 1
1 1 1 1
1 1 1 1
17
Après l'initiation à MATLAB, on
doit écrire des programmes plus
difficiles pour vérifier les
taches déffirent. Dans cette 2e
séance en va créer des scripts
pour trouver des solutions à des
problèmes simples .
18
➢ Exercice 1:
clear all
clc
x=input('entrer la valeur de x :');
y=input('entrer la valeur de y :');
z=input('entrer la valeur de z :');
if x>y && x>z
max=x;
elseif y>x && y>z
max=y;
else
max=z;
end
disp('la plus grand :')
disp(max)
entrer la valeur de x :6
entrer la valeur de y :9
entrer la valeur de z :1
la plus grand : 9
Interprétation des résultats :
Dans ce exercice en va écrire un programme pour trouver les nombres la
plus grand a couse de la instruction if .
➢ Exercice 2:
clear all
clc
x=input('la valeur de vecteur : ');
S=rand(1,x)
max=S(1);
for i=2:x
if S(i)>max
max=S(i);
end
end
disp('le plus grand élément :')
disp(max)
19
la valeur de vecteur : 7
S = 0.6160 0.4733 0.3517 0.8308 0.5853 0.5497 0.9172
le plus grand élément : 0.9172
Interprétation des résultats :
Dans ce exercice En utiliser la instruction for et if pour faire de test à
chaque élément d’un vecteur .
➢ Exercice 3:
clear all
clc
x=input('la valeur de de vecteur : ');
S=rand(1,x)
min=S(1);
for i=2:x
if S(i)<min
min=S(i);
end
end
disp('la plus petit élément : ')
disp(min);
la valeur de de vecteur : 5
S = 0.0759 0.0540 0.5308 0.7792 0.9340
la plus petit élément : 0.0540
➢ Exercice 4:
clear all
clc
L=input('entrer une valeur :');
E=rand(1,L)
s=0;
n=0;
20
for i=1:L
s=s+E(i);
end
disp('la somme des éléments d un vecteur')
disp(s)
entrer une valeur :4
E = 0.1524 0.8258 0.5383 0.9961
la somme des éléments d un vecteur 2.5127
Interprétation des résultats :
A cause de la instruction for en faire la somme des éléments d’un
vecteur.
➢ Exercice 6:
clear all
clc
L=input('entrer une valeur :');
E=rand(1,L)
s=0;
N=0;
for i=1:L
s=s+E(i);
end
disp('la somme des éléments d un vecteur')
disp(s)
N=s/L;
disp('la moyenne des éléments d un vecteur')
disp(N)
entrer une valeur :5
E = 0.2638 0.1455 0.1361 0.8693 0.5797
la somme des éléments d un vecteur 1.9944
la moyenne des éléments d un vecteur 0.3989
21
➢ Exercice 7:
clear all
clc
c10=0;
u=rand(1,30);
n=20*u
for i=1:30;
if n(i)>=10
c10=c10+1;
end;
end;
disp('le nombre des notes superieurs à 10 est : ')
disp(c10)
s=0;
disp('la moyenne des notes est : ')
m=mean(n)
for i=1:30;
if n(i)>=m
s=s+1;
end;
end;
disp('le nombre des notes superieurs à la moyenne est : ')
disp(s)
n =
Columns 1 through 13
11.5042 1.1956 4.6956 7.0632 16.4239
0.3081 0.8605 3.3798 12.9823 14.6344 12.9549
9.0185 10.9402
Columns 14 through 26
5.9264 14.8939 3.7791 13.7355 3.6702
7.3697 12.5124 15.6045 1.6225 18.5877 15.5143
9.7358 8.7172
Columns 27 through 30
8.9357 6.1270 10.1702 10.2154
le nombre des notes superieurs à 10 est :
14
la moyenne des notes est :
m =
9.1026
le nombre des notes superieurs à la moyenne est :
22
15
➢ Exercice 8:
clear all
clc
disp('1 en série 2 en parallèle')
R1=input('résistance 1 :');
R2=input('résistance 2 :');
R3=input('résistance 3 :');
A=input('le type de montage :');
if A==1
S=R1+R2+R3 ;
elseif A==2
S=1/R1+1/R2+1/R3 ;
Rt=1/S;
end
disp('résistance équivalence :')
disp(Rt)
1 en série 2 en parallèle 1 en série 2 en parallèle
résistance 1 :10 résistance 1 :1000
résistance 2 :20 résistance 2 :1000
résistance 3 :30 résistance 3 :1000
le type de montage :1 le type de montage :2
résistance équivalence : résistance équivalence :
60 333.3333
Interprétation des résultats :
Dans ce exercice en écrire un script pour donner la résistance total. En utiliser
instruction if pour choisi le type de montage .
23
Dans cette séance de TP en va essayer
d'écrire des scriptes Matlab pour
résoudre des fonctions mathématiques
et de tracer la solution et ensuite
tracer la série de Fourier
24
➢ Exercice 1:
clear all
clc
%% Question 1:
f1=@(x)((x.^5-3)./sqrt(x.^2+1))
%% Question 2:
f11=f1(1)
f10=f1(0)
%% Question 3:
x=[-5:0.1:5];
%% Question 4:
f2=((x.^5-3)./sqrt(x.^2+1));
plot(x,f2)
grid on
f1 = function_handle with value:
@(x)((x.^5-3)./sqrt(x.^2+1))
f11 = -1.4142
f10 = -3
Figure 1:
25
➢ Exercice 2:
clear all
clc
R=input('enter R : ');
P=2*pi*R;
disp('le perimetre du cerle de rayon R est: ');
disp(P)
S=P*R/2;
disp('la surface du cerle de rayon R est: ');
disp(S)
enter R : 5
le perimetre du cerle de rayon R est:
31.4159
la surface du cerle de rayon R est:
78.5398
➢ Exercice 3:
clear all
clc
N=input('veuillez entrer un entier N ');
S=0;
P=1;
for i=1:N
b=0; % nbr de division
for j=1:i
if mod(i,j)==0
b=b+1;
if b==2
S=S+i;
P=P*i;
end
if b>2
break
end
end
end
end
disp('la somme des N nombres premiers est: ');
26
disp(S)
disp('le produit des N nombres premiers est: ');
disp(P)
veuillez entrer un entier N 5
la somme des N nombres premiers est:
14
le produit des N nombres premiers est:
120
➢ Exercice 4:
clear all
clc
V=input('Veuillez entrer le volume en m^3 : ');
T=input('Veuillez entrer la temperature en cellesuis : ');
n=input('Veuillez entrer la quantit? de mati?re en mol:
');
t=T+273.15;
R=8.314;
disp(' ')
disp('la pression du gaz en Pa est : ')
P=(n*R*t)/V
Veuillez entrer le volume en m^3 : 4
Veuillez entrer la temperature en cellesuis : 5
Veuillez entrer la quantite de matiere en mol: 5
la pression du gaz en Pa est :
P=
2.8907e+03
27
Méthode Numérique
La méthode de dichotomie :
clear all
clc
R=vpa(sqrt(3),6)
x=0:0.1:5;
f=x.^2-3;
plot(x,f)
a=1;
b=2;
e=10^-6;
while b-a>e
m=(b+a)/2;
if (m.^2-3)*(a.^2-3)<0
b=m;
else
a=m;
end
end
disp(vpa(m),6)
Résultat:
28
R=
1.73205
1.73205089569091796875
Interprétation des résultats :
La méthode de dichotomie est une méthode pour trouver une solution
approchée à une équation .1ére en déclarer la fonction et intervalle et
erreur en suite en vérifier qui erreur inférieur de soustraction de
intervalle en fin en diviser intervalle en 2 .
Méthode de Newton :
clear all
clc
X0=input('entrer une valeur de depart X0');
tg=input('nombre de tg');
n=0;
x=[0:0.1:20];
f=x.^2-3;
e=10^-6;
plot(x,f)
title('courbe de x^2-3');
xlabel('x');
ylabel('f ( x )');
hold on
grid on
while (X0^2-3)>e && n<tg
X1=vpa(X0-(X0^2-3)/(2*X0),6);
y=2.*X0.*(x-X0)+X0.^2-3;
plot(x,y)
X0=X1;
n=n+1;
if n==tg
break
end
end
disp('____________')
disp('la valeur approche a sqrt(3) est: ')
disp(X1)
hold off
29
m=abs(vpa(sqrt(3)-X1,6))
Résultat :
entrer une valeur de depart X07
nombre de tg3
la valeur approche a sqrt(3) est:
1.79392
m=
0.061870186382503078448280575685203
Interprétation des résultats :
La méthode de newton et une méthode pour trouver numériquement
une approximation précise d’une fonction. Graphiquement, cela revient à
tracer la tangente à la courbe représentative de f et à chercher où elle
coupe l'axe des x. On considère donc la suite récurrente définie par une
valeur u0 proche de la racine
30
Intégration numérique
clc ;
clear;
a=0;
b=pi/2;
n=100;
I1=0;
I2=0;
I3=0;
I4=0;
x=a☹b-a)/n:b;
for i=1:n
I1=vpa(I1 +( x(i+1)- x(i)) * sin(x(i)),5);
I2=vpa(I2 +( x(i+1)- x(i)) * sin(x(i+1)),5);
I3=vpa(I3 +( x(i+1)- x(i)) * sin((x(i)+x(i+1))/2),5);
I4=vpa(I4 +( x(i+1)- x(i)) *
(sin(x(i))+sin(x(i+1)))/2,5);
end
disp(‘Methode des rectangles gauche :’) ;
disp(I1)
disp(‘Methode des rectangles droite :’) ;
disp(I2)
disp(‘Methode de point Milieu :’) ;
disp(I3)
disp(‘Methode des trapezes :’) ;
disp(I4)
Résultat :
Méthode des rectangles gauche :
0.99213
Méthode des rectangles droite :
1.0078
Méthode de point Milieu :
1.0
Méthode des trapezes :
0.99998
31
Interprétation des résultats :
Intégrale dans Matlab est plus facile .il y a 4 méthodes d’intégrale dans
Matlab. Méthode des rectangles gauche, Méthode des rectangles droite,
Méthode de point milieu et Méthode des trapez. La méthode la plus
efficace est la méthode de point milieu parce que il donne la valeur la plus
proche de intégrale.
32
Série de Fourier
scripte :
clear all
clc;
f=50;
T=1/f;
E=12; %amplitude du signale
N=100;% Le nombre de pas dans une periode T, pour
step=T/N;
syms t
x=0:step:T;
f1=(1/T)*int(E,t,0,T/2);
f2=(1/T)*int(-E,t,T/2,T);
a0=f1+f2;
k=5;
for n=1:k
f_b3=(2/T)*int(E*sin((2*pi*n*t)/T),t,0,T/2);
f_b4=(2/T)*int(-E*sin(2*pi*n*t/T),t,T/2,T);
bn=f_b3 + f_b4;
b(n)=double(bn);
end
for n=1:k
f_a3=(2/T)*int(E*cos((2*pi*n*t)/T),t,0,T/2);
f_a4=(2/T)*int(-E*cos(2*pi*n*t/T),t,T/2,T);
an = f_a3 + f_a4;
a(n)=double(an);
end
y=a0;
figure;
for n=1:k
y = a(n)*cos(2*pi*n/T*x) + b(n)*sin(2*pi*n/T*x);
plot(x,y);
hold all;
end
hold off
X=[Link]k*f;
Z= [a0, b];
figure;
bar(X,Z, 0.2)
33
Résultat :
Affichage du signale carré avec 100 harmoniques :
34
Ainsi, Matlab est un logiciel qui offre à l’utilisateur tant d’avantages
surtout au niveau du calcul matriciel, la manipulation des fonctions et
le traçage des courbes mathématiques.
En utiliser Matlab est un programme simple au niveau de résoudre
les fonctions mathématiques et l’intégrale.
En fin Matlab programme plus utiliser dans les industries parce que
rapide et simple d’utiliser est efficace.
35
