11.16~11.19绘制图表，导入EXCEL中数据，进行拟合

2023/12/4 8:49:30

这个错误通常是由于传递给curve_fit函数的数据类型不正确引起的。根据你提供的代码和错误信息，有几个可能的原因：

  1. 数据类型错误：请确保ce_data、lg_data和product_data是NumPy数组或类似的可迭代对象，且其元素的数据类型为浮点数。你可以使用print(type(ce_data))等语句检查数据类型。

  2. 数据维度不匹配：curve_fit函数的第一个参数是要拟合的模型函数，其后的参数应该是要拟合的数据。根据你的代码，模型函数是reaction_model，而数据是(ce_data, lg_data)和product_data。请确保数据的维度匹配，即ce_data、lg_data和product_data具有相同的长度或形状。

  3. 数据元素过于复杂：错误信息中的"object too deep for desired array"提示可能是由于数据元素过于复杂而导致的。在使用curve_fit函数之前，确保数据不包含任何复杂对象或嵌套结构。如果数据中包含复杂对象，可以尝试将其转换为简单的数值或处理数据以满足函数的要求。

检查以上可能的原因，并根据需要调整代码和数据类型，以解决这个错误。

绘制图表

  • fig = plt.figure() 创建一个新的图形窗口。
  • ax = fig.add_subplot(111, projection='3d') 在图形窗口中添加一个三维子图（子图索引为111），并指定投影方式为3D。
  • ax.scatter(x1, x2, y, c='b', marker='o') 在三维子图中绘制散点图，其中x1和x2表示两个自变量的值，y表示因变量的值。参数c设置散点的颜色为蓝色（'b'），marker设置散点的形状为圆形（'o'）。

多项式回归

将excel导入进python

只有excel和py的主文件在同一目录下时才可以只写文件名，不然就必须写地址，即使在桌面上也不可以

data = pd.read_excel('C:/Users/26861/Desktop/Annex I- Pyrolysis Product Yields of Three Pyrolysis Combinations.xlsx')

还有就是地址可以用/，用\可能会出现转义字符 

import pandas as pd
import numpy as np
from sklearn.preprocessing import PolynomialFeatures
from sklearn.linear_model import LinearRegression

# 读取Excel文件
data = pd.read_excel('data.xlsx')

# 提取自变量和因变量数据
X = data['自变量列名'].values.reshape(-1, 1)  # 将自变量数据转换为二维数组
Y = data['因变量列名'].values.reshape(-1, 1)  # 将因变量数据转换为二维数组

# 创建多项式特征
poly_features = PolynomialFeatures(degree=2)
X_poly = poly_features.fit_transform(X)

# 拟合模型
model = LinearRegression()
model.fit(X_poly, Y)

# 进行预测
Y_pred = model.predict(X_poly)

# 打印拟合结果
print(Y_pred)

# 可以继续进行可视化等操作...
# 生成一系列连续的自变量值，用于绘制曲线
X_plot = np.linspace(0, 1.2, 100).reshape(-1, 1)
X_plot_poly = poly_features.transform(X_plot)

# 进行预测
Y_pred = model.predict(X_plot_poly)

# 绘制原始数据点
plt.scatter(X, Y, color='blue', label='原始数据')

# 绘制拟合曲线
plt.plot(X_plot, Y_pred, color='red', label='拟合曲线')

# 设置图例和标签
plt.legend()
plt.xlabel('自变量')
plt.ylabel('因变量')
plt.title('多项式回归拟合')

# 显示图形
plt.show()
import pandas as pd
import numpy as np
from sklearn.preprocessing import PolynomialFeatures
from sklearn.linear_model import LinearRegression
import matplotlib.pyplot as plt
# 读取Excel文件
data = pd.read_excel('C:/Users/26861/Desktop/Annex I- Pyrolysis Product Yields of Three Pyrolysis Combinations.xlsx')

# 提取自变量和因变量数据
X = data['DFA/CS'].values.reshape(-1, 1)  # 将自变量数据转换为二维数组
Y = data['Tar yield'].values.reshape(-1, 1)  # 将因变量数据转换为二维数组

# 创建多项式特征
poly_features = PolynomialFeatures(degree=2)
X_poly = poly_features.fit_transform(X)

# 拟合模型
model = LinearRegression()
model.fit(X_poly, Y)

# 进行预测
Y_pred = model.predict(X_poly)

# 打印拟合结果
print(Y_pred)

# 可以继续进行可视化等操作...
# 生成一系列连续的自变量值，用于绘制曲线
X_plot = np.linspace(0, 1.2, 100).reshape(-1, 1)
X_plot_poly = poly_features.transform(X_plot)

# 进行预测
Y_pred = model.predict(X_plot_poly)

# 绘制原始数据点
plt.scatter(X, Y, color='blue', label='point')

# 绘制拟合曲线
plt.plot(X_plot, Y_pred, color='red', label='curve')

# 设置图例和标签
plt.legend()
plt.xlabel('DFA/CS')
plt.ylabel('Tar yelid')
plt.title('curve')

# 显示图形
plt.show()

在从Excel提取自变量和因变量后，将它们转换为二维数组的原因是为了与某些库和函数的要求相匹配。

许多机器学习和数据分析库（如scikit-learn、NumPy等）通常期望输入数据是二维数组的形式，其中每个样本位于行中，每个特征位于列中。通过将自变量和因变量转换为二维数组，可以将数据结构化并满足库和函数的要求。

对于自变量，如果你只有一个自变量列，你可以使用NumPy的`reshape`函数将其转换为二维数组。假设自变量数组为`X`，你可以使用以下代码转换为二维数组：

```python
import numpy as np
X = np.array(X).reshape(-1, 1)
```

这将将`X`数组从一维数组转换为二维数组，其中样本位于行中，特征位于列中。`-1`参数表示自动推断行数，`1`表示列数为1。

对于因变量，如果你有多个因变量列，你可以将它们合并为一个二维数组。假设因变量数组为`Y1`、`Y2`、`Y3`，你可以使用以下代码将它们合并为一个二维数组：

```python
import numpy as np
Y = np.column_stack((Y1, Y2, Y3))
```

这将使用`column_stack`函数将`Y1`、`Y2`和`Y3`列堆叠在一起，形成一个二维数组，其中每个列对应一个因变量。

通过将自变量和因变量转换为二维数组，你可以更方便地使用各种库和函数进行数据分析、建模和拟合操作。

希望这能解答你的问题！如果还有其他问题，请随时提问。

import pandas as pd
import numpy as np
from sklearn.preprocessing import PolynomialFeatures
from sklearn.linear_model import LinearRegression
import matplotlib.pyplot as plt
# 读取Excel文件
data = pd.read_excel('C:/Users/26861/Desktop/Annex I- Pyrolysis Product Yields of Three Pyrolysis Combinations.xlsx')

# 提取自变量和因变量数据
X = data['DFA/CS'].values. Reshape(-1, 1)  # 将自变量数据转换为二维数组
Y1 = data['Tar yield'].values. Reshape(-1, 1)  # 将因变量数据转换为二维数组
Y2 = data['Water yield'].values. Reshape(-1, 1)
Y3 = data['Char yield'].values. Reshape(-1, 1)
Y4 = data['Syngas yield'].values.reshape(-1, 1)
# 创建多项式特征
poly_features = PolynomialFeatures(degree=2)
X_poly = poly_features.fit_transform(X)

# 拟合模型
model = LinearRegression()
model.fit(X_poly, Y1)
model.fit(X_poly, Y2)
model.fit(X_poly, Y3)
model.fit(X_poly, Y4)
# 进行预测
Y_pred = model.predict(X_poly)

# 打印拟合结果
print(Y_pred)

# 可以继续进行可视化等操作...
# 生成一系列连续的自变量值，用于绘制曲线
X_plot = np.linspace(0, 1.2, 100).reshape(-1, 1)
X_plot_poly = poly_features.transform(X_plot)

# 进行预测
Y_pred = model.predict(X_plot_poly)

# 绘制原始数据点
plt.scatter(X, Y, color='blue', label='point')

# 绘制拟合曲线
plt.plot(X_plot, Y_pred, color='red', label='curve')

# 设置图例和标签
plt.legend()
plt.xlabel('DFA/CS')
plt.ylabel('Tar yelid')
plt.title('curve')

# 显示图形
plt.show()
import pandas as pd
import numpy as np
import matplotlib.pyplot as plt
from scipy.optimize import curve_fit

# 从Excel文件中读取数据
data = pd.read_excel('C:/Users/26861/Desktop/Annex I- Pyrolysis Product Yields of Three Pyrolysis Combinations.xlsx')

# 提取自变量数据列
x = data['DFA/CS'].values

# 提取多个因变量数据列
y_columns = ['Tar yield', 'Water yield', 'Char yield','Syngas yield']  # 替换为实际的数据列名称

# 自定义拟合函数
def func(x, a, b):
    return a * x + b

# 创建图形窗口和子图
fig, axs = plt.subplots(len(y_columns), 1, figsize=(8, 6), sharex=True)

# 遍历每个因变量数据列
for i, y_column in enumerate(y_columns):
    # 提取因变量数据列
    y = data[y_column].values

    # 执行拟合
    params, _ = curve_fit(func, x, y)

    # 生成拟合曲线
    fit = func(x, *params)

    # 绘制原始数据点和拟合曲线
    axs[i].scatter(x, y, label='Data')
    axs[i].plot(x, fit, label='Fit')

    # 添加图例
    axs[i].legend()

# 设置整体图形的标题和横轴标签
fig.suptitle('Fitted Curves')
plt.xlabel('X')

# 调整子图之间的间距
plt.tight_layout()

# 显示图形
plt.show()

得到EXCEL中，

每个单元格在其所在行中的占比

import pandas as pd
import numpy as np
import matplotlib.pyplot as plt
from scipy.optimize import curve_fit

# 从Excel文件中读取数据
data = pd.read_excel('C:/Users/26861/Desktop/Annex I- Pyrolysis Product Yields of Three Pyrolysis Combinations.xlsx')

# 提取自变量数据列
x = data['DFA/CS'].values

# 提取多个因变量数据列
y_columns = ['Tar yield', 'Water yield', 'Char yield','Syngas yield']  # 替换为实际的数据列名称

# 多项式拟合阶数
degree = 2

# 创建图形窗口和子图
fig, axs = plt.subplots(len(y_columns), 1, figsize=(8, 6), sharex=True)

# 遍历每个因变量数据列
for i, y_column in enumerate(y_columns):
    # 提取因变量数据列
    y = data[y_column].values

    # 执行拟合
    coeffs = np.polyfit(x, y, degree)
    poly = np.poly1d(coeffs)
    fit = poly(x)


    # 绘制原始数据点和拟合曲线
    axs[i].scatter(x, y, label='Data')
    axs[i].plot(x, fit, label='Fit')

    # 添加图例
    axs[i].legend()
 # 输出拟合函数的具体信息
    print(f"Fitted function for {y_column}:")
    print(poly)

# 设置整体图形的标题和横轴标签
fig.suptitle('Fitted Curves')
plt.xlabel('X')

# 调整子图之间的间距
plt.tight_layout()

# 显示图形
plt.show()
import pandas as pd
import numpy as np
from scipy.optimize import curve_fit
import matplotlib.pyplot as plt

# 自定义指数函数
def exponential_func(x, a, b, c):
    return a * np.exp(b * x) + c

# 从Excel文件中读取数据
df = pd.read_excel('C:/Users/26861/Desktop/Annex I- Pyrolysis Product Yields of Three Pyrolysis Combinations.xlsx')
x_data = df['DFA/CS'].values
y_data = df[['Tar yield', 'Water yield', 'Char yield','Syngas yield']].values.T


print(x_data,y_data)
# 创建子图
fig, axs = plt.subplots(2, 2, figsize=(10, 8))

# 针对每个因变量进行拟合和绘图
for i, ax in enumerate(axs.flat):
    y = y_data[i]  # 当前因变量的数据

    # 进行拟合
    popt, pcov = curve_fit(exponential_func, x_data, y, maxfev=10000)  # 增加maxfev的值

    # 绘制拟合曲线
    x_fit = np.linspace(min(x_data), max(x_data), 100)
    y_fit = exponential_func(x_fit, *popt)
    ax.plot(x_fit, y_fit, label='Fit Curve')

    # 绘制原始数据点
    ax.scatter(x_data, y, label='Data')

    # 显示方程信息
    equation_info = f'y = {popt[0]:.2f} * exp({popt[1]:.2f} * x) + {popt[2]:.2f}'
    ax. Text(2, max(y) / 2, equation_info)

    ax.set_xlabel('x')
    ax.set_ylabel('y')
    ax.set_title(f'Fit Curve for Variable {i+1}')
    ax.legend()

# 调整子图布局
plt.tight_layout()

# 显示图表
plt.show()
return a * np.exp(-(x - b)**2 / (2 * c**2)) + d
# 非线性最小二乘拟合
popt, pcov = curve_fit(exponential_func, x_data, y_data, maxfev=100000,p0=(1, -1, 1))

拟合时初始值很重要

curve_fit就是最小二乘拟合

对于单调减的初始值设为-1，开始时增加的，设为1

import numpy as np
import matplotlib.pyplot as plt
from scipy.optimize import curve_fit

x_data = np.array([0., 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.8, 1.])
y_data = np.array([19.46, 17.25, 15.43, 14.14, 13.89, 13.21, 12.84, 12.57, 12.13])

def exponential_func(x, a, b, c, d):
    return a * np.exp(b * x + d) + c

popt, pcov = curve_fit(exponential_func, x_data, y_data, p0=(1, -1, 1, 0))

a_fit, b_fit, c_fit, d_fit = popt
expression = f"{a_fit:.4f} * exp({b_fit:.4f} * x + {d_fit:.4f}) + {c_fit:.4f}"
print("expression：", expression)

x_fit = np.linspace(min(x_data), max(x_data), 100)
y_fit = exponential_func(x_fit, *popt)

plt.plot(x_fit, y_fit, label='Predicted value')
plt.scatter(x_data, y_data, label='Actual')
plt.xlabel('Mixing ratio of DFA/CS')
plt.ylabel('Tar yield')
plt.title('Kinetic model')
plt.legend()
plt.show()

y_pred = exponential_func(x_data, *popt)
comparison_table = np.column_stack((x_data, y_data, y_pred))
print("\nA table comparing predicted values to actual values：")
print("r     Actual      Predicted value")
print("-------------------------")
for row in comparison_table:
    print(f"{row[0]:.1f}    {row[1]:.2f}    {row[2]:.2f}")

多元线性回归拟合以及摘要 

import pandas as pd
import numpy as np
import statsmodels.api as sm
# 指定 Excel 文件路径
excel_file = 'C:/Users/26861/Desktop/Annex II-Pyrolysis Gas Yields of Three Pyrolysis Combinations.xlsx'
# 提取 X 列的数据作为自变量
df = pd.read_excel(excel_file)

X = df[['H2CE', 'H2LG']]
y = df['H2CS']

# 添加常数列
X = sm.add_constant(X)

# 拟合多元线性回归模型
model = sm.OLS(y, X)
results = model.fit()

# 打印回归结果摘要
print(results.summary())
equation = 'y = {:.4f} + {:.4f}*X1 + {:.4f}*X2 '.format(results.params['const'], results.params['H2CE'], results.params['H2LG'])
print(equation)


http://www.jnnr.cn/a/1081959.html

相关文章

【grafana ｜ clickhouse】实现展示多折线图

【grafana ｜ clickhouse】实现展示多折线图

说明： 采用的是 Visualizations 的 Time series，使用的 clickhouse 数据源 在工作中遇到了一个需求，写好了代码，需要在grafana上展示在一个项目中所有人的，随时间的代码提交量变化图 目前遇到的问题：展示…
阅读更多...
使用树莓派学习Linux系统编程的 --- 库编程（面试重点）

使用树莓派学习Linux系统编程的 --- 库编程（面试重点）

在之前的Linux系统编程中，学习了文件的打开；关闭；读写；进程；线程等概念.... 本节补充“Linux库概念 & 相关编程”，这是一个面试的重点！ 分文件编程 在之前的学习中，面对较大的…
阅读更多...
音视频同步笔记 - 以音频时间为基

音视频同步笔记 - 以音频时间为基

音视频同步 - 以音频时间为基 上图介绍： 该图是以音频的时间为基，对视频播放时间的延迟控制方案，只调整视频的播放延时。delayTime是视频播放的延迟时间，初始值是1 / FPS * 1000 (ms)，如果FPS为25帧率，初始…
阅读更多...
Linux进程通信——IPC、管道、FIFO的引入

Linux进程通信——IPC、管道、FIFO的引入

进程间的通信——IPC 进程间通信 (IPC，InterProcess Communication) 是指在不同进程之间传播或交换信息。 IPC的方式通常有管道 (包括无名管道和命名管道) 、消息队列、信号量、共享存储、Socket、Streams等。其中 Socket和Streams支持不同主机上的两个进程IPC。 …
阅读更多...
原理Redis-ZipList

原理Redis-ZipList

ZipList 1) ZipList的组成2) ZipList的连锁更新问题3) 总结 1) ZipList的组成 ZipList 是一种特殊的“双端链表” ，由一系列特殊编码的连续内存块组成。可以在任意一端进行压入/弹出操作, 并且该操作的时间复杂度为 O(1)。 ZipListEntry: ZipList 中的Entry并不像…
阅读更多...
python -opencv 边缘检测

python -opencv 边缘检测

python -opencv 边缘检测 边缘检测步骤: 第一步：读取图像为灰度图 第二步：进行二值化处理 第三步：使用cv2.findContours对二值化图像提取轮廓 第三步：将轮廓绘制到图中 代码如下： from ctypes.wintypes import SIZ…
阅读更多...
几个强力的nodejs库

几个强力的nodejs库

几个强力的nodejs库 nodejs被视为许多Web开发人员的理想运行时环境。 nodejs的设计是为了在运行时中使用JavaScript编写的代码，它是世界上最流行的编程语言之一，并允许广泛的开发者社区构建服务器端应用程序。 nodejs提供了通过JavaScript库重用代码的…
阅读更多...
C#开发的OpenRA游戏之属性QuantizeFacingsFromSequence（7）

C#开发的OpenRA游戏之属性QuantizeFacingsFromSequence（7）

C#开发的OpenRA游戏之属性QuantizeFacingsFromSequence（7） 前面分析了身体的方向，在这里继续QuantizeFacingsFromSequence属性，这个属性就是通过序列定义文件里获取身体的方向。 根据前面分析可知，同样有一个信息类QuantizeFacingsFromSequenceInfo： [Desc("Deriv…
阅读更多...
【2023云栖】陈守元：阿里云开源大数据产品年度发布

【2023云栖】陈守元：阿里云开源大数据产品年度发布

本文根据 2023 云栖大会演讲实录整理而成，演讲信息如下： 演讲人：陈守元 | 阿里云计算平台事业部开源大数据产品总监 演讲主题：阿里云开源大数据产品年度发布 随着云计算的不断发展，未来数据处理和应用的趋势将围绕C…
阅读更多...
基于安卓android微信小程序的好物分享系统

基于安卓android微信小程序的好物分享系统

运行环境 开发语言：Java 框架：ssm JDK版本：JDK1.8 服务器：tomcat7 数据库：mysql 5.7（一定要5.7版本） 数据库工具：Navicat11 开发软件：eclipse/myeclipse/idea Maven包&a…
阅读更多...
Springboot+vue的机动车号牌管理系统（有报告）。Javaee项目，springboot vue前后端分离项目

Springboot+vue的机动车号牌管理系统（有报告）。Javaee项目，springboot vue前后端分离项目

演示视频: Springbootvue的机动车号牌管理系统（有报告）。Javaee项目，springboot vue前后端分离项目 项目介绍： 本文设计了一个基于Springbootvue的前后端分离的机动车号牌管理系统，采用M（model&#xff09…
阅读更多...
LeetCode977.有序数组的平方（双指针法、暴力法、列表推导式）

LeetCode977.有序数组的平方（双指针法、暴力法、列表推导式）

LeetCode977.有序数组的平方 1.问题描述2.解题思路3.代码4.知识点 1.问题描述 给你一个按 非递减顺序 排序的整数数组 nums，返回 每个数字的平方 组成的新数组，要求也按 非递减顺序 排序。 示例 1： 输入：nums [-4,-1,0,3,10] …
阅读更多...
Node.js之TCP（net）

Node.js之TCP（net）

Hi I’m Shendi Node.js之TCP（net） 最近使用Nodejs编写程序，需要用到自己编写的分布式工具，于是需要将Java版的用NodeJs重新写一遍，需要使用到TCP通信，于是在这里记录下Node.js TCP 的使用方法 依赖 需要使…
阅读更多...
在UOS系统中编译CEF源码

在UOS系统中编译CEF源码

一、下载cef代码 git clone gitbitbucket.org:chromiumembedded/cef.git 二、执行自动下载代码 由于chromium的代码很大，至少需要准备大概80G的硬盘！！！整个代码量太大还是多准备一些空间吧（强烈建议使用固态硬盘保存否…
阅读更多...
指针变量和地址

指针变量和地址

A.指针变量和地址 理解了内存和地址的关系&#xff0c;我们再回到C语⾔&#xff0c;在C语⾔中创建变量其实就是向内存申请空间&#xff0c;比如&#xff1a; #include <stdio.h> int main() {int a 10;return 0; } ⽐如&#xff0c;上述的代码就是创建了整型变量a&…
阅读更多...
你好，我叫Python，欢迎你认识派森。（来自关于Python语言的全方位自我介绍。

你好，我叫Python，欢迎你认识派森。（来自关于Python语言的全方位自我介绍。

文章目录 自我简介一、Python的发展历程二、Python的特色1.语言特色2.语法特色 三、Python2与Python3的比较1.print 函数2.Unicode3.除法运算4.异常5.八进制字面量表示6.不等运算符7.python 3.0严格使用tab键进行缩进 四、Python适用开发场景及成果1.应用领域2.Python开发出的应…
阅读更多...
特效！视频里的特效在哪制作——Adobe After Effects

特效！视频里的特效在哪制作——Adobe After Effects

今天&#xff0c;我们来谈谈一款在Adobe系列中推出的一款图形视频处理软件&#xff0c;适用于从事设计和视频特技的机构&#xff0c;包括电视台、动画制作公司、个人后期制作工作室以及多媒体工作室的属于层类型后期软件——Adobe After Effects。 Adobe After Effects&#xf…
阅读更多...
OTP语音芯片 NV080D在智能空气检测仪的应用

OTP语音芯片 NV080D在智能空气检测仪的应用

随着人们对健康和环保的关注度不断提高&#xff0c;人们对看不见的家居环境也越来越重视。智能空气检测仪的市场需求也在不断增长中&#xff0c;呈现稳中向好的趋势。智能空气检测仪能够检测室内空气中的PM2.5、甲醛、TVOC等有害物质&#xff0c;同时还可以检测温湿度、空气质量…
阅读更多...
通过easyexcel实现数据导入功能

通过easyexcel实现数据导入功能

上一篇文章通过easyexcel导出数据到excel表格已经实现了简单的数据导出功能&#xff0c;这篇文章也介绍一下怎么通过easyexcel从excel表格中导入数据。 目录 一、前端代码 index.html index.js 二、后端代码 controller service SongServiceImpl 三、功能预览 四、后端…
阅读更多...
企业实现员工聊天和转账行为的实时监管

企业实现员工聊天和转账行为的实时监管

如何解决企业营销团队的管理问题&#xff1f; 在当今竞争激烈的市场环境中&#xff0c;企业营销团队的管理显得尤为重要。营销团队是企业发展的重要支柱&#xff0c;然而&#xff0c;一些常见的问题如员工飞单、私单、辱骂删除客户、离职带走公司客户以及工作不认真、工作量无…
阅读更多...
最新文章

Copyright @ 2022~2023