C#直线的最小二乘法线性回归运算实例

2024-12-02 04:35:39

本文实例讲述了C#直线的最小二乘法线性回归运算方法。分享给大家供大家参考。具体如下：

1.Point结构

在编写C#窗体应用程序时，因为引用了System.Drawing命名空间，其中自带了Point结构，本文中的例子是一个控制台应用程序，因此自己制作了一个Point结构

/// <summary>
/// 二维笛卡尔坐标系坐标
/// </summary>
public struct Point
{
  public double X;
  public double Y;
  public Point(double x = 0, double y = 0)
  {
    X = x;
    Y = y;
  }
}

2.线性回归

/// <summary>
/// 对一组点通过最小二乘法进行线性回归
/// </summary>
/// <param name="parray"></param>
public static void LinearRegression(Point[] parray)
{
  //点数不能小于2
  if (parray.Length < 2)
  {
    Console.WriteLine("点的数量小于2，无法进行线性回归");
    return;
  }
  //求出横纵坐标的平均值
  double averagex = 0, averagey = 0;
  foreach (Point p in parray)
  {
    averagex += p.X;
    averagey += p.Y;
  }
  averagex /= parray.Length;
  averagey /= parray.Length;
  //经验回归系数的分子与分母
  double numerator = 0;
  double denominator = 0;
  foreach (Point p in parray)
  {
    numerator += (p.X - averagex) * (p.Y - averagey);
    denominator += (p.X - averagex) * (p.X - averagex);
  }
  //回归系数b（Regression Coefficient）
  double RCB = numerator / denominator;
  //回归系数a
  double RCA = averagey - RCB * averagex;
  Console.WriteLine("回归系数A： " + RCA.ToString("0.0000"));
  Console.WriteLine("回归系数B： " + RCB.ToString("0.0000"));
  Console.WriteLine(string.Format("方程为： y = {0} + {1} * x",
    RCA.ToString("0.0000"), RCB.ToString("0.0000")));
  //剩余平方和与回归平方和
  double residualSS = 0;  //（Residual Sum of Squares）
  double regressionSS = 0; //（Regression Sum of Squares）
  foreach (Point p in parray)
  {
    residualSS +=
      (p.Y - RCA - RCB * p.X) *
      (p.Y - RCA - RCB * p.X);
    regressionSS +=
      (RCA + RCB * p.X - averagey) *
      (RCA + RCB * p.X - averagey);
  }
  Console.WriteLine("剩余平方和： " + residualSS.ToString("0.0000"));
  Console.WriteLine("回归平方和： " + regressionSS.ToString("0.0000"));
}

3.Main函数调用

static void Main(string[] args)
{
  //设置一个包含9个点的数组
  Point[] array = new Point[9];
  array[0] = new Point(0, 66.7);
  array[1] = new Point(4, 71.0);
  array[2] = new Point(10, 76.3);
  array[3] = new Point(15, 80.6);
  array[4] = new Point(21, 85.7);
  array[5] = new Point(29, 92.9);
  array[6] = new Point(36, 99.4);
  array[7] = new Point(51, 113.6);
  array[8] = new Point(68, 125.1);
  LinearRegression(array);
  Console.Read();
}

4.运行结果

希望本文所述对大家的C#程序设计有所帮助。

C# 基础之运算符

1.条件运算符条件运算符(?:)也称为三元(目)运算符,是if...else结构的简化形式,可以嵌套使用. 复制代码代码如下: int x = 1; string s = x + ""; ; s += (x == 1 ? "man" : "men"); Console.WriteLine(s);//输出1man 2.checked和unchecked 复制代码代码如下: byte b = 255; { b++; }
C#中的除法运算符与VB.NET中的除法运算符

VB.NET中的除法运算符有两个:/(浮点除法).\(整数除法) C#中的除法运算符只有一个:/(除法) VB.NET中的除法运算符与C#中的除法运算符存在很大的差异,使用时注意区分. 关于VB.NET中的除法运算符的介绍(摘自MSDN): /(浮点除法):将两个数相除并返回以浮点数表示的结果. 所得结果的数据类型取决于操作数的类型. 下表显示如何确定结果的数据类型. 操作数数据类型结果数据类型两个表达式都是整数数据类型(SByte.Byte.Short.UShort.Integer.UIn
C#数据结构与算法揭秘二线性结构

上文对数据结构与算法,有了一个简单的概述与介绍,这篇文章,我们介绍一中典型数据结构--线性结构. 什么是线性结构,线性结构是最简单.最基本.最常用的数据结构.线性表是线性结构的抽象(Abstract), 线性结构的特点是结构中的数据元素之间存在一对一的线性关系. 这种一对一的关系指的是数据元素之间的位置关系,即: (1)除第一个位置的数据元素外,其它数据元素位置的前面都只有一个数据元素: (2)除最后一个位置的数据元素外,其它数据元素位置的后面都只有一个元素.也就是说,数据元素是一个接一个的排
C#图像线性变换的方法

本文实例讲述了C#图像线性变换的方法.分享给大家供大家参考.具体如下: //定义图像线性运算函数(y=kx+v) private static Bitmap LinearOP(Bitmap a, double k, double v) { Rectangle rect = new Rectangle(0, 0, a.Width, a.Height); System.Drawing.Imaging.BitmapData srcData = a.LockBits(rect, System.Drawi
C#中矩阵运算方法实例分析

本文实例讲述了C#中矩阵运算方法.分享给大家供大家参考.具体分析如下: 一.测试环境: 主机:XP 开发环境:VS2008 二.功能: 在C#中实现矩阵运算三.源代码: using System; using System.Collections.Generic; using System.ComponentModel; using System.Data; using System.Drawing; using System.Linq; using System.Text; using Sy
C#线性渐变画刷LinearGradientBrush用法实例

本文实例讲述了C#线性渐变画刷LinearGradientBrush用法.分享给大家供大家参考.具体如下: using System; using System.Collections.Generic; using System.ComponentModel; using System.Data; using System.Drawing; using System.Text; using System.Windows.Forms; using System.Drawing.Drawing2D;
C#排序算法的比较分析

本文实例分析了C#的各种排序算法.分享给大家供大家参考.具体分析如下: 首先通过图表比较不同排序算法的时间复杂度和稳定性. 排序方法平均时间最坏情况最好情况辅助空间稳定性直接插入排序 O(n2) O(n2) O(n) O(1) 是冒泡排序 O(n2) O(n2) O(n) O(1) 是简单选择排序 O(n2) O(n2) O(n2) O(1) 是希尔排序 - O(nlog2n)~O(n2) O(nlog2n)~O(n2) O(1) 否快速排序 O(nlog2n) O(n2)
C#运算符重载用法实例分析

本文实例讲述了C#运算符重载用法.分享给大家供大家参考.具体分析如下: public class Plane { public virtual double TopSpeed() { return 300.0D;} public static bool operator>(Plane one, Plane two) { return one.TopSpeed() > two.TopSpeed(); } public static bool operator<(Plane one, Pla
c#泛型学习详解创建线性链表

术语表 generics:泛型type-safe:类型安全collection: 集合compiler:编译器run time:程序运行时object: 对象.NET library:.Net类库value type: 值类型box: 装箱unbox: 拆箱implicity: 隐式explicity: 显式linked list: 线性链表node: 结点indexer: 索引器泛型是什么? 很多人觉得泛型很难理解.我相信这是因为他们通常在了解泛型是用来解决什么问题之前,就被灌输了大量的理论
C#直线的最小二乘法线性回归运算实例

本文实例讲述了C#直线的最小二乘法线性回归运算方法.分享给大家供大家参考.具体如下: 1.Point结构在编写C#窗体应用程序时,因为引用了System.Drawing命名空间,其中自带了Point结构,本文中的例子是一个控制台应用程序,因此自己制作了一个Point结构 /// <summary> /// 二维笛卡尔坐标系坐标 /// </summary> public struct Point { public double X; public double Y; public
python 生成器协程运算实例

一.yield运行方式我们定义一个如下的生成器: def put_on(name): print("Hi {}, 货物来了,准备搬到仓库!".format(name)) while True: goods = yield print("货物[%s]已经被%s搬进仓库了."%(goods,name)) p = put_on("bigberg") #输出 G:\python\install\python.exe G:/python/untitled
C语言实现opencv提取直线、轮廓及ROI实例详解

一.Canny检测轮廓在上一篇文章中有提到sobel边缘检测,并重写了soble的C++代码让其与matlab中算法效果一致,而soble边缘检测是基于单一阈值的,我们不能兼顾到低阈值的丰富边缘和高阈值时的边缘缺失这两个问题.而canny算子则很好的弥补了这一不足,从目前看来,canny边缘检测在做图像轮廓提取方面是最优秀的边缘检测算法. canny边缘检测采用双阈值值法,高阈值用来检测图像中重要的.显著的线条.轮廓等,而低阈值用来保证不丢失细节部分,低阈值检测出来的边缘更丰富,但是很多边缘并
关于Python中定制类的比较运算实例

Python中的比较运算有几种:小于.小于等于.等于.大于等于.大于.不等于等.如果我们的数据对象具有明确的物理含义,比如说数值是带有数字与物理单位的字符串组合,那么进行大小比较的时候就可以做此定制. 下面针对6种比较方法中的一种来进行定制示范,选择小于判断来做尝试. 写如下示范代码: class myClass: def __init__(self,value): self.value = value def __lt__(self,other): return int(self.value[
python 负数取模运算实例

举例: 340%60 = 40 ,怎么来的? 340 - 60*5 = 40 340 - (比340小的那个可以被60整除的正整数) =. 40 如果是负数: -340%60 = -340 - (比-340小的那个可以被60整除的负整数) = -340 - (-360) = 20 如图: 也可以换个思路想, -340%60 =(比340大的那个可以被60整除的正整数) - 340 = 360 - 340 = 20 举例: 1%3 = 1 -1%3 = 2 补充知识:python中对负数求余的计算
Python-openCV开运算实例

我就废话不多说了,大家还是直接看代码吧~ #coding=utf-8 import cv2 import numpy as np img=cv2.imread('timg.jpeg',cv2.CV_LOAD_IMAGE_GRAYSCALE) #创建矩形结构单元 g=cv2.getStructuringElement(cv2.MORPH_RECT,(9,9)) #形态学处理,开运算 img_open=cv2.morphologyEx(img,cv2.MORPH_OPEN,g) img_hat=im
Python实现的简单线性回归算法实例分析

本文实例讲述了Python实现的简单线性回归算法.分享给大家供大家参考,具体如下: 用python实现R的线性模型(lm)中一元线性回归的简单方法,使用R的women示例数据,R的运行结果: > summary(fit) Call: lm(formula = weight ~ height, data = women) Residuals: Min 1Q Median 3Q Max -1.7333 -1.1333 -0.3833 0.7417 3.116
Go语言中的指针运算实例分析

本文实例分析了Go语言中的指针运算方法.分享给大家供大家参考.具体分析如下: Go语言的语法上是不支持指针运算的,所有指针都在可控的一个范围内使用,没有C语言的*void然后随意转换指针类型这样的东西.最近在思考Go如何操作共享内存,共享内存就需要把指针转成不同类型或者对指针进行运算再获取数据. 这里对Go语言内置的unsafe模块做了一个实验,发现通过unsafe模块,Go语言一样可以做指针运算,只是比C的方式繁琐一些,但是理解上是一样的. 下面是实验代码: 复制代码代码如下: packag
Python实现短网址ShortUrl的Hash运算实例讲解

本文实例讲述了Python实现短网址ShortUrl的Hash运算方法.分享给大家供大家参考.具体如下: shorturl实现常见的做法都是将原始Url存储到数据库,由数据库返回一个对应ID. 以下要实现的是不用数据库支持就对原始URL进行shorturl hash.说到这里我们很容易想到MD5,固定长度,冲突概率小,但是32个字符,太长?我们以MD5为基础,将其字符缩短,同时要保证一定数量范围内hash不会冲突. 我们分成两个步骤来实现. 第一步算法: ① 将长网址用md5算法生成32位签名串
Python3.5集合及其常见运算实例详解

本文实例讲述了Python3.5集合及其常见运算.分享给大家供大家参考,具体如下: 1.集合的定义:集合是一个无序的.无重复的数据的数据组合. 2.集合的特征: (1)去除重复元素:将一个列表变成集合就可实现去重. (2)测试关系:测试两组数据之间交集.并集.差集等关系. 3.集合常用的操作: (1)去重复:将列表变成集合,形式:集合=set(列表) list_1 = [1,3,5,3,6,8,9,6,8,1] list_1 = set(list_1) print(list_1,type(lis

C#直线的最小二乘法线性回归运算实例

相关推荐

随机推荐