在Mathematica中,如果我有一个符号代数方程,我如何获得乘以特定项的因子?
我的意思是,比如说,作为一些计算的结果,我得到了一个符号方程a x^2 + b x + c y,并说我对任何乘以x^2的东西感兴趣(即这里的a )。我应该编写什么语句,以便应用于a x^2 + b x + c y时返回a。
(请注意,在本例中有一个b x术语,因此我不能使用导数来提取因子)
我有两个红色和蓝色的(x,y)元组列表和一个ax+by+c形式的线方程列表,我的要求是从每条线方程中提取系数,并根据这2组点的图解确定这些点在线的两侧是否清楚地分开。挑战是我不能用矮胖。
我的方法是使用pyplot压缩2个列表中的红色和蓝色点。现在,我尝试使用正则表达式提取系数,如下所示。
lines = ["1x+1y+0","1x-1y+0","1x+0y-3","0x+1y-0.5"]
for i in lines:
z = re.match('(\d+)?(x)?\+(\d+)?(y)?\+(\d)?
假设我有下面的矩阵:
syms x y z
M = [x+y-z;2*x+3*y+-5*z;-x-y6*z];
我希望有一个矩阵,它包含变量x,y和z的系数:
CM = [1,1,-1;2,3,-5;-1,-1,6];
如果我把CM乘以[x;y;z],我希望得到M。
编辑
我有一个颂歌系统:
(d/dt)A = B
A和B是方阵。我想要解这组方程。我不想使用Matlab的ode解题命令。如果我将上述方程组转化为:
(d/dt)a = M*a
然后通过矩阵M的特征向量和值很容易地求解,这里a是包含变量的列向量,M是从B中提取的系数矩阵。
对于FIR系统的术语和模拟,我有一些困惑。我将感激帮助纠正我的错误和告诉什么是正确的。
假设FIR滤波器具有系数阵列A=[1,c2,c3,c4]。元素的数量是L,所以过滤器的长度是L,但是顺序是L-1。
Confusion1:被认为是一个系数吗?总是一个吗?
Confusion2:,我的理解是,对于给定的例子,长度L= 4和order=3是正确的吗?
Confusion3:从数学上讲,我可以把它写成:
其中u是输入数据,l从零开始。然后,为了模拟上述方程,我做了下面的卷积。对吗?:
N =100; %number of data
A = [1, 0.1, -0.5, 0.62];
u
我是通过给定的一组数据点生成三次样条图:
import matplotlib.pyplot as plt
import numpy as np
from scipy import interpolate
x = np.array([1, 2, 4, 5]) # sort data points by increasing x value
y = np.array([2, 1, 4, 3])
arr = np.arange(np.amin(x), np.amax(x), 0.01)
s = interpolate.CubicSpline(x, y)
plt.plot(x, y, 'b
如何提取R中正交多项式回归的系数?
它就像一个带有原始回归的护身符:
#create datas
set.seed(120)
x= 1:20
y=x^2 + rnorm(length(x),0,10)
datas = data.frame(x,y)
#compute raw model with function poly
mod= lm(data=datas, y~poly(x,2,raw=T))
#get coefficients with function coef()
coefficients = coef(mod)
#construct polynom and check fi
在发布实际代码之前,让我向您展示我的计算机的处理器和内存信息是很好的:
昨天,我发表了关于Lorenz方程(混沌理论中的经典方程)的文章,其中一位伟人帮助我给出了解,如下所示:
function f=lorenz(t,x,a,b,c)
% solve differential equation like this
%dx/dt=a*(y-x)
%dy/dt=-x*z+b*x-y
%dz/dt=xy-c*z/3
f=zeros(3,1);% preallocate result
f(1)=a