Scribd
Upload
29 views2 pages

Chapter5 Prob5

The document details the process of determining the principal moments of inertia of a cross-sectional area using the QR factorization and iteration method in MATLAB. It outlines the first four iterations, showing the progression of matrices and the approach towards an upper triangular matrix. The approximate eigenvalues obtained are λ1 = 7761.36 and λ2 = 1855.64, which align with previous results from Problem 5.4.

Uploaded by

skibidipapap
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
29 views2 pages

Chapter5 Prob5

The document details the process of determining the principal moments of inertia of a cross-sectional area using the QR factorization and iteration method in MATLAB. It outlines the first four iterations, showing the progression of matrices and the approach towards an upper triangular matrix. The approximate eigenvalues obtained are λ1 = 7761.36 and λ2 = 1855.64, which align with previous results from Problem 5.4.

Uploaded by

skibidipapap
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 2

1

5.5 Determine the principal moments of inertia of the cross-sectional area in Problem 5.4 by using the
QR factorization and iteration method. Carry out the first four iterations.
Solution
This is done in the MATLAB command window to preserve precision:
>> c=[5286;2914]; e=[1;0];
>> norm_c=sqrt((c(1)^2)+(c(2)^2)); v=c+(norm_c*e);
>> vT_v=v'*v; v_vT=v*v'; I=eye(2); H=I-((2/vT_v)*v_vT);
>> A=[5286 2914;2914 4331]; Q1=H; R1=H*A;
>> Q1*R1
ans =
1.0e+003 *
5.2860 2.9140
2.9140 4.3310
>> A2=R1*Q1
A2 =
1.0e+003 *
7.5274 -1.1519
-1.1519 2.0896
>> c=[7527.4;-1151.9];e=[1;0];norm_c=sqrt((c(1)^2)+(c(2)^2));
>> v=c+(norm_c*e);
>> vT_v=v'*v; v_vT=v*v'; I=eye(2); H2=I-((2/vT_v)*v_vT);
>> Q2=H2; R2=H2*A2;
>> Q2*R2
ans =
1.0e+003 *
7.5274 -1.1519
-1.1519 2.0896
>> A3=R2*Q2
A3 =
1.0e+003 *
7.7475 0.2861
0.2861 1.8695
>> c=[7747.5;286.1];e=[1;0];norm_c=sqrt((c(1)^2)+(c(2)^2));
>> v=c+(norm_c*e);
>> vT_v=v'*v; v_vT=v*v'; I=eye(2); H3=I-((2/vT_v)*v_vT);
>> Q3=H3; R3=H3*A3;
>> Q3*R3
ans =
1.0e+003 *
7.7475 0.2861
0.2861 1.8695
>> A4=R3*Q3

Excerpts from this work may be reproduced by instructors for distribution on a not-for-profit basis
for testing or instructional purposes only to students enrolled in courses for which the textbook
has been adopted. Any other reproduction or translation of this work beyond that permitted by
Sections 107 or 108 of the 1976 United States Copyright Act without the permission of the
copyright owner is unlawful.
2

A4 =
1.0e+003 *
7.7606 -0.0686
-0.0686 1.8564
>> c=[7760.6;-68.6]; e=[1;0]; norm_c=sqrt((c(1)^2)+(c(2)^2));
>> v=c+(norm_c*e);
>> vT_v=v'*v; v_vT=v*v'; I=eye(2); H4=I-((2/vT_v)*v_vT);
>> Q4=H4; R4=H4*A4;
>> Q4*R4
ans =
1.0e+003 *
7.7606 -0.0686
-0.0686 1.8564
>> A5=R4*Q4
A5 =
1.0e+003 *
7.7613 0.0164
0.0164 1.8557

From the above MATLAB session, it can be seen that the following sequence of matrices is generated, the
last one approaching an upper triangular matrix:
5286 2914 → 7527.4 – 1151.9 → 7747.5 286.1 → 7760.6 – 68.6 → 7761.3 16.4
2914 4331 – 1151.9 2089.6 286.1 1869.5 – 68.6 1856.4 16.4 1855.7
Thus, the approximate values of the eigenvalues are λ 1 = 7761.36 and λ 2 = 1855.64 , which agrees
well with the answers obtained in Problem 5.4.

Excerpts from this work may be reproduced by instructors for distribution on a not-for-profit basis
for testing or instructional purposes only to students enrolled in courses for which the textbook
has been adopted. Any other reproduction or translation of this work beyond that permitted by
Sections 107 or 108 of the 1976 United States Copyright Act without the permission of the
copyright owner is unlawful.

You might also like