python-插值多项式，其中系数是矩阵
插值多项式是一种数学方法，用于通过已知数据点的值来估计未知数据点的值。在Python中，可以使用多种方法来实现插值多项式，其中系数可以表示为矩阵。
插值多项式的主要分类包括拉格朗日插值、牛顿插值和埃尔米特插值。这些方法都基于不同的数学原理和算法，可以根据具体的需求选择合适的方法。
优势：
- 灵活性：插值多项式可以适用于各种不规则数据点的插值，可以处理任意数量的数据点。
- 高精度：通过插值多项式可以获得较高的插值精度，可以准确地估计未知数据点的值。
- 简单易用：Python提供了丰富的插值多项式库和函数，使得实现插值多项式变得简单易用。
应用场景：
- 数据重建：插值多项式可以用于数据重建，通过已知数据点的值来估计缺失数据点的值，常见于信号处理、图像处理等领域。
- 数据逼近：插值多项式可以用于逼近复杂函数，通过已知数据点的值来估计函数的近似值，常见于数值计算、科学计算等领域。
- 数据可视化：插值多项式可以用于生成平滑的曲线，使得数据在图形上更加直观和易于理解，常见于数据可视化、图表绘制等领域。
推荐的腾讯云相关产品和产品介绍链接地址：
- 腾讯云数学计算服务（https://cloud.tencent.com/product/scs）：提供了丰富的数学计算服务，包括插值多项式等数学方法的实现和应用。
- 腾讯云数据处理服务（https://cloud.tencent.com/product/dps）：提供了数据处理和分析的解决方案，可以用于插值多项式的数据处理和应用。
请注意，以上推荐的腾讯云产品仅为示例，实际选择产品应根据具体需求和情况进行评估和决策。
相关·内容
所以我们有一个矩阵，比如
12,32
24,12
...
长度为2xN和另一个
44,32
44,19
...
长度为2xN的函数f(x，y)返回z1，z2。我们得到的2个矩阵表示x，y和z1，z2的已知值对。在这种情况下，什么是有帮助的插值公式？
浏览 0提问于2011-10-17得票数 6
回答已采纳
2回答
下面的代码是一个简单的拉格朗日插值。我有两个向量k和kt，我想对它们进行计算。我把所有的答案放在最后的一个矩阵H中。但是，当我必须为许多k和kt执行此操作时，我的代码运行所用的时间比我预期的要多。
如何以更快的方式重写此代码？
N = 1:length(kt);
H = zeros(length(kt),length(k));
for m = 1:length(k)
for j = 1:length(kt)
h = 1;
index = find(N ~= j);
for jj = 1:length(index)
浏览 1提问于2015-08-25得票数 2
1回答
我试图用matlab ( windows 8上的R2014a版本)用样条插值来计算积分。
我有函数的3个值(对于x=0,0.5,1)。
所以我有两个向量-- x和y，它们包含函数的值，我正在执行
cube_spline = spline(x,y);
coefficients = qube_spline.coefs
我期望得到2个多项式，每个多项式为3次，也就是，我期望coefficients是一个大小为2*4的矩阵，但我得到的矩阵是1*4，这意味着只有2块的1多项式。
另一方面，如果我使用4个点(即3个面板)，那么coefficients的大小是3*4 和预期的，这意味着3个面板的3多项式。
浏览 0提问于2015-01-17得票数 0
回答已采纳
我正在研究一个函数，它将为给定的有序对集生成多项式插值。您当前在一个向量中输入节点的索引，在第二个向量中输入要插值的函数值。然后，我为插值这组点的拉格朗日多项式生成一个符号表达式。我希望能够从这个符号形式转到矢量形式，以便与测试函数等进行比较。也就是说，我有一些东西可以生成一些关于符号变量x的多项式P(x)，然后我想把这个多项式采样到一个向量，并得到(例如)线性空间(-1，1,1000)上的多项式的值。如果这是可能的，我该怎么做？
我想我将包括我到目前为止拥有的代码：
function l_poly = lpoly(x,f)
% Returns the polynomial interpo
浏览 2提问于2011-02-10得票数 0
回答已采纳
1回答
在MATLAB中实现二维自反卷积
、、、、
我有一些数据，我们称之为A的2D矩阵，我知道理论上可以用另一个称为B的矩阵的自卷积来描述：
A=conv2(B,B)
我正在尝试提取B。有没有一种方法可以在MATLAB中执行2D矩阵的自反卷积？有谁能给我指个方向吗？
浏览 1提问于2015-02-11得票数 5
我有下面的代码。它完美地计算了多项式的所有y点(并将它们打印到用gnuplot绘制的图中)，但是如何得到得到的多项式(在这种情况下是1-x 2)？
void twoDegreePoly() {
int n = 3;
double x[n],y[n];
printf ("#m=0,S=16\n");
for (int i=0; i<n ;i++) {
x[i] = ((double)2*i)/2 -1;
y[i] = f(x[i]);
printf ("%g %g\n", x[i
浏览 4提问于2014-11-18得票数 1
回答已采纳
我正在研究Shamir的秘密共享算法，我很清楚它是如何工作的，但是我不明白为什么我们需要找到一个素数，并使用这个素数做模运算。
在维基百科上，它说，如果不使用模块化算法，攻击者可以在没有足够共享的情况下获得有关该值的一些信息。
在喜悦的密码学中，似乎证明了模算法的必要性，认为多项式系数需要均匀分布在Z中，这是不可行的，所以我们使用Z_p，当一个均匀分布是可实现的。
在其他网站(和这里)，我看到一些人说，简单地说，模块化是必要的，以使价值观不会变得太大。
简而言之，我找不到一个明确的理由来解释为什么真正需要模块运算。当然，这可能是上述所有理由的综合，但我感到奇怪的是，看到所有这些资料来源给出了不
浏览 0提问于2021-12-12得票数 2
回答已采纳
1回答
我试图理解pandas中的插值，但我似乎不明白“立方”方法是3阶多项式插值还是样条插值。有人知道熊猫在这种方法背后用的是什么吗？
浏览 19提问于2019-09-05得票数 0
回答已采纳
1回答
哪些程序/算法在运行时改变其数据结构的表示以获得更好的性能？
上下文：数据结构“定义”如何在计算机内存中构造和表示真实世界的概念。对于不同类型的计算，应该/可以使用不同的数据结构来实现可接受的性能(例如，链接列表与数组实现)。
自适应的数据结构是根据具体的使用模式(例如，自平衡树)改变内部状态的数据结构。这些更改是内部的，即取决于数据。此外，这些变化是预期的设计。
其他算法可以从表示的外部变化中受益。例如，在矩阵乘法中，转置“第二矩阵”(使缓存被更有效地使用)是一种众所周知的性能技巧。这实际上是将矩阵表示从主要行更改为列的主要顺序。由于"A“与" transposed (A)”
浏览 4提问于2014-07-03得票数 11
1回答
我有一系列的点表示函数的值，下面是一个例子：
X和Y的值可以是实数(非整数)。该函数是单调的，非递减的。
我希望能够对任意X(例如1.5)插值/评估函数的值，这样连续的函数线将如下所示：
这是一个标准的插值问题，所以到目前为止我使用了拉格朗日插值。它非常简单，并且给出了足够好的结果。
插值的问题是，它还会插值作为输入给出的值，因此输入数据的最终结果将是不同的(例如x=1，x=2)
有没有一个算法可以保证所有的输入值在插值后都有相同的值？线性插值是一种解决方案，但它是线性的，X之间的距离不必是偶数(那么图形就很难看)。
请原谅我的英语/数学语言，我不是英语母语者。
浏览 3提问于2014-05-21得票数 1
我想把曲线拟合到一些数据上。为了获得最好的结果，我使用了PCHIP插值。此外，我想用ppval-function得到6个区间的系数。但是会弹出一个错误，如下所示：
Error using unmkpp (line 18)
The input array does not seem to describe a pp function.
Error in ppval (line 62)
[b,c,l,k,dd]=unmkpp(pp);
Error in SA (line 8)
v = ppval(p,xdata)
这是我的代码：
clear all
xdata = [0; 3.5; 6.8;
浏览 1提问于2017-09-10得票数 0
我有作业，我被要求建立牛顿和拉格朗日插值多项式。我对拉格朗日多项式没有问题，但牛顿多项式出现了一个问题:当拉格朗日插值多项式和原始函数完全匹配时，牛顿插值不会这样做。Here is the plot. 如果我没记错的话，牛顿和拉格朗日多项式插值是表示同一多项式的不同方法，它们应该在插值区间与原始函数完全匹配。 我认为牛顿系数的计算是错误的，所以我找到了另一个除差函数。我尝试了function和they gave me the same results. 我被困在这一刻了。我仍然认为计算除差函数有问题，但我看不出有什么错误。 有什么建议吗？ 代码如下： import numpy as np
f
浏览 65提问于2021-03-26得票数 1
2回答
我想计算直线和插值函数之间的交点。我可以绘制插值图以及图中的线：
import matplotlib.pyplot as plt
import numpy as np
import scipy.interpolate as scp
y = [ 815.97, 815.95, 815.98, ..., 815.9] #60 elements
x = [405383892, 405383894, 405383895, ..., 405383896] #60 elements
mediana = np.median(x)
f = scp.interp1d(x,y,bounds_error
浏览 1提问于2016-01-22得票数 3
我希望有人能帮我。我读了几篇论文，包括贝梅尔等人的原始论文“多线性秘密共享方案”，并观看了下面的视频https://www.youtube.com/watch?v=3GFdp7V_yog，但仍在为多重线性秘密共享方案的概念而奋斗。我知道，在LSSS中，秘密是一个领域的一个要素，经销商向各方分配股份。因此，各方可以独立地使用线性映射来恢复秘密。我将在视频中使用相同的示例；在字段F_2中共享一个秘密，矩阵如下所示：
📷
我的问题是：
矩阵的第一行和第二行被标记为同一方(即P_2 )，但如何将不同的份额分配给同一方？
只有P_2和P_4才能恢复秘密，因为当它们组合它们的共享时，r_2、r_3和r_4
浏览 0提问于2021-02-11得票数 0
回答已采纳
1回答
我有一个矩阵值函数，当x到1时，我试图求它的极限。
因此，在这个例子中，我有三个矩阵v1-3，分别表示在0.85，0.9，0.99的采样值。我现在所做的工作效率很低，主要有以下几点：
for i=1:101
for j = 1:160
v_splined = spline([0.85,0.9,0.99], [v1(i,j), v2(i,j), v3(i,j)], [1]);
end
end
必须有更好、更有效的办法来做到这一点。特别是当足够快的时候，我将面对v将是4-5维向量的情况。
谢谢!
浏览 0提问于2014-04-24得票数 0
2回答
我正在做的一些研究需要象征性地采用大型矩阵的行列式；矩阵从18x18到318x318不等。矩阵项要么是数值项，要么是在同一个变量omega中的二级多项式。
目前，我正试图在.det()中使用SymPy方法，但是它非常慢；一个18x18矩阵已经运行了超过45分钟，并且在我编写这篇文章时仍在计算。我意识到行列式的计算非常密集，但是我能做些什么来加速这个过程吗？我已经在读过这篇文章，但并没有从文章中拿走任何东西，也不知道如何才能加快这个过程。我能做什么?
浏览 15提问于2017-04-20得票数 1
回答已采纳
1回答
我有一个有两个数据集的图，它产生一个轻微的梯度，其中最适合的曲线可能会被过度绘制。
目前，我只得到了一条最合适的直线。我理解scipy.optimize.curve_fit应该能够帮助我，但这需要我知道我想要夸大的函数(我认为)。
下面是我的代码和情节。如何为这些数据集创建一个曲线图呢？
plt.figure(figsize=(15,6.6))
pl.subplot(1,2,1)
plt.plot(gg,AA, 'kx')
plt.xlabel('x')
plt.ylabel('y')
plt.gca().invert_yaxis()
y=AA
浏览 4提问于2013-12-11得票数 3
回答已采纳
1回答
给定一组点x和一组值y，我试图计算最适合于最小二乘意义的P(x) = y多项式。该函数应显示范德蒙矩阵，输出多项式c应以形式p (x) = c0*x^0 + c1*x^1 + c2*x^2 +.+ cn-1^(n-1)
我想清楚地看到在同一个图上的点(十一，一)，其中的函数是绘制的。
以下是我迄今所做的尝试：
function c = interpolation(x, y)
n = length(x);
V = ones(n);
for j = n:-1:2
V(:,j-1) = x.*V(:,j);
end
c = V \ y;
浏览 0提问于2015-11-08得票数 0
回答已采纳
1回答
我有一个多项式的形式：
p(y) = A + By + Cy^2 ... + Dy^n
这里，系数A,B,..,D中的每一个都是矩阵(因此p(y)也是矩阵)。假设我在n+1点处插值多项式。我现在应该能够解决这个系统了。我正试着在Numpy做这件事。我现在有以下代码：
a = np.vander([0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1.0,1.1,1.2]) #polynomial degree is 12, a -> (12x12)
b = np.random.rand(12,60,60) #p(x) is a 60x60 matrix that I
浏览 23提问于2018-09-10得票数 1
回答已采纳
2回答
我对编程很陌生，我想我会尝试写一个分段线性插值函数。(可能是用numpy.interp还是scipy.interpolate.interp1d完成的)
假设我得到的数据如下: x= 1，2.5，3.4，5.8，6 y=2，4，5.8，4.3，4
我想设计一个分段插值函数，用Python给出在1到2.5之间，2.5到3.4之间的所有线性多项式的系数。当然，matlab有这样做的interp1函数，但是我使用python，我想做和matlab完全一样的工作，但是python只给出了值，而不是线性多项式系数！(在matlab中，我们可以用pp.coefs实现这一点)。但是如何在python nump
浏览 7提问于2017-03-24得票数 4
回答已采纳
2回答
我正在寻找一种算法，可以将点与连续的曲线连接在一起。想象一下，从a点到b点再到c点，直到最后一个点，当你从一个点画到另一个点时，这条线必须是一条曲线，并且相对于上一个点和下一个点是连续的，就好像给定的点只是一个闭环的样本。如下图所示。
这样的事情有没有这样的算法呢？
*图中的圆圈是我的点列表。
浏览 2提问于2011-07-05得票数 4
回答已采纳
我已经完成了SSS方案的第一部分工作，所以我可以使用一些秘密数作为输入，生成一些随机多项式函数，并以对的形式创建简单的共享(xi，音译)。
任务是如何从股票中重建秘密？我们都知道，我们必须做一些聪明的数学猜测，以找到科弗斯。寻找合作伙伴的选择或算法/方法是什么？每种方法的优缺点是什么？在有限的领域里，它有什么不同？
浏览 0提问于2021-12-16得票数 -1
2回答
我试图取365x365矩阵的逆。有些值高达365**365，因此它们被转换成长数字。我不知道linalg.matrix_power()函数是否能处理长数字。我知道问题来自于这个问题(因为错误消息，因为我的程序对较小的矩阵很好地工作)，但我不确定是否有办法解决这个问题。代码需要为NxN矩阵工作。
这是我的密码：
item=0
for i in xlist:
xtotal.append(arrayit.arrayit(xlist[item],len(xlist)))
item=item+1
print xtotal
xinverted=numpy.linalg.matrix_pow
浏览 9提问于2015-04-25得票数 1
回答已采纳
3回答
嘿，伙计们，我在使用函数“根”方面有很多问题。
我必须找到's^1000 + 1‘的零。我做了Y=0(1,1000)，然后手动将第1000次矩阵更改为'1‘。但是，“根”函数不适用于它！
另一个问题是，我在矩阵乘法方面有困难。问题是找到(s^6 + 6*s ^5 + 15*s^4 + 20*s^3 +15*s^2 + 6*s +1)*(s^6 + 6s^5 + 15*s^4 + 15*s^2 +6*s +1)的零(根)，所以我这样做了：
a = [1 6 15 20 15 6 1]
b = [1 6 15 0 15 6 1]
y = a.*b;
roots(y)
但这给了我
浏览 2提问于2014-06-04得票数 0
回答已采纳
从这里开始- 的批开
一个可以证明多个评估(\phi(e_i) = y_i)_{i\in I}，对于任意点e_i使用一个恒定大小的KZG批处理证明，\pi_I = g^{q_I(\tau)}，其中
\begin{align} \label{eq:batch-proof-rel} q_I(X) &=\frac{\phi(X)-R_I(X)}{A_I(X)}\\ A_I(X) &=\prod_{i\in I} (X - e_i)\\ R_I(e_i) &= y_i,\forall i\in I\\ \end{align}
R_I(X)可以通过拉格朗日插值在O(\v
浏览 0提问于2023-01-12得票数 3
回答已采纳
在shamir秘密共享方案中，将秘密S设为方程中的常数。
y_p = s+ \sum_{i=0}^{i = t-1} a_i * x_p^i
S只能是常量项，否则最后的系数或秘密可能被泄露。以下是每一个问题的情况：
沙米尔秘密共享方案中的系数
我怎样才能用数学的方法来严格地展示这一点呢？在提供的答案中，只给出了一个例子。
更新：
我目前要展示的数学方法是(t，n)系统的情况，其中t= 3，两个人之间共享，y定义为
y_1 = a_0 + s*x_1 + a_1*x_1^2模P
y_2 = a_0 + s*x_2 + a_1*x_2^2模P
然后，通过重新排列(或取矩阵的逆)，我们有
s = (
浏览 0提问于2019-12-16得票数 3
回答已采纳
1回答
在Kate等人的这论文中，给出了一种常数大小的多项式组合方案。
承诺计划假定有一个公共参考字符串：
\begin{align*} \Big\{ \{ g^{ \tau^i }, g^{ \alpha \tau^i } \}_{i=-d}^{d} \}, \ \ \{ h^{ \tau^i }, h^{ \alpha \tau^i } \}_{i=-d}^{d}, \ \ e(g, h^{\alpha} ) \Big\} \end{align*}
其中\alpha, \tau \in \mathbb{F}是两个秘密值，d是一个大整数。承诺如下：
提交:让f \in \mathbb{F
浏览 0提问于2019-07-28得票数 1
1回答
这篇文章的目的主要是用于相关研究的关键词。
无约束N问题
用M (N<=M)棋盘数N个钩子在N上的所有可能排列，这样就不会有钩子互相攻击。
解是平凡的: C(M，N)N！
约束N问题
你不能把车放在棋盘的某些地方。
例如，如果棋盘是以0-1矩阵表示的，其中0是您不能将rook放在的位置。所以矩阵的解
1 1 1
1 1 1
0 1 1
是4：
R . . | R . . | . R . | . . R
. R . | . . R | R . . | R . .
. . R | . R . | . . R | . R .
相关问题
从问题中可以很容易地修改回溯算法。然而，现在我想解决一个
浏览 0提问于2019-01-17得票数 6
回答已采纳
嗨，大家好，我有两个关于Decasteljau算法的问题，它们更多的是一个一般性的问题，但是如果我是对的，它可以帮助解决许多problems.Here，它是：
我们有这样的条件：(i=0，n)Ʃ(j=o，m) Bi，n(U)，Bi，m(v) Pi，j分析，我已经发现，首先我们对一个参数U=Uo取一些值，然后我们给其他参数v -> 1取一组点，然后在代码language.My中用一个etc....for循环增加u，问题是我们能为任何值确定一个参数u=uo，然后计算参数v上的点吗?因为所有在一条曲线上的点都在表面上，如果曲线之间的距离接近于零(实际上是这样)，我们只能将DeCasteljau
浏览 1提问于2014-07-14得票数 0
回答已采纳
我对卫星轨道的纬度、经度和高度数据做了一个，并用插值(50个点)。
现在，我只想取一个窗口(经度0-4.5度)，然后做一个更高分辨率的插值(6,000点)。我认为我需要使用第一个低分辨率拟合的拟合系数来对我的经度窗口进行插值，而我不太确定如何做到这一点。
投入：
lat = [27.755611104020687, 22.50661883405905, 17.083576087905502, 11.53891099628959, 5.916633366002468, 0.2555772624429494, -5.407902834141322, -11.037514984810027, -16
浏览 4提问于2017-10-17得票数 0
回答已采纳
2回答
多项式的最大值
、、、、
我有一个N阶多项式(其中N是偶数)。这个多项式等于x减/加无穷大的负无穷大(因此它有一个最大值)。我现在做的是用polyder求多项式的导数，然后用Matlab语言中的roots函数求N-1阶多项式的根，得到N-1个解。然后我选择真正最大化多项式的实根。问题是，我经常更新我的多项式，并且在每个时间步，我使用上面的过程来寻找最大化。因此，root函数占用了太多的计算时间，使我的应用程序变得很慢。在Matlab或提出的算法中，有没有一种方法可以以计算高效的方式实现这种最大化(即只找到一个解，而不是N-1个解)？谢谢。
编辑:我还想知道Matlab中是否有一个例程只返回实数根，而不是返回所有实数/复
浏览 3提问于2012-10-25得票数 6
回答已采纳
1回答
我有一个2D2401*266矩阵，它对应于x值(t:存储在1*266数组中)和y值(z:存储在1*2401数组中)。
我想要外推矩阵K来预测一些未来值(对应于t(1,267:279)。到目前为止，我已经使用for循环扩展了t，使其成为一个1*279矩阵：
for tq = 267:279
t(1,tq) = t(1,tq-1)+0.0333333333;
end
然而，我对如何在不对每个单独的行拟合多项式的情况下推断K感到困惑？
我觉得一定有比这更有效的方法？？
浏览 4提问于2018-01-09得票数 0
回答已采纳
5回答
特征多项式
、、、、
平方矩阵特征多项式 A定义为多项式p_A(x) = \det(Ix-A)，其中I是单位矩阵，\det是行列式。注意，这个定义总是给我们一个单次多项式，这样解决方案是唯一的。
对于这个挑战，您的任务是为整数值矩阵计算特征多项式的系数，为此，您可以使用内建，但它是不鼓励的。
规则
输入是任何方便格式的N\times N (N \ge 1)整数矩阵。
您的程序/函数将以递增或递减的顺序输出/返回系数(请指定哪个)
对系数进行赋范，使x^N系数为1(见测试用例)。
不需要处理无效输入
测试案例
系数按递减的顺序给出。x^N, x^{N-1}, ..., x^2, x, 1)：
[0] -> [1
浏览 0提问于2017-12-14得票数 13
3回答
举个例子，我有5分，
(1,1) (2,-1) (3,2) (4,-2) (5,2)
现在,
( 1)在Matlab中需要一个函数对这些点进行插值。
2)我要绘制这个函数。
3)从输入读取一个数字，并将F(x)写到输出。
我该怎么做？？
浏览 0提问于2012-11-29得票数 6
回答已采纳
2回答
我在研究一种遗传算法。以下是它的工作原理：
输入: 2D点列表
输入:曲线的度
输出:通过点的最佳方式的曲线方程(尝试最小化从点的y到曲线的垂直距离之和)
该算法为简单直线和2次方程找到了较好的方程.
但对于4点和3度方程和更多，它变得更加复杂。我无法找到正确的参数组合:有时我不得不等待5分钟，而找到的曲线仍然非常糟糕。我试着修改很多参数，从人口规模到选择的父母数量……
GA编程中著名的组合/定理对我有帮助吗？
(谢谢！)
浏览 2提问于2017-04-15得票数 2
回答已采纳
几天来，我一直在用各种公式撞墙，似乎找不到一个对以下几点有效的公式：
100,0
35,10
25,37
15,56
10,75
5,90
1,99
0,100
有没有一个好的方法可以得到一条连接所有这些点的曲线，而不会有任何大的跳跃(IE -我所拥有的最好的公式是一条多项式曲线，它在x上的负数在35到100之间非常远。)
我需要一条平滑的渐变曲线，其中每个点之间的Y值都在相同的范围内(IE -从35到100，曲线需要从10到0平滑地移动)
浏览 0提问于2014-09-22得票数 1
2回答
这是我的问题:我的光学系统是由一个照相机和一个圆形有机玻璃“透镜”组成的，它根据压力(径向弯曲)改变其曲率。
这种曲率会导致摄像机捕捉到的图像发生变形。为了纠正这种变形，需要对图像进行校准。校准可以用栅格(棋盘、点、线)进行，压力范围必须用一定的步骤来离散。对于每个压力步骤，必须拍摄网格的图像。然后，必须将每幅图像与参考图像(P=0)进行比较，并计算和存储转换矩阵。最后，在实验过程中拍摄到的每一张特定压力的图像都需要通过变换矩阵进行校正。
变形是非线性的(不仅是旋转和平移的组合)，而且很可能是桶形变形。(同样不是由摄像机诱导)看起来是这样的：
我在ImageJ中找到了一个名为BunwarpJ，
浏览 6提问于2015-02-27得票数 2
回答已采纳
1回答
极大值的特征值计算
、、、
假设我有矩阵A的如下形式
然后，我得到了矩阵C的如下形式
和最后矩阵L的如下形式
我的目标是找出矩阵L的元素的公式，使矩阵A.的特征值大于矩阵A的特征值的"K“倍。
我首先介绍了矩阵的定义：
A: matrix(
[-a,0,b,c],
[0,a,-c,b],
[d,0,-e,-1],
[0,d,1,-e]
);
C: matrix(
[1,0,0,0],
[0,1,0,0]
);
L: matrix(
[l1,-l2],
[l2,l1],
[l3,-l4],
[l4,l3]
);
然后，我找到了矩阵A的特征多项式的公
浏览 2提问于2021-10-28得票数 2
1回答
我有一个数值问题，用普通的最小二乘法( LSM )求解方程组(10次多项式)。我得到的参数有很大的值和非常小的值-因此我不能用这种方法构造矩阵-精度很低，即使在扩展变量中也是如此。我试着用C++，Matlab，delphi语言来实现。有没有人知道应用仪器，我可以做到这一点的足够的准确性或数字提示，确实得到了很好的结果。不幸的是，矩阵的标准计算很难实现。
浏览 5提问于2013-01-11得票数 0
回答已采纳
我想将二次方程拟合到图中，如下所示 ? 显示。 我已经有了"From workspace“块和"Scope”块。 ? 。 现在我的问题是如何拟合二次图。我已经有了“常量”块，其中我有“汽车”矩阵。然后我有“坡道”区块，它创建1-365之间的斜率，如 ? 显示。 然后我有“最小二乘多项式拟合”和“多项式求值”模块 ? 。我做错了什么，因为它不工作？ 这是我的matlab代码： n = importdata('n.txt',';')
cars = n(:,2)
trucks = n(:,3)
bus = n(:,4)
t = linspace(1,
浏览 17提问于2019-03-02得票数 0
回答已采纳
8回答
重采样，插值矩阵
、、、、
我正在尝试插入一些数据，以用于绘图。例如，给定N个数据点，我希望能够生成一个由大约10*N个插值数据点组成的“平滑”图。
我的方法是生成一个N乘以10*N的矩阵，并计算原始向量和我生成的矩阵的内积，得到一个1乘以10*N的向量。我已经计算出了我想要用于插值的数学，但我的代码相当慢。我对Python非常陌生，所以我希望这里的一些专家能给我一些想法，让我可以尝试加速我的代码。
我认为部分问题是生成矩阵需要对以下函数进行10*N^2次调用：
def sinc(x):
import math
try:
return math.sin(math.pi * x) / (ma
浏览 8提问于2009-12-05得票数 9
回答已采纳
1回答
如何提取与构建在多项式模型上的cv.glmnet对象的特定lambda对应的系数？当我尝试使用可能用于二项式模型的语法时，coef函数返回一个系数稀疏矩阵列表，而不是一个特定的稀疏矩阵。
举个例子：
tempcv <- cv.glmnet(x=as.matrix(iris[,-5]), y=iris[,5], family="multinomial",
nfolds=20, alpha=0.5)
coefsMin <- coef(tempcv, s="lambda.min")
当我运行这个的时候，我得到了：
浏览 3提问于2016-04-29得票数 4
2回答
求一组整数的多项式
、、、、
我一直在研究Mathematica中的一个函数，它生成并修改一组整数，然后找到该集合的插值多项式。我可以很好地在C++中生成集合，但我不知道如何模拟Mathematica的Expand[InterpolatingPolynomial[]]命令。我知道这与多项式插值问题有关，但我甚至不知道从哪里开始为它编写C++代码。
我已经向Wolfram索要了一份MathCode C++的试用版，看看它是否能为我转换过来，但我认为我更愿意尝试自己解决这个问题，所以有人能给我指出正确的方向吗?我应该如何开始做这个工作？
浏览 3提问于2012-07-04得票数 2
1回答
熊猫数据的样条插值阶
、、、、
我有以下数据，它显示来自运动捕获的数据，其中每一列都是一个标记(即位置数据)，行是时间：
LTHMB X RTHMB X
0 932.109 872.921
1 934.605 873.798
2 932.383 873.998
3 940.946 875.609
4 941.549 875.875
... ... ...
14765 NaN 602.700
14766 562.350 NaN
14767 562.394 NaN
14768 562.421 NaN
14769 562
浏览 4提问于2020-08-28得票数 1
回答已采纳
1回答
我有一个系数到1次多项式的列表，用a[i][0]*x^1 + a[i][1]
a = np.array([[ 1. , 77.48514702],
[ 1. , 0. ],
[ 1. , 2.4239275 ],
[ 1. , 1.21848739],
[ 1. , 0. ],
[ 1. , 1.18181818],
[ 1.
浏览 7提问于2022-02-24得票数 0
我想根据以下关系来求解系数{c1，c2，...cn}： ? 其中N×N矩阵B是一组单项基{1，x，x^2，x^3，....x^n}，fa(x)是近似函数。如果我想遍历n= 3,10,20,50并给定xj = 1/2*5-1/2*3*cos(j*pi/n) (Chebyshev点分布)，我应该如何使用NumPy创建这样一个矩阵B？
浏览 15提问于2021-02-11得票数 0
回答已采纳
5回答
我正在编写一些代码，以不同的速度回放WAV文件，这样波形就会变慢和低音调，或者更快、更高。我目前正在使用简单的线性插值，如下所示：
int newlength = (int)Math.Round(rawdata.Length * lengthMultiplier);
float[] output = new float[newlength];
for (int i = 0; i < newlength; i++)
{
float realPos = i / l
浏览 0提问于2009-07-14得票数 22
回答已采纳
例如，我有一个20 x 20的矩阵。在矩阵中，背景值约为1，信号值约为3.5。我添加了一个图像来表示矩阵数据：
我想通过曲面拟合来对每个点(x和y)进行曲面拟合。我这样做的方法是创建一个权重矩阵，对于任何小于分界值(例如1.1)的值，它的权重为1，而任何更大的值的权重为0。我使用以下代码创建了权重矩阵：
def Weight_Values_Excluded(data,cutoff_val=1.1):
data_mat=np.copy(data)
weight_mat=np.zeros(data_mat.shape)
mask=(data_mat<=cutof
浏览 18提问于2019-08-14得票数 0
回答已采纳
2回答
我正在尝试创建一个点列表的插值。
我有一些坐标点( ti，xi)，其中ti是时间戳，xi是关联值。我想创建一个通过这些点的函数，并允许我找到与区间内的通用t相对应的x值。
我想用三阶插值对它们进行插值。我见过类似catmull-rom插值的东西，但它只有在点xi等距时才起作用。
例如，在这里，您可以发现时间戳点是等距的，就像这里的。
有什么方法可以应用于非正则点的三次插值吗？
浏览 3提问于2011-06-06得票数 3
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
即时通信 IM
实时音视频活动推荐
腾讯云开发者
