Добрый день.
Подскажите в чем проблема? Вариант для поверхности работает, а что здесь - не знаю.
Программа немного недоработана по алгоритму, но ошибка не в этом - просто не компилируется.
| Код |
//=========================================================================================================================
//--------------------------------------- Программа моделирования протекания токов ----------------------------------------
//=========================================================================================================================
#include <stdio.h>
#include <conio.h>
#include <math.h>
int main()
{
FILE *outfile;
outfile=fopen("1.txt","w");
short i,j,m,k,s,r,n,si,rj,h,hm,ii;
float Ps[15][15][31],Pr[15][16][30],Ph[16][15][30],s1,G1,I1;
float** G=new float*[6750];
for(i=0;i<6750;i++)
{
G[i]=new float[6750];
}
float* Fi=new float[6750];
float* F=new float[6750];
float* I=new float[6750];
float* S=new float[6750];
//Ввод параметров на моделирование------------------------------------------------------------------------------------------
s=15; //размеры тела
r=30;
h=15;
k=s*r*h-1;
for(i=1; i<=k; i++) //места подключения токовых присосок
I[i]=0;
I[3901]=1;
I[4000]=-1;
for(m=1; m<=h; m++)
for(i=1; i<=s; i++)
for(j=0; j<=r; j++)
Ps[m][i][j]=0;
for(m=1; m<=h; m++)
for(i=0; i<=s; i++)
for(j=1; j<=r; j++)
Pr[m][i][j]=0;
for(m=0; m<=h; m++)
for(i=1; i<=s; i++)
for(j=1; j<=r; j++)
Ph[m][i][j]=0;
for(m=1; m<=h; m++) //задание характеристики тела
for(i=1; i<=s; i++)
for(j=1; j<=r-1; j++)
Ps[m][i][j]=1;
for(m=1; m<=h; m++)
for(i=1; i<=s; i++)
for(j=1; j<=r-1; j++)
Pr[m][i][j]=1;
for(m=1; m<=h-1; m++)
for(i=1; i<=s; i++)
for(j=1; j<=r; j++)
Ph[m][i][j]=0.1;
//=========================================================================================================================
//Формирование матрицы проводимостей G-------------------------------------------------------------------------------------
for(i=1; i<=k; i++)
for(j=1; j<=k; j++)
G[i][j]=0;
for(i=1; i<=k; i++)
for(j=1; j<=k; j++)
if(i==j)
{
hm=(i-i%(r*s))/r*s;
if(hm==0)
hm=1;
ii=i-(hm-1)*r*s;
rj=ii%r;
if(rj==0)
rj=r;
if(ii==rj)
si=1;
else
si=(ii-rj)/r+1;
G[i][j]=Pr[hm][si-1][rj]+Pr[hm][si][rj]+Ps[hm][si][rj-1]+Ps[hm][si][rj]+Ph[hm][si][rj]+Ph[hm-1][si][rj];
}
/* for(i=1; i<=k; i++)
for(j=1; j<=k; j++)
{
if(j==i+1)
{
rj=i%r;
if(rj==0)
rj=r;
if(i==rj)
si=1;
else
si=(i-rj)/r+1;
G[i][j]=-Ps[hm][si][rj];
}
if(j==i-1)
{
rj=i%r;
if(rj==0)
rj=r;
if(i==rj)
si=1;
else
si=(i-rj)/r+1;
G[i][j]=-Ps[hm][si][rj-1];
}
if(j==i+r)
{
rj=i%r;
if(rj==0)
rj=r;
if(i==rj)
si=1;
else
si=(i-rj)/r+1;
G[i][j]=-Pr[hm][si][rj];
}
if(j==i-r)
{
rj=i%r;
if(rj==0)
rj=r;
if(i==rj)
si=1;
else
si=(i-rj)/r+1;
G[i][j]=-Pr[hm][si-1][rj];
}
if(j==i+1)
{
rj=i%r;
if(rj==0)
rj=r;
if(i==rj)
si=1;
else
si=(i-rj)/r+1;
G[i][j]=-Ph[hm][si][rj];
}
if(j==i+1)
{
rj=i%r;
if(rj==0)
rj=r;
if(i==rj)
si=1;
else
si=(i-rj)/r+1;
G[i][j]=-Ph[hm][si][rj];
}
} */
//==========================================================================================================================
// for(i=1; i<=k; i++)
// {
// printf("\n");
// for(j=1; j<=k; j++)
// printf("% 2.0f",G[i][j]);
// }
//Приведение матрицы G к треугольному виду----------------------------------------------------------------------------------
for(i=1; i<=k; i++)
{
for(j=1; j<=k; j++)
S[j]=G[i][j];
s1=S[i];
I1=I[i];
for(j=i+1; j<=k; j++)
if (G[j][i]!=0)
{
G1=G[j][i];
for(n=i; n<=k; n++)
G[j][n]=G[j][n]-S[n]*G1/s1;
I[j]=I[j]-I1*G1/s1;
}
}
//==========================================================================================================================
//Вычисление потенциалов----------------------------------------------------------------------------------------------------
Fi[k]=I[k]/G[k][k];
for(i=k-1; i>=1; i--)
{
F[i]=0;
for(j=i+1; j<=k; j++)
F[i]=F[i]+Fi[j]*G[i][j];
Fi[i]=(I[i]-F[i])/G[i][i];
}
//==========================================================================================================================
//Вывод массива потенциалов в файл------------------------------------------------------------------------------------------
for(i=1; i<=k; i++)
{
printf("\n");
// for(j=1; j<=k; j++)
// printf("% 3.2f",G[i][j]);
printf(" % 3.5f",I[i]);
}
printf("\n"); printf("\n");
Fi[k+1]=0;
n=1;
for(i=1; i<=s; i++)
{
for(j=1; j<=r; j++)
{
fprintf(outfile," % 7.6f",Fi[n],outfile);
n=n+1;
}
fprintf(outfile,"\n");
}
//==========================================================================================================================
fclose(outfile);
for(i=0;i<6750;i++)
{
delete [] G[i];
}
delete [] G;
getch();
return 0;
}
|
P.S. Как-то уж очень нецеремонно влезли... |
и особенно JackYF за добавление с примером.