solving a sequence in matlab

2015 3 6
1 个回答
回答已采纳
更新时间：2015 3 6
9 次查看（30 天）

0 个投票

Ax1 = b,
Ax2 = x1
Ax3 = x2,
Ax4 = x3,
Ax = x4,
how to solve this sequence A and b are defined
A = diag(2*ones(1,n)) - diag(ones(1,n-1),1) - diag(ones(1,n-1),-1)
b = [0:1:n/2-1 n/2-1:-1:0]’

1 个评论
显示 -1更早的评论 隐藏 -1更早的评论

James Tursa
James Tursa 2015-3-6
Please explain in more detail. Does Ax1 mean A*x1, so that the line Ax1 = b means A*x1 = b and you are solving for x1? Then you use that to solve for x1, x3, x4, and ultimately x?

请先登录，再进行评论。

请先登录，再回答此问题。

 采纳的回答

James Tursa
James Tursa 2015-3-6
编辑：James Tursa 2015-3-6

0 个投票

Plugging everything in I get
A^5 * x = b
So the solution would be
x = A^5 \ b
However, your expressions for A and b don't seem right since the dimensions are incompatible for the linear expression. E.g.,
>> n = 3
n =
3
>> A = diag(2*ones(1,n)) - diag(ones(1,n-1),1) - diag(ones(1,n-1),-1)
A =
2 -1 0
-1 2 -1
0 -1 2
>> b = [0:1:n/2-1 n/2-1:-1:0]'
b =
0
0.500000000000000
>> A^5\b
Error using \
Matrix dimensions must agree.

2 个评论
显示 无 隐藏 无

HIGH
HIGH 2015-3-6
编辑：HIGH 2015-3-6
Thanks for the reply. Actually A and b are correct and I figured the X's out. but I need to write a function for "upper bidiagonal system" to find these X'using the attached m files. I tried couple of ways but yet having trouble writing the function.
Ohh I am sorry, n should be equal to 10.. n=10
HIGH
HIGH 2015-3-6
n=10 sorry

请先登录，再进行评论。

更多回答（0 个）

类别

在 帮助中心 和 File Exchange 中查找有关 Data Type Conversion 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by