最小二乘法多项式拟合 C语言实现（引用）（附echarts画图代码）_c语言多项式拟合算法-程序员宅基地

技术标签：算法 C++

C++多阶拟合（附echarts画图代码）

细微修改，更通用

细微修改，更通用

原文章 https://blog.csdn.net/sunshineacm/article/details/79069561

/*
    本实验根据数组x[], y[]列出的一组数据，用最小二乘法求它的拟合曲线。
    近似解析表达式为y = a0 + a1 * x + a2 * x^2 + a3 * x^3;（三阶）
*/
#include <stdio.h>
#include <math.h>

#define maxn 5
#define rank_ 3 //阶数

int main()
{
    

    // double x[maxn] = {0, 0.25, 0, 5, 0.75};
    // double y[maxn] = {1, 1.283, 1.649, 2.212, 2.178};
    double x[maxn] = {
    0, 5, 10, 15, 20};//, 25, 30, 35, 40, 45, 50, 55};
    double y[maxn] = {
    0, 1.27, 2.16, 2.86, 3.44};//, 3.87, 4.15, 4.37, 4.51, 4.58, 4.02, 4.64};

    double atemp[2 * rank_ + 1] = {
    0}, b[rank_ + 1] = {
    0}, a[rank_ + 1][rank_ + 1];
    int i, j, k;

    for (i = 0; i < maxn; i++)
    {
     
        for (j = 0; j < rank_ + 1; j++)
        {
    
            b[j] += pow(x[i], j) * y[i];
        }

        for (j = 0; j < (2 * rank_ + 1); j++)
        {
    
            atemp[j] += pow(x[i], j);
        }
    }
    atemp[0] = maxn;
    for (i = 0; i < rank_ + 1; i++)
    {
     //构建线性方程组系数矩阵，b[]不变
        k = i;
        for (j = 0; j < rank_ + 1; j++)
            a[i][j] = atemp[k++];
    }
    //以下为高斯列主元消去法解线性方程组
    for (k = 0; k < rank_ + 1 - 1; k++)
    {
     //n - 1列
        int column = k;
        double mainelement = a[k][k];

        for (i = k; i < rank_ + 1; i++) //找主元素
            if (fabs(a[i][k]) > mainelement)
            {
    
                mainelement = fabs(a[i][k]);
                column = i;
            }
        for (j = k; j < rank_ + 1; j++)
        {
     //交换两行
            double atemp = a[k][j];
            a[k][j] = a[column][j];
            a[column][j] = atemp;
        }
        double btemp = b[k];
        b[k] = b[column];
        b[column] = btemp;

        for (i = k + 1; i < rank_ + 1; i++)
        {
     //消元过程
            double Mik = a[i][k] / a[k][k];
            for (j = k; j < rank_ + 1; j++)
                a[i][j] -= Mik * a[k][j];
            b[i] -= Mik * b[k];
        }
    }
    b[rank_ + 1 - 1] /= a[rank_ + 1 - 1][rank_ + 1 - 1]; //回代过程
    for (i = rank_ + 1 - 2; i >= 0; i--)
    {
    
        double sum = 0;
        for (j = i + 1; j < rank_ + 1; j++)
            sum += a[i][j] * b[j];
        b[i] = (b[i] - sum) / a[i][i];
    }
    //高斯列主元消去法结束，输出
    printf("P(x) = %f%+fx%+fx^2%+fx^3\n\n", b[0], b[1], b[2], b[3]);
    return 0;
}

#Echarts做的图（要先下echarts.min.js，记得勾选散点图)
echarts画的模拟图
html部分

<!DOCTYPE html>
<html lang="zh-CN">
	<head>
		<title>test</title>
		<meta charset="utf-8" />
		<meta name="renderer" content="webkit" />
	</head>
	<body>
		<div id="Graph" style="width: 800px;height: 500px;">123</div>
	</body>
	<script src="echarts.min.js?ver="></script>
	<script src="index.js?ver="></script>
</html>

js部分

draw();

function f1(x) {
    //二阶拟合
	let y = 0.230467 + 0.203691 * x - 0.002381 * x * x;
	return y;
}

function f3(x) {
    //三阶
	let y = 0.017839 + 0.263399 * x - 0.005216 * x * x + 0.000034 * x * x * x;
	return y;
}

function f4(x) {
    //二阶，数据减少为5个
	return (0.027714 + 0.259114 * x - 0.004486 * x * x) // + 0.000034 * x * x * x;
}

function draw() {
    

	let x = [0, 5, 10, 15, 20, 25, 30, 35, 40, 45, 50, 55];
	let y = [0, 1.27, 2.16, 2.86, 3.44, 3.87, 4.15, 4.37, 4.51, 4.58, 4.02, 4.64];

	let data = [];
	let data2 = []
	let data3 = []
	let data4 = []

	for (let idx in x) {
    
		let a = [];
		a.push(x[idx]);
		a.push(y[idx]);
		data.push(a);


		let b = [];
		b.push(x[idx]);
		b.push(f1(x[idx]))
		data2.push(b);

		let c = [];
		c.push(x[idx]);
		c.push(f3(x[idx]))
		data3.push(c);

		let c4 = [];
		c4.push(x[idx]);
		c4.push(f4(x[idx]))
		data4.push(c4);

	}

	let option = {
    
		tooltip: {
    
			trigger: 'axis',
			axisPointer: {
     // 坐标轴指示器，坐标轴触发有效
				type: 'shadow' // 默认为直线，可选为：'line' | 'shadow'
			}
		},
		legend: {
    
			data: ['原数据', '二阶', '三阶', '二阶5数据', "5"]
		},
		xAxis: {
    },
		yAxis: {
    },
		series: [{
    
				symbolSize: 10,
				data: data,
				type: 'scatter',
				name: '原数据'
			},

			{
    
				symbolSize: 10,
				data: data2,
				type: 'scatter',
				name: '二阶'
			},

			{
    
				symbolSize: 10,
				data: data3,
				type: 'scatter',
				name: '三阶'
			},
			{
    
				symbolSize: 10,
				data: data4,
				type: 'scatter',
				name: '二阶5数据'
			},

		]
	};

	// 基于准备好的dom，初始化echarts实例
	var myChart = echarts.init(document.getElementById('Graph'));
	myChart.setOption(option);
}

本文链接：https://blog.csdn.net/u012565113/article/details/106996382

原作者删帖不实内容删帖广告或垃圾文章投诉

智能推荐

《第一行代码》（第二版）广播的问题及其解决_代码里的广播错误-程序员宅基地

文章浏览阅读2.6k次，点赞5次，收藏13次。1）5.2.1弹出两次已连接或者未连接这是因为你同时打开了流量和WiFi，他就会发出两次广播。2）5.3.1中发送自定义广播问题标准广播未能弹出消息：Intent intent=new Intent("com.example.broadcasttest.MY_BROADCAST");sendBroadcast(intent);上述已经失效了。修改：Intent intent=new Intent("com.example.broadcasttest...._代码里的广播错误

K8s 学习者绝对不能错过的最全知识图谱（内含 58个知识点链接）-程序员宅基地

文章浏览阅读249次。作者 |平名阿里服务端开发技术专家导读：Kubernetes 作为云原生时代的“操作系统”，熟悉和使用它是每名用户的必备技能。本篇文章概述了容器服务 Kubernet..._k8知识库

TencentOS3.1安装PHP+Nginx+redis测试系统_tencentos-3.1-程序员宅基地

文章浏览阅读923次。分别是etc/pear.conf，etc/php-fpm.conf， etc/php-fpm.d/www.conf，lib/php.ini。php8安装基本一致，因为一个服务期内有2个版本，所以注意修改不同的安装目录和端口号。可以直接使用sbin下的nginx命令启动服务。完成编译安装需要gcc支持，如果没有，使用如下命令安装。安装过程基本一致，下面是安装7.1.33的步骤。执行如下命令，检查已经安装的包和可安装的包。执行如下命令，检查已经安装的包和可安装的包。执行如下命令，检查已经安装的包和可安装的包。_tencentos-3.1

urllib.request.urlopen()基本使用_urllib.request.urlopen(url)-程序员宅基地

文章浏览阅读3.1w次，点赞21次，收藏75次。import urllib.requesturl = 'https://www.python.org'# 方式一response = urllib.request.urlopen(url)print(type(response)) # <class 'http.client.HTTPResponse'># 方式二request = urllib.request.Req..._urllib.request.urlopen(url)

如何用ChatGPT+GEE+ENVI+Python进行高光谱，多光谱成像遥感数据处理？-程序员宅基地

文章浏览阅读1.5k次，点赞12次，收藏15次。如何用ChatGPT+GEE+ENVI+Python进行高光谱，多光谱成像遥感数据处理？

RS485总线常识_rs485 差分走綫間距-程序员宅基地

文章浏览阅读1.2k次。RS485总线常识 2010-10-12 15:56:36| 分类：知识储备 | 标签：rs485 总线传输差分 |字号大中小订阅RS485总线RS485采用平衡发送和差分接收方式实现通信：发送端将串行口的TTL电平信号转换成差分信号A,B两路输出，经过线缆传输之后在接收端将差分信号还原成TTL电平信号。由于传输线通常使用双绞线，又是差分传输，所_rs485 差分走綫間距

随便推点

移植、制作uboot、Linux（一）_uboot制作-程序员宅基地

文章浏览阅读621次。u-boot、linux烧录_uboot制作

windows下安装git和gitbash安装教程_64-bit git for windows setup.-程序员宅基地

文章浏览阅读1.2w次，点赞10次，收藏44次。windos上git安装，git bash安装_64-bit git for windows setup.

环形链表（算法java）_java 实现环形链表-程序员宅基地

文章浏览阅读196次。环形链表（算法java）的两种解决方法_java 实现环形链表

docker部署Airflow(修改URL-path、更换postgres --＞myslq数据库、LDAP登录)_airflow docker-程序员宅基地

文章浏览阅读5.7k次。Airflow什么是 Airflow？Airflow 的架构Airflow 解决哪些问题一、docker-compose 安装airflow（postgres）1、创建启动文件airflow-docker-compose.yml.1.1、添加挂载卷，需要修改airflow-docker-compose.yml的位置2、创建本地配置文件airflow.cfg2.1、如果想修改WEB URL地址，需要修改airflow.cfg中以下两个地方3、之后up -d直接启动即可web访问地址：二、存储数据库更换post_airflow docker

计算机毕业设计springboot高校教务管理系统532k79【附源码+数据库+部署+LW】-程序员宅基地

文章浏览阅读28次。选题背景：随着社会的发展和教育的普及，高校教务管理系统在现代高等教育中扮演着至关重要的角色。传统的手工管理方式已经无法满足高校日益增长的规模和复杂的管理需求。因此，开发一套高效、智能的教务管理系统成为了当今高校管理的迫切需求。选题意义：高校教务管理系统的开发具有重要的意义和价值。首先，它可以提高高校教务管理的效率和准确性。通过自动化处理学生选课、排课、考试安排等繁琐的事务，大大减轻了教务人员的工作负担，提高了工作效率。同时，系统可以实时更新学生信息和课程信息，减少了数据错误和冗余，保证了管理的准确性

javaint接收float_Java Integer转换double，float，int，long，string-程序员宅基地

文章浏览阅读132次。首页>基础教程>常用类>常用 Integer类Java Integer转换double，float，int，long，stringjava中Integer类可以很方便的转换成double，float，int，long，string等类型，都有固定的方法进行转换。方法double doubleValue() //以 double 类型返回该 Integer 的值。flo..._java integet接收float类型的参数