Creating a matrix using diag

2 views (last 30 days)
John Edward
John Edward on 27 Jan 2013
Hi,
Would you please help me to write a script with "diag" command to create the matrix as shown in image. http://s1.postimage.org/5i4gd1esv/matrix.jpg By they way, it's a 99 X 99 matrix.
Thank you, John
  2 Comments
Matt J
Matt J on 27 Jan 2013
Edited: Matt J on 27 Jan 2013
The toeplitz() command would be quicker.
John Edward
John Edward on 27 Jan 2013
Thanks for your reply and I forgot to mention that it's a 99 by 99 matrix. Can you be able to script it for me.
Thank you, John

Sign in to comment.

Accepted Answer

Matt J
Matt J on 27 Jan 2013
>> z=zeros(3,1); z(1)=2; z(2)=-1; A=toeplitz(z)
A =
2 -1 0
-1 2 -1
0 -1 2
  3 Comments
Cedric
Cedric on 27 Jan 2013
Edited: Cedric on 27 Jan 2013
Look at what happens if you replace the 3 in Matt's answer with e.g. 5.
But just to be sure, you want a tri-diagonal 99x99 matrix, right?
John Edward
John Edward on 27 Jan 2013
Thank YOU, YOU are a life saver. It worked.

Sign in to comment.

More Answers (2)

Walter Roberson
Walter Roberson on 27 Jan 2013
diag(this, -1) + diag(that) + diag(theother,1)
  1 Comment
John Edward
John Edward on 27 Jan 2013
Thanks for your reply but it's a 99 by 99 matrix.

Sign in to comment.


Matt J
Matt J on 27 Jan 2013
Edited: Matt J on 27 Jan 2013
A small example
>> A=diag([1,1,1])-diag([1,1],1); A=A+A.'
A =
2 -1 0
-1 2 -1
0 -1 2
  4 Comments
John Edward
John Edward on 27 Jan 2013
Honestly, I am new to matlab and I have been working on this for past 4 hours and I can't get the code right. Can you help me on that? Thank you.
Matt J
Matt J on 27 Jan 2013
I don't know what more help I can give without just doing the whole thing for you. Are you aware that
>> help diag
or
>> doc diag
will display an explanation of the diag command? If you've read that documentation, it should be really easy. Some more examples:
>> diag([1 2 3 4]), diag([10 20 30], 2), diag([5 7 6],-1)
ans =
1 0 0 0
0 2 0 0
0 0 3 0
0 0 0 4
ans =
0 0 10 0 0
0 0 0 20 0
0 0 0 0 30
0 0 0 0 0
0 0 0 0 0
ans =
0 0 0 0
5 0 0 0
0 7 0 0
0 0 6 0

Sign in to comment.

Community Treasure Hunt

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

Start Hunting!