Metodat numerike

[suli]

Pershendetje!!
kerkoj ndihme ne shkrimin e nje programi ne C per Metoden e Jacobit per te gjetur "vlerat e veta te matricave". Ku matrica eshte simetrike dhe me vlera reale.
do ju isha sh mirenjohes
flm - suli

[qoska]

Me fal qe po jua them po pse sprovoni me perpara te kerkoni njehere ne google, yahoo per te gjetur keto programe.
Per kete ne fjale ja kam gjetur nje personi para ca kohesh dhe ekzistonin versione me shumice.

[suli]

flm per keshillen....

[Ardi_Pg_ID]

Gjithashtu te sygjeroj googel po perseri ate script qe do gjesh e googel perdore si baze per te studjuar menyren edhe vijueshmerine e programit te rekomandoj edhe te sygjeroj qe te shkruash kodin tend personal do ndjehesh shume here me mire
Ardi

[suli]

nuk eshte me aktuale...
desha vetem te thoja se jam "i ri" ne kete forum, kur shkruajta ketu u nisa me besimin se ndihma do ishte me konkrete se berja moral dhe ca "keshilla"...
anyway

[Ardi_Pg_ID]

Ndihma kalon me konkrete ...
posto kodin qe ke shkruajtur duke ardhur me pyetje me konkrete edhe ne do te te ndihmojme ta permiresosh edhe ta regullosh ate kod
Ardi

[suli]

Ok.

problemi qendron per ta cuar matricen simetrike X->Ak=diag(lambdai);
A1=X;
Ak1 eshte e njejte me matricen Ak, por ndryshon ne elementet
Ak1.App=Ak.App+el.t*Ak.Apq;
Ak1.Aqq=Ak.Aqq-el.t*Ak.Apq;
Ak1.A[j][p]=Ak1.A[p][j]=Ak.A[p][j] - el.s*(Ak.A[q][j] + el.r*Ak.A[p][j]);
Ak1.A[j][q]=Ak1.A[q][j]=Ak.A[q][j] + el.s*(Ak.A[p][j] - el.r*Ak.A[q][j]);
per j!=p
ku t dhe r jane te dhena me formulat: t=c/s r=s/1+c
ku c dhe s jane vlerat e rrotulimeve givensa me kend fi;
kendi fi llogaritet ne baze elementit me te madh Apq (pervec diagonales) te matrices Ak
ndryshimet nga matrica Ak na Ak1 vazhdojne, deri sa elementet e matrices (pervec diagonales) te jene me te vogla se nje vlere eps.

problemi eshte se ky program nuk punon mire.elementet nuk zerohen.
do te isha mirenjohes po me ndihmove ketu.
te fala suli

struct Matrix{
    float A[MAX_SIZE][MAX_SIZE];
    int p;
    int q;
    float Apq;
    float App;
    float Aqq;
};
struct elem{
    float t;
    float r;
    float s;
};


//METODA JACOBIEGO

//funkcje potrzebne

/*
/**
1. funksioni per te gjetru vleren me te madhe te mactices
.kthen kete element, kooordinatat p, dhe q dhe vleren A[p][p] dhe A[q][q]
*/ 

Matrix MaxElement( float A[MAX_SIZE][MAX_SIZE], int n)
{
    Matrix wynik;
    wynik.Apq=0.0;
    for(int i=0;i<n;i++)
        for(int j=0;j<n;j++){
            if(A[i][j]>wynik.Apq && i!=j){
                wynik.Apq=A[i][j];
                wynik.p=i;
                wynik.q=j;
                wynik.App=A[i][i];
                wynik.Aqq=A[j][j];
                
            }
            wynik.A[i][j]=A[i][j];
        }
    return wynik;
}

//----------------------------------------
/*
/** funksjoni kthen vlrene t, r dhe s 
*/
elem DifferencElement(Matrix T, int n){
    
    //obliczenie kata
    float fi=0;
    if(T.App-T.Aqq==0){
        if(T.Apq > 0) fi = (float)PI/4;
        else if(T.Apq < 0) fi = (float) -PI/4;
    }else
        fi=0.5*atanf(2*T.Apq/(T.App-T.Aqq));

    //obliczenie wartosci t i r z korymi roznia sie macierz Ak+1
    elem wynik;
    wynik.t=(float)cos(fi)/sin(fi);
    wynik.r=(float)sin(fi)/(1+cos(fi));
    wynik.s=sin(fi);
    return wynik;
}

//----------------------------------------
/*
/**

*/
void overwrite(float A[MAX_SIZE][MAX_SIZE],float B[MAX_SIZE][MAX_SIZE],int n){
    for (int i=0;i<n;i++)
        for(int j=0;j<n;j++)
            A[i][j]=B[i][j];
}

//--------------------------------------
/*
/**
*/
bool IsEnd(Matrix Ak, int n){
    for(int i=0;i<n;i++)
        for(int j=0;j<n;j++){
            if(Ak.A[i][j] < eps && i!=j)
                return true;
            else
                return false;
        }
}

//--------------------------------------
/*
/**
funk kthen  Ak+1 ne baze te Ak
*/
Matrix Diagonalizacja(Matrix Ak, int n){
    Matrix Ak1;
    Ak=MaxElement(Ak.A,n);//znajdziemy najwiekszy element w macierzy C (oraz informacje
    int p=Ak.p;
    int q=Ak.q;
    elem el;
    el= DifferencElement(Ak,n);
    //float Ak1[MAX_SIZE][MAX_SIZE];
    overwrite(Ak1.A,Ak.A,n);
    Ak1.App=Ak.App+el.t*Ak.Apq;
    Ak1.Aqq=Ak.Aqq-el.t*Ak.Apq;
    for(int j=0;j<n;j++){
        if(p!=j){
            Ak1.A[j][p]=Ak1.A[p][j]=Ak.A[p][j] - el.s*(Ak.A[q][j] + el.r*Ak.A[p][j]);
            Ak1.A[j][q]=Ak1.A[q][j]=Ak.A[q][j] + el.s*(Ak.A[p][j] - el.r*Ak.A[q][j]);
        }
    }
    return Ak1;
}
//--------------------------------------
/*
/**
Metoda Jacobiego 
*/
void Jacobiego(float C[MAX_SIZE][MAX_SIZE], int n)
{
    Matrix A;
    overwrite(A.A,C,n); //przepisujemy macierz C na A
    elem wart;
    wart=DifferencElement(A,n);
    //sprawdzenie czy jest warunek stopu oraz interacja na macierz Ak+1
    while(IsEnd(A,n)==false){
        A=Diagonalizacja(A,n);
        
        printf("\n\n\n");
    }
    //return A;
} 

[Alket123]

ky nuk eshte forum matematike. Mburren moderatoret eshte reklame per forumin. Ne mezi sa kemi kaluar klasen.

ke tre shembuj nqs te ndihmojne sadopak.

[Alket123]

ose ka nje rruge me te thjeshte. Krijohe ne Matlab, pastaj exportohe/konvertoje ne C

help
http://www.mathworks.com/support/sol...lution=1-18L04

[suli]

flm prof.
e vleresoj ndihmen
respekte - suli