-
0 votesanswersviews
R中的Eta / Eta平方例程
除了在应用来自GLM族的一些技术之前使用的线性度的图形估计(凝视散点图方法)之外,有几种方法可以算术地进行这种估计(即没有图形) . 现在,如果两个变量之间的关系是线性的,那么我就得到了相关比:算术上,它是一个r(确定系数:r2) . 因此,您可以比较eta和r的值,并对关系类型(线性或非线性)进行评估 . 它提供了由自变量(线性或非线性)解释的因变量的方差百分比信息 . 因此,您可以在不满足线性... -
1 votesanswersviews
cplex中的分段非线性代价函数
是否可以在Cplex中建模非线性分段成本函数? 比如我在这里放的数字: non linear piece wise Cost function (black line) 我知道一种方法是将二次部分线性化为线性部分,但是,我想使用二次部分 . 您可以看到条件在决策变量本身上,成本函数可以表示如下: 如果x≲x0则成本是二次部分; 其他成本是线性部分 . 提前致谢 :) -
12 votesanswersviews
如何使用牛顿方法(代码非线性代数)找到最小的非线性,多元函数
我正在尝试进行一些参数估计,并希望选择最小化预测方程 over about 30 variables 中的平方误差的参数估计 . 如果方程是线性的,我只计算30个偏导数,将它们全部设为零,并使用线性方程求解器 . 但不幸的是 the equation is nonlinear 以及它的衍生品 . 如果方程式超过单个变量,我只会使用Newton's method(也称为Newton-Raphson)... -
0 votesanswersviews
用scipy的fsolve求解方程
我正试图解决方程 f(x) = x-sin(x) -n*t -m0 在这个等式中, n 和 m0 是我的类中定义的属性 . 此外, t 是等式中的常数整数,但每次都必须改变 . 我已经解决了这个等式,所以我得到了一个“新的等式” . 我已经导入了scipy.optimize def f(x, self): return (x - math.sin(x) -self.M0 - self.n*... -
0 votesanswersviews
在MATLAB中求解大量不明确的非线性方程
我将尝试描述我的问题,但最后还是有一些具体的问题 . 我在解决Matlab中的涡轮机械问题时遇到了一些麻烦 . 目标是用已知的物理测量和操作条件计算涡轮机的效率 . 要做到这一点,我使用科学论文中给出的方程式:http://www.sciencedirect.com/science/article/pii/S0196890409003124 在我的情况下的问题是,要知道效率,我需要出口气体参数(温... -
1 votesanswersviews
初始用户提供的目标函数评估失败 . FSOLVE无法继续
我试图在MATLAB 2013中使用fsolve解决一大组非线性方程 . 当我运行这段代码时: Unkn0=zeros(3249,1); % Guess [Unkn,feval,exitflag,output]=fsolve(Eqnh,Unkn0); 这里Eqnh是3249变量和3249方程的函数句柄 . 我收到以下错误: Error using Not enough input argumen... -
6 votesanswersviews
从python中的非线性方程组中寻找复杂的根
我一直在测试一种已经在文献中发表的算法,该算法涉及在Matlab和Python中求解一组“m”非线性方程 . 这组非线性方程涉及包含复数的输入变量,因此得到的解也应该是复杂的 . 截至目前,我已经能够通过使用以下代码行在Matlab中获得相当不错的结果: lambdas0 = ones(1,m)*1e-5; options = optimset('Algorithm','levenberg-mar... -
0 votesanswersviews
在matlab中求解大约束线性和非线性方程组的最佳方法是什么?
MATLAB中约束线性和非线性方程组的求解 我正在使用直接刚度方法在MATLAB中解决FEM问题 . 现在将问题表述为非线性方程组: KU=F or CF=U with C = K^{-1}. 问题是在3-Dd中制定的,有200多个节点 . 该模型用作变形估计器,力总是输入和变形输出 . 所以有很多线性方程 . 现在需要扩展这个问题,需要实现2种类型的约束 . 一些节点必须被约束,以便它们只能围... -
0 votesanswersviews
Scipy - 解包含CDF的非线性方程组
我有以下方程组(简化版) . y = x*scipy.special.gdtr(a, b, f(x,y)) x = y*scipy.speical.gdtr(c, d, f(x,y)) 这里 scipy.special.gdtr 指的是伽玛分布的CDF, a , b 是伽马CDF所采用的相应的两个参数 . 我的问题如下: 在这种情况下,是否有首选求解器? fsolve 或 brentq 或其... -
1 votesanswersviews
在MATLAB中求解两个困难的非线性方程组
我有两个非线性方程,有三个变量,这些变量是在论文中得出的 . 在本文中,第三个变量从一定值变化到另一个值,另外两个变量相对于它绘制 . 这两个方程是(1) abs(C(z,theta1,theta2))=0.1 和(2) diff(abs(C(z,theta1,theta2)),w)=0 其中 C(z,theta1,theta2)=A/B 与 A=j*(2*(z-1/z)*tan(theta1)+... -
0 votesanswersviews
在matlab中求解非线性方程组
我有一个我无法回答的作业问题 . 这是问题提示: De fi ne Eq . 8.3和Eq . 8.4在一个函数中 . 此函数应采用关节角度矢量(和)作为输入,并应返回包含在这些角度评估的两个函数(和)的列向量 . 该函数应包含和的链接长度和值以进行计算 . 在解决方案中包含您的代码 . 公式8.3和8.4是: (8.3) (8.4) 哪里 这是我写的函数: function F... -
1 votesanswersviews
Matlab中3变量非线性系统的Newton-Raphson方法
我试图在MATLAB中使用Newton-Raphson方法求解3个变量的3个非线性系统 . 以下是3个非线性方程: c * (alpha*I + k_f + k_d + k_n * s + k_p*(1-q))-I *alpha = 0 s * (lambda_b * c* P_C + lambda_r *(1-q))- lambda_b* c * P_C = 0 q * ( ga... -
3 votesanswersviews
求解python中的colebrook(非线性)方程
我想在MATLAB中的python what this guy did中做 . 我已经安装了anaconda,所以我有numpy和sympy库 . 到目前为止,我已尝试使用numpy nsolve,但这不起作用 . 我应该说我是python的新手,而且我知道如何在MATLAB中做到这一点:P . The equation: -2*log(( 2.51/(331428*sqrt(x)) ) + ( ... -
2 votesanswersviews
Matlab的fsolve收敛*但*似乎给出了错误的解决方案
我试图用fsolve求解一个非线性方程组;我们说吧 F(x; lambda)= 0,其中lambda是参数的向量,x是我想要求解的向量 . 我正在使用Matlab的fsolve . 我有2个参数lambda值,我想解决系统问题 . 对于lambda的一个值,我得到一个解决方案,这似乎没问题 . 对于lambda的另一个值,我再次得到一个解决方案(matlab以1的标志退出 . 但是我知道这不是一个... -
0 votesanswersviews
在MATLAB中求解非线性有限元
我试图用MATLAB中的节点热源解决四面体有限元上的热扩散问题,这取决于解向量 . 非线性方程系统如下所示: BU'AU = q(T) B是热容矩阵,A是电导率矩阵,q是源项,U是温度 . 我使用Adams-Bashforth / Trapezoid Rule预测器 - 校正器方案,Picard迭代后跟时间步长控制 . 源项的温度精确地在最后一个时间步的温度和预测器的温度之间进行评估 . 这是预测... -
0 votesanswersviews
fsolve / fzero:找不到解决方案,看似常规
我正在尝试使用fsolve或fzero执行以下算法: K5=8.37e-2 P=1 Choose an A S2=(4*K5/A)^(2/3) S6=3*S2 S8=4*S2 SO2 = (5*P)/149 - (101*S2)/149 - (293*S6)/149 - (389*S8)/149 H2O = (40*P)/149 + (556*S2)/4... -
1 votesanswersviews
使用MATLAB 's ' decic'函数实现一致初始条件的收敛失败
我有一组高度非线性的隐式ODE,我需要使用ode15i等隐式求解器来积分方程 . 为了使用ode15i,您必须获得衍生物和因变量的一致初始条件 . MATLAB有一个内置函数'decic',可以计算一致的IC,但是,当我执行该函数时,我得到以下错误: “DECIC的融合失败 . ” 现在,一般来说,如果你的方程式不一致,'decic'会给你一个警告,要求你释放一些Y0的变量 . 但在我的情况下,... -
0 votesanswersviews
在python中遇到ODEINT问题
我是Python的新手,并试图用它来解决二阶非线性微分方程,特别是电解质中的Poisson-Boltzmann方程 . phi''(r) + (2/r)*phi'(r) = (k^2)*sinh(phi(r)) 基本上它描述了静电势(phi)远离电解质中的带电表面的衰减,其中衰减速率由参数k控制 . phi(r) - r的潜力 dphi(r) - r的潜在导数 r - 距离曲面... -
0 votesanswersviews
在求解方程时它意味着什么,只有括号以最大值回归
在求解方程时它意味着什么,只有括号以最大值回归? 我试图解决phi,它保持双方的phi . 我试着解决A并用括号回来任何想法为什么? 这是等式 eq2:acos(A*cos((2*%pi*f*t+phi)/A) -2*%pi*f*t)=phi; wxWidgets:3.0.0 Unicode支持:是Maxima版本:5.32.1 即使我在其中添加值仍然拒绝将phi放在一边 -
1 votesanswersviews
使用scipy.optimize非线性求解器时出错
我试图解决一组M个同时具有M个变量的eqns . 我输入一个M X 2矩阵作为我的函数的初始猜测,它返回一个M X 2矩阵,如果我的猜测正确,每个条目将等于零 . 因此,对于 k=1,2,...N ,我的函数可以表示为 f_k(u1,u2,...uN) = 0 . 下面是我的函数的代码,(为了简单起见,我遗漏了与这段代码一起使用的模块,例如p . 或phi . 我更想知道是否有其他人之前有这个错... -
0 votesanswersviews
非线性方程MATLAB
我不知道如何用matlab解决这个等式: F(1)= 0.098253 ×(1)-atan((TAN(X(1)) - 棕褐色(X(2)))/ 2) * 0.531268-0.433015 X(2)-0.27994 f(2)= 0.9951 * x(1)-atan((tan(x(1)) - tan(x(2)))/ 2)* 0.12909 0.866022 * x(2)-0.350005; 我尝试使... -
4 votesanswersviews
比较python和matlab中的fsolve结果
我对几天前写的帖子有一个跟进问题,谢谢你以前的反馈: Finding complex roots from set of non-linear equations in python 我现在已经在python中设置了设置的非线性方程,因此fsolve将独立处理实部和虚部 . 但是,仍然存在python“fsolve”收敛到正确解决方案的问题 . 我有完全相同的输入在Matlab中使用,经过双重检查... -
1 votesanswersviews
非线性拟合
我有一些实验数据,我想使用最小二乘法(Levenberg-Marquardt)来拟合我的参数 . 我正在使用两个非线性方程式,我正在使用一些计算程序(Origin和Matlab) . 首先是: y=A+B*(((2*pi*x)^2+Alfa4^2)*((2*pi*x)^2+Alfa5^2))/(((2*pi*x)^2+Alfa1^2)*((2*pi*x)^2+Alfa2^2)*((2*pi*x)^... -
1 votesanswersviews
matlab:耦合非线性方程组的大系统
我正试图解决一个非常大的耦合非线性方程组 . 在此thread以及Matalb的相关帮助(第一个示例)后,我尝试编写以下代码: %% FSOLVE TEST #2 clc; clear; close all %% global a0 a1 a2 a3 a4 h0 TM JA JB a0 = 2.0377638272727268; a1 = -7.105521894545453; a2 = 9.23... -
1 votesanswersviews
Scipy - 非线性方程组的所有解
我有一个非线性方程组,其中任何n都可以选择,因此向量x =(x1,...,xn)的长度可以不同 . 例如,系统可以是这样的: f1(x1,...,xn) = sum( xi + xi^2 ) = 0, i={1,n} f2(x1,...,xn) = sum( e^xi + xi + sin(xi*pi) ) = 0, i={1,n} 根据this example,我使用scipy库的fs... -
1 votesanswersviews
Scipy的解决方案?
我试图在数值上解决一个非线性方程组: def func(p): x, f = p return (math.exp(-x/O)-f, L - L*((1 - math.exp(-x/O))**W) - x*math.exp(-x/O)) 我正在以下列方式使用scipy.fsolve: x, f = fsolve(func, (10, 0.2)) 我确定我...