How to create a loop on a matrix whose elements are also a row/column matrix.
조회 수: 1 (최근 30일)
이전 댓글 표시
x= [ 0:pi/18:pi/2];
A=11;
B=12;
C=9;
D=8;
K1=12;
K2=18;
syms V
E = A.*sin(x)+B.*cos(x);
F = C.*sin(x)+D.*cos(x);
Z = [ E-V K1; K2 F-V];
DT=det(Z); %%%%%determint of Matrix Z which comes a polynomial of order 2 in V
cx = coeffs(DT,V);
S0 = collect(simplify(cx(:,3)),k) %%%%%coefficient of V.^2
S1 = collect(simplify(cx(:,2)),k) %%%%%coefficient of V
S2 = collect(simplify(cx(:,1)),k) %%%%constant term
v1 = zeros(1,length(x));
for ii = 1:numel(x)
V1 = roots([S0(ii),-S1(ii),S2(ii)]);
v1(ii) = sqrt((V1(1,1))./p);
end
Firstly,I want to calculate determint of matrix Z which comes a polynomial of degree 2 in V for all value of x. Then I want to find two values of V from that polynomial....Plz Help
댓글 수: 2
Sargondjani
2020년 6월 30일
Please be more specific which part of the code is not working as you expected.
채택된 답변
KALYAN ACHARJYA
2020년 7월 1일
편집: KALYAN ACHARJYA
2020년 7월 1일
%Partial Code
DT=zeros(1,length(x));
for i=1:length(x)
Z=[E(i)-V K1; K2 F(i)-V];
DT(i)=det(Z);
end
추가 답변 (0개)
참고 항목
카테고리
Help Center 및 File Exchange에서 Polynomials에 대해 자세히 알아보기
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!