闲来无事,自学了一点python,刚好最近要准备数学建模比赛,就学了一点用python建模。这里只有一个最小线性二乘的一个简单的例子,因为刚学,好不容易跑出来了图,就发篇博客纪念一下。
假定模型为y=2x+1,我们一共有100组数据(xi,yi),我们需要在这100组数据的基础上找出xi和yi的线性关系方程y=2x+1。以下是代码实现:
import scipy.linalg as la
import numpy as np
import matplotlib.pyplot as plt
m = 100
x = np.linspace(-1, 1, m)
y_exact = 1 + 2 * x
xi = x + np.random