ORTAGONAL MATRIS

Yazan ve Cevaplayan:Fatih KABAKCI

Bu soru Chip Online sayfasında 30.12.2010 tarihinde Matella adlı kullanıcı tarafından sorulmuştur.

Gelen Soru: 3.Dışarıdan girilen bir matrisin bir ortogonal matris olup olmadığını bulan c prograamını yazınız.

Verilen Cevap:

M11 M12 M13       *        M11 M21 M31       =      1  0  0
M21 M22 M23                M12 M22 M32              0  1  0 
M31 M32 M33                M13 M23 M33              0  0  1

#include <stdio.h>
#define row 2
#define field 2

int birimMatris(int sonucMatris[][field]) {

int flag=1,i,j;
for(i=0;i<row;i++) {
for(j=0;j<field;j++) {
if(i==j && sonucMatris[i][j]!=1) 
flag=0;
else if(i!=j && sonucMatris[i][j]!=0) 
flag=0;
}
if(!flag)
break;
}
if(flag)
printf("Matris ortogonal bir matristir\n");
else
printf("Matris ortogonal bir matris degildir\n");
}
void tranzpozAl(int tranzpozMatris[][field],int ilkMatris[][field]) {

int i,j;
for(i=0;i<row;i++)
for(j=0;j<field;j++)
tranzpozMatris[i][j]=ilkMatris[j][i]; // transpose alinir
}

void matrisCarp(int resultMatrix[][field],int ilkMatris[][field],int transpozMatris[][field]) {
int i,j,k;
for(i=0;i<row;i++)
for(j=0;j<field;j++)
for(k=0;k<2;k++) 
resultMatrix[i][j]+=ilkMatris[i][k]*transpozMatris[k][j]; // iki matris carpilir
}

int main()  {

int sonucMatris[row][field]={0,0};
int ilkMatris[row][field]={{1,0},{0,-1}};
int transpozMatris[field][row];
tranzpozAl(transpozMatris,ilkMatris);
matrisCarp(sonucMatris,ilkMatris,transpozMatris);
birimMatris(sonucMatris);
getchar();
return 0;
}