수다닷컴

고조파의 웨이브폼 데이터를 집어넣어 차수별 고조파의 크기를 나타내도록 만들었습니다.

역행렬만드는게 젤 어려웠는데 가우스 조던으로 해결했네여.

#include stdio.h
#include math.h
#define M_PI 3.14
#define n (31)
double sqr_X[n][n] = {0}; // S_n 입력
double mat[n][2*n] = {0}; // S_n에 정방행렬 붙인 형태
double inv_mat[n][n] = {0}; // S_n의 역행렬

double table1 [n][1] = {0}; //차수별 고조파의 크기 입력 A_n
double table3 [n][1] = {0}; // 차수별 고조파의 크기 출력V_t
double table5 [n][1] = {0}; // A_n result

int main (void)
{

int n;
int i , j, k;
double mul_cont, div_cont;

double fGet;

for ( i = 0 ; i n ; i++ )
{
for ( j = 0 ; j n ; j++ )
{
fGet = sinl (((double)i+1) * (2 * M_PI) * (((double)j+1) / 32));
sqr_X[i][j] = fGet;

}
}
printf(\n 입력한 행렬 S_n = \n);
for ( i = 0 ; i n ; i++ )
{
;
for ( j = 0 ; j n ; j++ )
{
printf(%3.2lf\t, sqr_X[i][j]);
}
printf(\n);
}
printf(\n);

for (i=0;i n;i++)
{
for (j=0;j n;j++)
{
mat[i][j] = sqr_X[i][j];
}
for (j=n;jn*2;j++)
{
mat[i][j] = 0;
if (i+n == j)
mat[i][j]=1;
}
}

for (j=0;j n;j++)
{
div_cont = mat[j][j];
for (k = 0;k n*2; k++)
{
mat[j][k] = mat[j][k] / div_cont;
}
if (j != 0)
{
for (i=0;i j;i++)
{
mul_cont = mat[i][j];
for (k=0;k n*2;k++)
{
mat[i][k] = mat[i][k] - mat[j][k] * mul_cont;
}
}
}

if (j != n-1 )
{
for (i=j+1;i n;i++)
{
mul_cont = mat[i][j];
for (k=0;k n*2;k++)
&nbsr / {
mat[i][k] = mat[i][k] - mat[j][k] * mul_cont;
}
}
}
}

for (i=0;i n;i++)
{
for (j=n;j n*2;j++)
{
inv_mat[i][j-n] = mat[i][j];
}
}
printf(\n 입력한 행렬 S_n의 역행렬 = \n);
for (i =0 ; i n; i++)
{
for (j=0 ; j n; j++)
{
printf(%8.2lf\t , inv_mat[i][j]);
}
printf(\n);
}

printf(\n---------------------------------------------------------------\n);
printf(각 차수별 고조파의 크기를 넣으세여 \n ); /* A_n 값 입력 */
printf(\n);

for( j = 0 ; j n ; j++)
{
for( k = 0 ; k 1 ; k++)
{
printf(%d 차 고조파의 크기 : ,j+1);
scanf(%lf, &table1[j][k]);
}
}
printf(\n);
printf(\n);
printf(**********************\n);
printf( V_t \n);
printf(**********************\n);

for( i = 0 ; i n ; i++)
{
for( j= 0 ; j n ; j++)
{
for( k = 0 ; k 1 ; k++)
{
table3[i][k] += sqr_X[i][j] * table1[j][k];
}
}
}
for( i = 0 ; i n; i++)
{
for( k= 0 ; k 1 ; k++)
{
printf(V_t_%d : %lf \n,i+1, table3[i][k], \n );
}
}
printf(\n);
printf(**********************\n);
printf( A_t result \n);
printf(**********************\n);
for( i = 0 ; i n ; i++)
{
for( j= 0 ; j n ; j++)
{
for( k=0 ; k 1 ; k++)
{
table5[i][k] += inv_mat[i][j] * table3[j][k];
}
}
}
for( i = 0 ; i n; i++)
{
for( k= 0 ; k 1 ; k++)
{
printf(A_t_%d : %lf \n,i+1, table5[i][k], \n );
}
}

return 0;

}

신청하기