在R中为L求解方程
的一种方法是使用R的内置函数solve()。solve()函数可以用于求解线性方程组或矩阵的逆。
假设我们要求解以下线性方程组： Lx = b
其中，L是一个n x n的矩阵，x是一个n维向量，b是一个n维向量。
首先，我们需要定义矩阵L和向量b。在R中，可以使用matrix()函数定义矩阵，使用c()函数定义向量。例如：
L <- matrix(c(1, 2, 3, 4), nrow = 2, ncol = 2) # 定义一个2x2的矩阵L
b <- c(5, 6) # 定义一个2维向量b接下来，我们可以使用solve()函数来求解方程。例如：
x <- solve(L, b) # 求解方程Lx = b，得到向量x最后，我们可以打印求解得到的向量x的值。例如：
print(x) # 打印向量x的值关于R中求解线性方程的更多信息，你可以参考腾讯云的产品文档：R语言求解线性方程组。
相关·内容
1回答
我读过，还很难理解Matlab求解方程中多变量的顺序。
例如
syms f1(R,L) f2(R,L);
f1 = (80+R)^2 + (120*pi*L)^2 - 232^2;
f2 = R^2 + (120*pi*L)^2 - 176^2;
[R,L] = solve (f1,f2);
fprintf('Resistor is %s Ω\n',round(R,3));
fprintf('Inductance is %s H\n',round(L,3));
如果将变量L替换为Z，则结果将是相反的。这似乎是字母顺序，我不知道确切的规则是什么。
浏览 9提问于2021-03-18得票数 2
回答已采纳
5回答
我有六个参数方程，使用了18个不同的变量(实际上不是26个)，其中6个是未知的。
我可以坐下来几页纸，计算出每个未知数的方程式是什么，但有没有一个简单的程序解(我正在考虑在Matlab中)，会吐出我正在寻找的六个方程式？
编辑：耻辱这已经关闭了，但我想我能理解为什么。如果有人还感兴趣，方程(我相信)是非线性的：
r11^2 = (l_x1*s_x + m_x)^2 + (l_y1*s_y + m_y)^2
r12^2 = (l_x2*s_x + m_x)^2 + (l_y2*s_y + m_y)^2
r13^2 = (l_x3*s_x + m_x)^2 + (l_y3*s_y + m_y)^2
浏览 3提问于2009-12-11得票数 3
回答已采纳
我有一个类似下面的matlab格式的微分方程。
syms x y m g r l J
% x,y are variables, the others are constant
1: 0.5*m*(r^2*x^2+l^2*(Dx-Dy)^2+2*r*l*Dx*(Dx-Dy)*cos(y))+0.5*J*(Dx-Dy)^2=m*g*
(l*sin(x-y)-r*(1-cos(x)));
2: J*(D2x-D2y)+l^2*(D2x-D2y)-r*l*(Dx)^2*sin(y)+r*l*D2x*cos(y)-m*g*l*cos(x-y)=0,
3: x(0)=pi/2,y(0)=pi/2,
浏览 8提问于2016-09-02得票数 0
我用matlab编写了这段代码
syms L1 L0 P1 P0 s
L = expand(s*(s^2+L1*s+L0)*(s-75)+1400*(P1*s+P0))
R = expand((s+75)*(s+150)^3)
我想解决P1、P0、L1和L0的问题
L = R
1400*P0 - 75*L0*s + 1400*P1*s + L0*s^2 - 75*L1*s^2 + L1*s^3 - 75*s^3 + s^4 = s^4 + 525*s^3 + 101250*s^2 + 8437500*s + 253125000
通过观察：
1400*P0 = 253125000
等等。
如何从
浏览 0提问于2014-09-10得票数 0
我试着解决最大化函数的问题。14个约束，15个变量。有些是线性的，有些则不是。它们都是方程(没有不等式)。
我试过使用'fsolve‘和’求解‘，使用拉格朗日(最终得到29个方程和30个变量)
我转到了fmincon。我在一个名为objectfun.m的文件中设置了一个带有目标函数的脚本：
function f = objectfun(x,I,rho)
% SWAPPING VARIABLE NAMES FOR READABILITY:
w = x(1);
t = x(2);
beta = x(3);
r = x(4);
% VALUE FUNCTION TO BE MINIMIZ
浏览 0提问于2016-02-27得票数 0
谁能解释一下如何在SparkR操作glm中指定优化方法？当我尝试用glm拟合OLS模型时，我只能指定"normal"或"auto"作为求解器类型。SparkR不能解释求解器规范"l-bfgs"，这使我相信当我指定"auto"时，SparkR只是假定"normal“，然后使用LS法方程解析地估计模型系数。
用随机梯度下降和L-BFGS拟合GLM在SparkR中是不可用的，还是我写错了下面的评估？
m <- SparkR::glm(y ~ x1 + x2 + x3, data = df, solver = "l-
浏览 1提问于2016-07-12得票数 2
1回答
我有这样的方程式：
wl[n]= w[n] + Δx*v[n]
ϕl[n]= ϕ[n] + Δx*ρ[n]
ρl[n] = ρ[n] - Δt*fρ(ρ,v,w,n)
vl[n] = v[n] - Δt*fv(ρ,v,w,Δx,n)
模拟这些方程：
我可以使用ODE解算器在时间上，然后在空间中进行解算吗？或者，有没有Julia软件包可以帮助我使这些方程无条件稳定？
浏览 6提问于2021-11-19得票数 1
我想知道在mathematica中是否有一种方法，可以在给定其他变量和一组方程的情况下，解决特定变量的问题。基本上有6个变量，我得到了其中的3个，必须使用这些方程来计算其他的-
变量-
B,Qs,f0,R,c,L
方程式-
f0=1/(2*Pi*Sqrt[L*c])
Qs=(w*L)/R
w=2*Pi*f0
B=f0/Qs
我们得到了这些变量中任何3个变量的值，并必须使用这些值计算出其余的变量。
我想也许可以使用消除，但我不确定这将是如何构造的，因为我以前只使用过设置不变的变量和单个输出。
浏览 1提问于2013-06-02得票数 3
回答已采纳
Introduction
注:代码和指示：
单d是一阶导数双d是二阶导数
通过数值求解控制LaGrange方程，用Matlab对一些动态系统进行了数值模拟。基本上是一组二阶常微分方程。我正在使用ODE45。我从Mathworks (下面的教程链接)上找到了一个关于如何求解基本的二阶常微分方程组的很好的教程。
根据本教程，我通过获得两个二阶常微分方程(一个用于角θ，另一个用于弹簧伸长)来模拟弹性弹簧摆的运动，如下所示：
theta双素数方程：
M*thetadd*(L + del)^2 + M*g*sin(theta)*(L + del) + M*deld*thetad*(2*L + 2*de
浏览 4提问于2017-01-25得票数 0
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
我有以下问题：
最大CEQ(w) s.t.W在(0,1)中，除了由形式为CEQ(w) = F( CEQ(w) )的不动点方程给出外，我对CEQ(W)一无所知。如果我固定一个w，我可以用fzero函数求解不动点方程，并得到一个CEQ的值。如果我选择一个不同的w，我会得到另一个CEQ的值。因此，我可以遍历w的所有可能值，然后选择CEQ最高的值。这似乎是糟糕的编程，我想知道我是否可以在MATLAB中更有效地实现这一点:我想将不动点方程的解建模为w的函数，但我不知道如何实现它。
更准确地说，下面是一个示例代码：
clear all
clc
NoDraws = 1000000;
T_hat = 12;
浏览 5提问于2015-04-30得票数 2
1回答
我对Python比较陌生，我试图用它来解一个二阶非线性微分方程，特别是电解质中的Poisson-Boltzmann方程。
phi''(r) + (2/r)*phi'(r) = (k^2)*sinh(phi(r))
从本质上讲，它描述了静电势(phi)在电解质中远离带电表面的衰减，衰减速率受参数k的控制。
phi(r) -r处的电位
酸酐(R)-r处电位的导数
距离曲面的r-距离(在这种情况下，我正在求解r=1到r=R)。
和边界条件
phi(1) =5
伤寒杆菌(R)=0
代码的问题位如下
from scipy.integrate i
浏览 4提问于2014-12-12得票数 0
1回答
当我运行这个程序时，我在最后没有得到解决方案，但是应该有一个解决方案(我相信)。知道我做错了什么吗？如果你把Q从e2方程中去掉，它看起来是正确的。
#!/usr/bin/python
from sympy import *
a,b,w,r = symbols('a b w r',real=True,positive=True)
L,K,Q = symbols('L K Q',real=True,positive=True)
e1=K
e2=(K*Q/2)**(a)
print solve(e1-e2,K)
如果我们做以下工作的话，这是可行的：
设置Q=
浏览 4提问于2013-10-03得票数 2
2回答
如何使点轨道成一条线，三维
、、、、
本质上，我想找出环的方程，它在空间中绕点旋转，这个环垂直于离这个点的法线。
我有一条线，以2分的形式，L1，L2；
L1->L2，N正常；
我有一个与L1->L2正常的平面，L1就在这个平面上。ax + by + cz = d；
我的半径远离L1，R；
->，我想做一个点，，V，，绕点L1转这条线；
我想我必须在这个平面上建立一个以L1为原点的圆形方程。我不知道如何在三维平面上绘制二维方程。
或者有人知道如何用另一种方式来做，交叉产品或其他什么？
浏览 0提问于2014-04-08得票数 0
回答已采纳
1回答
我有以下方程，并希望R为L求解。 有什么想法吗？ Average = 370.4
m = 2
p = 0.2
n = 5
#L = ?
log10(Average) = 0.379933834 -0.107509315* m + 0.104445717 * p + 0.016517169 * n -0.025566689* L + 0.014393465 * m * p + 0.001601271 * m * n - 0.014250365 * n * L + 0.002523518 * m^2 + 0.237090759 * L^2
浏览 15提问于2021-11-10得票数 0
2回答
给定L和U 以及常数向量b，例如LU*x=b，是否有内置的函数可以找到x？意思是-
X = functionName(L,U,b)
请注意，在L和U中，我们都在处理三角形矩阵，这些矩阵可以通过直接求解，而无需使用过程。
编辑:
求解这个线性方程组的步骤如下-
1. define y - s.t Ux=y
2. solve Ly=b by forward substitution
3. solve Ux=y by backward substitution
4. return y
编辑2 :
我找到了，但我没有尝试，因为我的版本太旧了(R2010a)。对任何人都有用吗？
浏览 1提问于2013-05-28得票数 6
回答已采纳
1回答
我试图在SageMath中求解以下两个微分方程(数值)：
我的目标是获得M(r)-r.的图
我尝试了以下代码：
sage: r = var('r')
sage: M = function('M')(r)
sage: a = function('a')(r)
sage: de1 = (M*a*a*diff(M,r) + (M*M*a+6*a)*diff(a,r) + 1/(r*r) == 0)
sage: de2 = (a*r*diff(M,r) + 7*M*r*diff(a,r) + 2
浏览 3提问于2020-09-07得票数 0
回答已采纳
我想建立一个线性规划(如果可能的话，我不确定)在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
回答已采纳
我很难把数学公式转换成代码。我可以很容易地在纸上解决它，但是我很难把它变成代码形式。一般情况下，是否有可能将数学公式直接引入程序，或者，如果不能，如何将其传输？
例如，毕达哥拉斯定理：a² + b² = c²。在代码中定义了a和b，直接应用这个公式。
c = sqrt(pow(a, 2) + pow(b, 2))
目前，我很难把我在纸上解决过的直线交汇处放到我的程序中。假设一条线被定义为L: O + a * D，则可以通过将两行设置为相等来找到它们的交集。
L1: O1 + a * D1
L2: O2 + b * D2
=> O1 + a * D1 = O2 + b * D2
用笔一张
浏览 0提问于2013-04-25得票数 3
回答已采纳
1回答
数字数学
导论
这是“你好，世界！”偏微分方程(偏微分方程)。拉普拉斯或扩散方程在物理学中经常出现，如热方程、变形、流体力学等。现实生活是3D的，但我们想说“你好，世界！”不唱"99瓶啤酒“这项任务是在1D中完成的。你可以把它解释为用某种力量绑在两端墙壁上的橡胶长袍。
在[0,1]域中，为给定的源函数f和边界值u_L和u_R查找函数u，以便：
-u'' = f
u(0) = u_L
u(1) = u_R
u''表示u的二阶导数
这完全可以从理论上解决，但您的任务是在离散域x上对N点进行数值求解：
x = \{\frac i {N-1} : 0
浏览 0提问于2016-11-06得票数 15
回答已采纳
1回答
我已经得到了基本的SIR模型以下代码。我需要在没有ode()函数或任何其他函数的情况下求解微分方程。你能帮我在不使用ode()函数的情况下计算这个微分方程吗？
library(deSolve)
library(shiny)
ui <- fluidPage(
titlePanel("Choose amounts"),
fluidRow(column(width = 4,
sliderInput(inputId = "time_values", label = "Days", value = 10, min = 1, m
浏览 2提问于2020-02-18得票数 0
1回答
例如，矩形框的体积可以计算为V=L*W*H 假设我们知道V，L，W，那么我们可以求解H。 假设我们知道V，H，W，那么我们就可以求解V。 假设我们知道L，W，H，那么我们可以求解V.和e.t.c 有没有办法在python中解决(我现在正在尝试Sympy )来根据给定的输入来解决它？当然我可以使用if的情况，但我需要写4个方程来求解，这是一个简短的方程。 任何建议都是值得感谢的。 致以亲切的问候，伊万
浏览 4提问于2020-01-12得票数 1
回答已采纳
1回答
鲤鱼C++ LU分解
、、、、
我正在使用Armadillo C++库来求解中、大尺寸线性系统(1000-5000方程)。
因为我要解不同的线性系统
AX=b
在A总是一样，b总是变化的情况下，我只想对A进行LU分解一次，并在不同的b中重用LU分解。不幸的是，我不知道如何在Armadillo中执行这种操作。
我所做的只是A矩阵的LU分解：
arma::mat A;
// ... fill the A matrix ...
arma::mat P,L,U;
arma::lu(L, U, P, A);
但是现在我想用矩阵P，L和U来解几个不同b向量的线性系统。
你能帮帮我吗？
浏览 1提问于2013-07-30得票数 2
1回答
Java中的微分方程
、、、
我正在尝试用java创建一个简单的SIR-流行病模型的模拟程序。
基本上，SIR是由三个微分方程系统定义的：
S'(t) =- l(t) * S(t)
I'(t) = l(t) * S(t) - g(t) * I(t)
R'(t) = g(t) * I(t)
S易感人群，I感染人群，R恢复人群。
l(t) =c*x* I(t) / N(T)
C-接触数，x-传染性(接触病人后发病的概率)，N(t) -总人口(不变)。
我如何在Java中求解这样的微分方程？我不认为我知道任何有用的方法来做到这一点，所以我的实现产生了垃圾。
public class Main {
publ
浏览 0提问于2010-12-05得票数 11
回答已采纳
2回答
参数方程中的直线相交
、、、
我相信这个问题以前也有人问过。然而，我试图通过把纸上的方程转换成一个实际的函数来连接这些点。我想在这里问一问数学网站会很有趣(因为它无论如何都会被用于游戏)。
假设我们有向量方程：
x = s + Lr;
其中x是得到的向量，S是我们的起点/向量。L我们的参数和我们的方向向量。
(不确定它的名称是这样的，请纠正我)正态方程是：
x.n = c;
如果我们用向量方程来代替，我们得到：
(s+Lr).n = c.
我们现在需要孤立L，这导致了
L = (c - s.n) / (r.n);
L需要0
我的问题是:我想知道L是什么，所以如果我用L代替两个向量方程(或两条线)，他们应该给我相同的交点坐
浏览 0提问于2012-11-30得票数 7
回答已采纳
下面是矩阵方程：
[M] = [A][R] + [B][L]
其中：A和B =2×2矩阵
M，R，L和=2×1矩阵。
为了估计A和B的一个特定值，我有9种不同的[M]、[R]和[L]数据。假设这9个数据是一个组，我有数百个组数据。
我需要找到与特定组数据相关的A和B的值。
那么，有人是否对matlab有或知道如何通过简单地将[A]、[R]和[L]数据输入到matlab中来获取每个组的[M]和[L]值？
浏览 5提问于2009-12-10得票数 0
5回答
数值微分方程求解器算法意外分段
、、、、
我试着用八度解微分方程，但用我选择的微分单位永远用不了多久，所以我决定用C语言编码。
#include <stdio.h>
double J = 5.78e-5; // (N.m)/(rad/s^2)
double bo = 6.75e-4; // (N.m)/(rad/s)
double ko = 5.95e-4; // (N.m)/rad
double Ka = 1.45e-3; // (N.m)/A
double Kb = 1.69e-3; // V/(rad/s)
double L = 0.311e-3; // mH
double R = 150; // ohms
dou
浏览 2提问于2013-05-17得票数 1
回答已采纳
1回答
FFT系数线性系统的求解
、、、、
如果我想在(2d)空间中用FT积分周期域中的热方程，我会用来写：
(1-kappa*dt*L/2)ftT(now)=(1+kappa*dt*L/2)ftT(last)
其中L是一个数列，它将解ftT缩放到它的Laplacian。如果ftT阵列是方形的，就像在Matlab/Octave中那样，用LAPACK很容易解决这个问题。然而，在r2c、c2r模式下使用FFTW，复杂的转换数组并不是正方形的。很容易看到如何修改L以处理这部分非正方形，并且很清楚如何填充ftT数组以使其成为正方形，但这似乎是浪费的。我还没有想出一种解决非平方线性系统的方法。对于这样一个系统，方程的数量必须等于变量的数量，我对此
浏览 2提问于2017-06-23得票数 1
回答已采纳
目前，我使用具有已知/未知参数的不同组合的单个方程。因为我没有任何花哨的计算器，所以在Mathematica中定义方程并传递已知参数来计算未知值会容易得多。
如果你们中的任何人能给出一个示例解决方案(如果可能的话，使用给定的方程)，我将非常感谢。
假设我们有一个椭圆轨道上给定点的卫星速度方程：
V= sqrt(u(2/r - 1/a))
哪里
V=速度u=常数3.986 * 10^14 m^3/s^2 r=半径(距地球中心的距离)a=椭圆的半长轴
这个方程可以用来计算速度，例如，我们知道一次操纵将货物移动到其他轨道所需的速度，并必须模拟给定半径(r)的轨道(a)。
谢谢!
浏览 0提问于2009-09-13得票数 3
回答已采纳
2回答
我想用MATLAB消除方程中的变量。例如，让我们考虑以下方程：
p = (m + n)
q = (m - n)
r = (m^3 - n^3)
现在，r可以用p和q来表示，方法是完全消除m和n，比如r = (3*p^2*q + q^3)/4。
可以使用以下方法在Mathematica中实现这一点：
Eliminate[{p == (m + n), q == (m - n), r == (m^3 - n^3)}, {m, n}]
如果有可能的话，如何在MATLAB中得到同样的结果。为此，在不同的应用程序之间切换是非常不方便的。
浏览 5提问于2018-07-13得票数 5
回答已采纳
我需要解一个线性方程(Ax = b)，它有一个很大的稀疏的A和多次不同的b。因此，LU (如果A是对称的，则是Cholesky )，因式分解是非常可取的。
我正在使用armadillo库，我听说可以使用下面的函数：
spsolve(x, A, b, "superlu");
为了解决这样一个系统。我不太关心检索L和U矩阵。然而，最重要的是，每次我调用spsolve时，L和U都不会被重新计算。
spsolve(x, A, b, "superlu")是否存储LU分解?如果不是，是否有方法检索所述矩阵？
浏览 0提问于2018-10-20得票数 1
我正在尝试解决一个非线性优化问题，其中目标函数是非线性的，约束是线性的。我在R中读了一点ROI包，我决定使用同样的包。然而，我在解决优化问题时遇到了一个问题。
我本质上是在试图最小化供需曲线下的面积。供应曲线和需求曲线的方程在代码中定义：
目标函数:最小化(供给曲线积分+需求曲线积分)，
受到大于或等于34155 (存储在称为ICR的变量中)的约束Q的约束，
Q大于或等于0
并且Q小于或等于40000。
我曾尝试在RStudio中运行ROI包，但一直收到错误消息，告诉我找不到解算器。
library(tidyverse)
library(ROI)
library
浏览 4提问于2018-12-21得票数 0
所以，我有一个压力损失的管道模型，其中未知的是质量流量。一般情况下，在该问题的大多数模型中，都采用守恒方程来计算质量流量，但这类模型存在大量的收敛问题(因为管端的阻塞流导致了末端的无限压力导数)。关于左边和右边问题的表示，见下图，这是显示无限压力导数的图。
正因为如此，我使用了一个更稳健的模型，虽然它输出的不是质量流量，而是管道长度，这是已知的。因此，需要一个迭代循环来确定质量流量。然后，我编写了一个函数length，给定管的几何形状、质量流量和边界条件，它输出计算出的管长，并使方程如下所示：
parameter Real L;
Real m_flow;
...
eq
浏览 0提问于2018-07-10得票数 2
回答已采纳
1回答
我需要解决一个非线性方程系统，如下所示：
def trySolveEquation(V, L):
#The equation to solve is:
#{ (V . Ct) ^ 2 = 1
#{ (L + u) . Ct = 0
#C and u are the unknowns, C is a vector and u is a scalar, Ct is a vector transposed from C.
#V is a vector with dimension equal to Ct, L is a square matri
浏览 1提问于2021-05-06得票数 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
2回答
我正在寻找一个第三方自定义控件，可以处理编辑和解决简单的数学方程(与图像生成)。
我了解LaTeX并使用它的语法来生成数学方程式的图像(我在这里使用了这个实现LaTeX：的项目)，并且我找到了几个可以将方程式作为输入字符串、编译、解析和求解方程式的项目(这里有一个：)。
但是，我还没有找到任何方法，允许您使用输出的预览图像构建方程(就像LaTeX一样)，然后允许您为方程中的变量指定值，并使用输出的结果求解方程。
有没有人知道是否存在可以编辑和求解数学方程的第三方控件？我使用的是C# .NET 4.0，所以希望我可以使用.NET中的库/控件。
任何帮助或链接到类似的东西将非常感谢。
谢谢，马克
浏览 2提问于2011-04-13得票数 9
我正在尝试用MatLab求解一些常微分方程，因为方程中的变量是总体，它们需要被约束为正数。因此，在调用方程求解器之前，我尝试使用odeset()将它们设为非负数，但在之后绘制这些值时，它们有时实际上是负数(在下面的代码中，它是洋红色的线)。我做错了什么？
下面是一些代码：
%Lots of variables
includeJ=1;
cullLIRate=1/2000;
cullDIRate=1/2000;
N = 16800;
beta = 2e-7;
delta = 0.5;
gamma = 1/50;
sigma = 1/400;
mu = 1/365;
maxTime = 30*365
浏览 0提问于2013-02-25得票数 0
回答已采纳
2回答
我正在用C语言编写一个光线跟踪器，为了绘制一个球体，我使用了笛卡尔方程：
x^2 + y^2 + z^2 = R^2.
我有我的眼睛位置(x_eye，y_eye，z_eye)和我的眼向量(Vx，Vy，Vz)。我这条线的参数方程是：
x = x_eye + k * Vx
y = y_eye + k * Vy
z = z_eye + k * Vz
我把我的直线的参数方程放在球的笛卡尔方程中，以便求解它。
(x_eye + k * Vx)^2 + (y_eye + k * Vy)^2 + (z_eye + k * Vz)^2 = R^2
(Vx^2 + Vy^2 + Vz^2) * k
浏览 2提问于2011-03-12得票数 3
回答已采纳
1回答
我试图在数值上找到一个问题的解，即求圆的中心和半径，它与给定的直线相切，并通过另一个给定的点。问题如下图所示：
圆必须与直线Y=mx +c切线，并且必须通过点。由于问题是对称的，圆的中心总是位于x轴(y=0)上，x轴是这两条线的平分线。A点一直是已知的，角θ也是如此。点P指定为输入
我可以写下直角三角形的一些方程：Sin theta = BC / AC和我知道BC == PC == R，但是我不知道如何用数值来解决这个问题，因为每次调用代码时，直线和点的方程都会不同。用手求解将需要将A和P替换回我认为的线的方程中。
我想我需要用一些scipy模块来解一些联立方程，但是我不知道我到底该
浏览 15提问于2022-10-04得票数 0
回答已采纳
1回答
这是我的载体。注意，向量m由未知值x组成。
r=c(1, 3, 5)
p=c(5, 3, 1)
s=c(1, 3, 5)
m=c(x, x, x)
如何在R中解出以下条件向量方程？我想找出x的值，它使p*(s-m)*(ifelse((s-m)<0, r, 6-r))向量中的元素之和为零。
sum(p*(s-m)*(ifelse((s-m)<0, r, 6-r)))=0
我被告知rSymPy可能会处理这个问题，但我不认为这适用于向量方程。有什么想法或建议吗？
事实上，我能够在Excel中使用“目标寻求”得到答案，但我想得到解决这个问题的R命令。
浏览 1提问于2016-04-22得票数 1
回答已采纳
1回答
机器学习中损失函数的图形来自哪里？我正在学习关于机器学习的知识。我有时不理解使用正则化项进行优化的模型。在正则化的解释中，可能会出现如下图。
下面是L1正则化项的一个示例。我假设模型有两个权重参数w1，w2。也就是说，模型y的方程由下面的方程表示。
y = w1x1 + w2x2
为了简单起见，我忽略了偏置项。
红色方块表示正则化项。蓝色椭圆表示不含正则项的损失函数。正则化项由下式给出
| w1 | ^ q + | w2 | ^ q = r ^ q (r is const.)
因此，在w1> 0和w2> 0处的图的方程表示如下。
w2 = (r ^ q-| w1 | ^ q) ^
浏览 24提问于2020-02-26得票数 1
1回答
R二次规划
、、、
我有一个问题想在R中解决。
我看到，我可以使用包lsei中的函数limSolve来最小化以矩阵形式写成的线性方程组( Ax=b )，但条件是等式约束Ex=f和不等式约束Gx>=h。
然而，不是线性方程组，我现在有一个二次方程组，它可以写成t(x)Ax=b。
我看到了二次型情况下的包quadprog，但是它似乎不允许一组二次方程，只有一个方程。
有人知道我能用什么来最小化等式和不等式约束下的二次方程组吗？
这是我的例子。我试图结合三个概率P(A)，P(B)，P(C) -这就创建了7个片段v1到v7，其中v1是P(A solus)等等。v4是P((A和B)不是C)等。v7为P(A、B、C)。
浏览 2提问于2016-01-13得票数 0
我正在研究如何使用Nlopt来解决Python中的优化问题。
我有一系列形式为Ax =b的联立方程
其中A是一个NxM矩阵，有x的解。另一种思考方法是，我有N个形式的联立方程，形式是x_1c_1m + x_2c_2m +……+ x_Nc_Nm = k_M，其中x_i是要求解的变量，c_im是方程M=m中与x_i相关的常数，k_m是方程M=m中的常数，c_im和k_m都是已知的。
让我困惑的是，在Nlopt如何处理这个问题。Nlopt要求您有实际的可调用函数，而我没有吗？我想我可以把上面那个矩阵方程中的每一个方程推广到这样的情况：
def fn(x,c_m,k_m):
val = 0
浏览 25提问于2022-01-26得票数 0
我知道如何用numpy来求解基本的线性矩阵方程。
但是，我有一个矩阵A，方程A^2 + xA + yI = 0，其中x和y不是向量，而是标量。I是单位矩阵，0是维数匹配A的零矩阵。
对于小型矩阵来说，这是一个非常简单的方法(当然，假设存在一个解决方案)，但是我正在练习一个编码面试，并希望用python来解决类似的问题。也许给出的矩阵很大..。
下面是一个样本矩阵A，它得到解x=-2，y=1：
np.array([[1,1,0],
[0,1,0]
[0,0,1]]
在纸上，这与解线性方程组x= -2和x+y=-1一样容易。我所面临的问题是将上述形式的方程解
浏览 4提问于2020-02-22得票数 1
回答已采纳
1回答
以下是计算结果
L1 =线连接点A(x1，y1)和B(x2，y2) L2 =线连接点c(x3，y3)和D(x4，y4)
For Line L1
直线方程:y= m1*x + c1
斜率m1：(y2-y1)/(x2-x1)
Y截距: c1 = (y1 - m1*x1)
For Line L2
直线方程:y= m2*x + c2
斜坡m2：(y4-y3)/(x4-x3)
Y截距: c2 = (y3 - m2*x3)
求交点
解上述方程我们得到
X= (c2 -c1)/(m1-m2)
Y= (c1*m2 - c2*m1)/(m2-m1)
这些计算用于计算我的java程序中的交点。
问题
当这两条线分
浏览 0提问于2014-11-06得票数 0
回答已采纳
我正在尝试解决R中的一个数值方程，但我希望有一种方法，它的性能类似于Matlab中的vpasolve。我有一个非线性方程(涉及许多对数函数)，当用单根在R中求解时，与matlab中的vpasolve相比，它给出了完全不同的答案。
浏览 0提问于2020-09-01得票数 0
1回答
我试图用numpy来求解如下形式的非线性方程组：
a(y-2.7)(1-exp(-a*z)) = (x-2.7)(1-exp(-z))
b(w-2.7)(1-exp(-b*z)) = (x-2.7)(1-exp(-z))
c(w-2.7)(1-exp(-b*z)) = (y-2.7)(1-exp(-a*z)
d([y+w]/2-2.7)(1-exp(-d*z)) = (x-2.7)(1-exp(-z))
很明显，系统中的方程和未知数一样多。a，b，c，d是上述系统的常量。这是最简单的系统，在其他情况下会有更多的方程。这些方程的解具有相似的数量级，因此我知道Levenberg-Marquar
浏览 4提问于2014-06-09得票数 0
回答已采纳
在求解如下形式的非线性抛物型偏微分方程时，我遇到了一个问题。(将u_x视为u w.r.t.的区别。x)
U(1)_t = a1(U(1)，x)*U(2)_t + b1(U(1)，x，U(3)) =( D1(x，U(1))*U(1)_x )_x + c1( U(1)，U(3))*U(1)_x -(1)
和它的耦合方程
U(3)_t = a2(U(1)，U(3))*U(2)_t + b2(U(1)，x，U(3)) =( D3*U(3)_x )_x -(2)
如你所见，U(1)，U(2)，U(3)只有两个偏微分方程。这是因为U(2)不是x的函数，它有自己的pde作为
U(2)_t = (1/r)(
浏览 3提问于2012-10-23得票数 1
1回答
解复矩阵方程
、、、
我有一个涉及矩阵的复杂方程式：
R = expm(X)*A + (expm(X)-I)*inv(X)*B*U;
其中R、B和U是已知矩阵。
I是一个单位矩阵。
我需要解决X的问题。有没有办法在MATLAB中解决这个问题？
浏览 0提问于2019-04-28得票数 1
1回答
我在下面的代码中定义了两个方程f1(b，bb) =0和f2(b，bb) =0。我试图找到满足这两个条件的(b，bb)集。 import sympy as sp
from sympy import symbols, simplify, factor
a, b, aa, bb, q, l, h = symbols('a b aa bb q l h')
pb = l*q+h*(1-q)
pb0 = (l*(1-l)*q+h*(1-h)*(1-q))/((1-l)*q+(1-h)*(1-q))
pb1 = (l**2*q+h**2*(1-q))/(l*q+h*(1-q))
a =
浏览 20提问于2019-09-07得票数 2
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
实时音视频
即时通信 IM活动推荐
腾讯云开发者
