How can I generate a matrix as in the picture?

A is the matrix and n is an arbitrary number that larger than 3,thank you very much!

Answer 1

[Andrei Bobrov]

n = 6; % Let n = 6

z = [-1 3 -3 1];
m = numel(z);
k = n-m+1;
A = full(spdiags(ones(k,1)*z,0:m-1,k,n));

other way with Communications System Toolbox:

n = 6; % Let n = 6

z = [-1 3 -3 1];
A = convmtx(z,n-numel(z)+1);

[xiyou fu]

great!thank you!

Answer 2

[Walter Roberson]

You can add together diag() calls

[diag(-1 * ones(1, n-3)), zeros(n-3,3)] + ...
[zeros(n-3, 1), diag(3 * ones(1, n-3)), zeros(n-3,2)] + ...
[zeros(n-3, 2), diag(-3 * ones(1, n-3)), zeros(n-3,1)] + ...
[zeros(n-3, 3), diag(1 * ones(1, n-3))]