Последние ответы: SadZhaba

Вычисляю параметры полиномиальной аппроксимации методом МНК, вот формулы:

Вот сам код, ничего сложного, но при выводе b1 считается как 29.75, что довольно странно, хотя в экселе получается 0.04.

Код:

double b1 = sumYT/sumT2;
double b0 = (sumT4*sumY - sumT2*sumYT2) / (50*sumT4 - sumT2*sumT2);
double b2 = (50*sumYT2 - sumT2*sumY)/(50*sumT4 - sumT2*sumT2);

Вот так вычисляю суммы рядов:

Код:

for(int i =0 ;i< 50;++i){
        sumY+=x[i];
        sumT+=i+1;
        sumT2+=(i+1)*(i+1);
        sumT4+=std::pow(y[i], 4.);
        sumYT2+=y[i]*y[i]*x[i];
        sumYT+=x[i]*(i+1);
}

y в формулах это t из изображения, x это y из изображения.
Где x вводится из файла через QDataStream, y просто равен (i+1), так как это просто значения времени(шкала Х). В чем может быть проблема неправильного рассчета? Мне кажется это связано с экспоненциальной формой записи больших сумм.
Вывод qDebug()
sumY 50123.9 sumT 1275 sumT2 42925 sumT4 6.56667e+07 sumYT2 4.2983e+07 sumYT 1.27728e+06
b0 1003.92 b1 29.7562 b2 -0.00167683

Весь исходник кнопки:

Код:

void MainWindow::on_pushButton_2_clicked()
{
    double b1 = sumYT/sumT2;
    double b0 = (sumT4*sumY - sumT2*sumYT2) / (50*sumT4 - sumT2*sumT2);
    double b2 = (50*sumYT2 - sumT2*sumY)/(50*sumT4 - sumT2*sumT2);
    qDebug() << "sumY " << sumY<< " sumT "<<sumT<<" sumT2 " <<sumT2<<" sumT4 "<<sumT4<<" sumYT2 "<<sumYT2<<" sumYT "<<sumYT;
    qDebug() << "b0 "<<b0<<" b1 "<<b1<<" b2 "<<b2;
    QVector<double> approx;
    for(int i =0;i<50;++i){
        approx.push_back(b0 + b1* i + b2*i*i);
    }
    ui->widget->addGraph();
    ui->widget->graph()->setPen(QPen(QColor("#00FF00")));
    ui->widget->graph(isPressed)->addData(y, approx);
    //ui->widget->rescaleAxes();
    ui->widget->replot();
    QString text_b0 = QString::number(b0, 'g', 4);
    QString text_b1 = QString::number(b1, 'g', 4);
    QString text_b2 = QString::number(b2, 'g', 4);
    ui->textBrowser->setText("y^2(t) = " + text_b0 + " + " + text_b1+"*t "+text_b2+"*t^2");
    isPressed++;
}

Данные x:

Код: