-
4 votesanswersviews
使用BLAS或LAPACK与Xcode
我正在寻找一个很好的指南，介绍如何将BLAS或LAPACK函数合并到我通过Xcode开发的Objective C程序中 . 我可以在线找到BLAS / LAPACK中程序的唯一来源是用Fortran编写的 - 这很令人困惑 . 如何在Objective C中使用BLAS / LAPACK和代码？有没有在线示例/指南？我希望能够执行矩阵运算，其中矩阵的元素具有复杂的值 . 干杯， Shaddy -
2 votesanswersviews
Sympy复杂矩阵求逆/运算
使用来自Continuum Anaconda发行版的Sympy 0.7.2，（懒得构建我自己的堆栈，不讨厌），在做复杂矩阵代数时我一直遇到这个问题，其中一些原生函数取代了典型的'1j'使用'I' . 这很烦人，因为它不能很好地与其他算术运算一起使用，并且基本上保持其形式并增加因子而不是计算浓缩结果 . 我将给出一个例子： >>> T [ 1, 1.0 - 1.0... -
0 votesanswersviews
如何绘制复杂传递函数的根轨迹？
我在尝试获取以下复杂传递函数（H）的根轨迹图时遇到问题： % Ax = b Form A = 1.0e+06 *[0.0713 - 0.0073i,... -1.8793 - 0.0357i]; b = -32.4298; %TF G = A\b; %PI controller gain Kp = 100; %proportional gain ... -
6 votesanswersviews
从python中的非线性方程组中寻找复杂的根
我一直在测试一种已经在文献中发表的算法，该算法涉及在Matlab和Python中求解一组“m”非线性方程 . 这组非线性方程涉及包含复数的输入变量，因此得到的解也应该是复杂的 . 截至目前，我已经能够通过使用以下代码行在Matlab中获得相当不错的结果： lambdas0 = ones(1,m)*1e-5; options = optimset('Algorithm','levenberg-mar... -
-2 votesanswersviews
使用Python查找Laurent系列函数
我已经找到了一个名为SymPy的库，用于Python中的符号计算，但问题是我不知道如何在程序中生成Laurent系列 . 当然，我总是使用泰勒系列以临时方式计算Laurent系列，从未使用算法方法 . 如果有人帮助我找到一个使用下面问题中提到的输入生成Laurent系列的算法，我将不胜感激:-) Given a fractional function containing polynomials ... -
2 votesanswersviews
Tensorflow - 从complex64转换为2x float32
我正在尝试使用Tensorflow LSTM RNN进行一些音频处理 . 我正在使用tf.contrib.signal.stft来帮助神经网络更容易理解我的数据，但它返回了类型为complex64的Tensor . 如果我尝试将其提供给dynamic_rnn，我会收到以下错误： ValueError：需要变量rnn / basic_lstm_cell / kernel的初始值设定项 所以我需要... -
1 votesanswersviews
Python 3.6：复数的分层
我需要对一个复数进行分区，但似乎在Python 3.x中删除了复数的分区 . 这有什么变通方法吗？任何想法，将不胜感激 . 例如 . (2 + 3i) // (complex anything) will throw an error. -
4 votesanswersviews
比较python和matlab中的fsolve结果
我对几天前写的帖子有一个跟进问题，谢谢你以前的反馈： Finding complex roots from set of non-linear equations in python 我现在已经在python中设置了设置的非线性方程，因此fsolve将独立处理实部和虚部 . 但是，仍然存在python“fsolve”收敛到正确解决方案的问题 . 我有完全相同的输入在Matlab中使用，经过双重检查... -
0 votesanswersviews
具有复杂方程组的lsqnonlin
我有两个复杂的非线性方程，有两个复变量（然后是四个实变量） . 我试图使用lsqnonlin而不是fsolve来定义我的未知变量的边界（lb，ub） . options=optimset('Display','iter','MaxFunEvals',1e6,'TolX',1e-16); fun = @myfun; x0 = [0,0,1000/2,100/2]; lb = [0,0,0,0]; u... -
0 votesanswersviews
如何检查复数是否为整数
给定一个值数值 x ，你可以这样做 float(x).is_integer() 来检查它是否是一个整数 . 对复杂的数字有没有办法做到这一点？ 我正在尝试使用列表推导来仅获取有限域上的多项式的整数根，这些整数是整数 . [r for r in solve(f,domain=FiniteField(p)) if float(r).is_integer()] 但如果solve函数返回复杂的根，则不起... -
1 votesanswersviews
在C和Fortran代码之间传递复数数组
我正在使用C和Fortran代码 . C代码需要调用Fortran子例程 . Fortran子例程的一个参数具有Fortran类型 complex . C代码在一个文件中，Fortran代码子例程在另一个文件中 . 我在64位 GNU/Linux 系统上使用 gcc 和 gfortran . 这是一个显示Fortran子例程声明的片段（以及一些额外的行）： SUBROUTINE minp... -
2 votesanswersviews
最小二乘拟合python中具有错误的复数
这里抱歉是代码 . 我试图将贝塞尔函数拟合到某些2D数据，贝塞尔函数很复杂 . 但是，当我适应数据时，我总是最终得到一个数组错误，我无法排序（我是python的新手） . 最后，代码理想情况下会给出错误来解释导入数据的分散 . from numpy import * from pylab import * from scitools.std import * import csv import m... -
1 votesanswersviews
Gfortran复杂的实际到真正的虚拟参数
我试图使用fftpack与gfortran，但我得到的错误，我认为与虚拟参数被声明为真实时传递复杂数组的一些例程有关 . 我在intel fortran页面上阅读了一条评论，可以禁用“ check routine interface ” . 有谁知道gfortran是否有类似的选择？ 我想不必编辑fftpack ...（我想这是因为内存中的复杂由两个实数表示，数组参数作为引用传递但如果我错了请纠正... -
2 votesanswersviews
tensorflow CNN具有复杂的功能和标签？
我最近发现了一篇论文，他们使用带有复杂2D特征图的CNN作为输入 . 但是，网络也输出复杂的输出向量 . 他们使用了具有张量流后端的Keras . 这是链接：https://arxiv.org/pdf/1802.04479.pdf 我问自己是否有可能构建复杂的深度神经网络，如具有张量流的CNN . 据我所知，这是不可能的 . 我错过了什么吗？ 还有其他相关问题解决了相同的问题，没有答案：Compl... -
0 votesanswersviews
使用sympy解决三次方程：错误的答案
任务是解决qubic方程 . 这个方程的系数是依赖于x的函数，所以我在每个点求解方程并得到一个根数组 . 问题在于有些点，其中由sympy给出的所有三个根都是虚构的，但是具有实系数的qubic方程是不可能的 . 此外，不同的sympy版本给出了不同的答案 . import numpy as np import pylab import sympy as sp from sympy import s... -
0 votesanswersviews
如果需要复杂数组，则传递一个真实数组
我有一个Fortran子程序，它需要一个复杂的数组 subroutine foo(cnumbers, n) integer :: n complex :: cnumbers(n) ... end subroutine foo 后来我想称之为 real :: rnumbers(40) ... call foo(rnumbers, 20) 但是，我得到编译器错误： 错误＃... -
1 votesanswersviews
MATLAB angle（）到C＃的转换
我想转移到C＃，这是一个计算MATLAB， angle() 表达式的相量角的函数 . 我找到 angle(x+yi)=atan2(y,x) 但是这里出现了我的问题，我有一个平方根，根据我给它的值是 positive 或 negative . 但是，在MATLAB中，如果sqrt函数变为负数，则返回 imaginary ，与C＃不同，它返回 NaN . 那么，我怎样才能使两个代码给出相同的结果呢...