수다닷컴

질문 제목 : 구조체 인자전달질문 요약 :제가 만든 코드에 책에있는 고유값고유벡터 구하는 코드 붙이려고하는대요질문 내용 : 이 책에있는 코드 함수의 매개변수가 무슨뜻인지 이해가 안가네요ㅠㅠ
빨간색이 책에있는 코드고요;; 녹색부분이 이해가안갑니다ㅠ

#include stdio.h
#include math.h

const float epsilon = 1.0e-10f;
const int maxsweeps = 32;
void calculateeigensystem(const matrix3d& m, float *lambda, matrix3d& r);
struct arr getdata(void);

struct matrix3d
{
float n[3][3];
float& operator()(int i,int j)
{
return(n[j][i]);
}
const float& operator()(int i,int j)const
{
return(n[j][i]);
}
void setidentity(void)
{
n[0][0]=n[1][1]=n[2][2]=1.0f;
n[0][1]=n[0][2]=n[1][0]=n[1][2]=
n[2][0]=n[2][1]=0.0f;
}
};
struct arr{
float a[4][3];//벡터 성분변수 선언
};
int main()
{
struct arr s=getdata();
float p[3];//평균위치 m 연산
p[0]=(s.a[0][0]+s.a[1][0]+s.a[2][0]+s.a[3][0])/4;
p[1]=(s.a[0][1]+s.a[1][1]+s.a[2][1]+s.a[3][1])/4;
p[2]=(s.a[0][2]+s.a[1][2]+s.a[2][2]+s.a[3][2])/4;
printf(평균 위치 m=(%.2lf %.2lf %.2lf)\n,p[0],p[1],p[2]);

int i,j;
float c[3][3];
c[0][0]=(pow(s.a[0][0]-p[0],2)+pow(s.a[1][0]-p[0],2)+pow(s.a[2][0]-p[0],2)+pow(s.a[3][0]-p[0],2))/4;//공분산 행렬 연산
c[1][1]=(pow(s.a[0][1]-p[1],2)+pow(s.a[1][1]-p[1],2)+pow(s.a[2][1]-p[1],2)+pow(s.a[3][1]-p[1],2))/4;
c[2][2]=(pow(s.a[0][2]-p[2],2)+pow(s.a[1][2]-p[2],2)+pow(s.a[2][2]-p[2],2)+pow(s.a[3][2]-p[2],2))/4;

c[0][1]=c[1][0]=((s.a[0][0]-p[0])*(s.a[0][1]-p[1])+(s.a[1][0]-p[0])*(s.a[1][1]-p[1])
+(s.a[2][0]-p[0])*(s.a[2][1]-p[1])+(s.a[3][0]-p[0])*(s.a[3][1]-p[1]))/4;

c[0][2]=c[2][0]=((s.a[0][0]-p[0])*(s.a[0][2]-p[2])+(s.a[1][0]-p[0])*(s.a[1][2]-p[2])
+(s.a[2][0]-p[0])*(s.a[2][2]-p[2])+(s.a[3][0]-p[0])*(s.a[3][2]-p[2]))/4;

c[1][2]=c[2][1]=((s.a[0][1]-p[1])*(s.a[0][2]-p[2])+(s.a[1][1]-p[1])*(s.a[1][2]-p[2])
+(s.a[2][1]-p[1])*(s.a[2][2]-p[2])+(s.a[3][1]-p[1])*(s.a[3][2]-p[2]))/4;

printf(공분산 행렬 c는 다음과 같습니다\n);
for(i=0;i3;i++)//공분산 행렬 표시
{
for(j=0;j3;j++)
{
printf(%8.2lf,c[i][j]);
}printf(\n);
}
float lamb[3];
float matrix3d te;
calculateeigensystem(c, lamb, te);------------요부분을 어떻게 해야할까요?
return 0;

}
void calculateeigensystem(const matrix3d& m, float *lambda , matrix3d& r)
{
float m11=m(0,0);
float m12=m(0,1);
float m13=m(0,2);
float m22=m(1,1);
float m23=m(1,2);
float m33=m(2,2);
r.setidentity();
for(int a=0;amaxsweeps;a++)
{
if((fabs(m12)epsilon)&&(fabs(m13)epsilon)&&
(fabs(m23)epsilon)) break;
if(m12 !=0.0f)
{
float u = (m22-m11)*0.5f/m12;
float u2=u*u;
float u2p1=u2+1.0f;
float t= (u2p1 != u2) ?
((u0.0f)?-1.0f:1.0f)*(sqrt(u2p1)-fabs(u))
:0.5f/u;
float c=1.0f/sqrt(t*t+1.0f);
float s=c*t;
m11 -= t*m12;
m22 += t*m12;
m12=0.0f;
float temp = c*m13-s*m23;
m23=s*m13+c*m23;
m13=temp;
for(int i=0; i3; i++)
{
float temp = c*r(i,0)-s*r(i,1);
r(i,1) = s*r(i,0)+c*r(i,1);
r(i,0)=temp;
&n}
}

if(m13 !=0.0f)
{
float u = (m33-m11)*0.5f/m13;
float u2=u*u;
float u2p1=u2+1.0f;
float t= (u2p1 != u2) ?
((u0.0f)?-1.0f:1.0f)*(sqrt(u2p1)-fabs(u))
:0.5f/u;
float c=1.0f/sqrt(t*t+1.0f);
float s=c*t;
m11 -= t*m13;
m33 += t*m13;
m13=0.0f;
float temp = c*m12-s*m23;
m23=s*m12+c*m23;
m12=temp;
for(int i=0; i3; i++)
{
float temp = c*r(i,0)-s*r(i,2);
r(i,2) = s*r(i,0)+c*r(i,2);
r(i,0)=temp;
}
}
if(m23 !=0.0f)
{
float u = (m33-m22)*0.5f/m23;
float u2=u*u;
float u2p1=u2+1.0f;
float t= (u2p1 != u2) ?
((u0.0f)?-1.0f:1.0f)*(sqrt(u2p1)-fabs(u))
:0.5f/u;
float c=1.0f/sqrt(t*t+1.0f);
float s=c*t;
m22 -= t*m23;
m33 += t*m23;
m23=0.0f;
float temp = c*m12-s*m13;
m13=s*m12+c*m13;
m12=temp;
for(int i=0; i3; i++)
{
float temp = c*r(i,1)-s*r(i,2);
r(i,2) = s*r(i,1)+c*r(i,2);
r(i,1)=temp;
}
}
}
lambda[0]=m11;
lambda[1]=m22;
lambda[2]=m33;

int i;
for(i=0;i3;i++)
{
printf(%f,lambda[i]);
}
}
struct arr getdata(void)
{
struct arr t;
int i,j;
for(i=0;i4;i++)
{
printf(%d번째 벡터 \n,i+1);
for(j=0;j3;j++)
{
printf( %d번째 성분입력:,j+1);
scanf_s(%f,&t.a[i][j]);
}
}
return t;
}

신청하기