1

Здесь я хочу сравнить кривую мощности двух разных турбин. Предположим, у вас есть две таблицы tb1 для турбины 1 (есть X = скорость И X1 = значения мощности) и tb2 для турбины 2 (имеют y = скорость и y1 = значения мощности). Теперь с помощью данных tb1; Я подсчитал binEdge (с помощью X), а затем сделали биннинг (используя x1), так что с помощью станда и средних значений мы знаем распределение, как вы можете видеть код ниже:Сравнение и оценка новых данных кривых через использование контрольные кривые средние и значение STD?

binEdges = linspace(min(X),max(X),numBins); 
[N,binEdges] = histcounts(X, binEdges); 

for i = 1:numBins 
flagBinMembers = (binEdges == i); 
binMembers = X1(flagBinMembers); 
BinMean(i) = mean(binMembers); 
BinStd(i)  = std(binMembers); 
end 

Сейчас в ТВ2, я хотите проверить значения скорости или y распределения выше вычисленных std и средних значений. Как только я узнаю местоположение или в котором распределение y или скорость принадлежит, я хочу оценить его значение мощности, выполнив некоторый вероятностный тест. Эта мощность может быть проверена с доступной мощностью y1, чтобы увидеть, что оценочные значения мощности для погоды верны или нет.

+0

Что вы подразумеваете под «вероятностным тестом», в частности? –

+0

@OliverCharlesworth chi square test с функцией распределения вероятности –

+0

Вы пытались найти «Matlab chi squared test»? –

ответ

1

Если у вас есть статистика и Machine Learning Toolbox, и все вы хотите оценить параметры binMembers, а затем посмотреть, если X1 происходит от того же распределения:

numBins = 5; 
binEdges = linspace(min(X),max(X),numBins); 
N = histcounts(X, binEdges); 

for k = 1:numBins-1 
    binMembers = X(X>=binEdges(k) & X<binEdges(k+1)); 
    if ~isempty(binMembers) 
     is_X1_fit = ~chi2gof(X1,'CDF',makedist('Normal',... 
      'mu',mean(binMembers),'sigma',std(binMembers))) 
    else 
      disp('bin is empty') 
    end 
end 

is_X1_fit возвратит 1, если X1 принадлежит binMembers и 0 в противном случае.

+0

«CDF» - это аргумент ' chi2gof', и он используется для сравнения дистрибутивов, так же хорошо, как вы можете сделать это с помощью PDF, но PDF не является аргументом 'chi2gof'. В любом случае, это просто, чтобы «объект распределения вероятности» перешел на 'chi2gof'. – EBH

+0

Это еще один вопрос? ваши «данные» - это X или некоторые другие данные? Что вы хотите сравнить с X1? – EBH

+0

Мои извинения, если вы не смогли понять проблему. Поэтому я написал образец кода с более подробным описанием моей проблемы. –