zuoyefengqi 2022-10-12 19:04 采纳率: 85.7%
浏览 130
已结题

python拟合曲线为什么画的图中一条曲线一条直线

python拟合曲线为什么画的图中一条曲线一条直线?给定了函数形式让程序拟合的。

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


#去掉了最后一个水平


# x就是等效应变
def dfin(x, a, b, c):
    return (a*z**b-30.0)*(1.0-np.exp(-np.log(2.0*(x/ee)**c)))+30.0

# D初始晶粒尺寸
D = 30.0
# z不同温度需要重新输入
z = 1.92*10**20
ee = 0.0016*D**0.2709*z**0.1147
# 给出x、y对应坐标,即给出等效应变和对应的真实最终晶粒尺寸
x = [0.74,0.69,0.56,0.41,0.3,0.21,0.15,0.12,
    1.38,1.33,1.18,1.01,0.82,0.64,0.48,0.35]
x = np.array(x)
num = [26,29,28,30,31,40,39,41,
    21,22,23,24,23,27,28,30]
y = np.array(num)

# 非线性最小二乘法拟合
parameter_bounds = ([15*10**5,-0.25,0],[np.inf,0,3])
popt, pcov = curve_fit(dfin, x, y, bounds=parameter_bounds)
# 获取popt里面拟合系数
print(popt)
a = popt[0]
b = popt[1]
c = popt[2]
yvals = dfin(x, a, b, c)
print('系数a:', a)
print('系数b:', b)
print('系数c:', c)
print('系数yvals:', yvals)

# 绘图
plot1 = plt.plot(x, y, 's')
plot2 = plt.plot(x, yvals, 'r')
plt.xlabel('ε')
plt.ylabel('d')
plt.title('1120,v=14,zhui')
plt.show()


图片如下:

img

  • 写回答

2条回答 默认 最新

  • 天元浪子 Python领域优质创作者 2022-10-12 22:39
    关注

    x数组没有排序,最小值索引7,最大值索引8,直线就是这两点的连线。建议将 yvals = dfin(x, a, b, c) 改为 yvals = dfin(x[np.argsort(x)], a, b, c)

    本回答被题主选为最佳回答 , 对您是否有帮助呢?
    评论
查看更多回答(1条)

报告相同问题?

问题事件

  • 系统已结题 10月21日
  • 已采纳回答 10月13日
  • 创建了问题 10月12日

悬赏问题

  • ¥15 易优eyoucms关于二级栏目调用的问题
  • ¥15 寻找公式识别开发,自动识别整页文档、图像公式的软件
  • ¥15 为什么eclipse不能再下载了?
  • ¥15 编辑cmake lists 明明写了project项目名,但是还是报错怎么回事
  • ¥15 关于#计算机视觉#的问题:求一份高质量桥梁多病害数据集
  • ¥15 特定网页无法访问,已排除网页问题
  • ¥50 如何将脑的图像投影到颅骨上
  • ¥15 提问一个关于vscode相关的环境配置问题,就是输入中文但是显示不出来,代码在idea可以显示中文,但在vscode不行,不知道怎么配置环境
  • ¥15 netcore使用PuppeteerSharp截图
  • ¥20 这张图页头,页脚具体代码该怎么写?