-
-1 votesanswersviews
如何评估或绘制SYM对象？
我正在使用一个脚本来解决matlab中的拉格朗日插值问题 . 我是新手，我不知道它是如何工作的 . 插值数据是下一个 X = [0 1 2 4 8] Y = [1 5 10 24 50] 然后我以这种方式构建多项式od lagrange： syms x l1 = ( (x-X(2)) * (x-X(3)) * (x-X(4)) * (x-X(5)) ) / ( (X(1)-X(2)) * (X... -
1 votesanswersviews
MATLAB中的符号数学，解决了简单的集成问题
我有一个问题，通过MATLAB解决一个简单的集成 . 我想解决这个象征性问题，并且通过其他程序执行此操作时没有任何问题 . 好吧，我有这个等式： syms k x 对于2 <= x <= 6，fX（x）= k * e ^（ - 3 * x） 我希望从区间2到6进行积分 . 然后我将求解方程，使得fX（x）= 1，并求解k的方程 . 我键入： S = solve（'int（k * ex... -
1 votesanswersviews
MATLAB：symfun的符号积分
我希望随着时间的推移进行符号整合;代码如下 . syms x1 u1 t x1 = symfun(sym('x1(t)'), [t]); x1dot = p1 + p4*p8 - p13*x1; int(x1dot,t) 答案应该是： e^(-p13*t)*x1(0)+(p1 + p4*p8)/(-p13)*[1-e^(-p13*t)] 我得到的是： Warning: Explicit int... -
0 votesanswersviews
MATLAB中学生密度的符号整合
我想在MATLAB中集成学生的t密度，但MATLAB似乎失败了，只返回相同的表达式而不进行计算 . 关于如何进行的任何想法 . 我知道它应该评估为1，但是迫使MATLAB去做的诀窍是什么 . syms x s mu nu g = @(x) gamma(1/2*(nu+1))/(gamma(nu/2)*sqrt(pi*nu*s^2))*(1 + 1/nu*(x-mu)^2/s^2)^(-1/2... -
2 votesanswersviews
Matlab中的双重符号积分
我正在尝试计算一些积分，例如： a = 1/sqrt(2); b = -5; c = 62; d = 1; f = exp(-x^2-y^2)*(erfc((sym(a) + 1/(x^2+y^2)*(sym(b)*x+sym(d)*y))*sqrt((x^2+y^2)*sym(10.^(c/10))))... + erfc((sym(a) - 1/(x^2+y... -
0 votesanswersviews
Matlab中使用梯形规则的符号积分
我想用梯形规则象征性地整合一个函数 . 问题是我不想最终表达集成的样子 . 我想要所有的步骤 . 例如，对于一个函数 function y = myfun(x, y) y = 2*x+y.^2; 从t = 0到t = 2进行积分，使用梯形法则得到t = 0.5的时间步长，我想要每个时间步长的求和表达式： t = 0.0 to t = 0.5 t = 0.5 to t = 1.0 t = ... -
2 votesanswersviews
如何在Matlab中象征性地获得复杂表达式的虚构和实部？
我试图获得以下复杂表达的虚构和真实部分 . syms a b F = imag((cos(a)-j*sin(a))/(1+j*a*b-cos(a)+j*sin(a))) simplify(F) 无论是否简化，输出都相同 . 有没有办法用标准的复数格式（x jy）表达这个表达式？ -
3 votesanswersviews
使用sympy和Mathematica对Modified Bessel函数进行符号积分
我想用 sympy 解决一个我知道Mathematica可以解决的定积分 . 在Mathematica中，以下行 Integrate[z^2 (BesselI[0, z^2] - BesselI[1, z^2]) Exp[-z^2], {z, 0, x}] 产量 1/3 x^3 HypergeometricPFQ[{1/2,3/2},{1,5/2},-2 x^2]-1/10 x^5 Hyperg... -
5 votesanswersviews
Sympy“全球”替代
我有很多符号表达式，我可能会发现其中一个系数为零 . 我想，也许是因为我习惯了mathematica，以下是有道理的： from sympy import Symbol x = Symbol('x') y = Symbol('y') f = x + y x = 0 f 令人惊讶的是，返回的是 x + y . 除了在每个等式上明确调用"subs"之外，有没有办法让 f 只返回... -
4 votesanswersviews
如何在同情中采用矩阵的元素对数？
使用同情矩阵或numpy数组的同情符号，如何采用元素对数？ 例如，如果我有： m=sympy.Matrix(sympy.symbols('a b c d')) 然后 np.abs(m) 工作正常，但 np.log(m) 不起作用（"AttributeError: log"） . 有解决方案吗 -
3 votesanswersviews
SymPy无法计算此矩阵的特征值
我想计算拉普拉斯矩阵的第二个特征值，以检查相应的图形是否连接，但是当我尝试使用SymPy的 eigenvals 时，很多时候会发生抛出错误 MatrixError: Could not compute eigenvalues for Matrix([[1.00000000000000, 0.0, 0.0, 0.0, -1.00000000000000, 0.0, 0.0, 0.0, 0.0, 0... -
0 votesanswersviews
matlab与真实的符号整合给出了一个复杂的答案
我希望得到Matlab中以下积分的解析（封闭形式）解决方案 . 但是，Matlab给出了一个实部和虚部的答案 . 如何通过“真实”部分来获得答案 . 这是完整的代码 . close all; clear all; clc; syms t real; syms thetak real; syms sik real; syms tbar real; syms sjk real; expr = ex... -
3 votesanswersviews
如何在Simulink Matlab功能块中使用syms
我想解决 a 和 b+tau 之间 expm(A*s) 的积分，其中 tau 是时变延迟 . 我在Simulink中创建了一个带有 tau 作为输入的Matlab功能块，如下所示： function y = compute_int(u, tau) syms s gamma=double(int(expm(A*s),s,a,b+tau)); B = [gamma; 1] y = B*u; 之前定... -
6 votesanswersviews
使用matlab求和系列
当我在matlab中写这个 syms x; f=x^3-cos(x); g=diff(f) 它放弃了 g = 3 * x ^ 2 sin（x） 现在我想生成求和系列 我google并发现“symsum”命令，但是当我编写以下命令时，它不执行我所需的任务 syms k symsum(k^2, 0, 10) symsum(1/k^2,1,Inf) 它给出了输出 ans = 385 ans ... -
1 votesanswersviews
用符号表达式求解线性方程组
嗨，我试图用mathematica解决线性方程组 . 我有18个方程和18个未知数，系数矩阵具有满秩 . 所有条目都是象征性的，因为我试图通过分析解决问题 . 不幸的是，Mathematica永远不会停止评估 . 我准备了一个最小的工作示例： n = 18 A = Table[AA[i, j], {i, 1, n}, {j, 1, n}]; A // MatrixForm x = Table[xx... -
3 votesanswersviews
使用Matlab中的标准绘图函数绘制符号方程
为了获得流体行为的图形表示，通常的做法是绘制其流线 . 对于具有速度分量u = Kx和v = -Ky的给定二维流体（其中K是常数，例如：K = 5），可以获得流线方程组件的流线方程，如下： 流线方程：∫dx/ u =∫dy/ v 求解的等式如下所示：A = B C（其中A是第一积分的解，B是第二积分的解，C是积分常数） . 一旦我们实现了这一点，我们就可以通过简单地为C赋值来开始绘制流线，例如：C... -
1 votesanswersviews
如何解决极大的方程组
我需要解决一个庞大的方程组 . 我现在需要解决的系统有30个方程和30个未知数，但我需要能够使用同一个系统来解决1,000个方程和1,000个未知数 . 我不需要知道所有未知变量的答案，我只需要知道其中一些变量的答案（30个方程式中的2个，以及1,000个方程式中的10个） . 我需要解决的方程式各自为： (c0 * x0) + (c1 * x1) + ... + (c_n * x_n) 其中 ... -
0 votesanswersviews
如何在Simulink Matlab函数中使用符号变量和函数（Syms）？
我想在Simulink Matlab函数中创建一个符号函数来解决变量h和t1 . Matlab产生错误“独立代码生成不支持函数'syms' . 请参阅coder.extrinsic文档，了解如何在模拟中使用此函数 . ”当我尝试使用以下代码编译Simulink Matlab函数时 . syms Eq1(h,t1); Eq1(h,t1) = h*t1; 我尝试在顶部添加“coder.extrin... -
3 votesanswersviews
JModelica结合CasADi在Python中获得线性化模型
我最近发现了一篇论文“制作可用于PYTHON控制系统库分析的MODELICA模型”，解释了如何将JModelica与CasADi结合使用以获得Python中的线性化模型 . 这正是我正在寻找的，但显然该论文是从2014年开始使用JModelica 1.12 . 不推荐使用所描述的方法，并且不再使用已使用的功能 . 关于此主题，JModelica的文档非常有限，我不知道如何使用最新版本的JModel... -
2 votesanswersviews
MATLAB的符号数学无法解决简单的不等式
我有一个简单的不等式，MATLAB的符号数学工具箱正在做一些非常奇怪的事情 . 以下是变量： >> syms X ndot4B xiA ndot4A xiB 我正试图解决以下不平等问题（请在“纸上”自己动手）： >> solve(X*ndot4B*xiA - ndot4B*xiA + X*ndot4A*xiB > 0, xiA) 答案是： ans = (X*nd... -
2 votesanswersviews
在MATLAB / Maple /等上有效地解决符号线性方程组
基本上我有一个线性方程组，其中每个系数都是一个变量 . 例如，对于2x2系统： a11 * x1 a12 * x2 = f1 a21 * x1 a22 * x2 = f2 x1，x2 = ....（带有a11，a12，a21，a22，f1，f2的符号表达式） 我尝试通过Cramer的规则在MATLAB中解决它，以符号形式编写系统的矩阵并计算适当的决定因素 . 问题是，所需的时间随着系统... -
1 votesanswersviews
在Maple中进行高效编程
我需要一些关于在Maple中有效编程的技巧（如在计算机代数软件中） . 通过高效，我并不是说我的程序/代码更高效（性能），而是有效地使用Maple代码 . 这是具有讽刺意味的，但是我为用C / C等语言编程开发的软件和工具而烦恼，例如： ：IDE，debuging工具（ Watch ，断点等），编译器和语言结构本身（面向对象，程序更灵活，重载，通过引用或值传递等） . 从这个意义上来说，我觉得在C... -
0 votesanswersviews
MATLAB：在symsum中使用匿名函数，然后使用fsolve
在MATLAB中，我希望定义一个匿名函数，其中包含一个明确的和，以及另一个匿名函数 . 这是一个MWE，希望能够描述我想要做的事情： clear; n=1; syms j; a=0; b=sqrt(0.5); Finv = @(x) logninv(x,a,b); fun = @(x) 0.5-symsum(Finv(j*x), j, 1, n+1); fsolve(fun,0.1) 返回的错误... -
7 votesanswersviews
如何在SymPy中扩展矩阵表达式？
在SymPy中，我试图执行矩阵乘法并在之后展开它 . 但是，SymPy似乎不支持矩阵表达式的扩展 . 例如，这是矩阵的4阶Runge-Kutta（RK4）： from sympy import init_session init_session() from sympy import * A = MatrixSymbol('A', 3, 3) x = MatrixSymbol('x', 3, 1... -
3 votesanswersviews
动态形状和带Tensorflow的索引
我正在为Tensorflow改编this MLP tutorial . 但是，我无法弄清楚如何正确有效地计算负对数似然 . 当我尝试构建以下图表时，会发生错误： x = tf.placeholder(tf.float32, shape=(None, 784), name='x') y = tf.placeholder(tf.int32, shape=(None,), name='y') # [.... -
0 votesanswersviews
Matlab中隐式符号表达的梯度
我的情况如下：我有一个象征性的表达式： syms X Y Z K Ra=51.7; P=[0 0 200]; Sa=sym('Ra^2==(Z-P(3))^2+(Y-P(2))^2') 其中 Y 和 Z 被定义为符号 . Ra 和 P 是向量 . 我需要获得 Sa 的渐变但是我收到一个错误： G=gradient(Sa,[Y Z]); Error using symengine (lin... -
8 votesanswersviews
sympy：'Transpose'对象没有属性tolist
我正试着用sympy做一些符号矩阵计算 . 我的目标是获得一些矩阵计算结果的符号表示 . 我遇到了一些问题，我已经归结为这个简单的例子，其中我尝试评估指数矩阵的取幂结果并将其乘以任意向量 . >>> import sympy >>> v = sympy.MatrixSymbol('v', 2, 1) >>> Z = sympy.zeros(2,... -
2 votesanswersviews
sympy - 象征性数量的元素的象征性和
在'sympy'中表达以下内容的最恰当方式是什么： 样本'x [i]'与“i”从混凝土0到符号'N'的总和 . 'x [i]'本身应具有象征性，即始终显示为变量 . 目标是在线性方程组中使用这些表达式 . 示例（平凡的最小二乘法）： 给定一组样本（x [i]，y [i]），它们应该位于由'y = m * x a'给出的行上 . 也就是说，估计的行由'm'和'a'确定 . 样本和估计线之间的误差可... -
1 votesanswersviews
Sympy中带符号的布尔运算
符号上布尔变量的布尔运算产生 TypeError ，但反过来没有问题： >>> from sympy import * >>> x = Symbol('x', bool=True) >>> x ^ True Not(x) >>> True ^ x Traceback (most recent call last): Fi... -
2 votesanswersviews
计算机代数系统支持可变大小的矩阵
我熟悉sympy，matlab符号包，减少，并尝试过其他一些计算机代数系统 . 但是，据我所知，它们似乎都不能在可变大小的矩阵上进行代数 - 它们只能用于固定大小的矩阵 . 是否可以为可变大小的矩阵做代数？我知道会有相当多的重大案例，但我觉得有很大一部分是可行的，因为很容易进行许多简化/代数，可以用R ^ nxn中的矩阵手工完成 . 可以在许多这些中使用非共同代数元素，因此涵盖了矩阵的加法和Had...