How to get numerical answers for solving determinant is 0
조회 수: 19 (최근 30일)
이전 댓글 표시
I want to solve determinant of an 8x8 matrix is zero. Some of the element of the matrix has symbolic "p" in it. I want to know what value of "p" will make the determinat of the matrix to be 0.
The result returns a different super long functions with Z. Each funtion is displayed on a line.
I am not really sure if this is the correct result and I want a numerical answer for "p".
format shortE
E = 210 * 10^9;
k1 = zeros(8, 8);
k2 = zeros(8, 8);
syms p
k1(1, 1) = 15.24;
k1(1, 2) = 2.28;
k1(1, 3) = -15.24;
k1(1, 4) = 2.28;
k1(2, 1) = 2.28;
k1(2, 2) = 0.46;
k1(2, 3) = -2.28;
k1(2, 4) = 0.23;
k1(3, 1) = -15.24;
k1(3, 2) = -2.28;
k1(3, 3) = 22.86;
k1(3, 4) = 0;
k1(3, 5) = -7.63;
k1(3, 6) = 2.28;
k1(4, 1) = 2.28;
k1(4, 2) = -0.23;
k1(4, 3) = 0;
k1(4, 4) = 1.37;
k1(4, 5) = -2.28;
k1(4, 6) = -0.46;
k1(5, 3) = -7.62;
k1(5, 4) = -2.28;
k1(5, 5) = 22.86;
k1(5, 6) = 0;
k1(5, 7) = -15.24;
k1(5, 8) = 2.28;
k1(6, 3) = 2.28;
k1(6, 4) = 0.46;
k1(6, 5) = 0;
k1(6, 6) = 1.37;
k1(6, 7) = -2.28;
k1(6, 8) = 0.23;
k1(7, 5) = -15.24;
k1(7, 6) = -2.28;
k1(7, 7) = 15.24;
k1(7, 8) = -2.28;
k1(8, 5) = 2.28;
k1(8, 6) = -0.23;
k1(8, 7) = -2.28;
k1(8, 8) = 0.46;
k2(1, 1) = 4;
k2(1, 2) = 0.1;
k2(1, 3) = -4;
k2(1, 4) = 0.1;
k2(2, 1) = 0.1;
k2(2, 2) = 0.04;
k2(2, 3) = -0.1;
k2(2, 4) = -0.01;
k2(3, 1) = -4;
k2(3, 2) = -0.1;
k2(3, 3) = 6;
k2(3, 4) = 0;
k2(3, 5) = -2;
k2(3, 6) = 0.1;
k2(4, 1) = 0.1;
k2(4, 2) = -0.01;
k2(4, 3) = 0;
k2(4, 4) = 0.12;
k2(4, 5) = -0.1;
k2(4, 6) = -0.02;
k2(5, 3) = -2;
k2(5, 4) = -0.1;
k2(5, 5) = 6;
k2(5, 6) = 0;
k2(5, 7) = -4;
k2(5, 8) = 0.1;
k2(6, 3) = 0.1;
k2(6, 4) = -0.02;
k2(6, 5) = 0;
k2(6, 6) = 0.12;
k2(6, 7) = -0.1;
k2(6, 8) = -0.01;
k2(7, 5) = -4;
k2(7, 6) = -0.1;
k2(7, 7) = 4;
k2(7, 8) = -0.1;
k2(8, 5) = 0.1;
k2(8, 6) = -0.01;
k2(8, 7) = -0.1;
k2(8, 8) = 0.04;
k = k1 .* E .* (10^-5) - p .* k2;
result = det(k);
solve(result == 0)
댓글 수: 0
채택된 답변
Torsten
2022년 11월 22일
편집: Torsten
2022년 11월 22일
This is a generalized eigenvalue problem. Use "eig" to solve for p:
E = 210 * 10^9;
k1 = zeros(8, 8);
k2 = zeros(8, 8);
k1(1, 1) = 15.24;
k1(1, 2) = 2.28;
k1(1, 3) = -15.24;
k1(1, 4) = 2.28;
k1(2, 1) = 2.28;
k1(2, 2) = 0.46;
k1(2, 3) = -2.28;
k1(2, 4) = 0.23;
k1(3, 1) = -15.24;
k1(3, 2) = -2.28;
k1(3, 3) = 22.86;
k1(3, 4) = 0;
k1(3, 5) = -7.63;
k1(3, 6) = 2.28;
k1(4, 1) = 2.28;
k1(4, 2) = -0.23;
k1(4, 3) = 0;
k1(4, 4) = 1.37;
k1(4, 5) = -2.28;
k1(4, 6) = -0.46;
k1(5, 3) = -7.62;
k1(5, 4) = -2.28;
k1(5, 5) = 22.86;
k1(5, 6) = 0;
k1(5, 7) = -15.24;
k1(5, 8) = 2.28;
k1(6, 3) = 2.28;
k1(6, 4) = 0.46;
k1(6, 5) = 0;
k1(6, 6) = 1.37;
k1(6, 7) = -2.28;
k1(6, 8) = 0.23;
k1(7, 5) = -15.24;
k1(7, 6) = -2.28;
k1(7, 7) = 15.24;
k1(7, 8) = -2.28;
k1(8, 5) = 2.28;
k1(8, 6) = -0.23;
k1(8, 7) = -2.28;
k1(8, 8) = 0.46;
k2(1, 1) = 4;
k2(1, 2) = 0.1;
k2(1, 3) = -4;
k2(1, 4) = 0.1;
k2(2, 1) = 0.1;
k2(2, 2) = 0.04;
k2(2, 3) = -0.1;
k2(2, 4) = -0.01;
k2(3, 1) = -4;
k2(3, 2) = -0.1;
k2(3, 3) = 6;
k2(3, 4) = 0;
k2(3, 5) = -2;
k2(3, 6) = 0.1;
k2(4, 1) = 0.1;
k2(4, 2) = -0.01;
k2(4, 3) = 0;
k2(4, 4) = 0.12;
k2(4, 5) = -0.1;
k2(4, 6) = -0.02;
k2(5, 3) = -2;
k2(5, 4) = -0.1;
k2(5, 5) = 6;
k2(5, 6) = 0;
k2(5, 7) = -4;
k2(5, 8) = 0.1;
k2(6, 3) = 0.1;
k2(6, 4) = -0.02;
k2(6, 5) = 0;
k2(6, 6) = 0.12;
k2(6, 7) = -0.1;
k2(6, 8) = -0.01;
k2(7, 5) = -4;
k2(7, 6) = -0.1;
k2(7, 7) = 4;
k2(7, 8) = -0.1;
k2(8, 5) = 0.1;
k2(8, 6) = -0.01;
k2(8, 7) = -0.1;
k2(8, 8) = 0.04;
A = k1 .* E .* (10^-5);
B = k2;
format long
p = eig(A,B)
댓글 수: 0
추가 답변 (0개)
참고 항목
카테고리
Help Center 및 File Exchange에서 Linear Algebra에 대해 자세히 알아보기
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!