how to remove conj from equations? AND SOLVE THEM ?

7 views (last 30 days)
hii i am getting this kind of equations can someone please help me to solve these equations
here "t" is a variable i need values of A1,A2,A3,A4,A5,A6,Q1,Q2,Q3,Q4,Q5,Q6.
conj((A4*sin(Q4 + (2^(1/2)*t)/5))/10 == 0)
conj((A5*sin(Q5 + (2^(1/2)*t)/5))/10 == 0)
conj((A6*sin(Q6 + (2^(1/2)*t)/5))/10 == 0)
conj((10^(1/2)*A1*sin(Q1 + (2*2^(1/2)*5^(1/2)*t)/25))/100 == 1)
conj((10^(1/2)*A2*sin(Q2 + (2*2^(1/2)*5^(1/2)*t)/25))/100 == 1)
conj((10^(1/2)*A3*sin(Q3 + (2*2^(1/2)*5^(1/2)*t)/25))/100 == 0)
conj((2^(1/2)*A4*cos(Q4 + (2^(1/2)*t)/5))/50 == 0)
conj((2^(1/2)*A5*cos(Q5 + (2^(1/2)*t)/5))/50 == 0)
conj((2^(1/2)*A6*cos(Q6 + (2^(1/2)*t)/5))/50 == 0)
conj((2^(1/2)*5^(1/2)*10^(1/2)*A1*cos(Q1 + (2*2^(1/2)*5^(1/2)*t)/25))/1250 == 0)
conj((2^(1/2)*5^(1/2)*10^(1/2)*A2*cos(Q2 + (2*2^(1/2)*5^(1/2)*t)/25))/1250 == 0)
conj((2^(1/2)*5^(1/2)*10^(1/2)*A3*cos(Q3 + (2*2^(1/2)*5^(1/2)*t)/25))/1250 == 0)
THANKS IN ADVANCE

Accepted Answer

Walter Roberson
Walter Roberson on 4 Aug 2022
At some point in your code, you construct a row vector of equations, and then you use the ' operator thinking that it will transpose the row vector into a column vector. However, ' is not the transpose operator: it is the conjugate transpose operator. You need to switch from using ' to using .' which is the transpose operator

More Answers (1)

Rahul Jangid
Rahul Jangid on 4 Aug 2022
clc
clear all
close all
%% Intoroduction of variables
syms m k K1 K2 K3 K4 K5 K6 K7 K8 C1 C2 C3 C4 C5 C6 C7 C8 t
% sample data
a=4;
a1=2;
a2=(a-a1);
b=4;
b1=2;
b2=(b-b1);
c=4;
c1=2;
c2=(c-c1);
m=100;
J1=1000;
J2=1000;
J3=1000;
K=1;
C=1;
K1=K;K2=K;K3=K;K4=K;K5=K;K6=K;K7=K;K8=K;
C1=C;C2=C;C3=C;C4=C;C5=C;C6=C;C7=C;C8=C;
Xo = [0 0 0 1 1 0];
Vo = [0 0 0 0 0 0];
I = [Xo Vo].';
Mg=[m 0 0 0 0 0
0 m 0 0 0 0
0 0 m 0 0 0
0 0 0 J1 0 0
0 0 0 0 J2 0
0 0 0 0 0 J3];
%% STIFFNESS MATRIX 1
Kg = [K1+K2+K3+K4+K5+K6+K7+K8 0 0 0 -(K1+K2+K5+K6)*c1+(K3+K4+K7+K8)*c2 (K1+K2+K3+K4)*b1-(K5+K6+K7+K8)*b2
0 K1+K2+K3+K4+K5+K6+K7+K8 0 (K1+K2+K5+K6)*c1-(K3+K4+K7+K8)*c2 0 -(K1+K5+K4+K8)*a1+(K2+K3+K6+K8)*a2
0 0 K1+K2+K3+K4+K5+K6+K7+K8 -(K1+K2+K3+K4)*b1+(K5+K6+K7+K8)*b2 (K1+K5+K4+K8)*a1-(K2+K3+K6+K7)*a2 0
0 (K1+K2+K5+K6)*c1-(K3+K4+K7+K8)*c2 -(K1+K2+K3+K4)*b1+(K5+K6+K7+K8)*b2 (K1+K2+K5+K6)*c1^2+(K3+K4+K7+K8)*c2^2+(K1+K2+K3+K4)*b1^2+(K5+K6+K7+K8)*b2^2 -(K1+K4)*b1*a1+(K2+K3)*b1*a2+(K5+K8)*b2*a1-(K6+K7)*b2*a2 -(K1+K5)*c1*a1+(K2+K6)*c1*a2-(K3+K7)*c2*a2+(K4+K8)*c2*a1
-(K1+K2+K5+K6)*c1+(K3+K4+K7+K8)*c2 0 (K1+K4+K5+K8)*a1-(K2+K3+K6+K7)*a2 -(K1+K4)*b1*a1+(K2+K3)*b1*a2+(K5+K8)*b2*a1-(K6+K7)*b2*a2 (K1+K2+K5+K6)*c1^2+(K3+K4+K7+K8)*c2^2+(K1+K4+K5+K8)*a1^2+(K2+K3+K6+K7)*a2^2 -(K1+K2)*b1*c1+(K4+K3)*b1*c2+(K5+K6)*b2*c1-(K8+K7)*b2*c2
(K1+K2+K3+K4)*b1-(K5+K6+K7+K8)*b2 -(K1+K5+K4+K8)*a1+(K2+K6+K7+K3)*a2 0 -(K1+K5)*c1*a1+(K2+K6)*c1*a2-(K3+K7)*c2*a2+(K4+K8)*c2*a1 -(K1+K2)*c1*b1+(K3+K4)*b1*c2+(K5+K6)*c1*b2-(K7+K8)*c2*b2 (K1+K2+K3+K4)*b1^2+(K5+K6+K7+K8)*b2^2+(K1+K4+K5+K8)*a1^2+(K2+K3+K6+K7)*a2^2];
%% Equation of motion
format short
[U,w2] = eig(Kg,Mg);
%% Modal displacements
syms A1 A2 A3 A4 A5 A6 Q1 Q2 Q3 Q4 Q5 Q6 t real
for i = 1:6 % here [1 for x dir , 2for y dir, 3 for z dir, 4for rotation about x(i.e. alpha)]
x(i) = U(i,1)*A1*sin(sqrt(w2(1,1))*t+Q1) + U(i,2)*A2*sin(sqrt(w2(2,2))*t+Q2) + U(i,3)*A3*sin(sqrt(w2(3,3))*t+Q3) + U(i,4)*A4*sin(sqrt(w2(4,4))*t+Q4) + U(i,5)*A5*sin(sqrt(w2(5,5))*t+Q5) + U(i,6)*A6*sin(sqrt(w2(6,6))*t+Q6);
v= diff(x);
end
x;
v;
F = [x v].'
equ = F-I==0
S = solve(equ,[A1 A2 A3 A4 A5 A6 Q1 Q2 Q3 Q4 Q5 Q6])
can u please help me out .....
  1 Comment
Walter Roberson
Walter Roberson on 4 Aug 2022
What difficulty are you encountering?
When I execute on my system, it thinks for a while, and eventually responds with 625 solutions.
The A* variables are all constant for the solutions, but there are 5 unique solutions for each of Q3, Q4, Q5, Q6, each of which appears in every combination, for a total of 5^4 = 625 solutions.
Some of the solutions can be simplified if you can assume that

Sign in to comment.

Categories

Find more on Numerical Integration and Differential Equations in Help Center and File Exchange

Community Treasure Hunt

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

Start Hunting!