矩阵的运算

链接: B站原视频

下面展示一些 内联代码片。

需要安装Eigen库

// An highlighted block
//
// Created by dennis on 23-2-27.
//
#include <iostream>
#include <chrono>
#include <Eigen/Core>
#include <Eigen/Dense>

using namespace std;
using namespace Eigen;


int main()
{

    Eigen::Matrix<double, 15, 2> x_vector;
    x_vector.setOnes();
    x_vector.col(1) << 75, 71, 83, 74, 73, 67, 79, 73, 88, 80, 81, 78, 73, 68, 71;
    // cout << x_vector << endl;


    Eigen::VectorXd z_vector{15, 1};
    z_vector.setZero();
    z_vector << 183, 175, 187, 185, 176, 176, 185, 191, 195, 185, 174, 180, 178, 170, 184;
    //cout << z_vector << endl;

    Matrix<double, 2, 15> x_vecotr_T = x_vector.transpose();
    Matrix<double, 2, 1> Y_Star;
    Y_Star << 100, 10;
    //cout << Y_Star << endl;
    Matrix<double, 2, 2> rate;
    rate.setZero();
    rate(0, 0) = 0.001;
    rate(1, 1) = 0.00001;
    for (int i = 0; i < 100000; ++i)
    {
        Y_Star = Y_Star - rate * x_vecotr_T * (-z_vector + x_vector * Y_Star);
        if (i %1000 ==0) cout << Y_Star  << endl << " rate: " << (rate * x_vecotr_T * (-z_vector + x_vector * Y_Star)).transpose()<<endl << endl;
    }
    //Y_Star = Y_Star - 0.0002 * (x_vecotr_T * z_vector - x_vecotr_T * x_vector * Y_Star);
    cout << Y_Star << endl;
    return 0;
}