Двоичный циклический код Хэмминга

M[i+10]=(Pdop-Pls), i=i+1;

13) Выводим Poo, Pno, Pls, lgPls, переходим к шагу 8;

14) i=17;

15) Если i<=31, то Pbo=(1-P)^i, Poo=0, k=1, иначе переходим к шагу 35;

16) Выводим Pbo;

17) Если k<=2, то Poo=, иначе переходим к шагу 19;

18) k=k+1, переходим к шагу 17;

19) Pno=1-Pbo-Poo, Pls=Pno/(1-Poo),

lgPls=log10(Pls),

M[i+10]=(Pdop-Pls), i=i+1;

20) Выводим Poo, Pno, Pls, lgPls, переходим к шагу 15;

21) Если d=4, то i=11, иначе переходим к шагу 35;

22) Если i<=15, то Pbo=(1-P)^i, Poo=0, k=1, иначе переходим к шагу 28;

23) Выводим Pbo;

24) Если k<=3, то Poo=, иначе переходим к шагу 26;

25) k=k+1, переходим к шагу 24;

26) Pno=1-Pbo-Poo, Pls=Pno/(1-Poo), lgPls=log10(Pls),

M[i+10]=(Pdop-Pls), i=i+1;

27) Выводим Poo, Pno, Pls, lgPls, переходим к шагу 22;

28) i=17;

29) Если i<=31, то Pbo=(1-P)^i, Poo=0, k=1, иначе переходим к шагу 35;

30) Выводим Pbo;

31) Если k<=3, то Poo=, иначе переходим к шагу 33;

32) k=k+1, переходим к шагу 31;

33) Pno=1-Pbo-Poo, Pls=Pno/(1-Poo), lgPls=log10(Pls),

M[i+10]=(Pdop-Pls), i=i+1;

34) Выводим Poo, Pno, Pls, lgPls, переходим к шагу 29;

35) h=0, i=0;

36) Если i<=60, то переходим к шагу 37, иначе переходим к шагу 38;

37) Если M[i]>0, то h=h+1, i=i+1, иначе i=i+1 и переходим к шагу 36;

38) Выделяем память под массив Н из h элементов.

39) Если i<=60, то переходим к шагу 40, иначе переходим к шагу 41;

40) Если M[i]>0, то H[k]=M[i], k=k+1, i=i+1, иначе i=i+1 и переходим к шагу 39;

41) i=0;

42) Ищем минимальный элемент в массиве Н;

43) Если i<=60, то переходим к шагу 44, иначе переходим к шагу 50;

44) Если M[i]=минимальному элементу, то и переходим к шагу 45, иначе i=i+1 и переходим к шагу 43;

45) Если i>=0 и i<=20, то выводим (i+11,i+10)-код, иначе переходим к шагу 46;

46) Если i>=21 и i<=25, то выводим (i-10,i-14)-код, иначе переходим к шагу 47;

47) Если i>=26 и i<=40, то выводим (i-9,i-14)-код, иначе переходим к шагу 48;

48) Если i>=41 и i<=45, то выводим (i-30,i-35)-код, иначе переходим к шагу 49;

49) Если i>=46 и i<=60, то выводим (i-29,i-35)-код, иначе i=i+1 и переходим к шагу 39;

50) Выводим минимальный элемент из массива Н, как минимум разницы Рдоп-Рлс;

51) Конец.

Распечатка программы

Программа написана на языке С++.

#include <vcl.h>

#include <math.h>

#include <stdio.h>

#include <vector>

#include <algorith>

#pragma hdrstop

#include "Unit1.h"

//---------------------------------------------------------------------------

#pragma package(smart_init)

#pragma resource "*.dfm"

float P = 0.0006;

float Pdop = 0.0000002;

using namespace std;

float M[61];

vector<float>H;

char B[128];

TForm1 *Form1;

//---------------------------------------------------------------------------

fastcall TForm1::TForm1(TComponent* Owner)

: TForm(Owner)

{

}

//---------------------------------------------------------------------------

float C(int n,int m)

{float c=1.0;

for(int i=n;i>=n-m+1;i--)c*=i;

for(int i=1;i<=m;i++)c/=i;

return (int)c;

}

void fastcall TForm1::ComboBox1Select(TObject *Sender)

{int i=0, k;

double Pbo,Poo,Pno,Pls,lgPls;

AnsiString s;

ListBox1->Clear();

ListBox2->Clear();

ListBox3->Clear();

ListBox4->Clear();

ListBox5->Clear();

ListBox6->Clear();

ListBox7->Clear();

//d=2

if(ComboBox1->ItemIndex==0)

for(i=11;i<=31;i++)

{s="("+IntToStr(i)+","+IntToStr(i-1)+")";

ListBox1->Items->Add(s);

Pbo=pow(1-P,i);

sprintf(B,"%.8f",Pbo);

ListBox2->Items->Add(B);

Poo=0;

Poo=C(i,1)*pow(P,1)*pow(1-P,i-1);

sprintf(B,"%.8f",Poo);

ListBox3->Items->Add(B);

Pno=1-Pbo-Poo;

sprintf(B,"%.8f",Pno);

ListBox4->Items->Add(B);

Pls=Pno/(1-Poo);

sprintf(B,"%.8f",Pls);

ListBox5->Items->Add(B);

lgPls=log10(Pls);

sprintf(B,"%.2f",lgPls);

ListBox6->Items->Add(B);

Series1->AddXY(i,lgPls,s,clRed);

M[i-11]=(Pdop-Pls);

}

//d=3

if(ComboBox1->ItemIndex==1)

{for(i=11;i<=15;i++)

{s="("+IntToStr(i)+","+IntToStr(i-4)+")";

ListBox1->Items->Add(s);

Pbo=pow(1-P,i);

sprintf(B,"%.8f",Pbo);

ListBox2->Items->Add(B);

Poo=0;

for(k=1;k<=2;k++)

Poo+=C(i,k)*pow(P,k)*pow(1-P,i-k);

sprintf(B,"%.8f",Poo);

ListBox3->Items->Add(B);

Pno=1-Pbo-Poo;

sprintf(B,"%.8f",Pno);

ListBox4->Items->Add(B);

Pls=Pno/(1-Poo);

sprintf(B,"%.8f",Pls);

ListBox5->Items->Add(B);

lgPls=log10(Pls);

sprintf(B,"%.2f",lgPls);

ListBox6->Items->Add(B);

Series2->AddXY(i,lgPls,s,clLime);

M[i+10]=(Pdop-Pls);

}

for(i=17;i<=31;i++)

{s="("+IntToStr(i)+","+IntToStr(i-5)+")";

ListBox1->Items->Add(s);

Pbo=pow(1-P,i);

sprintf(B,"%.8f",Pbo);

ListBox2->Items->Add(B);

Poo=0;

for(k=1;k<=2;k++)

Poo+=C(i,k)*pow(P,k)*pow(1-P,i-k);

sprintf(B,"%.8f",Poo);

ListBox3->Items->Add(B);

Pno=1-Pbo-Poo;

sprintf(B,"%.8f",Pno);

ListBox4->Items->Add(B);

Pls=Pno/(1-Poo);

sprintf(B,"%.8f",Pls);

ListBox5->Items->Add(B);

lgPls=log10(Pls);

sprintf(B,"%.2f",lgPls);

ListBox6->Items->Add(B);

Series2->AddXY(i,lgPls,s,clLime);

M[i+9]=(Pdop-Pls);

}

}

//d=4

if(ComboBox1->ItemIndex==2)

{for(i=11;i<=15;i++)

{s="("+IntToStr(i)+","+IntToStr(i-5)+")";

ListBox1->Items->Add(s);

Pbo=pow(1-P,i);

sprintf(B,"%.8f",Pbo);

ListBox2->Items->Add(B);

Poo=0;

for(k=1;k<=3;k++)

Poo+=C(i,k)*pow(P,k)*pow(1-P,i-k);

sprintf(B,"%.8f",Poo);

ListBox3->Items->Add(B);

Pno=1-Pbo-Poo;

sprintf(B,"%.8f",Pno);

ListBox4->Items->Add(B);

Pls=Pno/(1-Poo);

sprintf(B,"%.8f",Pls);

ListBox5->Items->Add(B);

lgPls=log10(Pls);

sprintf(B,"%.2f",lgPls);

ListBox6->Items->Add(B);

Series3->AddXY(i,lgPls,s,clYellow);

M[i+30]=(Pdop-Pls);

}

for(i=17;i<=31;i++)

{s="("+IntToStr(i)+","+IntToStr(i-6)+")";

ListBox1->Items->Add(s);

 Скачать реферат