How to automate the matrix-vector multiplication using the index of a vector?

1 次查看（过去 30 天）
Romio
Romio 2022-7-29
评论： Matt J 2022-7-29
I have the following
m = 50;
a = 5
b = 2
c = 2
J = @(t,y) [diag(a*ones(1,m)) + diag(b*ones(1,m-1),1) + diag(c*ones(1,m-1),-1)];
I want to obtain the function/function handle (which takes (t,y) as argument) that gives essentially
A = @(t,y) [J*y]
but I want the output to be in terms of the elements of the vector y, these are y(1), y(2), ..., y(m). That is
[5*y(1) + 2*y(2);
2*y(1) + 5*y(5) + 2*y(3);
.
.
.
2*y(m-1) + 5*y(m)]

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

回答（1 个）

Matt J
Matt J 2022-7-29
Ran in:
m = 50;
a = 5;
b = 2;
c = 2;
J = [diag(a*ones(1,m)) + diag(b*ones(1,m-1),1) + diag(c*ones(1,m-1),-1)];
syms y [m,1]
expression = J*y
expression = 
  2 个评论
Romio
Romio 2022-7-29
Thanks a lot Matt. The problem is that I can't use the sympolic toolbox in my application (ode solver)
Matt J
Matt J 2022-7-29
If you can't use the symbolic toolbox, your question doesn't make sense because you've asked for a symbolic result. If you don't need a symbolic result, then J*y already gives you that.

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Number Theory 的更多信息

标签

产品


版本

R2021b

Community Treasure Hunt

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

Start Hunting!

Translated by