MATRIS DE OKLIT,HILBERT ve MUTLAK NORM

Yazan ve Cevaplayan:Fatih KABAKCI

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

Gelen Soru:Bir matrisin max mutlak norm, hilbert normu ve euclidiean normunu veren programları yazabilir misiniz? yardımlarınız için teşekkürler.
1.=max |aij| mutlak norm
2.= xi burada xi, B=AA^T matrisinin en büyük öz değeridir. hilbert normu
3.=[ijE ai,j'2]'1/2 =[i=1E xii]'1/2 euclidiean normu (E toplam sembolü)

Verilen Cevap:

#include <iostream>
#include <cmath>

using namespace std;

const int SIZE=2;

double mutlakNorm(double matrix[][SIZE])
{
     double max=matrix[0][0];
     int i,j;
     for(i=0;i<SIZE;i++)
     for(j=0;j<SIZE;j++)
     if(max < matrix[i][j])
     max=matrix[i][j];
     return max;
     }

double hilbertNorm(double matrix[][SIZE])
{

double resultMatrix[SIZE][SIZE]={0,0};
int i,j,k;

for(i=0;i<SIZE;i++)
for(j=0;j<SIZE;j++)
for(k=0;k<2;k++)
resultMatrix[i][j]+=matrix[i][k]*matrix[k][j];

for(int i=0;i<SIZE;i++)
for(int j=0;j<SIZE;j++)
matrix[i][j]=resultMatrix[j][i];

double b=matrix[0][0]+matrix[1][1]; // warning:! only for two-two matrix
double c=matrix[0][0]*matrix[1][1]-matrix[1][0]*matrix[0][1];
double delta=pow(b,2)-4*1*c;
double root1=(-b+sqrt(delta))/2;
double root2=(-b-sqrt(delta))/2;
if(root1 > root2)
return root1;
return root2;
}          

double euclidNorm(double matrix[][SIZE])
{
       double total=0.0;
       int i,j;
       for(i=0;i<SIZE;i++)
       for(j=0;j<SIZE;j++)
       total+=pow(matrix[i][j],2);
       return sqrt(total);
       }