Exemple de Programme Matlab
Transféré par
ScribdTranslationsExemple de Programme Matlab
Transféré par
ScribdTranslationsExemples de base et mise en œuvre
Utilisation du programme MatLab
% 1. Exercice de création de graphiques
tout effacer;
format court;
x=[0 1 2 3 4 5];
y=[0 20 60 68 77 110];
tracer(x,y,'b');
titre('Exercice de Création de Graphiques');
xlabel('nombre entier');
nombre de dizaines
grille
% 2. Programme graphique de la rotation du moteur à induction
format court;
effacer tout;
Le texte fourni ne contient pas de phrases à traduire, mais plutôt un exemple de données ou une déclaration incomplète.
45, 50, 55, 60];
y=[0, 0.46, 0.91, 1.39, 1.88, 2.36, 2.89, 3.38,...
3.87, 4.36, 4.86, 5.35,5.86];
axe([0,60,0,6]);
tracer(x,y,'r-');
titre('Graphique Fréquence vs Tension tacho');
xlabel('Fréquence (Hz)');
ylabel('Tension du tachymètre (volt)');
x:fréquence
tension
grille;
% 3. Onde de phase différente
format court;
effacer tout;
wt=0:0.05:3*pi;
v1=10*sin(wt);
v2=10*sin(wt-120);
v3=10*sin(wt+240);
tracer(wt,v1,'r',wt,v2,'y',wt,v3,'b');
Graphique de décalage de phase
xlabel('Angle (degrés)');
ylabel('Tension (volt)');
grille
% 4. Graphique d'onde triphasée
format court;
tout effacer;
wt=(0:0.05:3*pi);
v1=10*sin(wt);
v2=10*sin(wt-180);
v3=10*sin(wt-90);
tracer(wt,v1,'r',wt,v2,'g',wt,v3,'b');
Titre('Graphique de sortie d'onde à trois phases');
xlabel('wt (drajat)');
ylabel('Tension (volt)');
text(1.2,8.5,'V1');
V2
texte(4.3,8.5,'V3');
grille
1
% 5. Onde du Système Triphasé
format court;
effacer tout;
V=12;
Vmag=abs(V);
theta=linspace(0,6*pi);
Va=12*sin(thêta+120);
Vb=12*sin(theta+240);
Vc=12*sin(theta+0);
% Créer un graphique d'onde
sdt=theta*180/pi;
tracer(sdt,Va,'r',sdt,Vb,'g',sdt,Vc,'b');
Tension de ligne
xlabel('Angle (degrés)');
ylabel('Tension (kV)');
texte(30,12.5,'Va');
text(395,12.5,'Vb');
texte(810,12.5,'Vc');
grille
% 6. Graphique d'onde de différence de phase
format court;
effacer tout;
wt=(0:0,05:2*pi);
v1=1*sin(wt);
v2=1*-cos(wt);
tracer(wt,v1,'r',wt,v2,'b');
title('Graphique de décalage de phase');
xlabel('wt (radian)');
ylabel('Tension (volt)');
grille
% 7. matrice e et y
effacer tout;
format court;
A=[8 10.3426 27.5197
10.3426 27.5197 96.1149
27.5197 96.1149140 375.5433];
b=[9.4; 12.4276; 31.2741];
x=inv(A)*b
% 8. Exercice de création d'un programme de nombres complexes
format court;
effacer tout;
% 9. De polaire à rectangulaire
z=5*(cos(53.13*pi/180)+sin(53.13*pi/180*j))
% 10. Méthode de division des nombres complexes
i=(1+0j)/(0.5j+0.2j)
% 11. De Rectangulaire à Polaire
zmag=|3+4j|
% 12. Multiplication of Polar Numbers
z1=abs(5);
z2=abs(2);
z=z1*z2
z=z1*z2,(angle(4+3j)*180/pi+angle(8+6j)*180/pi)
2
% 13. Exercice de création de graphiques
effacer tout;
format court;
x=[0 1 2 3 4 5];
y=[0 20 60 68 77 110];
tracer(x,y,'b');
titre('Exercice de création de graphique');
xlabel('nombre entier');
ylabel('nombre de dizaines');
grille
% 14. Remplissage du Capaciteur
format court;
tout effacer;
disp(' ')
disp('_____________________________________________')
disp(' ')
disp(' CALCULER LE COURANT DE TRANSITION ')
disp(' AVEC LA MÉTHODE DE RÉSOLUTION ')
disp(' D'UNE ÉQUATION DIFFÉRENTIELLE ')
disp(' ORDE PREMIER ')
disp('_____________________________________________')
afficher(' ')
epsilon=0.0001;
x=1;
i=input ('Courant initial = ');
t=input ('Heure de début = ');
h=input('Prix h = ');
disp('_____________________________________________')
disp(' ')
disp(['Iter t je abs(i1 -i) '])
disp('_____________________________________________')
disp(' ')
format court g
iter=0;
tant que x >= epsilon
iter=iter + 1;
i1=2*(1-exp(-5*t));
x=abs(i1-i);
i=i1;
t=t+h;
fprintf('%i',iter),disp([t, i, x])
fin
disp('_____________________________________________')
3
% 15. Comparaison
format court
effacer tout
num=[4.5];
den1=[0.35 1];
t=0:0.05:2.5;
Réponse plante
[yp,x1,t]=étape(num,den1,t);
den2=[0.1 1];
[ym,x2,t]=étape(num,den2,t);
n=50;
ys=4.5;
ys
pour i=2:51
e(i)=ys-ym(i);%ym ou yp
de(i)=e(i)-e(i-1);
u(i)=ym(i)/yp(i)*4.5;
fin;
(1:51)' t' e' de' u'
tracer(t,yp,'r.-',t,yp,'b-',t,ym,'.',t,ym,'r-')
Modèle de réponse
xlabel('Temps (seconde)')
ylabel('Signal de contrôle (sortie)')
Comparaison de la réponse du modèle avec la réponse de la plante
%yp:plante
%ym:modèle
grille
% 16. Tension de travail du diode
format court;
effacer tout;
disp(' ')
disp('_____________________________________________')
disp(' ')
disp(' CALCULER LA TENSION DE MARCHE D'UN DIODE ')
disp(' AVEC LA MÉTHODE DIRECTE ')
disp('_____________________________________________')
disp(' ')
0,0001
x=1;
Est=input ('Est= ');
R= input ('R= ');
Vb=
v=input ('Tension initiale Tension v= ');
disp('_____________________________________________')
disp(' ')
afficher('Itération Tension f(v) abs(v1 - v) ')
disp(' ke : Fonctionnement d'une diode ')
disp(' (volt) ')
disp('_____________________________________________')
disp(' ')
format court g
iter=0;
tant que x >= epsilon
iter=iter + 1;
fv=R*Is*(exp(40*v)-1)+v-Vb;
v1=(log(((Vb-v)/(R*Is))+1))/40;
x=abs(v1-v);
v=v1;
fprintf('%i',iter),disp([v1, fv, x])
4
fin
disp('_____________________________________________')
% 17. Résolution avec Matlab;
format court
tout effacer;
1,15 %réactance synchrone
Xe=0.15; %rektansi équivalente
1,2 %daya sortie
V=1.0; %tension de bus
%Réactance totale
Xt=Xs+Xe;
A. Calculez la tension d'excitation
%D'après l'équation Pm=E*V/Xt
E=Pm*Xt/V
%B. Si la puissance de sortie diminue de 0,7 pu
P=0,7;
%Nouveau flux et décalage d'angle
D'après l'équation P=Pm*sin(theta)
theta=asin(P/Pm);
theta*180/pi
E*(cos(theta)+i*sin(theta));
I=(E_kompleks-V)/Xt*i;
abs(I)
eta=atan(imag(I)/real(I));
eta * 180 / pi
Pergeseran_sudut_I=eta*180/pi
% 18. Exemple de question
Deux résistances R1 = 10 ohms, R2 = 20 ohms sont connectées en série à une source de 24
%volt. Calculez :
%a. Tahanan total (Rt)
%b. Courant de circuit (I)
%a. Énergie absorbée (P)
%Résolution
formatcourt;
effacer tout;
10
R2=20;
Vs=24;
Rt=R1+R2
I=Vs/Rt
P=Vs*I
P=I^2*Rt
% 19. Exemple de sous-graphe
formatcourt;
effacer tout;
x=linspace(0,2*pi,30);
y=sin(x);
z=cos(x);
a=2*sin(x).*cos(x);
b=sin(x)./(cos(x)+eps);
sous-graphe(2,2,1) %Configurer le sous-graphique supérieur gauche de l'array de sous-graphiques 2x2
tracer(x,y), axe([0 2*pi -1 1]), titre('sin(x)')
grille
sous-graphique(2,2,2) %Configurer le sous-graphique en haut à droite de 4 sous-graphiques
tracer(x,z), axe([0 2*pi -1 1]), titre('cos(x)')
grille
5
sous-graphe(2,2,3) %Configurer le sous-graphique inférieur gauche de 4 sous-graphes
tracer(x,a), axe([0 2*pi -1 1]), titre('2sin(x)cos(x)')
grille
sous-graphique(2,2,4) Configurer le sous-graphique en bas à droite de 4 sous-graphiques
tracer(x,b), axe([0 2*pi -1 1]), titre('sin(x)/cos(x)')
grille;
Canal de transmission
% 20. Exemple de Canal Court
formatcourt
tout effacer;
% ==============================================================
%
% TRANSMISSION ÉLECTRIQUE CANAL COURT
%
% ==============================================================
Vrll= 23000;
f= 60;
Z= 2,48+6,57i;
magZ= abs(Z);
sudutZ= angle(Z)*180/pi;
9000000;
0,85
acospi = acos(kospi) * 180 / pi;
sinpi = sin(acospi * pi / 180)
%Solution
Tension ligne à neutre
Affichage('Tension Ligne au Neutre :')
Vrln= Vrll/sqrt(3);
magVrln= abs(Vrln);
angle(Vrln) * 180 / pi;
disp(['Vrln = ',num2str(Vrln),' volt']);
fprintf('Tension Ln = %g<%g \n', magVrln,sudutVrln);
%Ligne Arus
Affichage('Ligne de courant :')
Iln=(Charge/(sqrt(3)*Vrll*kospi))*(kospi-sinpi*i);
magIln = abs(Iln);
sudutIln=angle(Iln)*180/pi;
disp(['Iln = ',num2str(Iln),' A']);
fprintf('Arus Ln = %g<%g ', magIln,sudutIln);
%Chute de tension
Afficher 'Chute de Tension :'
Iln*Z;
abs(Vdrop);
angle(Vdrop)*180/pi;
disp(['Vdrop = ',num2str(Vdrop),' volt']);
fprintf('Baisse de Tension = %g<%g \n', magVdrop,sudutVdrop);
%Tension du côté d'envoi
Tension ligne à neutre
Afficher('Tension Côté Envoyer :')
Affichage('Tension Ligne à Neutre :')
Vsln= Vrln+Vdrop;
Forme rectangulaire
magVsln= abs(Vsln);
angle(Vsln)*180/pi;
disp('Tension du côté d'envoi')
disp(['Vsln = ',num2str(Vsln),' volts']);
6
fprintf('Tension Ln = %g<%g \n', magVsln,sudutVsln);
Tension entre les lignes
disp('Tension Ligne à Ligne :')
Vsll = racine(3) * Vsln;
forme rectangulaire
abs(Vsll);
angle(Vsll)*180/pi;
disp(['Vsll = ',num2str(Vsll),' volts']);
fprintf('Tension ll = %g<%g\n', magVsll,sudutVsll);
%Angle de phase
disp('Angle de Phase :')
sudutVsll
disp(['teta = ',num2str(teta),'']);
