Дисперсия

[Vago]

Да, там опечатка, конечно, во второй формуле.

...

159.6, 181.6, 151.0, 143.6,..
Равномерное распределение питоновский uniform() в теории должен давать...

[justso]

Сообщение от Vago

159.6, 181.6, 151.0, 143.6,..

Вот и у меня как то так. А говорят надо 300-500 примерно.

[Vago]

Ну, давайте тех, "кто говорит", сюда

...

Генератор имитирует равномерное распределение. Возведите 11.25 в квадрат. Сколько получится? Неужели 500?..

[CesaR_JC]

Ну может я чтото не так сделал.

[justso]

Ну тогда сейчас в С++ сделаю функцию вывода х и у в текстовый файл, потом на Matlab сделаю запрос этого файла для графиков. Единственное что могу сообразить.

Хотя может можно как-то проще?

Дополнено:

Такс..... Вроде разобралась.

Вот код для Dy и вывода Х У на экран, написан на С++ (не компилируется - попробуйте на другом компайлере)

#include <iostream>
#include <ctime>
#include <fstream>
#include <cmath>

using namespace std;

int main ()
{
    
    
double X, Y, c, d, k, f;
int i;
double  N, Xrand[2][130], Yrand[2][130], b, a, M1, My, D1, Dy ;
    srand(time(NULL));
    ofstream outfile;
    outfile.open("XY.txt");
    
X = 324;
Y = 225;
a = X*0.2;
b = Y*0.2;
N = 128;

for (i=1; i<=128; i++)
{
    k = (rand()%9);
    f = (rand()%11);
c = (rand()%64) + k/10.0000;
d = rand()%44 + f/10.0000;
Xrand[1][i] = X - (a/2.00) + c;
Yrand[1][i] = Y - (b/2.00) + d;
}

M1 = 0;
for (i=1; i<=128; i++)
{
M1 = M1 + Yrand[1][i];

}

My = M1/N;
D1 = 0;
/*cout << M1 <<endl;
cout << My;*/

for (i=1; i<=128; i++)
{
D1 = D1 + (Yrand[1][i] - My)*(Yrand[1][i] - My);

}

Dy = D1/(127.0000);

cout << "Dy = " << Dy;
cout << "\n\n X         Y " << endl;

for (i=1; i<=128; i++)
{
cout << Xrand[1][i] << "      " << Yrand[1][i] << endl;
outfile << Xrand[1][i] << " " << Yrand[1][i] << endl;
}
    
system ("pause");
return 0;
}

Он большой и неаккуратный, но правильный.

Я не волшебник, я только учусь. (С)

Этот код выдаст вам текстовый файл со значениями Х У. В Матлабе запустите это:

fid = fopen('XY.txt');
XY_plot = fscanf(fid, '%g %g', [2 inf]);
fclose(fid);
XY_plot = XY_plot';


plot(XY_plot(:,1),XY_plot(:,2),'b') ;
Title('X vs. Y');
Xlabel('X');
Ylabel('Y');


!!! Убедитесь что папка в матлабе указывает на папку с файлом ХУ.txt


Если что то все же не так, стучитесь в асю.

[CesaR_JC]

За программу в С++ миллион спасибо.

Как в МАТЛАБ сформировать выборку из значений в диапазоне 324+-10%?