Pivot de gauss partiel en langage c

Fermé
Leila - 6 mai 2004 à 18:38
 Wolf - 15 mai 2010 à 17:18
je souhaite avoir la programmation du pivot de gauss partiel en langage c. C'est une méthode de résolution des matrices
merci d'avance

11 réponses

si tu as trouvé la resolution de systeme d'equation par le pivaot de gauss veux tu bien me l'envoyer a mon mail
merci.
5
j'attend vos reponses..et merci!!
1
juanpablo Messages postés 142 Date d'inscription mercredi 30 mai 2007 Statut Membre Dernière intervention 16 juin 2008 10
10 janv. 2008 à 23:52
salut

exemple de mise en oeuvre de la methode de gauss . attention les indices des tableaux vont de 1 à n (comme en math)

a : la matrice du système b : le second membre du système n : l'ordre de la matrice.

renvoie la valeur 1 si la matrice est numériquement inversible et la valeur 0 sinon ;
quand la valeur renvoyée est 1, la solution du système se trouve dans le tableau b.

NMAX est égale à la dimension maximale de la matrice +1 ; N2MAX =NMAX+1

int sl_gauss(double a[NMAX][NMAX],double b[NMAX],int n)
{
int i,j,k,l,err;
double max,pivot,coef,s;
double t[NMAX][N2MAX];
for(i=1;i<=n;i++)
{
for(j=1;j<=n;j++)
t[i][j]=a[i][j];
t[i][n+1]=b[i];
}
err=1;
k=1;
while (err==1 && k<n)
{
max=fabs(t[k][k]);
l=k;
for(i=k+1;i<=n;i++)
if(max<fabs(t[i][k]))
{
max=fabs(t[i][k]);
l=i;
}
if(max!=0)
{
if(l!=k)
for(j=k;j<=n+1;j++)
{
pivot=t[k][j];
t[k][j]=t[l][j];
t[l][j]=pivot;
}
pivot=t[k][k];
for(i=k+1;i<=n;i++)
{
coef=t[i][k]/pivot;
for(j=k+1;j<=n+1;j++)
t[i][j] -= coef*t[k][j];
}
}
else err=0;
k++;
}
if(t[n][n]==0) err=0;
if(err==1)
{
b[n]=t[n][n+1]/t[n][n];
for(i=n-1;i>=1;i--)
{
s=t[i][n+1];
for(j=i+1;j<=n;j++)
s-=t[i][j]*b[j];
b[i]=s/t[i][i];
}
}
return(err);
}
//------------------------------------------------------------------------------
0
je vs remercie infiniment pour votre aide ...merci!!
1
je souhaite avoir la programmation du pivot de gauss partiel en langage c. C'est une méthode de résolution des matrices
merci d'avance.
1

Vous n’avez pas trouvé la réponse que vous recherchez ?

Posez votre question
Amoureux3 Messages postés 940 Date d'inscription vendredi 9 avril 2004 Statut Membre Dernière intervention 5 mai 2007 9
18 mai 2004 à 01:11
suffit de déterminer les étapes de la solution mathématique avec ses cas particuliers (exemple le calcul du déterminant), puis préciser les structures que tu dois utiliser (tableau 2d, matrice, variables, indices...)

te reste ensuite que le plus simple, l'interprétation de la solution mathématique en programme informatique..

[ Ton Doux Visage ]
0
Merci pour ta reponce
mais le probleme et que je ne c pas le resoudre mathematiquement
je c qu'il y a le choix du pivot mais comment le determiner?
la trangulation de la matrice mais qu'elle sont les etapes?
et enfin la resolution.

en realité mon projet est a faire ezn ADA et donc si j'avais un algo ou un cour de maths assez bien expliqué je commencerai sans pb.

merci.
0
bonjour s.v.p si tu as trouvé la resolution de systeme d'equation par la methode de gauss si le pivot est null veux tu bien me l'envoyer a mon mail
merci.
0
JSS Messages postés 3745 Date d'inscription jeudi 22 avril 2004 Statut Contributeur Dernière intervention 17 décembre 2004 32
18 mai 2004 à 17:31
Bonjour ou Bonsoir,


Recherche google : pivot gauss explication :

http://www.dil.univ-mrs.fr/~garreta/generique/sujets/PB01.pdf

++

JSS
0
je cherche comment effectuer un programme en langage c pour la methode pivot de gauss
0
bonjour juanpablo!j'ai regardé ton programme et je ne comprends pas comment fonctionne ta boucle "tant que" while.est ce que ce serait pour proceder a l'echange entre les equations pour la suite des calculs?et a quoi correspond "err"?
0
Il y'a un problème des pivots dans les système matricielle quelle est la meilleure méthode pour résoudre ce problème
0
Salut, ça fait longtemps que j'ai travaillé la dessus , j'espere que cela t'aidra
bonne chance!!





#include <iostream.h>
int main(){
int n;
double e[11][10];
double s[10];
cout<<"programme du pivot de gauss\nCombien dequations ?\nN= ";
cin>>n;
cout<<"\n";
for (int i=0;i<n;i++){
cout<<"equation "<<i<<"\n";
for (int p=0;p<n;p++){
cout<<"v"<<p<<" = ";
cin>>e[p][i];
}
cout<<"\n";
cout<<"equation "<<i<<" = ";
cin>>e[n][i];
cout<<"\n";
}
// on a saisi les facteurs des equations ds e[][]
int y=0;
double var1=0,var2=0;
double temp;
int a,t;
for(int x=0;x<n-1;x++){
for(a=1+x;a<n;a++){
temp=e[x][a];
for (t=x;t<n+1;t++){
e[t][a]=e[t][a]*e[x][x]-e[t][x]*temp;

/// triangulation du systeme.
}
}
}
//Passons au remplacage..
int af;
s[n-1]=e[n][n-1]/e[n-1][n-1];
e[n][n-1]=0;
e[n-1][n-1]=0;
for (int ligne=1;ligne<=n;ligne++){
for (int sol=2;sol<=n;sol++){
e[n-ligne][n-sol]*=s[n-ligne];
e[n][n-sol]-=e[n-ligne][n-sol];
e[n-ligne][n-sol]=0;

}

s[n-(ligne+1)]=e[n][n-(ligne+1)]/e[n-(ligne+1)][n-(ligne+1)];

}

// Resolution..

// Affichage..
for (af=0;af<n;af++){cout<<" var "<<af<<" = "<<s[af];}
cin>>n;

return 0;}
0