Вот такая программа:
| Код |
#include <vcl.h>
#pragma hdrstop
#include<iostream.h>
#include<math.h>
int maximal(int n,double R0[]);
void main(){
int i,j,n,f,k,iter;
double S,det;
cout<<"N = ";cin>>n;
double *x=new double [n];
double **b=new double *[n];
for(i=0;i<n;i++)
b[i]=new double[n+1];
double **a=new double *[n];
for(i=0;i<n;i++)
a[i]=new double[n+1];
cout<<"N Interation:";
cin>>iter;
cout<<"Full Matrix:\n";
for(i=0;i<n;i++){
for(j=0;j<=n;j++)
cin>>b[i][j];
}
cout<<"Matrix for relaxion\n";
for(i=0;i<n;i++){
for(j=0;j<n;j++)
a[i][j]=-b[i][j]/b[i][i];
a[i][n]=b[i][n]/b[i][i];
}
for(i=0;i<n;i++){
for(j=0;j<n+1;j++)
cout<<" "<<a[i][j]<<" || ";
cout<<"\n";
}
double *x0=new double [n];
for(i=0;i<n;i++)
x[i]=0.00;
double *R0=new double [n];
cout<<"X1(0)..Xn(0):\n";
for(i=0;i<n;i++)
cin>>x0[i];
S=0.00;
for(i=0;i<n;i++){
for(j=0;j<n;j++)
S=S+a[i][j]*x0[i];
}
for(i=0;i<n;i++){
R0[i]=a[i][n]-x0[i]+S;
cout<<"R("<<i<<")="<<R0[i]<<" | ";
}
f=maximal(n,R0);
det=R0[f];
for(k=0;k<iter;k++){
cout<<"det{"<<k<<"}="<<det<<"\n";
for(i=0;i<n;i++){
if(i!=f) R0[i]=R0[i]+a[i][f]*det;
else R0[i]=R0[i]-det;
}
for(i=0;i<n;i++)
cout<<"R["<<i+1<<"]="<<R0[i]<<" ";
x[f]=x[f]+det;
f=maximal(n,R0);
det=R0[f];
}
cout<<"\n";
for(i=0;i<n;i++)
cout<<"X{"<<i+1<<"}="<<x[i]<<"\n";
cin.sync();
getchar();
delete []x;
delete []R0;
delete []x0;
delete []a;
}
int maximal(int n,double R0[]){
int i,f;
f=0.0;
for(i=0;i<n-1;i++){
if(R0[i+1]>R0[i]) f=i+1;
}
return f;
}
|
Решает пример:
10x1 - 2x2 - 2x3 = 6
-x1 + 10x2 - 2x3 = 7
-x1 - x2 + 10x3 = 8
При семи итерациях навязки обращаются в ноль и получаем три корня {1;1;1}
Считаю в маткаде: тоже самое.
К примеру для этой системы:
4x1 + 3x2 + 3x3 = 4
2.5x1 + 4x2 + 6x3 = 5
6x1 + 2x2 + 4x3 = 7
по маткаду должны быть следующие ответы {0.757;-.581;0.905}
считаю программой: немогу получить нулевые навязки (не считает короче)
В основном я код этой проги понимаю, но не весь (: Может кто подскажет где загвоздка? И как сделать чтобы учитывалась точность при доведении до нуля навязок?
|
)