ActiveUser

490 Reputation

10 Badges

7 years, 115 days
i would also not like to ask, but if not ask, what should i do?

MaplePrimes Activity


These are replies submitted by ActiveUser

@Kitonum 

i edited your function for all elements in matrix, but not all cases considered, such as transpose or partly transpose

with(LinearAlgebra);
IsEquivalent:=proc(obj1::Matrix, obj2::Matrix, var::list)
local P, S;
uses combinat, LinearAlgebra;
P:=permute(var);
#print(P):
S:=[seq(subs(zip(`=`,var,p),obj2), p=P)];
#print(S):
SS := choose(S,2);
result := false:
for i from 1 to nops(SS) do
if `or`(seq(Equal(obj1,s), s=SS[i])) then 
 result := true:
end if:
od:
return result:
end proc:
elementss := proc(AA)
local elements1:
elements1 := [op(indets(AA))];
elements1 := [op(elements1), op(indets(AA,float))];
elements1 := [op(elements1),op(indets(AA,constant))];
return elements1:
end proc:
A1 := Matrix([[777,x,x],[77,7,7],[2,y,y]]);
B1 := Matrix([[777,y,y],[77,7,7],[2,x,x]]);
IsEquivalent(A1, B1, elementss(A1)); # correct
A1 := Matrix([[777,x,x],[77,7,7],[2,y,y]]);
B1 := Matrix([[777,y,y],[77,7,7],[2,x,y]]);
IsEquivalent(A1, B1, elementss(A1)); # correct

A1 := Matrix([[x,x,x],[77,7,7],[z,y,y]]);
B1 := Matrix([[z,y,y],[77,7,7],[x,x,x]]);
IsEquivalent(A1, B1, elementss(A1)); # incorrect

A1 := Matrix([[x,x,x],[77,7,7],[z,y,y]]);
B1 := Matrix([[z,y,y],[77,7,7],[x,x,y]]);
IsEquivalent(A1, B1, elementss(A1)); # correct

A1 := Matrix([[777,x,x],[77,7,7],[2,y,y]]);
B1 := Transpose(A1);
IsEquivalent(A1, B1, elementss(A1)); # incorrect

A1 := Matrix([[777,x,x],[77,7,7],[2,y,y]]);
B1 := Matrix([[777,x,2],[77,7,7],[x,y,y]]);
IsEquivalent(A1, B1, elementss(A1)); # incorrect

A1 := Matrix([[777,x,x],[77,7,7],[2,y,y]]);
B1 := Matrix([[777,77,2],[x,7,7],[x,y,y]]);
IsEquivalent(A1, B1, elementss(A1)); # incorrect
 

@Kitonum 

i got error when input a list having 3 elements x, y, z

it expects 2 arguments

how to edit the function for more than 2 elements?

@Kitonum 

is there any clues to write the extension?

what things need to consider?

@Kitonum 

 

in maple 12 , it return not true, can it change to run in maple 12?

 

second question is about two operations, such as symmetric flip and then rotate

for example:

IsEquivalent(Matrix([[777,x,x],[77,7,7],[2,y,y]]), Matrix([[777,y,y],[77,7,x],[2,7,x]]), [x,y]);

@Kitonum 

can this express in print or printf to show the operation it recognised ?

because expect to count how many same operations in series.

I have some understanding not clear is that

if there are two or three operations to achieve between two matrix , is it one equivalent can be considered as one operations for two or more operations?

@Joe Riel 

i mean self modifying maple code application

and need a database to record the change

it write maple code itself

@Kitonum 

Not use procedure 

and run directly , it can show animation now

@Kitonum 

do not know why i run in maple 2015 , no result.

 

@Kitonum 

i run your code in maple 12 and maple 2015 , no graph and animation was shown.

@Kitonum 

assume rotation matrix is (Matrix([[cos(t),sin(t)],[sin(t),cos(t)]])

how to animate a particle rotate around another particle?

how to plot the track of x and y coordinate with the change of t from -Pi to Pi?

would like to draw this track in 2D diagram but there are 3 variable?

how to set an initial point for x and y and then change t to see how x and y change in 2D coordinate?

@vv 

i understand now,

by the way,

do you know how to predict a future matrix 

belong to which one of 3 matrix?

actually eigenvector has very small difference in 

Decimal value , can linear regression help to predict the exact eigenvector values?

i find book mention magnetron , neucleon these words, but I use Time series only to find Pauli equation. Is it the end after find Pauli equation ?

how to use Pauli equation for Time series?

@asa12 

i succeed to do and completed today

@Kitonum 

@Kitonum 

 

with(LinearAlgebra):
A := Matrix(2, 2, {(1, 1) = (1/6)*sqrt(3)+(1/2)*I, (1, 2) = (1/6)*sqrt(3)-(1/2)*I, (2, 1) = (1/6)*sqrt(3)-(1/2)*I, (2, 2) = (1/6)*sqrt(3)+(1/2)*I}):
with(MTM):
factor(Re(A));
A3 := Matrix([[0,0],[0,0]]):
for i from 1 to 2 do
 for j from 1 to 2 do
  A3[i,j] := Im(A[i,j]);
 od:
od:
counterd := 0:
for i from 1 to 2 do
 for j from 1 to 2 do
  if A2[i,j] mod ((1/6)*sqrt(3)) = True then
   counterd := counterd + 1:
  end if:
 od:
od:
counterd;
if counterd = 4 then
 for i from 1 to 2 do
 for j from 1 to 2 do
  A2[i,j] := A2[i,j]/((1/6)*sqrt(3)):
 od:
od:
end if:
A2;
A2 := Matrix([[0,0],[0,0]]):
for i from 1 to 2 do
 for j from 1 to 2 do
  A2[i,j] := Im(A[i,j]);
 od:
od:
counterd := 0:
for i from 1 to 2 do
 for j from 1 to 2 do
  if A2[i,j] mod (1/2) = True then
   counterd := counterd + 1:
  end if:
 od:
od:
counterd;
if counterd = 4 then
 for i from 1 to 2 do
 for j from 1 to 2 do
  A2[i,j] := A2[i,j]/(1/2):
 od:
od:
end if:
Error, invalid input: modp received 1/2, which is not valid for its 2nd argument, m;
 
i do it in another method
mod has error in maple 12
 
how to mod in maple 12 for fraction ?

@Kitonum 

moreover , is it possible sqrt(3)/6 and 1/2 combine into one term and the rest is one matrix result?

@Kitonum 

Your doing is correct.

A := Matrix(2, 2, {(1, 1) = (1/6)*sqrt(3)+(1/2)*I, (1, 2) = (1/6)*sqrt(3)-(1/2)*I, (2, 1) = (1/6)*sqrt(3)-(1/2)*I, (2, 2) = (1/6)*sqrt(3)+(1/2)*I}):
A1:=(sqrt(3)/6)%*(Matrix(2,(i,j)->Re(A[i,j]))/(sqrt(3)/6))%+((1/2)%*(Matrix(2,(i,j)->A[i,j]-Re(A[i,j]))/(1/2)));
 
Error, (in rtable/Product) invalid arguments
in maple 12
 
how to run in maple 12?
 
 
5 6 7 8 9 10 11 Last Page 7 of 44