가우스 조던으로 짠 31차 역행렬
희나리
고조파의 웨이브폼 데이터를 집어넣어 차수별 고조파의 크기를 나타내도록 만들었습니다.
역행렬만드는게 젤 어려웠는데 가우스 조던으로 해결했네여.
#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;
}
번호 | 제 목 | 글쓴이 | 날짜 |
---|---|---|---|
2655471 | 링 버퍼 구현된 소스 인데요.. 이해가 잘안되서요.. | 소심한녀자 | 2024-05-20 |
2655440 | test파일이 안 되요. (2) | 리라 | 2024-05-19 |
2655416 | 고수님들 도움이 절실합니다. (2) | 세리 | 2024-05-19 |
2655361 | 엔터 입력받을때까지 기다리는 getchar (1) | 큰모음 | 2024-05-19 |
2655328 | 함수에서 포인터 (3) | 모아 | 2024-05-18 |
2655299 | 아이콘바꾸는법..? | 소 | 2024-05-18 |
2655271 | scanf 입력받으면 쭉 실행되는.. (5) | ComeOn | 2024-05-18 |
2655184 | 이 세개 프로그램좀 해결해주실.... 고수님들???? (3) | 말달리자 | 2024-05-17 |
2655129 | 다차원 배열 원소들 주소계산하는거 질문이요~ | 터전 | 2024-05-17 |
2655103 | AVRstudio MoveAxis. MoveLine 함수질문 | 앨버트 | 2024-05-16 |
2655072 | 콘솔창의 위치와 크기를 정하고 싶습니다. | 딥레드 | 2024-05-16 |
2655044 | 자료구조랑 online judge site에 대해서 질문있습니다... (1) | 뽀송뽀송 | 2024-05-16 |
2655016 | 별찍기 문제 중 이해안되는 부분이 있어요... | 도전중 | 2024-05-16 |
2654986 | 입력받은 성적을 배열에 저장하기 (3) | 안녕나야 | 2024-05-15 |
2654959 | ===========-이런거 미리 지정하는법 (1) | 김자영 | 2024-05-15 |
2654930 | 도와주십시요 ㅜㅜ 자료구조 ㅜㅜ | 도널드 | 2024-05-15 |
2654881 | PID 모듈 검색 부분 질문합니다 . | 이름없음 | 2024-05-14 |
2654850 | 조건에 만족하는 다음 수를 찾고 싶습니다. | 초엘 | 2024-05-14 |
2654790 | 두파일간 byte단위 txt파일 복사 中 | 슬찬 | 2024-05-14 |
2654768 | 사전만들기 입니다. | 여우By | 2024-05-13 |