[已解决] 第二差分矩阵求解为什么无法得出抛物线?

[复制链接]
qiuhuizuo 发表于 2021-9-18 14:22:59
题目:假设 K 是 9x9 的第二差分矩阵(2’s 在对角线,-1 在对角线的上方及下方),求 解方程式 Kx = b = (10, ..., 10)。如果在 x 轴的点 1, ..., 9 的上方画出对应的点画 出 x1 , ...., x9 ,我想这 9 个点会是一个抛物线。

答案:For the −1, 2, −1 matrix K(9 by 9) and constant right side b = (10, · · · , 10), the solution x = K-1 b = (45, 80, 105, 120, 125, 120, 105, 80, 45) rises and falls along the parabola x i = 50i − 5i 2 . (A formula for K −1 is later in the text.)


可是我无论是手算,还是用matlab计算结构都是一致的,都无法求出类似抛物线的结果。我的计算结果如下。请老师看看错在哪里了?
  1. A=[1 -1 0 0 0 0 0 0 0;-1 2 -1 0 0 0 0 0 0;0 -1 2 -1 0 0 0 0 0;0 0 -1 2 -1 0 0 0 0;0 0 0 -1 2 -1 0 0 0;0 0 0 0 -1 2 -1 0 0;0 0 0 0 0 -1 2 -1 0;0 0 0 0 0 0 -1 2 -1;0 0 0 0 0 0 0 -1 2];

  2. E=[10 10 10 10 10 10 10 10 10]
  3. E=E'
  4. F=A\E

  5. F =

  6.    450
  7.    440
  8.    420
  9.    390
  10.    350
  11.    300
  12.    240
  13.    170
  14.     90
复制代码


最佳答案


15631202458 发表于 2021-9-18 15:44:55
A中第一行第一列等于2,你写错了
  1. A=[2 -1 0 0 0 0 0 0 0;-1 2 -1 0 0 0 0 0 0;0 -1 2 -1 0 0 0 0 0;0 0 -1 2 -1 0 0 0 0;0 0 0 -1 2 -1 0 0 0;0 0 0 0 -1 2 -1 0 0;0 0 0 0 0 -1 2 -1 0;0 0 0 0 0 0 -1 2 -1;0 0 0 0 0 0 0 -1 2];

  2. E=[10 10 10 10 10 10 10 10 10]
  3. E=E'
  4. F=A\E
复制代码
回复此楼

2 条回复


15631202458 发表于 2021-9-18 15:44:55
A中第一行第一列等于2,你写错了
  1. A=[2 -1 0 0 0 0 0 0 0;-1 2 -1 0 0 0 0 0 0;0 -1 2 -1 0 0 0 0 0;0 0 -1 2 -1 0 0 0 0;0 0 0 -1 2 -1 0 0 0;0 0 0 0 -1 2 -1 0 0;0 0 0 0 0 -1 2 -1 0;0 0 0 0 0 0 -1 2 -1;0 0 0 0 0 0 0 -1 2];

  2. E=[10 10 10 10 10 10 10 10 10]
  3. E=E'
  4. F=A\E
复制代码
回复此楼

qiuhuizuo 发表于 2021-9-18 19:28:35
15631202458 发表于 2021-9-18 15:44
A中第一行第一列等于2,你写错了

谢谢你的回复。问题已经解决。

我看的书前几章中给出一个3x3的第二差分矩阵。第一个元素是1。而不是2。请问所谓“第二差分矩阵”是否有严格的定义比如第一个元素。在度娘上搜了一圈没有结果。

我第一次看到书上3x3的第二差分矩阵的例子,觉得很奇怪。为什么对角线上只有第一个元素不是2。后来想是不是如此定义。既对角线上元素等于其(上方元素+左边元素)*-1 。因为对角线第一个元素没有上方元素+左边元素。所以被定义为1.
您需要登录后才可以回帖 登录 | 注册

本版积分规则

相关帖子
相关文章
热门教程
站长推荐
快速回复 返回顶部 返回列表