R函数输入值的求解
是指在R语言中，通过给定函数的输入值，计算出函数的输出值的过程。R是一种用于统计分析和数据可视化的编程语言，具有丰富的数学和统计函数库，可以进行各种数值计算和数据处理操作。
在R中，可以使用函数调用的方式来求解函数的输入值。函数调用的一般形式是函数名后跟一对圆括号，括号中包含函数的输入值。例如，对于一个简单的函数f(x)，可以通过f(2)来求解在输入值为2时的函数输出值。
R函数的输入值可以是单个数值，也可以是向量、矩阵或数据框等数据结构。对于向量或矩阵等多个输入值的情况，R会自动对每个输入值进行相应的计算，并返回相应的输出值。
R函数的求解可以通过内置函数来实现，也可以通过自定义函数来实现。内置函数是R语言提供的已经实现好的函数，可以直接调用使用。自定义函数是根据用户的需求编写的函数，可以根据具体的问题来定义函数的输入和输出。
在R中，有许多用于求解函数输入值的函数和技术。例如，可以使用optimize()函数来求解函数的最优输入值，使用uniroot()函数来求解函数的根，使用integrate()函数来求解函数的积分值等。
对于R函数输入值的求解，可以应用于各种领域和场景。例如，在统计分析中，可以通过求解函数的输入值来进行参数估计和模型拟合；在数据可视化中，可以通过求解函数的输入值来生成图形和图表；在机器学习和数据挖掘中，可以通过求解函数的输入值来进行模型训练和预测等。
腾讯云提供了一系列与云计算相关的产品和服务，可以帮助用户进行云计算的开发和运维工作。具体推荐的产品和产品介绍链接地址可以根据具体的需求和场景来选择。
相关·内容
1回答
我在上遇到了一个奇怪的问题。下面是一个最低限度的工作示例：
import sys
import time
import numpy as np
import matplotlib.pyplot as plt
import matplotlib
from scipy.integrate import ode, complex_ode
def fun(t):
return np.exp( -t**2 / 2. )
def ode_fun(t, y):
a, b = y
f = fun(t)
c = np.conjugate(b)
dt_a = -2j*f
浏览 2提问于2021-03-13得票数 0
回答已采纳
1回答
我试图用Python来求解一个耦合的一阶ODEs系统。我是新手，但到目前为止，来自的SciPy.org一直是一个很大的帮助。
在我的例子中，一个重要的区别是用于“驱动”我的ODEs系统的输入数据在不同的时间点突然改变了，我不知道如何最好地处理这个问题。下面的代码是我能想到的最简单的例子来说明我的问题。我欣赏这个例子有一个直截了当的解析解，但是我的实际ODEs系统更复杂，这就是为什么我试图理解数值方法的基础。
简化示例
假设水桶底部有一个洞(这种“线性水库”是许多水文模型的基本构件)。桶的输入流量为R，孔的输出为Q。假设Q与桶中的水量V成正比。比例的常数通常被写成
，其中T是商店的“停留时间”
浏览 3提问于2015-08-26得票数 4
回答已采纳
1回答
我正在使用Pyomo求解MINLP，它采用整数输入的向量(图的链路容量)，并在edge_enum字典所做的network.What中返回丢包率，它枚举图的所有边。例如， {0: (0, 4),
1: (1, 4),
2: (1, 5), ... a=0
for e in G.edges():
edge_enum[a]=e
a+=1
model=ConcreteModel
key_list=list(edge_enum.keys())
edge_list=list(edge_enum.values())
model.A=Set(initialize=edge_list
浏览 60提问于2020-03-25得票数 0
这个问题涉及到Python和Mathematica代码。
我想用数值方法解决一个复杂的常微分方程系统。最初，我使用NDSolve，但是因为我最终需要进行大量的参数搜索，所以我决定改用Python，特别是SciPy的solve_ivp函数。然而，我很快就遇到了以下问题：
A、B、C和D是我想要求解的四个耦合的复常微分方程。值得注意的是，abs(A)^2 + abs(B)^2 + abs(C)^2 + abs(D)^2 =1总是成立的(由于问题的性质)。实际上，在NDSolve的解决方案中总是如此
。
但是，当我将代码转移到Python时，使用solve_ivp得到的结果是不正确的。
浏览 3提问于2020-07-12得票数 1
我想建立一个线性规划(如果可能的话，我不确定)在R中解决这个问题：
我想把功能最小化：
abs(x1) + abs(x2) + abs(x3) + abs(x4)
制约因素：
x1 + x2 + x3 + x4 = 0
0.2x1 + 0.3x2 + 0.5x3 + 0.1x4 = 0.2
约束通常被写成不等式，所以我假设每个方程都有两个方程，一个是<= 0和<= 0.2，另两个是>= 0和>= 0.2。决策变量是连续的，可以是正的，也可以是负的。
我可以为约束设置使用lpSolve的矩阵，但对于如何在目标函数中实现绝对值函数，我感到非常困惑。这个站点()是一个很好的读
浏览 3提问于2021-01-16得票数 1
回答已采纳
R语言能找到一阶微分方程的一般解吗？
例如：
(5x-6)^2 y' = 5(5x-6) y - 2
PS：
可以用手解决的问题，即特定的解决办法是：
y = 1/(5(5x-6))
和泛型
C*(5x-6)
我需要弄清楚R能不能做到这一点？
浏览 2提问于2018-05-08得票数 3
回答已采纳
我有一个2d网格，其中某些单元格值产生不同的分数，我想通过给单元格分配值来最大化分数。在求解过程中，是否有可能跟踪求解者尝试过的每个网格？下面是模型的一个片段，只是为了演示这个想法。
int: i_count;
int: e_count;
set of 1..3: ROW;
set of 1..3: COL;
enum CELL = { U, I, E };
array[ROW,COL] of var CELL: grid;
constraint sum(r in ROW, c in COL)(grid[r, c] = I) <= i_count;
constraint sum(
浏览 3提问于2020-11-15得票数 3
回答已采纳
1回答
由于R对函数参数的求值，可以指定一组一致的输入参数，并自动计算其他参数。
考虑下面的函数，在化学中将浓度、质量、体积和摩尔重量联系起来，
concentration <- function(c = m / (M*V), m = c*M*V, V = m / (M*c), M = 417.84){
cat(c("c=", c*1e6, "micro.mol/L\n",
"m=", m*1e3, "mg\n",
"M=", M, "g/mol\n",
浏览 2提问于2011-08-26得票数 1
我的问题是，我没有看到unsat_core使用api跟踪在整个块中提供的任何断言。
f = Z3_parse_smtlib2_string(c, "unsat_core_example1.smt2",0,0,0,0,0,0);
params p(c);
p.set(":unsat-core", true);
s.set(p);
// enabling unsat core tracking
expr r = to_expr(c, f);
unsat_core_example1.smt2：
(declare-fun p () Bool)
(declare-fun
浏览 0提问于2014-04-11得票数 1
所以我这里有个代码(针对一个电气问题)
R_1= input('Enter the value of R_1 ');
R_2= input('Enter the value of R_2 ');
R_3= input('Enter the value of R_3 ');
V_1= input('Enter the value of V_1 ');
V_2= input('Enter the value of V_2 ');
V_3= input('Enter the value of V_3 ')
浏览 1提问于2017-06-30得票数 0
我正在使用poplib从POP3服务器获取电子邮件。
但是发生了这个错误：
*cmd* 'USER xxx@gmail.com'
*put* b'USER xxx@gmail.com'
*get* b'+OK send PASS\r\n'
*resp* b'+OK send PASS'
*cmd* 'PASS asdasd\n'
*put* b'PASS asdasd\n'
*get* b'-ERR bad command qj1pf12899100pbb.19\r\n'
*resp
浏览 0提问于2012-12-06得票数 1
回答已采纳
我想使用最优brent寻根算法计算隐含波动率：
def calcimpliedvol(S,K,T,r,marketoptionPrice):
d1=(np.log(S/K)+(r-0.5*sigma**2)*T)/(sigma*np.sqrt(T))
d2=d1-(sigma*np.sqrt(T))
BSprice_call=S*si.norm.cdf(d1,0,1)-K*np.exp(-r*T)*si.norm.cdf(d2,0,1)
fx=BSprice_call-marketoptionPrice
return optimize.brentq(fx,0,1,maxi
浏览 3提问于2021-03-15得票数 1
回答已采纳
我有一个Simulink R2015b模型包含SimElectronics和SimMechanics块。这种模式是可行的。我的一位同事要求将模型导出到R2015a。但是，当我执行这个模型时，我会得到以下错误消息：
组件: Simulink类:模型错误
*“型号名称”：没有为一个或多个Simulink转换器块提供足够的输入导数，以供所选择的求解器使用。隐式求解器(ode23t、ode15s和ode14x)通常比显式求解器需要更少的输入导数，而本地求解器不需要任何输入导数。您可以通过在Simulink-PS转换器块的输入处理选项卡上选择不同的选项(例如，打开输入过滤)来提供更多的输入衍生物：.‘
浏览 2提问于2016-03-11得票数 2
回答已采纳
1回答
我试图用Z3py优化一个基于最小化的集覆盖问题(SCP41)实例。
结果如下：
使用
(1)我知道Z3支持优化()。很多时候，我在SCP41和其他实例中达到了最佳状态，但也有一些没有达到最佳效果。
(2)我理解，如果我使用Z3py API而不使用优化模块，我将不得不执行@Leonardo在()中描述的典型的顺序搜索。它给不了我任何结果。
我的方法
(3)我试图通过实现类似于解释@Philippe ()的二进制搜索来改进顺序搜索方法，当我运行我的算法时，它等待并且没有得到任何结果。
我知道二进制搜索应该更快，在这种情况下可以工作吗？我还知道实例SCP41很大，产生了许多限制，并且变得极其组合，这
浏览 0提问于2019-06-25得票数 1
回答已采纳
2回答
如何将蓝色曲线值转换为线性(红色曲线)？我正在excel中做一些测试，但基本上，我在3D应用程序中有那些蓝线值，我想用python来操作这些值，这样我就可以将这些值线性化。有什么数学方法我错过了吗？
X轴从0到90，y轴从0到1。
例如:在图的中间，蓝线给我一个"0,70711“的值，我知道在线性上它是"0,5”。我想知道是否有一个简单的公式可以把所有输入的非线性值转换成线性的。
我不知道什么是“公式”创造的非线性蓝线，也忽略了黄线，因为我只是试图“反向工程师”，看看是否会引导我得出任何结论。
谢谢
浏览 5提问于2014-06-23得票数 1
回答已采纳
我正在尝试使用python为lambdaλ求解这个方程。 ? 因为，对于索引列中的i，数据帧中给出了P_0、c、t、D、F，如下所示： ? 在给定r = constant = 0.6的情况下，Excel Solver可以方便地求解λ方程。但我是python的新手，我一直在努力解决python症状，以找到解决它的方法。请帮帮我！非常感谢！ 编辑:这是我的数据帧dataframe 这是我的python代码。它需要永远运行，而Excel Solver只用了不到1秒，所以我认为它一定是出了什么问题。 import sympy as sy
x = sy.Symbol('x'
浏览 14提问于2020-07-24得票数 0
2回答
给定柯西问题
Y '= (y-1) ^2 (y-2) (y-a)且y(1) = y0；
其中a= 0、a=1或a= 2：
写一个程序，在不同的情况下绘制问题的解决方案。它必须遵循以下步骤: a)程序将输入变量。
b)将视口划分为三个子地块，位于两个...第一个和第二个中的两个。如果第一个子图选择a= 0，第二个子图选择a=1，第三个子图选择a= 2，则使用命令。
c)创建固定大小为1，10 x -1，3的图形轴。使用命令ode45和plot draw in red fixed points。放置轴名称和标题图，指示我们(a = 0，a=1或a= 2)的情况。
这是我的决心：
funct
浏览 1提问于2013-04-06得票数 1
1回答
我想用R模拟一个\周期的ARMA(1,1)，我找到了R包perARMA，但我不知道如何使用它。有允许模拟parma(1,1)的函数makeparma。但我不明白输入参数和用来模拟周期性物质的模型。
这是由试图模拟parma(2，1)的包提供的源代码：
T=12
nlen=480
p=1
a=matrix(0,T,p)
q=1
b=matrix(0,T,q)
a[1,1]=.8
a[2,1]=.3
phia<-ab2phth(a)
phi0=phia$phi
p
浏览 1提问于2018-12-19得票数 1
1回答
带数值的线性方程的非零解
、、、
如何用Numpy求这类方程的非平凡解？
r1 = r1 * 0.03 + r2 * 0.88 + r3 * 0.2425 + r4 * 0.03 + r5 * 0.03
r2 = r1 * 0.455 + r2 * 0.03 + r3 * 0.2425 + r4 * 0.03 + r5 * 0.88
r3 = r1 * 0.455 + r2 * 0.03 + r3 * 0.03 + r4 * 0.03 + r5 * 0.03
r4 = r1 * 0.03 + r2 * 0.03 + r3 * 0.2425 + r4 * 0.03 + r5 * 0.03
r5 = r1 * 0.03
浏览 0提问于2012-11-09得票数 1
回答已采纳
我尝试模拟未知积分中的值(创建一个气候预报器)
我的函数是：$\int_{x = 0}^{x = 0.25} 4*y^(-1/x) dx$
通常情况下，输入变量y并获得一个值作为输出。但是，我想输入此积分等于的值，并获得y的值作为输出。我有3个长度分别为1000、10000和100000(值在0和1之间)的runif向量，我将它们用作输入值。假设第一个值是0.3，第二个值是0.78
我想计算哪个y，上面的积分等于0.3 (或者第二个值等于0.78 )。
我怎么才能在R中做到这一点呢？我已经尝试了一些关于集成函数的东西，但是我需要一个y值才能使其工作
浏览 10提问于2020-05-18得票数 0
1回答
我有一个代码来执行优化来推断一个参数：
import numpy as np
from scipy.integrate import odeint
import matplotlib.pyplot as plt
from scipy.optimize import root
from scipy.optimize import minimize
import pandas as pd
d = {'Week': [1, 2,3,4,5,6,7,8,9,10,11], 'incidence': [206.1705794,2813.420201,1
浏览 10提问于2022-02-03得票数 0
回答已采纳
因此，我试图找出如何使用R中的qr.solve命令来根据特定的数据集求解最佳匹配线(斜率和截距)的组件。虽然我已经知道如何使用lm函数来完成这个任务，但据我所知，矩阵的qr分解可能产生同样的结果。有人能解释一下这样的事情是如何正确的吗?也许它比一个简单的线性模型命令更有好处呢？qr.solve函数的输入可以在R中实现，还是必须先自己求解，然后再插入R中呢？
我尝试将我的数据集(每行2列表示图上的点)输入到一个矩阵中，并在函数qr.solve中使用它作为x参数。但是，我不太确定该为b插入什么。
浏览 6提问于2022-11-26得票数 1
回答已采纳
#define _CRT_SECURE_NO_WARNINGS
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#define SIZE 11
double cubicspline(double val);
void TDMA(void);
double x[SIZE] = { -1.0, -0.8, -0.6, -0.4, -0.2, 0.0, 0.2, 0.4, 0.6, 0.8, 1.0 };
double y[SIZE] = { 0.038, 0.058, 0.1, 0.2, 0.
浏览 0提问于2018-04-23得票数 0
为了解决SAT问题，我决定使用微软和Python3的Z3-solver。目的是使用一个很长的模型(多达500,000个特征)并找到所有可能的解决方案。为了找到它们，我想将第一个解S1添加到初始方程中，并排除它，以此类推。我将使用while循环来完成此操作。解决SAT问题对我来说很重要，因为我想分析特征模型。
但是我遇到了一个问题，就是把某物加到初始方程式中。我将分享一个最小的例子：
# Import statements
import sys
sys.path.insert(0,'/.../z3/bin')
from z3 import * # https://
浏览 3提问于2016-09-15得票数 0
我有一个非线性系统，它有二次输入(操纵变量，MV)。为了使用GEKKO对该系统进行MPC模拟，我定义了以下函数来处理二次输入： def NTVin(T, u):
n = T.shape[0]
Tntv = zeros(n)
u2 = zeros_like(u, dtype = float64)
for i in range(len(u2)):
u2[i] = copy(u[i].VALUE)
for i in range(n):
Tntv[i] = dot(dot(u2, T[i,:,:]), u2)
retu
浏览 8提问于2019-10-02得票数 3
回答已采纳
以下是问题所在：
常数：
受抚养价值函数：
def v_ego(self,ego_vel_x,ego_a_ini_x,t):
v_ego = ego_vel_x + ego_a_ini_x * t
return v_ego
def x_obj(self,x_obj_ini,obj_vel_x,obj_a_ini_x,t):
x_obj = x_obj_ini + obj_vel_x * t + 0.5 * obj_a_ini_x * t ** 2
return x_obj
def x_ego(self,x_ego_ini,ego_vel_x,ego_a_i
浏览 7提问于2022-11-22得票数 0
1回答
我正在尝试使用Pyomo解决一个优化问题。为此，我需要声明一个二维向量和一个三维向量。这两个向量只能存储0,1。
Sri和Xrij
R=3 V=8 1<=i，j<=V 1<=r<=R
我尝试在pyomo中使用range
model.IDXV = range(v+1)
model.IDXR = range(r+1)
model.x=Var(model.IDXR,model.IDXV,model.IDXV,within=Binary,initialize=0)
model.s=Var(model.IDXR,model.IDXV,within=Binary,initiali
浏览 83提问于2018-08-12得票数 0
回答已采纳
我面临着一个在目标函数中乘以一个项的问题。我的目标函数是： \begin{equation}
min \sum_{t=1}^T p_{t} * (1+r) * y_{t}
\end{equation} 其中$p$是价格，$y$是数量(决策变量)，$r$是价格附加费。$t$是时间段的索引。 问题是$r$的值取决于$y$ (决策变量)的值位于哪个区间之间，如下面的分段函数所示。 \begin{equation}
\r_t (Y_t) =
\begin{cases}
% s_0
0 & \mathrm{if }\; Y_t
浏览 83提问于2021-09-20得票数 0
1回答
我试图通过输入以下命令来寻找多项式P的根
solx = solve(P)
哪里
P = p(p(p(p(x)))
使用p(x) = r*x*(1-x)
但是，我所有的root都被返回为
RootOf(Q)
其中Q是P的复多项式。
我如何返回完成分析所需的确切符号解(我需要找到这些根的导数)？
完整代码：
syms r x;
p = symfun(r*x*(1-x),x);
P = p(p(p(p(x))));
[Q,R] = quorem(P-x, p(p(x))-x)
solve(Q)
浏览 4提问于2016-03-07得票数 1
我正在尝试用我的研究数据进行非线性回归。因为它是非线性的，所以我不能使用Simplex LP。相反，我做的是GRG非线性，所有参数都有上下界。
这很奇怪，因为Excel有时会产生正确的答案，而有时会产生错误的答案。我必须手动将参数更改为任意数字，然后再次运行Solver，希望它是正确的。让我给你看看我的excel表格。
基本上，我希望使用两条线性线来拟合我拥有的数据。原始数据可以分为两部分，都是线性的。这两条线相交的点是临界值。
我的原始数据的正确输出是临界值= 0.006707。我最后一次运行它时，在excel工作表上，您可以看到该临界值= 3.36E-06。
如果它能帮助你更好地理解这一
浏览 1提问于2014-06-27得票数 0
1回答
下面是我的合并排序实现，它所做的是接受一个整数对数组，并根据对中的第二个元素对它们进行排序(第一个元素在配对中断开了关系，所有元素都是不同的)
void mer(vector<pair<int, int>> a, int l, int m, int r, vector<pair<int, int>> res) {
int i = l;
int j = m;
int k = l;
while (i < m && j < r) {
if (a[i].second < a
浏览 2提问于2020-10-11得票数 1
回答已采纳
1回答
我正在尝试解决一组模拟淀粉酶(酶)水解淀粉的常微分方程。当我试着解这组方程时，我得到一个
lsoda-- at current t (=r1), mxstep (=i1) steps
taken on this call before reaching tout
In above message, I1 = 500
In above message, R1 = 0.6333483931400E+00
Excess work done on this call (perhaps wrong Dfun type).
Run
浏览 0提问于2014-07-17得票数 2
1回答
在sas 9.4中是否有类似于matlab中求解函数的函数，我可以用它来求解方程：
syms x sig1 sig2 mu1 mu2;
solve(1/sig1/sqrt(2*pi) * exp(-1/2*((x-mu1)/sig1)^2) == 1/sig2/sqrt(2*pi) * exp(-1/2*((x-mu2)/sig2)^2), x)
或者在sas中寻找2条或2条以上直线的交点最有效的方法是什么。
谢谢
浏览 7提问于2014-09-15得票数 0
回答已采纳
2回答
我试图在MATLAB中找到使用fmincon的最小值，并面临以下问题：
优化完成是因为初始点梯度的大小小于函数公差的默认值。
我的目标函数的曲面显示“步骤”，因此它在某些输入变量范围内具有相同的值(如果我正确的话，梯度的大小是零)：
当从初始点移动时，求解器不会看到目标函数的值发生任何变化，并完成了优化：
Iteration Func-count f(x) Step-size optimality
0 3 581.542 0
Initial poi
浏览 3提问于2014-11-13得票数 3
回答已采纳
1回答
正如标题所示，我正在寻找一个类似的函数，或者是关于如何实现deval在MATLAB中使用Python所做的相同事情的建议。
代码：
import numpy as np
from scipy.integrate import odeint
import matplotlib.pyplot as plt
#define model function that returns ds/dt, di/dt/ and dr/dt
def epidemic_model(state, t):
s, i, r = state
d_s = -0.06 * s * i
d_i =
浏览 1提问于2018-05-27得票数 1
本文尝试用牛顿法求方程=tanh(zJ/k_B _ T)。当我运行程序时，我得到了一个错误：
图:输入参数#2和#3:不兼容维度的错误大小。
还有一张空白图。请帮我解决我代码中的问题。
z=4\\ no. of nearest neighbours
m=1;//value of J/K
T=[0.1:0.1:8]\\value of temerature
s(1)=-0.5;//initial value
n=100;\\no. of iterations
for i=1:n \\ running of for loop
f= s(i)-tanh(z*m*s(i)./
浏览 1提问于2020-05-24得票数 0
我正在使用scipy.optimize.minimize来解决一个有效的投资组合。
使用默认设置时，我经常在使用某些输入时遇到"BaseException:用于行搜索的正向导数“错误。我注意到，如果我将容差设置得足够高，问题会变得更少，但不会消失。有什么建议吗？
import numpy as np
import pandas as pd
import scipy.optimize
def fx(TOLERANCE):
#TOLERANCE = 1.5
def solve_weights(R, C, rf, b_):
def port_mean_v
浏览 26提问于2014-06-03得票数 2
此Meteor客户端代码需要按字段“类别”对mongo集合中的文档进行分组，并返回字段“类别”的唯一值的排序数组。我使用的是Blaze模板帮助程序，它返回一个由collection.find.foreach填充的数组，因为我无法很容易地找到关于如何聚合客户端的mongo集合的示例，如何让该助手等待收集操作准备就绪？谢谢
'jobs': function(){
let categories = []
Jobs.find({}).forEach((doc) => {
if(!categories.includes(doc.category
浏览 1提问于2021-05-15得票数 0
对于输入，我已经设置了
设置对象单元格：$A$1 (即。可能是18)
致: Min
通过改变变量$B$1
选中“使无约束变量为非负变量”：
求解方法:非线性
我不是这里的数学专家，但我试图计算我应该使用哪一种无约束最小化函数？看起来像fminunc --问题是：
我只是在objfun中使用Matlab min函数吗？
(如所解释)
将会
fminunc(@objfun，x0)
如果xo相当于上面的Excel中的$B$1？
我会假设，没有任何额外的选择，以及获得预期的结果。至少我没有看到任何控制负输出的选项。
这听起来对吗？谢谢
浏览 1提问于2014-10-07得票数 1
回答已采纳
我有四个号码: 4，3.6，1.8，1.2，
有没有办法找出每一个的权重，这样当我把(数字*它们的权重)相加时，我会得到一个特定的数字，例如: 2。
因此，对于数字2，结果将是：
2是输入，输出是权重。我知道有几个权重可以回答这个问题，但我想知道我是否能得到答案。
有方法在Excel中计算吗？
谢谢
浏览 2提问于2020-06-08得票数 0
回答已采纳
3回答
假设我有一个在空间需求中是nlogn的函数，我想为给定的可用空间求出该函数的最大输入大小。也就是说，我想找n个地方，nlogn=c。
我跟随计算n，这在R中是这样的：
step = function(R, z) { log(log(R)-z)}
guess = function(R) log(log(R))
inverse_nlogn = function(R, accuracy=1e-10) {
zi_1 = 0
z = guess(R)
while(abs(z - zi_1)>accuracy) {
zi_1 = z
z = step(R, z)
}
exp(
浏览 2提问于2010-12-09得票数 2
1回答
时变输入ODE求解
、、、、
我试图用ode45来求解一组ODE方程。我的一些参数已经是时间的函数，但我总是会出错。
function odo
dx(1,1) = (vl+vr)/2*cos(x(3));
dx(2,1) = (vl+vr)/2*sin(x(3));
dx(3,1) = obz
其中obz、vr和vl都是向量，例如：
obz = sin(t), t = 0:dt:tf;
我使用以下语法：
[t, x1] = ode45(@(t,x) odo(t,x,b,obz,vr,vl), 0:dt:tf, [0;0;0]);
使用R15，但继续获取错误：
赋值的非单例rhs维度比非单例rhs的下标要多。
浏览 0提问于2019-04-13得票数 0
我正在尝试实现本参考文献中描述的数学模型：
刘查德，袁力，范远。凤凰崛起路径的数学建模。Plos One，计算生物学，第10卷，2014年。
我希望输出变量dPLAdt的数值计算符合Heaviside函数，它计算变量C3和C7。从下面的代码中可以看到，Heaviside函数将在1440年之后被激活。该程序正在运行，但没有产生与纸所附图片相匹配的结果。我想和报纸的结果相匹配。
这是我的代码：
1)功能
function dPLAdt = PLAprod(t,PLA)
dPLAdt = zeros(2,1);
k1 = 2.8*10^(-5);
k3 = 3.24*10^(-6);
浏览 16提问于2014-08-26得票数 2
Parameter st(t) starting time of curfew on airport t
/11 540
22 540
33 540
44 540
55 540/;
Table arr(i,t) arrival
11 22 33 44 55
101 0 1 0 0
浏览 0提问于2020-07-09得票数 0
1回答
我需要解决一个问题:假设湖中的鱼类种群P(t)在时间t= 0被疾病攻击，导致鱼类停止繁殖(因此出生率为β= 0)，死亡率δ(每周每条鱼的死亡人数)与1/√P成正比。如果最初湖中有900条鱼，6周后剩下441条鱼，那么湖中所有的鱼需要多长时间才会死亡？我不明白如何添加条件P (6) = 441？ r=np.array([6])
p=np.array([441])
def fun (y , t, r, k):
P=y
dydt = -k*math.sqrt(P)
return dydt
y0 = np.array([900])
t = np.linspace(1, 9,
浏览 8提问于2019-08-29得票数 1
1回答
我想要解一个由七个耦合二阶微分方程组成的边值问题。有7个函数，y1(x),...y7(x)，每个函数都由一个形式的微分方程来描述。
d^2yi/dx^2 = -(1/x)*dyi/dx - Li(y1,...,y7) for 0 < a <= x <= b,
其中Li是一个给出y1,...,y7线性组合的函数。我们给出了一阶导数dyi/dx在x=a和函数yi在x=b的边界条件。
dyi/dx(a) = Ai,
yi(b) = Bi.
因此，我们可以把它重写为一个由14个耦合的一阶ODE组成的系统：
dyi/dx = zi,
dzi/dx = -(1/x)*zi - Li(y1
浏览 6提问于2022-07-15得票数 1
回答已采纳
给出了由径向畸变和切向畸变定义的针孔模型的畸变系数D = k1，k2，p1，p2，k3。如何用这些系数来解构图像呢？我知道cv2.undistort()函数。有人能解释我如何实现非失真函数来执行迭代最小化吗？
# TODO: Iterative minimization required
r2 = x_ * x_ + y_ * y_
# Remove tangential distortion
x_tangential = x_ - (2 * distort_coeffs[2] * x_ * y_ + distort_coeffs[3] * (r2 + 2 *
浏览 3提问于2019-02-08得票数 1
回答已采纳
这是我定义的函数的全部代码 def solve(eq, var=('x', 'y')):
import re
var_re = re.compile(r'(\+|\-)\s*(\d*)\s*\*?\s*(x|y)')
const_re = re.compile(r'(\+|\-)\s*(\-?\d+)$')
constants, eqns, coeffs, default = [],[], {'x': [], 'y': []}, {'': &#
浏览 28提问于2020-09-25得票数 0
回答已采纳
1回答
应该有一个与每个约束相关联的对偶解决方案。
我有len(E)+len(D)约束，它等于84+20 = 104
但我只有23个二重值？
我的输入数据如下：
#Input data
from pyscipopt import Model, quicksum
def parameter_a(E,R,D):
a = [[[0 for d in range(len(D))] for r in range(len(R))] for e in range(len(E))]
for e in range(len(E)):
for d in range(len(D)):
浏览 19提问于2020-03-05得票数 0
回答已采纳
1回答
我正在尝试z3中的递归函数，我很好奇是否存在模型构造中的一个bug。考虑：
(define-fun-rec f ((x Int)) Int
(ite (> x 1)
(f (- x 1))
1))
(check-sat)
(get-value ((f 0)))
这里，f实际上是常量函数1，只是以一种愚蠢的方式定义的。对于此输入，z3打印：
sat
(((f 0) 0))
这似乎是不正确的，因为f 0应该等于1。
有趣的是，如果我断言z3的结果是什么，那么我就得到
浏览 0提问于2017-02-11得票数 0
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
实时音视频
对象存储
云直播活动推荐
腾讯云开发者
