python拟合曲线求方程,曲线拟合方程python

我试图用一个包含6个元件的rlc电路的导纳方程来拟合数据。我按照给定的he[fit]1re的例子，插入我的方程。该方程是用Mathcad简化的六元件电路导纳的实部。在附图中，x轴是ω(w=2*pi*f)，y是以milli Siemens为单位的导纳。在

该程序运行，但它没有做拟合，尽管有一个良好的试验功能。我感谢任何帮助，为什么适合是一条直线。我还附上了一个高斯拟合的例子。在

这就是当我试图拟合方程时得到的结果。数据为左侧峰值较小的数据，试验函数为虚线。拟合是一条直线

from numpy import sqrt, pi, exp, linspace, loadtxt

from lmfit import Model

import matplotlib.pyplot as plt

data = loadtxt("C:/Users/susu/circuit_eq_real5.dat")

x = data[:, 0]

y = data[:, 1]

def circuit(x,C0,Cm,Lm,Rm,R0,Rs):

return ((C0**2*Cm**2*Lm**2*R0*x**4)+(Rs*C0**2*Cm**2*Lm**2*x**4)+(C0**2*Cm**2*R0**2*Rm*x**2)+(Rs*C0**2*Cm**2*R0**2*x**2)+(C0**2*Cm**2*R0*Rm**2*x**2)+(2*Rs*C0**2*Cm**2*R0*Rm*x**2)+(Rs*C0**2*Cm**2*Rm**2*x**2)-(2*C0**2*Cm*Lm*R0*x**2)-(2*Rs*C0**2*Cm*Lm*x**2)+(C0**2*R0)+(Rs*C0**2)-(2*Rs*C0*Cm**2*Lm*x**2)+(2*Rs*C0*Cm)+(Cm**2*Rm)+(Rs*Cm**2))/((C0**2*Cm**2*Lm**2*x**4)+(C0**2*Cm**2*R0**2*x**2)+(2*C0**2*Cm**2*R0*Rm*x**2)+(C0**2*Cm**2*Rm**2*x**2)-(2*C0**2*Cm*Lm*x**2)+(C0**2)-(2*C0*Cm**2*Lm*x**2)+(2*C0*Cm)+(Cm**2))

gmodel = Model(circuit)

result = gmodel.fit(y, x=x, C0=1.0408*10**(-12), Cm=5.953*10**(-14),

Lm=1.475*10**(-7), Rm=1.571, R0=2.44088, Rs=0.42)

print(result.fit_report())

plt.plot(x, y, 'bo')

plt.plot(x, result.init_fit, 'k--')

plt.plot(x, result.best_fit, 'r-')

plt.show()

下面是Fit报告

^{pr2}$

非常感谢M Newville和Mikuszefski以及其他人提供的见解和反馈。我同意我所说的也许是一个混乱的程序。从python代码中可以明显看出，我并不精通python或编程。在

Mikuszefsky，感谢您发布rlc示例代码。你的方法既简洁又有趣。我不知道Python是直接进行复杂拟合的，我会试试你的方法，看看能不能做到。我要同时拟合Y(准入)的实部和虚部。我肯定会被困在某个地方，并将我的进步张贴在这里。

最好的，

苏苏

python拟合曲线求方程,曲线拟合方程python相关推荐

python拟合曲线误差分析_python曲线拟合不能给出合理的拟合结果
您错误地调用了curve_fit,下面是用法curve_fit(f, xdata, ydata, p0=None, sigma=None, absolute_sigma=False, check_fi ...
python惰性求值_让Python中类的属性具有惰性求值的能力
起步我们希望将一个只读的属性定义为 property 属性方法,只有在访问它时才进行计算,但是,又希望把计算出的值缓存起来,不要每次访问它时都重新计算. 解决方案定义一个惰性属性最有效的方法就是利 ...
python迭代法求极值_用Python实现最速下降法求极值的方法
对于一个多元函数 ,用最速下降法(又称梯度下降法)求其极小值的迭代格式为其中为负梯度方向,即最速下降方向,αkαk为搜索步长. 一般情况下,最优步长αkαk的确定要用到线性搜索技术,比如精确线性搜 ...
python编程求导数_用python怎么计算导数最简单？
谢邀,请恕我微积分学得不扎实,我记得常数的一阶导数均为0. 如果列表中传入的为含变量x的式子,代码可能如下. from sympy import Symbol, diff x = Symbol('x' ...
python 多项式求系数_在Python中用于计算“多项式系数”的numpy / scipy函数
是否有任何 python函数(可能来自numpy或scipy)计算扩展中x ** r的系数(1 xx ** 2 x ** 3 - x **(k-1))** n ,其中k> = 1,n> = ...
python递归求5!_用Python解数独[6]：递归获得最终答案
import copy import time # 求每一行单元格行值域 def valueRange(row): temp = copy.deepcopy(row) row_value_range ...
python编程求导数_在python中如何计算导数？
这个问题的答案很简单.当然,在另一个答案中给出的subs选项可用于计算某个数的导数,但如果要绘制导数,则不起作用.有一种方法可以解决这个问题:lambdify,如下所述. 使用lambdify将所有s ...
maya python 创建求_如何使用python在Maya中创建列表
预期输出我想在Maya 2014中创建一个用户界面,其中包含图像中给定的图层和相机.我不知道该使用什么小部件.我尝试创建,代码如下所示.在import maya.cmds as cmds windo ...
python如何求列表的平均值_python如何求列表平均值？
推荐教程:<python视频教程> python如何求列表平均值? python函数求列表平均值的方法: 用法:mean(matrix,axis=0)其中matrix为一个矩阵,axis为 ...
python中的sympy模块求导和解方程
利用python中的sympy模块求导和解方程对sympy这个模块我不是很了解,我学习python主要是想往机器学习和深度学习方向发展.这次了解到有这个模块,是因为有同学在写论文中要用到求导,并且计 ...

python拟合曲线求方程,曲线拟合方程python

python拟合曲线求方程,曲线拟合方程python相关推荐

最新文章

热门文章