我要评论最小二乘法介绍✨
最小二乘法（least squares method）是一种常见的数学优化技术，广泛应用于数据拟合、回归分析和参数估计等领域。其目标是通过最小化残差平方和来找到一组参数，使得模型预测值与观测值之间的差异最小化。
最小二乘法的原理✨
线性回归模型将因变量 （y） 与至少一个自变量 （x） 之间的关系建立为：
在 ols 方法中，我们必须选择一个b1和b0的值，以便将 y 的实际值和拟合值之间的差值的平方和最小化。
平方和的公式如下：
我们可以把它看成是一个关于b1和b0的函数，分别对b1和b0求偏导，然后让偏导等于0，就可以得到最小平方和对应的b1和b0的值。
先说结果，斜率最后推导出来如下所示：
截距推导出来结果如下：
don't worry about that，慢慢推导总是可以弄明白的（不感兴趣可以直接略过）：
用c#实现最小二乘法✨
创建数据点✨
首先创建想要拟合的数据点：
ndarray? x, y;
x,y为全局变量。
//使用numsharp创建线性回归的数据集
x = np.arange(0, 10, 0.2);
y = 2 * x + 3 + np.random.normal(0, 3, x.size);
使用到了numsharp，需要为项目添加numsharp包：
x = np.arange(0, 10, 0.2);
的意思是x从0增加到10（不包含10），步长为0.2：
np.random.normal(0, 3, x.size);
的意思是生成了一个均值为0，标准差为3，数量与x数组长度相同的正态分布随机数数组。这个数组被用作线性回归数据的噪声。
使用oxyplot画散点图✨
oxyplot是一个用于在.net应用程序中创建数据可视化图表的开源图表库。它提供了丰富的功能和灵活性，使开发者能够轻松地在其应用程序中集成各种类型的图表，包括折线图、柱状图、饼图等。
添加oxyplot.windowsforms包：
将plotview控件添加到窗体设计器上：
// 初始化散点图数据
var scatterseries = new scatterseries
{
markertype = markertype.circle,
markersize = 5,
markerfill = oxycolors.blue
};
表示标志为圆形，标志用蓝色填充，标志的大小为5。
for (int i = 0; i < x.size; i++)
{
scatterseries.points.add(new scatterpoint(x[i], y[i]));
}
添加数据点。
plotmodel? plotmodel;
将plotmodel设置为全局变量。
// 创建 plotmodel
plotmodel = new plotmodel()
{
title = "散点图"
};
plotmodel.series.add(scatterseries);
// 将 plotmodel 设置到 plotview
plotview1.model = plotmodel;
这样就成功绘制了散点图，效果如下所示：
使用最小二乘法拟合数据点✨
double a = 0;
double c = 0;
double x_mean = x?.mean();
double y_mean = y?.mean();
//计算a和c
for(int i = 0; i < x?.size; i++)
{
a += (x[i] - x_mean) * (y?[i] - y_mean);
c += (x[i] - x_mean) * (x[i] - x_mean);
}
//计算斜率和截距
double m = a / c;
double b = y_mean - m * x_mean;
//拟合的直线
var y2 = m * x + b;
套用公式就可以，a表示上面斜率公式的上面那部分，c表示上面斜率公式的下面那部分。
double x_mean = x?.mean();
double y_mean = y?.mean();
计算x与y的平均值。
使用oxyplot画拟合出来的直线✨
//画这条直线
var lineseries = new lineseries
{
points = { new datapoint(x?[0], y2[0]), new datapoint(x?[-1], y2[-1]) },
color = oxycolors.red
};
// 创建 plotmodel
plotmodel?.series.add(lineseries);
// 为图表添加标题
if (plotmodel != null)
{
plotmodel.title = $"拟合的直线 y = {m:0.00}x + {b:0.00}";
}
// 刷新 plotview
plotview1.invalidateplot(true);
points = { new datapoint(x?[0], y2[0]), new datapoint(x?[-1], y2[-1]) },
画直线只要添加两个点就好了x?[0], y2[0]表示x和y的第一个点，x?[-1], y2[-1])表示x和y的最后一个点，使用了numsharp的切片语法。
画出来的效果如下所示：
c#实现的全部代码：
using numsharp;
using oxyplot.series;
using oxyplot;
namespace olsregressiondemousingwinform
{
public partial class form1 : form
{
ndarray? x, y;
plotmodel? plotmodel;
public form1()
{
initializecomponent();
}
private void button1_click(object sender, eventargs e)
{
//使用numsharp创建线性回归的数据集
x = np.arange(0, 10, 0.2);
y = 2 * x + 3 + np.random.normal(0, 3, x.size);
// 初始化散点图数据
var scatterseries = new scatterseries
{
markertype = markertype.circle,
markersize = 5,
markerfill = oxycolors.blue
};
for (int i = 0; i < x.size; i++)
{
scatterseries.points.add(new scatterpoint(x[i], y[i]));
}
// 创建 plotmodel
plotmodel = new plotmodel()
{
title = "散点图"
};
plotmodel.series.add(scatterseries);
// 将 plotmodel 设置到 plotview
plotview1.model = plotmodel;
}
private void button2_click(object sender, eventargs e)
{
double a = 0;
double c = 0;
double x_mean = x?.mean();
double y_mean = y?.mean();
//计算a和c
for(int i = 0; i < x?.size; i++)
{
a += (x[i] - x_mean) * (y?[i] - y_mean);
c += (x[i] - x_mean) * (x[i] - x_mean);
}
//计算斜率和截距
double m = a / c;
double b = y_mean - m * x_mean;
//拟合的直线
var y2 = m * x + b;
//画这条直线
var lineseries = new lineseries
{
points = { new datapoint(x?[0], y2[0]), new datapoint(x?[-1], y2[-1]) },
color = oxycolors.red
};
// 创建 plotmodel
plotmodel?.series.add(lineseries);
// 为图表添加标题
if (plotmodel != null)
{
plotmodel.title = $"拟合的直线 y = {m:0.00}x + {b:0.00}";
}
// 刷新 plotview
plotview1.invalidateplot(true);
}
}
}
用python实现最小二乘法✨
import numpy as np
import matplotlib.pyplot as plt
# 用最小二乘法拟合 y = mx + b
# 设置随机数种子以保证结果的可复现性
np.random.seed(0)
# 生成一个在[0, 10]区间内均匀分布的100个数作为x
x = np.linspace(0, 10, 100)
# 生成y，y = 2x + 噪声，其中噪声是[0, 10)之间的随机整数
y = 2 * x + 5 + np.random.randint(0, 10, size=100)
# 计算x和y的均值
x_mean = np.mean(x)
y_mean = np.mean(y)
a = 0
c = 0
for i in range(x.shape[0]):
a += (x[i] - x_mean) * (y[i] - y_mean)
c += (x[i] - x_mean) ** 2
# 计算斜率和截距
m = a / c
b = y_mean - m * x_mean
# 画这条直线
y2 = m * x + b
plt.plot(x, y2, color='red')
# 画数据点
plt.scatter(x, y)
plt.xlabel('x')
plt.ylabel('y')
plt.title(f'y = {m:.2f}x + {b:.2f}')
plt.show()
运行效果如下所示：
总结✨
参考✨
1、
2、
发表评论