vs.blog :)

Zeki Optimizasyon Teknikleri - 2 (Tavlama Benzetimi)

2009-10-22T21:16:00.023+03:00

Optimizasyon algoritmalarımızdan ikincisi Simulated Annealing olarak ta bilinen tavlama benzetimi algoritması. İsmi Metalurji biliminden gelmektedir. Metallerin tavlanması işleminden esinlenerek ortaya çıktığı için bu ismi almıştır. Genellikle ayrık optimizasyon problemleri için kullanılır.

Yöntemden kabaca bahsetmek için Hill Climbing yöntemimizi hatırlayalım. Hill-Climbing yönteminde ilk çözüm üretmiş, bir komşu üreteci fonksiyonu ile ilk çözümümüze komşu olan yeni bir çözüm kümesi üretmiştik.*Bu yeni kümeyi kendi içerisinde değerlendirip en iyi elemanını seçmiştik. Eğer bu seçtiğimiz en iyi eleman bizim eski çözümümüzden daha iyi bir çözümse yeni çözümümüz olarak bu değeri kabul etmiştik. Ayrıca Hill-Climbing yönteminin kolay uygulanmasının yanında kötü bir özelliği olarak local çözümlere takılma olarak belirtmiştik.
Algoritmanın bu lokal çözümlere takılmasının nedeni bulduğumuz bir eğim doğrultusunda ilerlememiz, bu yönde çözümleri kısıtlamamızdır. Lokal çözüme doğru hızla ilerlerken çözüm kümesinin diğer kısımlarını göz ardı ederiz. Eğer şans eseri ilk çözüm üretecimiz başlangıç çözümünü global çözüme yakın bir bölgede üretirse muhtemelen global çözüme ulaşırız. Bundandır ki, Hill Climbing algoritması başlangıç noktasına çok bağımlıdır demiştik.

Hill Climbing algoritmamız üzerine bir takım modifikasyonlar yaptığımızı düşünelim. Yine ilk çözüm üretecimiz bir çözüm oluştursun, komşu çözüm üretecimiz komşu çözümler kümesini oluştursun, yine üretilen komşu çözüm kümesini kendi içerisinde değerlendirerek buradan en iyi çözümü seçelim. Klasik Hill Climbing algoritmamızda komşu çözüm kümesinden seçilen en iyi çözüm eğer eski çözümümüzden daha kötü ise bu değeri kullanmayıp elimizdeki çözümü saklıyorduk. Şimdi biraz daha insaflı davranıp komşu kümeden seçilen kötü çözüme bir şans verelim. Bu şans P olsun. Belkide verdiğimiz bu şans hatırına kötü çözüm bizi komşusu olan daha iyi çözümlere götürebilir. Hem böylece Hill-Climbing in local çözümlere takılıp kalma özelliğinden kurtulabiliriz. Eğer kötü çözümlere verdiğimiz seçilme şansı(olasılığı) P çalışma boyunca sabit kalırsa bu sefer tam güzel çözümlere doğru ilerlerken oradan oraya oradan oraya atlarız. Bu nedenle bu P olasılığının(iyi çözümü feda ederek yerine kötü çözümü kabul etme olasılığı) dinamik olması gereklidir. Bu dinamiklik iterasyonla P nin azaltılması şeklinde(ve daha pek çok şekilde) yapılablir. Böylece problem çözümünün ilk kısımlarında çözüm bölgeleri arasında sıçrayış daha fazla olurken iterasyon sayımız artıp elde ettiğimiz çözümler iyice güzelleştikçe 0 a yaklaşır. (Arama bölgemiz daralır.)

İşte yukarıda modifiye edilmiş bir Hill-Climbing algoritması gördük. Tavlama Benzetimi algoritmasının temel prensibi tam olarak budur. Kötü çözümü seçme olasılığı sistemli bir şekilde sıcaklıkla azaltılır. Sıcaklık iterasyona bağlı(genellikle düzgün veya logaritmik azalan) bir ifadedir.

Tavlama benzetimi algoritmasının temel amacı çözüm uzayında aranmadık bölge bırakmamaktır.

Ayrıca Tavlama Benzetimi algoritmasında komşu çözüm kümesinin büyüklüğüde dinamik şekilde değiştirilebilir. (Probleme göre değişse de genel olarak çözümü çok fazla etkilememektedir) Bu şekilde bir kullanım benimsenirse yaygın kullanım küçük bir kümeyle başlayıp iterasyonla büyütmek şeklinde olur.

Sıcaklığın azaltılması, olasılık değerinin sıcaklıkla azaltılması, çözüm kümesinin değiştirilmesi için pek çok yöntem kullanılır.
Bunlara girmeden basit bir problemi tavlama benzetimi algoritması ile çözelim. Aşağıda C++ kodu çözüm olarak string şeklinde "VOLKANSALMA@BLOGSPOT@COM" dizisini oluşturmaya çalışmaktadır.
Örnek bu sitedeki genetik algoritma ile Hello! World yazısını oluşturma probleminden esinlenilmiştir. Çözümün kalitesini değerlendirmek için gerekli fonksiyon (Fitness function) linki verilen siteden direk alınmıştır. .

Bu yazının hazırlanmasında hocam Doç.Dr. M.Ali Akcayol' un Zeki Optimizasyon Teknikleri Dersi ders notları kullanılmıştır.

*Hill-Climbing için yazmış olduğum C++ kod örneğinde bir yanlışlık yaptığımın farkına vardım. Aslında çözüm kümesi kendi içerisinde değerlendirilip en iyi eleman seçilmesi gerekirken, ben ilk bulduğum iyi çözüm ile mevcut çözümü değiştirdim. Sonrada diğer iterasyona geçtim. Bu yanlışlık daha güzel çözümleri kullanamamıza neden oldu. Böyle bir hata çözüme yakınsama hızımızı düşürür.

*Biçimlendirme güçlükleri nedeni ile kodda "<" ve ">" işaretleri yerine "-" işareti kullanılmıştır.

//Simulated Annealing C++ Sample Code
//generates "VOLKANSALMA@BLOGSPOT@COM" string.
//inspired from http://www.generation5.org/content/2003/gahelloworld.asp
//A "Hello World!" Genetic Algorithm Example
//by volkan salma
//http://volkansalma.blogspot.com

#include -iostream-     // for cout etc.
#include -vector-     // for vector class
#include -string-     // for string class
#include -algorithm-    // for sort algorithm
#include -time.h-     // for random seed
#include -math.h-     // for abs()

#define SA_POPSIZE  50  // sa population size
#define SA_MAXITER  300  // maximum iterations
#define COEF            0.01
#define SA_TARGET  std::string("VOLKANSALMA@BLOGSPOT@COM")

using namespace std;

struct sa_struct
{
 string str;      // the string
 unsigned int fitness;   // its fitness
};

typedef vector-sa_struct- sa_vector;// for brevity

void generateNeighbourSolutionsToMySolution(sa_struct myCurrentSolution,sa_vector &neighbourSolutionSet);
void calcFitnessOfEveryElementInNeighbourSet(sa_vector &neighbourSolutionSet);
void init_population(sa_vector &population, sa_vector &buffer);
void calc_fitness(sa_vector &population);
bool fitness_sort(sa_struct x, sa_struct y);
void sort_by_fitness(sa_vector &population);
sa_struct generateFirstSolution();
sa_struct getNeighbourSolutionSetsBestElement(sa_vector &population);
void print_best(sa_vector &sav);


int main()
{
    srand(unsigned(time(NULL)));
    sa_struct myCurrentSolution;
    sa_struct buffer;
 vector -sa_struct- neighbourSolutionSet;
 double myRandomNumber;

 myCurrentSolution = generateFirstSolution();

 double Temperature = 100;
 double boltzman;


 for (int i=0; i - SA_MAXITER; i++)
 {
     generateNeighbourSolutionsToMySolution(myCurrentSolution,neighbourSolutionSet);

  calcFitnessOfEveryElementInNeighbourSet(neighbourSolutionSet);

  buffer = getNeighbourSolutionSetsBestElement(neighbourSolutionSet);

  if( buffer.fitness < myCurrentSolution.fitness) myCurrentSolution = buffer;

  else ///here is critical section. ;)
                     // we can change or not change the ourSolution with worst one.
  {
            myRandomNumber = rand()%10;
            myRandomNumber /= 10.0;   // i want to generate a rondom value in [0,1] range.
            boltzman = COEF * (1.0/ (exp(buffer.fitness - myCurrentSolution.fitness)/Temperature));

            if( boltzman > myRandomNumber)
            {
                cout--"Boltzman: "--boltzman--" > MyRandom: "--myRandomNumber--endl;
                myCurrentSolution = buffer;
                cout--"**********************************************************************"--endl;
            }
  }

  Temperature -= (100.0 / SA_MAXITER);

  cout --i--" My Solution: " -- myCurrentSolution.str -- " (" -- myCurrentSolution.fitness << ") Temp:"--Temperature-- endl;

  if(myCurrentSolution.fitness == 0) break;

 }

    cout -- "My Best Solution Is: " -- myCurrentSolution.str -- " (" -- myCurrentSolution.fitness -- ")" -- endl;
    cin--boltzman;

    return 0;


}

void calcFitnessOfEveryElementInNeighbourSet(sa_vector &population)
{
 string target = SA_TARGET;
 int tsize = target.size();
 unsigned int fitness;

 for (int i=0; i - SA_POPSIZE; i++)
 {
  fitness = 0;
  for (int j=0; j-tsize; j++)
  {
   fitness += abs(int(population[i].str[j] - target[j]));
  }

  population[i].fitness = fitness;
 }
}

bool fitness_sort(sa_struct x, sa_struct y)
{
    return (x.fitness < y.fitness);
}

void sort_by_fitness(sa_vector &population)
{
    sort(population.begin(), population.end(), fitness_sort);
}

sa_struct getNeighbourSolutionSetsBestElement(sa_vector &population)
{
    sa_struct bestMemberOfSolutionSet;
    sort_by_fitness(population);
    bestMemberOfSolutionSet = population[0];

    return bestMemberOfSolutionSet;
}

sa_struct generateFirstSolution()
{
        int tsize = SA_TARGET.size();
        sa_struct result;

  result.fitness = 0;
  result.str.erase();

  for (int j=0; j-tsize; j++)
  {
   result.str += (rand() % 26) + 64;
  }

  string target = SA_TARGET;

  for (int j=0; j-tsize; j++)
  {
   result.fitness += abs(int(result.str[j] - target[j]));
  }


  return result;
}

void generateNeighbourSolutionsToMySolution(sa_struct myCurrentSolution,sa_vector &neighbourSolutionSet)
{
    int tsize = SA_TARGET.size();
    sa_struct buffer;
    neighbourSolutionSet.clear();

    for (int i=0; i-SA_POPSIZE; i++)
 {
     buffer.str = myCurrentSolution.str;
     buffer.str[i%tsize] = (rand() % 26) + 64;
     neighbourSolutionSet.push_back(buffer);
 }
}

Programın şu şekilde bir çıktısı oluyor: (Yıldızlarla gösterilen bölgelerde bulunan çözüm daha kötü olmasına rağmen olasılık değerlendirmesini geçmiş, çözüm olarak seçilmiştir. )

0 My Solution: SGGYVNFHLND@HEL@RFOYMGEO (143) Temp:99.6667
1 My Solution: SGGYHNFHLND@HEL@RFOYMGEO (129) Temp:99.3333
2 My Solution: SGGYHNFHLND@HEL@RFOYAGEO (117) Temp:99
3 My Solution: SGGYHNRHLND@HEL@RFOYAGEO (105) Temp:98.6667
4 My Solution: SGGLHNRHLND@HEL@RFOYAGEO (92) Temp:98.3333
5 My Solution: SGGLHNRHLND@HEL@RFOYAGRO (85) Temp:98
6 My Solution: SGGLHNRHLND@HEL@RMOYAGRO (78) Temp:97.6667
7 My Solution: SGGLHNRHLND@HML@RMOYAGRO (72) Temp:97.3333
8 My Solution: SNGLHNRHLND@HML@RMOYAGRO (65) Temp:97
9 My Solution: SNGLBNRHLND@HML@RMOYAGRO (59) Temp:96.6667
10 My Solution: SNLLBNRHLND@HML@RMOYAGRO (54) Temp:96.3333
11 My Solution: SNLLBNRHLND@HMLERMOYAGRO (49) Temp:96
12 My Solution: VNLLBNRHLND@HMLERMOYAGRO (46) Temp:95.6667
13 My Solution: VNLLBNRHLND@BMLERMOYAGRO (40) Temp:95.3333
14 My Solution: VNLLBNRHLND@BMLERMORAGRO (37) Temp:95
15 My Solution: VNLLBNRHLND@BMLERMORAARO (35) Temp:94.6667
16 My Solution: VNLLBNRALND@BMLERMORAARO (28) Temp:94.3333
17 My Solution: VNLLBNRALNB@BMLERMORAARO (26) Temp:94
18 My Solution: VNLLANRALNB@BMLERMORAARO (25) Temp:93.6667
19 My Solution: VNLLANRALNA@BMLERMORAARO (24) Temp:93.3333
20 My Solution: VNLLANRALNA@BMOERMORAARO (21) Temp:93
21 My Solution: VNLLANRALNA@BMOERMORAAOO (18) Temp:92.6667
22 My Solution: VNLLANRALNA@BMOERMORAAOM (16) Temp:92.3333
23 My Solution: VNLLANRALNA@BMOERMORACOM (14) Temp:92
24 My Solution: VNLLANRALNA@BMOERMOTACOM (12) Temp:91.6667
25 My Solution: VNLLANRALNA@BMOESMOTACOM (11) Temp:91.3333
26 My Solution: VNLLANRALNA@BMOESPOTACOM (8) Temp:91
Boltzman: 0.91 > MyRandom: 0.8
**********************************************************************
27 My Solution: VNLLANRALNA@BMOESPOTACOM (8) Temp:90.6667
28 My Solution: VOLLANRALNA@BMOESPOTACOM (7) Temp:90.3333
29 My Solution: VOLLANRALNA@BMOGSPOTACOM (5) Temp:90
Boltzman: 0.9 > MyRandom: 0.7
**********************************************************************
30 My Solution: VOLLANRALNA@BMOGSPOTACOM (5) Temp:89.6667
Boltzman: 0.896667 > MyRandom: 0
**********************************************************************
31 My Solution: VOLLANRALNA@BKOGSPOTACOM (5) Temp:89.3333
Boltzman: 0.893333 > MyRandom: 0.5
**********************************************************************
32 My Solution: VOLLANRALNA@BKOGSPOTACOM (5) Temp:89
33 My Solution: VOLLANRALNA@BKOGSPOTACOM (5) Temp:88.6667
Boltzman: 0.326186 > MyRandom: 0.3
**********************************************************************
34 My Solution: VNLLANRALNA@BKOGSPOTACOM (6) Temp:88.3333
Boltzman: 0.883333 > MyRandom: 0.2
**********************************************************************
35 My Solution: VNLLANRALNA@BKOGSPOTACOM (6) Temp:88
36 My Solution: VNLKANRALNA@BKOGSPOTACOM (5) Temp:87.6667
Boltzman: 0.876667 > MyRandom: 0.6
**********************************************************************
37 My Solution: VNLKANRALNA@BKOGSPOTACOM (5) Temp:87.3333
38 My Solution: VOLKANRALNA@BKOGSPOTACOM (4) Temp:87
39 My Solution: VOLKANSALNA@BKOGSPOTACOM (3) Temp:86.6667
40 My Solution: VOLKANSALNA@BKOGSPOTACOM (3) Temp:86.3333
Boltzman: 0.863333 > MyRandom: 0.8
**********************************************************************
41 My Solution: VOLKANSALNA@BKOGSPOTACOM (3) Temp:86
42 My Solution: VOLKANSALMA@BKOGSPOTACOM (2) Temp:85.6667
Boltzman: 0.31515 > MyRandom: 0.3
**********************************************************************
43 My Solution: VPLKANSALMA@BKOGSPOTACOM (3) Temp:85.3333
Boltzman: 0.313924 > MyRandom: 0
**********************************************************************
44 My Solution: VQLKANSALMA@BKOGSPOTACOM (4) Temp:85
45 My Solution: VQLKANSALMA@BKOGSPOTACOM (4) Temp:84.6667
46 My Solution: VOLKANSALMA@BKOGSPOTACOM (2) Temp:84.3333
Boltzman: 0.843333 > MyRandom: 0.1
**********************************************************************
47 My Solution: VOLKANSALMA@BKOGSPOTACOM (2) Temp:84
Boltzman: 0.309019 > MyRandom: 0
**********************************************************************
48 My Solution: VPLKANSALMA@BKOGSPOTACOM (3) Temp:83.6667
49 My Solution: VOLKANSALMA@BKOGSPOTACOM (2) Temp:83.3333
Boltzman: 0.833333 > MyRandom: 0.3
**********************************************************************
50 My Solution: VOLKANSALMA@BKOGSPOTACOM (2) Temp:83
51 My Solution: VOLKANSALMA@BKOGSPOTACOM (2) Temp:82.6667
Boltzman: 0.826667 > MyRandom: 0.2
**********************************************************************
52 My Solution: VOLKANSALMA@BKOGSPOTACOM (2) Temp:82.3333
Boltzman: 0.823333 > MyRandom: 0.5
**********************************************************************
53 My Solution: VOLKANSALMA@BKOGSPOTACOM (2) Temp:82
54 My Solution: VOLKANSALMA@BLOGSPOTACOM (1) Temp:81.6667
Boltzman: 0.816667 > MyRandom: 0.7
**********************************************************************
55 My Solution: VOLKANSALMA@BLOGSPOTACOM (1) Temp:81.3333
Boltzman: 0.813333 > MyRandom: 0.3
**********************************************************************
56 My Solution: VOLKANSALMA@BLOGSPOTACOM (1) Temp:81
57 My Solution: VOLKANSALMA@BLOGSPOT@COM (0) Temp:80.6667
My Best Solution Is: VOLKANSALMA@BLOGSPOT@COM (0)

OpenCV 2.0 Windows .lib dosyaları?

2009-10-18T15:05:00.003+03:00

OpenCV v2.0 kararlı sürüm yayınlandı. Pek çok güncelleme ve yeni fonksiyonlar içeriyor.(Neler Mesela?)

OpenCV yi windows ortamında visual studio ile kullanıyorsanız bir süprizle karşılaşabilirsiniz. Kurulumla birlikte lib dosyaları gelmiyor. Bu dosyaları elde etmek için kaynak kodun isteğinize uygun konfigürasyon parametreleri ile derlenmesi gerekiyor.

Derleme işlemine girmeden hazır olarak kullanmak isterseniz. Kendim için derlediğim libleri sizinle paylaşabilirim :)

Download OpenCV 2.0 lib,dll,include files for Visual Studio.

Derleme ve konfigürasyon işlemin nasıl yapıldığını öğrenmek isterseniz bluekid abimize bakalım.
Bu dosyaları kullanarak nasıl opencv projesi oluşturacağımızı bilmiyorsak buraya bakalım.(Visual Studio,Code::Blocks)
Baktınız ama işin içinden çıkamadıysanız yorum şeklinde sorun. Cevap vermeye çalışayım.

*Kütüphaneler Visual Studio Prof. Editionda OpenMP desteği ile derlenmiştir fakat test edilmemiştir.

Zeki Optimizasyon Teknikleri - 1 (Hill Climbing)

2009-10-08T00:09:00.005+03:00

Optimizasyon Algoritmalarımızdan ilki Stochastic(içerisinde rastgelelik bulunan) tekniklerden birisi olan Hill Climbing algoritması. Bu yöntem pek çok uygulamada uygulama kolaylığı nedeniyle seçilebilir.

Optimizasyonu istenen problemin gösterimi, elde olan çözüme göre komşu çözüm üreteci, ilk çözüm üreteci ve çözüm değerlendirme fonksiyonu uygulama için yeterli olmaktadır.

Algoritmanın dezavantajları olarak yerel çözümlere takılabilmesini ve çözümün seçilen başlangıç noktasına çok bağlı olmasını sayabiliriz.

Hill Climbing yöntemden bahsedersek:

İlk olarak başlangıç çözümü üretilir ve bu çözüm eniyi çözüm olarak kabul edilir. (Bu çözümün üretimi probleme göre rastgele olabileceği gibi problem konusunda uzman görüşü veya daha önce yapılmış çalışmalar yardımıyla da olabilir.)

Daha sonra bu probleme komşu bir çözüm kümesi üretilir.(Değişik sayıda olabilir.) Çözüm kümesi elemanları değerlendirilir. Çözüm kümesinin en iyi elemanı eğer bir önceki adımdaki eniyi çözümden daha iyi ise, eniyi çözüm güncellenir. Sonrasında en iyi çözüme göre yeni komşular üretilir ve iteratif olarak devam edilir. Belirli iterasyon sayısına ulaşıldığında, yeterli bir çözüm elde edildiğinde veya artık çözümler değişmediğinde iterasyon sonlandırılır.

Hill Climbing yöntemini kullanarak derste geçen bir basit problem için C++ uygulama kodu hazırladım. Problem: bd 30 karakterlik bir binary dizisi olmak üzere f(bd) = |20 * dizideki_bir_sayisi(bd) - 100| fonksiyonunu max yapan bd değeri? (Problemimizin çözümünün "111111111111111111111111111111" olduğu rahatça görülüyor.)

Bu yazının hazırlanmasında hocam Doç.Dr. M.Ali Akcayol' un Zeki Optimizasyon Teknikleri Dersi ders notları kullanılmıştır.

*Biçimlendirme güçlükleri nedeni ile kodda "<" ve ">" işaretleri yerine "-" işareti kullanılmıştır.

//f(bd) = |20 * bir_sayisi(bd) - 100| fonksiyonunu max yapan
// 30 karakterlik binary dizisinin(bd) hill climbing algoritmasi ile
//bulunmasi.
//Volkan SALMA. http://volkansalma.blogspot.com

#include -iostream-
#include -vector-
#include -stdlib.h- //for rand funct.
#include -time.h- //time for rand function seed.
#include -math.h- // for abs function


using namespace std;

#define BINARY_ARRAY_SIZE 30
#define ITERATION_NUM 15
#define SOLUTION_SET_ELEMENT_CNT 30


struct Solution
{
    bool value[BINARY_ARRAY_SIZE];
};

vector-solution- generate40NeighbourSolutionToInput(Solution); //generate new solution depending old one.
Solution generateFirstSolutionRandomly(void);
int countOnes(bool); //number of 1's
long evaluteMyFunctionAndGiveResult(Solution); // f(bd) = |20 * bir_sayisi(bd) - 100|
void printSolutionResult(Solution,int,bool);


int main()
{
    Solution myBestSolution; //holds best solution.
    vector-solution- testSolutionSet;

    myBestSolution = generateFirstSolutionRandomly();

    for(int i = 0; i < ITERATION_NUM; i++)
    {
        testSolutionSet = generate40NeighbourSolutionToInput(myBestSolution);

        for(int j = 0 ; j < SOLUTION_SET_ELEMENT_CNT; j++)
        {
            if(evaluteMyFunctionAndGiveResult(testSolutionSet[j]) > evaluteMyFunctionAndGiveResult(myBestSolution) )
            {
                myBestSolution = testSolutionSet[j];
                break;
            }
        }

         printSolutionResult(myBestSolution,i+1,false);
    }

    printSolutionResult(myBestSolution,-1,true);

    return 0;
}


vector-solution- generate40NeighbourSolutionToInput(Solution in)
{
    Solution newSln;
    vector-solution- result;

    for(int i = 0; i < SOLUTION_SET_ELEMENT_CNT; i++)
    {
            for(int j = 0; j < BINARY_ARRAY_SIZE; j++ )
            {
                if(i != j) newSln.value[j] = in.value[j];
                else newSln.value[j] = !in.value[j];
            }

            result.push_back(newSln);
    }

    return result;
}

Solution generateFirstSolutionRandomly(void)
{
    Solution result;

    srand((unsigned)time(NULL));

    for(int i = 0; i < BINARY_ARRAY_SIZE; i++)
    {
        result.value[i] = rand()%2;
    }

    return result;
}

int countOnes(Solution in)
{
    int cnt = 0;

    for(int i= 0; i < BINARY_ARRAY_SIZE; i++)
    {
        if(in.value[i] == 1) cnt++;
    }

    return cnt;
}

long evaluteMyFunctionAndGiveResult(Solution solution)
{
     // f(bd) = |20 * bir_sayisi(bd) - 100|
     return abs(20*countOnes(solution)-100);
}
void printSolutionResult(Solution bestSolution,int iteration, bool isBest)
{
    if (isBest)cout--"Best Solution is:"--endl;
    else cout--iteration--". iteration:"--endl;
    for(int i = 0; i < BINARY_ARRAY_SIZE; i++)
    {
        cout--bestSolution.value[i];
    }

    cout--endl;
}

OpenCV altında SVM kullanmak. (OpenCV SVM Sample Code)

2009-08-23T20:46:00.001+03:00

Bilindiği gibi OpenCV kütüphanesinin machine learning modülünde libsvm wraplenmiştir. Böylece biz OpenCV kullanıcıları harici bir kütüphane kullanmadan SVM(Support Vector Machines) i uygulamalarımızda kullanabiliriz. Aşağıda bu konuda bir örnek hazırlamaya çalıştım. Yorumlarda açıklamalar mevcuttur.

#include
#include 
#include 
#include 

using namespace std;

int main()
{
 CvSVM SVM;
 CvMat trainData;
 CvMat labels;
 CvTermCriteria criteria = cvTermCriteria(CV_TERMCRIT_EPS,100,0.001); //eğitim için term kriterimiz.

 int array_labels[5]= {1,1,1,-1,-1}; //5 eğitim verisinin bildiğimiz etiketleri.(labels, responses).
           //ilk 3 veri 1 sınıfına, diğer iki veri ise  -1 sınıfına ait.

 //eğitim verilerini diziye doldurdum. Gerçek uygulamalarda filedan okuma veya hesaplama şeklinde olur.
 //kodu basit tutmak için bu şekilde gösterdim.
 #pragma region data
 //5 egitim verisi olsun herbiri 40 tane feature içersin.
 
 float array_trainData[5][40] = { 1,0.443137,0.333333,0.392157,0.0313725,
          0.270588,0.870588,0.964706,1,0.439216,
          0.337255,0.380392,0.0235294,0.266667,0.862745,
          0.964706,1,0.431373,0.337255,0.368627,
          0.0156863, 0.262745, 0.854902, 0.964706, 1,
          0.423529, 0.341176, 0.360784, 0.0156863, 0.258824,
          0.847059, 0.964706, 1, 0.419608, 0.341176,
          0.34902, 0.0117647, 0.254902, 0.839216, 0.968627,

          0.572549,0.4,0.113725,0.0784314,0.137255,
          0.160784,0.00784314,0.337255,0.580392,0.407843,
          0.121569, 0.0745098, 0.133333, 0.156863, 0.00392157,
          0.341176, 0.584314, 0.411765, 0.12549, 0.0666667 ,
          0.129412, 0.152941, 0, 0.34902, 0.592157,
          0.419608, 0.129412, 0.0588235, 0.12549, 0.14902,
          0.00392157, 0.352941, 0.596078, 0.423529, 0.137255,
          0.054902, 0.121569, 0.145098, 0.00784314, 0.356863,

          0.882353, 0.741176, 0.305882, 0.305882, 0.278431,
          0.141176, 0.156863, 0.509804, 0.882353, 0.741176,
          0.305882, 0.305882, 0.27451, 0.141176, 0.14902,
          0.501961, 0.886275, 0.741176, 0.305882, 0.305882,
          0.270588, 0.141176, 0.141176, 0.498039, 0.886275,
          0.741176, 0.305882, 0.305882, 0.266667, 0.141176,
          0.137255, 0.490196, 0.886275, 0.737255, 0.305882,
          0.305882, 0.266667, 0.137255, 0.129412, 0.486275,

          0.929412, 0.917647, 0.717647, 0.356863, 0.133333,
          0.156863, 0.231373, 0.709804, 0.929412, 0.913725,
          0.721569, 0.356863, 0.137255, 0.164706, 0.219608,
          0.694118, 0.92549, 0.905882, 0.729412, 0.356863,
          0.141176, 0.168627, 0.207843,0.67451, 0.917647,
          0.901961, 0.733333, 0.360784, 0.141176, 0.172549,
          0.196078, 0.658824, 0.909804, 0.894118, 0.737255,
          0.360784, 0.141176, 0.180392,0.180392, 0.639216,

          0.760784, 0.372549, 0.407843, 0.298039, 0.160784,
          0.72549, 0.988235, 0.960784, 0.772549, 0.372549,
          0.411765, 0.305882, 0.160784, 0.709804, 0.992157,
          0.960784, 0.780392, 0.376471, 0.419608, 0.309804,
          0.160784, 0.698039, 0.992157, 0.960784, 0.784314,
          0.376471, 0.419608, 0.317647, 0.164706, 0.686275,
          0.996078, 0.964706, 0.792157, 0.376471, 0.423529,
          0.32549, 0.160784, 0.670588, 1, 0.964706
 };

#pragma endregion data

 cvInitMatHeader(&trainData,5,40,CV_32FC1,array_trainData); //eğitim verimizi tutacak matris.
 cvInitMatHeader(&labels,1,5,CV_32SC1,array_labels); //etiketlerimizi tutacak matris.
            //classification olduğu için CV_32SC1 kullandık.


 SVM.train(&trainData,&labels,0,0,CvSVMParams(CvSVM::C_SVC,CvSVM::RBF,3,0.025,0,1,0.5,0.1,NULL,
  criteria));  //train fonksiyonumuz (libsvm ile kontrol etmek için libsvm ile
      //gelen svm-train.exe nin default parametreleri kullanıldı.)
      //Parametrelerin detaylı açıklaması için opencv docs.

 //SVM.save("test.volkan"); //eğittiğimiz modeli kaydedip daha sonra kullanabiliriz.
for(int i = 0; i < 5; i++ )
 {
  double r;
  CvMat sample;
  cvGetRow(&trainData, &sample, i );

  r = SVM.predict( &sample );  //predict edip sonuçları kontrol edelim.
 }
}

Programınız için performans ipuçları.(x86)

2009-05-01T21:48:00.013+03:00

Genel amaçlı üretilmiş mimariler üzerinde sinyal işleme, görüntü işleme gibi ağır matematiksel işlemler gerektiren uygulamalar yapıyorsanız optimizasyon çok önemli bir hale gelir. Bilgisayar tasarımı çalışırken bir zamanlar Türkiye Görüntü İşleme mail grubuna gönderilen bir yazı aklıma geldi. Faydalı bir yazı olduğunu düşündüm ve blogumda paylaşmak istedim. Yazı 17 Şubat 2008 tarihinde Mustafa Sakar tarafından AMD x86 Code Optimization Guide tan yararlanılarak yazılmış.
Özeti biraz daha özetleyip bilgimiz olan yerlere burnumuzu sokarsak ;) 1- 32-bit data tiplerini kullanın. (Neden çünkü mimari bunun üzerine. Bütün aritmatik işlem birimleri 32 bit için design edilmiş. Daha büyük veri tipleri bölünüp yapılır. Küçük olanlarda 32 bitlik registerlara taşınır işaretli ise işaret genişletilir ekstra işlemler...) 2- Integer işlemlerinin işaretine önceden karar verin. Bölüm ve kalan hesaplarken, döngü sayaçları, array indexleri için unsigned tipler, Integer-to-float dönüşümleri için signed tipler kullanın. Kötü örnek:

double x;        // MOV [tmep+4],0
unsigned int i;  // MOV EAX, i
        // MOV [temp], EAX
x = i;           // FILD QWORD PTR [temp]
        // FSTP QWORD PTR [x}

İyi örnek:

double x;        // FILD DWORD PTR [i]
int i;           // FSTP QWORD PTR [x]
x = i;

3- Bölme işlemlerini işaretsiz tiplerle yapmak daha hızlıdır. Kötü örnek:

int i;
i = i / 4;         // MOV EAX,i
        // CDQ
        // AND EDX,3
        // ADD EAX, EDX
        // SAR EAX, 2
        // MOV i, EAX

İyi Örnek:

unsigned i;
i = i / 4;         // SHR i,2

4- Pointer-style kod yerine Array-style kod kullan. //Bu çok önemli.. Kötü örnek:

typedef struct {
float x,y,z,w;
} VERTEX;
typedef struct {
float m[4][4];
} MATRIX;

void XForm (float *res, const float *v, const float *m, int
numverts) {
float dp;
int i;
const VERTEX* vv = (VERTEX *)v;
for (i = 0; i < dp =" vv-">x * *m++;
     dp += vv->y * *m++;
     dp += vv->z * *m++;
     dp += vv->w * *m++;
     *res++ = dp; /* write transformed x */
     dp = vv->x * *m++;
     dp += vv->y * *m++;
     dp += vv->z * *m++;
     dp += vv->w * *m++;
     *res++ = dp; /* write transformed y */
     dp = vv->x * *m++;
     dp += vv->y * *m++;
     dp += vv->z * *m++;
     dp += vv->w * *m++;
     *res++ = dp; /* write transformed z */
     dp = vv->x * *m++;
     dp += vv->y * *m++;
     dp += vv->z * *m++;
     dp += vv->w * *m++;
     *res++ = dp; /* write transformed w */
     ++vv; /* next input vertex */
     m -= 16; /* reset to start of transform matrix */
}
}

iyi Örnek:

typedef struct {
float x,y,z,w;
} VERTEX;
typedef struct {
float m[4][4];
} MATRIX;

void XForm (float *res, const float *v, const float *m, int
numverts) {
int i;
const VERTEX* vv = (VERTEX *)v;
const MATRIX* mm = (MATRIX *)m;
VERTEX* rr = (VERTEX *)res;
for (i = 0; i <>x = vv->x*mm->m[0][0] + vv->y*mm->m[0][1] +
             vv->z*<div id=":1l" class="ii gt"><wbr>mm->m[0][2] + vv->w*mm->m[0][3];
     rr->y = vv->x*mm->m[1][0] + vv->y*mm->m[1][1] +
             vv->z*<wbr>mm->m[1][2] + vv->w*mm->m[1][3];
     rr->z = vv->x*mm->m[2][0] + vv->y*mm->m[2][1] +
             vv->z*<wbr>mm->m[2][2] + vv->w*mm->m[2][3];
     rr->w = vv->x*mm->m[3][0] + vv->y*mm->m[3][1] +
             vv->z*<wbr>mm->m[3][2] + vv->w*mm->m[3][3];
}
}

5- Kısa donguleri tamamen kaldırmalısınız. Kötü örnek:

// 3D-transform: multiply vector V by 4x4 transform matrix M
for (i=0; i<4; i++) {
        r[i] = 0;
        for (j=0; j<4; j++) {
                r[i] += M[j][i]*V[j];
        }
}

İyi Örnek:

// 3D-transform: multiply vector V by 4x4 transform matrix M
r[0] =  M[0][0]*V[0] + M[1][0]*V[1] + M[2][0]*V[2] +
      M[3][0]*V[3];
r[1] =  M[0][1]*V[0] + M[1][1]*V[1] + M[2][1]*V[2] +
      M[3][1]*V[3];
r[2] =  M[0][2]*V[0] + M[1][2]*V[1] + M[2][2]*V[2] +
      M[3][2]*V[3];
r[3] =  M[0][3]*V[0] + M[1][3]*V[1] + M[2][3]*V[2] +
      M[3][3]*v[3];

6- Gereksiz Store-to-Load hesaplamalardan kacınmalısınız. Kötü Örnek: for (k = 1; k < VECLEN; k++) { x[k] = x[k-1] + y[k]; //Buralar facia durmadan memory e başvur. registera load et //sonucu geri registerdan memory'e store et.. } for (k = 1; k < VECLEN; k++) { x[k] = z[k] * (y[k] - x[k-1]); } İyi örnek:

double x[VECLEN], y[VECLEN], z[VECLEN];
unsigned int k;
double t;
t = x[0];

for (k = 1; k < VECLEN; k++) {
t = t + y[k];
x[k] = t;
}
t = x[0];
for (k = 1; k < VECLEN; k++) {
t = z[k] * (y[k] - t);
x[k] = t;
}

7 - Çarpma işlemi bölmeden hızlıdır. Kötü örnek:

double a,b,c,e,f;
e = a/c;
f = b/c;

İyi örnek:

double a,b,c,e,f,t;
t = 1/c;
e = a*t
f = b*t;

Kötü örnek:

int i,j,k,m;
m = i / j / k;

İyi örnek:

int i,j,k,l;
m = i / (j * k);

8 - Hızlı Floating-Point-to-Integer cevrimi kullanın. Yavaş Örnek:

double x;
int i;
i = x;

Hızlı Örnek:

#define DOUBLE2INT(i,d) \
{double t = ((d)+6755399441055744.0); i=*((int *)(&t));}

double x;
int i;

DOUBLE2INT(i,x);

***Bazı maddeleri atladım. En önemli gördüklerimi seçmek istedim. Görüntü işleme ile uğraşıyor performans sıkıntısı yaşıyorsanız hem de çift çekirdekli işlemciniz varsa ilk işiniz bazı işleri paralel yaptırmak olsun. Bluekid abimize link verelim. OpenMP genel. OpenMP uygulama. İmkanınız varsa güncel compiler lar kullanın. Yeni compilerla yeni donanımlar için en iyi optimizasyon seçenekleri sağlarlar. Performans bakımından aynı mfc uygulamasının MS un VC++ 9 ile derlenmesi ile VC++ 6 ile derlenmesi arasında inanılmaz farklar gözlenmiştir. $$$* Ayrıca intel işlemci kullanıyorsanız intelin kendi compiler ı için vaatleri büyük. (Denenmemiştir.)

minix3 macerası - 2 (kurulum)

2009-04-25T14:37:00.011+03:00

Minix3 için gerekli sistem gereksinimleri şöyle.

386, 486, or Pentium CPU.
16 veya 8 MB ram.
Kernel için 50 mB. CD deki diğer paketlerin hepsini kuracaksanın 600 mB harddisk alanı.
Boş bir CD iso yun yazmak için. Flashtan boot edebilen bir makinanız varsa flash diskten de kurulum yapabilirsiniz(Ben CD den yaptım.)

VMware, VirtualPC gibi simulatörler üzerinede kurulum yapılabiliyor ama biliyorsunuz bizi bozar onlar. Kurulum Adımları:

Maceraya 00.40 civarında başladım. Minix3 ün resmi download adresinden IDE-3.1.2a.zip (298 MB) zip dosyasını indirdim. Bu dosya indiriledururken içerisinde kurulumla alakalı bilgiler içeren pdf'i okumakta yarar var.
Belgede kurulum sırasında ethernet kartının chip seti ile alakalı bilgiler sorulacağı yazıyor bu nedenle ethernet kartımın modelini SURECOM EP-320 X-V/VI 32 Bit olarak not ediyorum bi kenara. (Kartınızın modelini bilmiyorsanız belirtmeden de geçebiliyorsunuz. )
Bu sırada iso muş inmiş. Sıkıştırılmış olduğundan windows ortamındaki favori open source zip açıcımız 7zip ile bir güzel açıyoruz.
Açılan dosyayı windows ortamının vazgeçilmez cd yazma programlarından ImgBurn programı(2 mb lık efsane) ile cd ye yazıyoruz. 2$ bağışta bulunulabilir ImgBurn un sitesine girmişken ;)
Sonra cd mizi sisteme kuracağımız makinaya takıyoruz. Bios tan cd den boot etmesi için ayarı yapıyoruz. (Saat ~02.00)

Ben P2 400, 128 MB makina üzerine üzerine kurulum yaptım. Hardiskte yedeklemem gereken hiç bir şey olmadığından kurulumun otomatik olarak yaptığı partition yapısını kullandım. Bu nedenle partition lama işlemlerini çok iyi bilemeyeceğim. Siz eğer sistemininiz yanına kurulum yapmak isterseniz bu dökümanı okuyun derim.
Kurulum başladığında önce bir açıklama geliyor. Sonra klavye tipini soruyor. Enter a basarak default değeri seçebiliyorsunuz. (ing klavye diziliminde kullanacağız.) Ethernet kartınızın chipseti hakkında bilgi istiyor. 0 seçerek direk geçebilirsiniz.(no networking) Veya bunlar dışında (7) diyip daha sonra configure etmek için bırakabilirsiniz. Ben burada kendi ethernet kartı bilgilerimi göremediğime şaşırmadım çünkü windows bile tanımıyor. Her kurulumda driver arıyorum. Realtek vs çok kullanılan chipsetler mevcut.
Nasıl bir kurulum yapmak istediğimizi soruyor. Minimal kurulum veya tam kurulum yapabiliriz. Ben tam kurulum için "F" seçeneğini kullandım.
Partitioning işlemlerine sıra geldi ben direk tek disk üzerine tek partionda kurulum yapmak istediğimden default seçeneklerle ilerliyorum. disk number "0" enter region number "0" enter.
home klasoru için kaç mb yer ayırmak istediğimizi soruyor. 200 mb değer verdim.
Kilobyte cinsinden blocksize ı soruyor. Default olan 4Kb ı seçtim. Daha kısıtlı kaynaklar için değişik kb lar seçilebiliyor.
Daha sonra disk yüzeyinde bad sector taramasına başlıyor. İkinci adımda dayanamadım(epey uzun sürecekti. ilerleme grafiğine bakılırsa) ctrl+c ile işlem tamamlanmadan devam ettim.
Tüm işlemler 10 dk civarı sürdü. Ve sistemimiz hazır oldu. shutdown diyerek sistemden logout oluyoruz.(çünkü live cd üzerinden çalışıyor) ve boot komutu ile sistemi yüklediğimiz yerden boot ediyoruz. Komut "boot dXp0" şeklinde. Burada X biosta görülen sürücü numarası.(Kaçıncı disk olduğu.) p den sonrasıda partition. Bizim kurulumumuz için komut "boot d0p0" şeklinde...
shutdowndan sonra konsole tekrar düşüyoruz yukarıdaki şekilde boot komutunu verdiğimizde login: root, pass: ""(hiç bir şey yazmadan enter) şeklinde login oluyoruz. Sistemi unix komutları ile gezebiliriz. Bu esnada. Eğer paket kurmak istersek packman isminde paket yönetici var. packman komutu ile başlatabiliyoruz.
packman ile internetten veya cd den paket kurulumu yapabiliriz. cd cd sürücüde iken packman komurunu verdiğimizde cd deki paketleri tarıyor 1 seçeneği ile tüm paketleri kurabiliriz. 2 ile tüm paketler + source kodlarını kurabiliriz. 3 ile istediğimiz paketi seçebiliriz. 4 seçeneği ile ise packmanden çıkabiliyoruz.
Benim ethernet kartım tanınmadığı için kurulum cdsi ile gelen paketleri kurdum. Bu aşama biraz uzun sürebiliyor. (İlk aşamada sadece binary leri kurmanızı tavsiye ederim.) 30-45 dakika civarı sürebilir. Kurulum sırasında sistem tepkisiz görünebilir hemen resetlemeyin.(nede olsa yüksek kararlı sistem kullanıyoruz ;) ) Kurulum tamamlansın. Devamı inşallah sonraki yazıya.

*Biraz aksilikler çıktı ama sabah 06.00 civarında xwindowu gördüm. Yazının devamı gelecek. It's time to study. ;)

minix3 macerası - 1

2009-04-25T06:16:00.012+03:00

Uzun bir zamandır görüntü işleme uygulamaları geliştirip çalıştırabileceğim hızlı bir şekilde açılıp kapanabilen ama üzerinde arayüz uygulayabileceğim gömülü bir dağıtım bakıyordum.

Kendisinden çok şey öğrendiğimiz Bluekid abimiz sağolsun minix3 ü önerdi. Minix3 tam aradığım özelliklere sahip ayrıca üzerinde FLTK çalışabiliyor. ;)

Minix3 ün sitesinde sistemin hedef alanları olarak şunlar gösteriliyor.

Yüksek derecede kararlılık (reliability) istenen uygulamalarda.
Tek-chip, küçük ram e sahip, düşük güç tüketen cihazlar vs.
Embedded sistemler. (camera, cep telefonu vs.).
GPL lisansının çok kısıtlayıcı olduğu uygulamalar.(minix3 BSD lisansına sahiptir.)
Eğitim.(Üniversitelerde Operating System dersleri için.)

Linuxun gelişimini takip edememiş birisi olarak hep linux kerneli çok karmaşık, linux configurasyonu zor gelir(ya da gözümde büyütüyorum.) Uzun zamandır değişik linux dağıtımlar kullansamda alt seviyelere hiç inemedim.

Minix3 dağıtımının bir güzel özelliği de bu. Henüz bakmadım ama ~5000 satırlık bir kernel kodu varmış. Üzerinde çalıştırılabilinecek kütüphaneler ise bir hayli fazla. Bknz.

Buradan screenshot lara göz atabilirsiniz. Bir ayda yaklaşık 12000 kişi iso dosyasını siteden indirip kurulum yapıyormuş.Buranın yalancısıyım.

Ben de evdeki emektar makinayı(Intel P2 400,256 Mb Ram) bu sistemden mahrum bırakmak istemedim. Bir gece uğraşarak xwindow u görmeyi başardım. Gerçekten sisteme çok içim ısındığını söylemeliyim. Keşke derslerden projelerden biraz daha vakit olsada doya doya kurcalayabilsem. Uzun bir süre makinada kalacak gibi görülüyor. Üzerinde önce fltk yı çalıştırıp basit uygulamalar yapmak sonrada opencv çalıştırmak gibi düşüncelerim var.

* Yalnız sistemde floating point donanımsal olarak desteklenmiyormuş. Bu iş nasıl etkiler bilemiyorum. Beraber görelim..

x86 --> MIPS Cross Derleme

2009-03-13T15:30:00.005+02:00

MIPS RISC yapısına sahip bir mikroişlemci mimarisidir.Karşılaştırmak ne kadar uygun olur bilmiyorum ama komut seti Intel x86 komut setine göre oldukça temiz ve basit. Bu tasarım avantajından dolayı üniversitelerdeki bilgisayar mimarisi derslerinde genellikle MIPS mimarisi okutulur.(vikipedi). Benim yeni okulumda da öyle..

Intel x86 işlemcili makinam üzerinde bir C kodunu derleyip MIPS assembly si elde etmem gerekiyordu. Gcc ile target mimari belirtilebildiğini biliyordum, deneyince bu kadar kolay olmadığını gördüm. Bu işlemi yapabilmek için gcc nin kaynak kodunu tekrardan özel parametrelerle derlemem gerekiyormuş. (ubuntu'da başarılı olamadım...)

Daha sonra aramalarıma devam ederken princeton uni. den "lcc" isimli bir cross compiler buldum. (http://www.cs.princeton.edu/software/lcc/) Windows ve linux versiyonları mevcut. Ben windows versiyonunu kullandım. Gayet kullanışlı. MIPSR3000 için little endian ve big endian derleme yapabiliyor. (ALPHA, SPARC ve Intel x86 için de yapabiliyormuş.)

Microsoft VC++ 4 embedded da Microsoftun artık ücretsiz olarak dağıttığı çok değişik mimariler için derleme yapabileceğimiz bir araç. MIPS asm side üretiyor.

*SPIM ide atlamayalım. MIPS32 assembly sini simule etmek için harika bir araç. Windows,Linux, Unix,MacOS altında çalışıyor. (http://pages.cs.wisc.edu/~larus/spim.html)

Dalgacık.

2008-07-10T00:36:00.002+03:00

...Foton, kütlesi 0 olan; saniyede 300 000 km hızla giden en küçük dalga/parçacık tır. Gerçek cisim tanecikleri de ışık gibi aynı temele dayanıyor. Gezegenleri, yıldızları ve galaksileri oluşturan madde, enerji dalgalarından üretilen taneciklerdir.

Durağan görünen bedenimizin temelinde, sayısız foton süper hızla akış halindedir. Bedenimiz ışık hızında gerçekleşen bir dalga akışıdır.

Vücudumuz kuantum dalgalarından örülen taneciklerden, onlardan örülü atomlardan, onlardan örülü moleküllerden veonlardan örülü hücrelerden yapılmıştır. Işık hızında fotondan örülen tanecik yavaşlar; onlardan örülen atom daha da yavaşlar. Böylece örgü elementleri kayalara, kıtalara, gezegenlere, galaksilere ulaşır.

Yerinde duruyor gibi görülen her madde, bedenin derinlerinde baş döndürücü titreşimleri, fırtınaları, uçuşmaları gizlemektedir.

Maddenin aslı, süper hızlı enerji dalgalarıdır. Madde aslında enerjidir; ama, ilginçtir ki enerjinin maddi vücudu yoktur. Sobanızda yaktığınız odunlar kül olup havaya uçar. Kütlesi enerjiye yani evinizdeki maddenin atomlarında yaşanan titreşimine/ısıya dönüşür.

Odun yanarak yok olur ; karşılığında sadece çevresindeki maddeleri titreştirir.
... Evrenin tüm enerjisi vücutsuz bir zerreye sıkışabilir.

Dr. Muhammed Bozdağ. 2005

Anten Temel - 1 Empedans Uygunlaştırma

2008-06-06T16:04:00.001+03:00

Anten teorisi finalini de atlattığımıza göre işin temel pratiğe dönük eğlenceli kısımları ile biraz ilgilenelim. Eski dergileri karıştırırken Circuit Cellar Şubat 2008 sayısında Antenna Technology 101 başlıklı bir yazı ile karşılaştım. Oldukça güzel bir makale. Bahsettiğim sayıyı ayrıca Elektronik müh. ile alakalı onlarca seçme pdf e buradan ulaşabilirsiniz.

Şimdi biraz çeviri biraz yorum bazı kavramları açıklamaya çalışalım.

Empedans Uygunlaştırma: Empedans uygunlaştırma pek çok elektronik uygulamada olduğu gibi anten sistemlerinde de temel kavramlardan biridir. Basit bir şekilde açıklayacak olursak; Pil ile çalışan bir dirençli ısıtıcımız olsun :) Bu ısıtıcımızın direnci ısıtıcıdan maksimum güç alabilmek için acaba nasıl seçilmelidir. Çok küçük bir direnç üzerinden çok büyük bir akım geçer ancak üzerinde düşen gerilim az olacağından alabileceğimiz güç V*I dan oldukça küçük olacaktır. Büyük bir direnç seçildiğinde üzerinde güşen gerilim oldukça büyük olmasına rağmen akım küçülecek. En uygunu optimum noktada bir direnç seçmektir. Böyle bir uygulama için seçilmesi gereken direnç pilimizin iç direncine eşit bir ısıtıcı direnci olmalıdır. Devre teorisi dersi alan arkadaşların maximum güç transferi teoremi bunu ifade ediyor.

Yukarıda DC için açıkladığımız kurallar bazı değişikliklerle AC için de geçerlidir. AC için empedans reel ohmik direnç ve imajiner reaktans şeklinde gösterilir. Eğer devremiz kapasitif ise reaktans - değer indüktif ise + değerler alacaktır.

Kaynak empedansımızı Zs = Rs + j Xs şeklinde gösterebiliriz. Yukarıda DC ısıtıcı örneğinde olduğu gibi bu şekilde bir kaynaktan maksimum gücü transfer edebilmemiz için yükümüzün empedansınında sadece ohmik kısımdan oluşmaması, sanal ksmının da olması gerekir. Değeri ise kaynak empedansının kompleks eşleniği şeklindedir. Yani reaktans kısım ters işaretli olmalıdır. Kaynak endüktif ise yük kapasitif, kapasitif ise yük endüktif olmalıdır. Burada kaynağımızın kompleks bir değerde olduğunu bu şekilde en yüksek gücü alabilmek için bağlanabilecek yükü belirledik fakat pratikte yükümüz ve kaynak sabit bir karakteristik gösterir en uygun verimi alabilmek için bağlantı noktasında eşleştirme yaparız. Endüktif bir yük için (ör motor) hatta bir karşı reaktans yani kapasite, kapasitif bir yük için ise hatta bir endüktans bağlayarak reaktif kısmı eşleştiririz. Ohmik kısım için ise genellikle empedans uygunlaştırıcı transformatör kullanılır. Bu işlem empedans uygunlaştırma oalrak bilinir. Anten sistemlerinde empedans frekansada bağlı olarak değiştiğinden alıancak veya yayın yapılacak frekansa uygun olarak dizayn yapılması da gereklidir.

İnşallah sonraki yazımızda bunları grafik üzerinde inceler ve temel anten kavramlarına devam ederiz. Belkide yazıların sonunda wireless ağlara uzaktan bağlanabileceğimiz bir el yapımı anten yaparız :)

C++ ile curl e giriş.

2008-03-12T02:48:00.009+02:00

Curl u uzun zamandır denemek istiyorum. Ancak bugün bir şeyler elde edebildim. Konuya yabancı olanlar için curl http requestler gönderebildiğimiz bir komut satırı programı. libCurl ise curl un platformdan bağımsız kütüphanesi. Çok daha güzel anlatımı ise şöyle :)

libcurl is a free and easy-to-use client-side URL transfer library, supporting FTP, FTPS, HTTP, HTTPS, SCP, SFTP, TFTP, TELNET, DICT, LDAP, LDAPS and FILE. libcurl supports SSL certificates, HTTP POST, HTTP PUT, FTP uploading, HTTP form based upload, proxies, cookies, user+password authentication (Basic, Digest, NTLM, Negotiate, Kerberos4), file transfer resume, http proxy tunneling and more!

Aklıma gelen bazı uygulama alanları.. Console uygulamamızdan bir sayfadan veri çekebiliriz, kendi rss reader ımızı yazabiliriz. Mailimize login olup yeni ileti var mı kontrol edebiliriz.Programımızı weble ilişkilendirebiliz. Yazılım güvenliği konusunda kullanılabilir. Tabi yaramaz işlerde kullanmaya da oldukça uygun... Burada anlatıldığına göre MIT/X türevi bir lisansa sahip. açık kaynak kodlu bir proje ticari uygulamalarda da gayet güzel kullanılabiliyor. Yukarıda saydığım uygulamalar belki .net veya başka RAD(delphi, cbuilder, vb) ürünlerini kullanarak çok daha kolay şekilde gerçekleştirilebilir ama bildiğim kadarıyla kütüphane c/c++ ve php dünyasında pek yaygın.

UYGULAMA

Kütüphanenin resmi kaynağı burası. Pek çok işletim sistemi için değişik opsiyonlarla derlenmiş kütüphane dosyaları bulunabilir. Ben ms vc++ 6 ile denediğim derleme denemelerimde başarılı olamadım. Hem de hazır proje dosyası şeklinde olmasına rağmen!! Uygulamanın pek çok bağımlılığı var ve çoğu ağ teknolojileri ile alakalı benim yabancı olduğum konulular. Sizde aynı durumdaysanız libcurl-7.15.1-msvc-win32-ssl-0.9.8a-zlib-1.2.3.zip paketini önerebilirim. Paket içerisinde 32 bit win sistemler için derlenmiş ssl desteği sağlayan lib ve dll ler var. Bunu uygun yere(benim sistemimde c:\curl) kaydettikten sonra derleyicimizde gerekli linker ve header directory ayarlarını yapmamız gerekiyor. Bunlar Ms VC++ 6 için şöyle:

-Project>Settings>Link>Object/Library Modules e "libcurl.lib" ifadesini ekleyin.

-Tools>Options>Directories>Show Directory for Include Files seçili iken C:\curl\include satırını ekleyin.

-Tools>Options>Directories>Show Directory for Library Files seçili iken C:\curl\ satırını ekleyin.

Diğer derleyicilerde de menüler benzer şekilde olacaktır. Ayrıca c:/curl de libcurl_imp.lib dosyasını libcurl.lib yapıp curl klasorunde lib içine kaydetmeyi unutmayın. Dll ler ise libcurl.dll ssleay32.dll zlib1.dll libeay32.dll bunlarıda uygulamamızın exesinin bulunduğu klasore kopyalayın. simple.c Burada örnek bir uygulama var. İlk önce bu kodu denemenizi tavsiye ederim. Yaptığınız ayarlarla bu programı derlediğinizde haxx.se adresini bir http request gönderilecek ve gelen yanıt console ekranınıza yansıyacaktır. http://curl.haxx.se/lxr/source/docs/examples/ Adresinde daha gelişmiş uygulamalar bulunmakta. Kütüphane oldukça kullanışlı ve kolay. Hiç deneyimim olmamasına rağmen bir saat kadar bir süre içinde okulumuz öğrenci sistemine programla login olabildim. Belki bir sonraki yazıya konu olabilir. :)

Bismillah..

2008-02-29T02:43:00.014+02:00

...

Besmele ile başlanmayan her önemli iş noksan kalır.

Besmele ile yazı yazanın haceti kolaylaşır, Allahü teâlâ da razı olur.