Items tagged with matrices matrices Tagged Items Feed

Hi all.

Assume that we have partitioned [0,a], into N equidistant subintervals and in each subinterval we have M sets of poly nomials of the following form:

where Tm(t)=tm( namely Taylor Series) and tf is a(final point)
for Example with N=4, M=3 we have:

now we want to approximate a function, asy f(t), in this interval with following form:

If we have:

(Tau is a constant number)
then: can we find L and Z matrices using maple? I have found it manualy but for large M and N this procedure is longsome.

I used Inverse command in maple but the vector's inverse is wrong(operation not defined for Vectors)

Best Wishes

Mahmood   Dadkhah

Ph.D Candidate

Applied Mathematics Department

Hey Guys

I'm designing some stuff for controll theory

There I met following problem:

I want calculate the matrix A which is the product of X and G.

A is 6x6 matrix and all values in the diagonal  are ones, the rest is undefined and should be calculated

X is 6x6 diagonal matrix with the values X1...X6

and G is also 6x6 and have all these values.

 

now I'm stuck 'cause I haven't used Maple frequently (only some simple stuff up to now, nothing with matrices)

I'd really appreciate it, if someone can tell how to tackle this problem

Cheers

Let A and B be matrices of dimensions 3,2 and 2,3 corresondingly. It is
known Multiply(A,B)=Matrix([[8, 2, -2], [2, 5, 4], [-2, 4, 5]])


 What is Multiply(B,A)?

I want to introduce a matrix of order M by M as for any m, M, pl help as show in file

 

Help.mw

 

 

Say I have a 3x3 matrix with real entries (like 3.4) and want to compute, say, A^(11). How can I do it and only have 5 digits in the output? I set Digits:=5; and used MatrixPower -- but got 10 digits, which is way too many for the classroom presentation.

I'm trying to create and populate the Array `A` with a sequence of lower triangular matrices based on the following relation:

A[k+1][i,j] = a (1+k) A[k][i,j] + b ((1+j) A[k][i-1,j] - j A[k][i-1,j-1])

The rows and columns of each matrix are given by i and j; a and b are scalars. The structure of A is not important provided a matrix can be accessed for a given k. The highest k value will not be more than about 100.

The starting matrix for k=1 is

[-1  0]
[ 1 -1]

Any guidance on where to start with setting up the appropriate procedure would be greatly appreciated.

Thanks.

Greetings to all.

It is a new year (for some time now) and I am writing to indicate that the mathematical adventures with cycle index computations and Maple continue!

Here are the previous installments:

My purpose this time is to alert readers who might be interested to a new cycle index computation that is neither an application of the classical form of the Polya Enumeration Theorem (PET) nor of Power Group Enumeration. The former counts objects being distributed into slots with a group acting on the slots and the latter objects going into slots with a second group which permutes the objects in addition to the slots being permuted. What I am about to present treats a third possible case: when the slot permutation group and the object permutation group are one and the same and act simultaneously (not exactly the same but induced by the action of a single group).

This requires quite radical proceedings in the etymological sense of the word, which is to go back to the roots of a problem. It seems that after working with the PET sooner or later one is confronted with enumeration problems that demand the original unmitigated power of Burnside's lemma, sometimes called the lemma that is not Burnside's. This is the case with the following problem. Suppose you have an N-by-N matrix whose entries are values from 1 to N, with all assignments allowed and the symmetric group on N elements acts on the row and column indices permuting rows and columns as well as the entries simultaneously. We ask how many such matrices there are taking these double symmetries into account. This also counts the number of closed binary operations on a set of N elemnents and there is a discussion as well as the Maple code (quite simple in my opinion and no more than a few lines) that solves this problem at the following Math Stackexchange link, which uses Lovasz Formula for the cycle index of the symmetric group which some readers may remember.

In continuing the saga of Polya and Burnside exploration I have often reflected on how best to encapsulate these techniques in a Maple package. With this latest installment it would appear that a command to do Burnside enumeration probably ought to be part of such a package.

Best regards,

Marko Riedel

Is it just me that is being daft, or am I correct in ascerting that there does still not exist (even in Maple 17) a command for taking the commutator of two matrices?

As far as I can tell from the help pages, there exists only a commutator of vector fields (in the tensor package), and a commutator af algebraic quantities (in the Physics package).

Not having such a command for matrices is, of course, not critical, as it is very easily programmed, but I am just surprised, for something as fundamental as checking some representation of a Lie algebra, say, requires taking such commutators.

Hi, I am trying to write a procedure that inputs two matrices (A,B) and vector (C), it then multiplies C by A, and thenC by B and adds these to the array, it then multiplies each other them by A and B and adds them to the array etc up to (N-1)^2 (N will be given), this is what I hve so far, I am trying to write an if loop inside to say if it's already in the array then do not add it again but i'm unsure. Any help would be appreciated thanks!

 

proc_cerny1:=proc(A::Matrix,B::Matrix,C::Vector,N)
local x, S, i, j, T, R, y;
x:=(N-1)^2;
S:=Array([]);
S[0]:=C;
i:=0;
j:=0;
while (i<(2^N)) do
T=S[i].A;
S[j]=T;
j=j+1;
R=S[i].B;
S[j]=R;
i=i+1;
for y from 0 to j do
if (S[y]=T) then
S[i]=S[i]-T;
fi:
od:
print(S);
end proc:

Hi, I am currently creating a procedure so that i can input two matrices and a vector and multiply them together in different patterns, I want the result to be either one of the vectors {4,0,0,0} or {0,4,0,0} or {0,0,4,0} or {0,0,0,4}, I am trying to create an if loop that says is one of the S's is equal to one of these print ("Synchronising"), I have decalred these vectors inside the procedure and all the S's but my if statement is still not working. I think S needs to be decalred more clearly inside but i'm unsure what to do, any help would be much appreciaited! Here is my procedure:

 

proc_cerny1:=proc(A::Matrix,B::Matrix,C::Vector)
local s, s1, s2, s3, s4, s5, s6, s7, s8, s9, s10, s11, s12, s13, s14, e, f, g, h;
s:=C.A;
s1:=C.B;
s2:=C.A.B;
s3:=C.B.A;
s4:=C.A.A;
s5:=C.B.B;
s6:=C.A.A.A;
s7:=C.B.B.B;
s8:=C.A.B.B;
s9:=C.A.B.A;
s10:=C.B.B.A;
s11:=C.B.A.A;
s12:=C.A.A.B;
s12:=C.B.A.B;
s13:=C.B.A.A;
s14:=C.B.A.A.A;
e:=Vector([4,0,0,0],orientation=row);
f:=Vector([0,4,0,0],orientation=row);
g:=Vector([0,0,4,0],orientation=row);
h:=Vector([0,0,0,4],orientation=row);
g:=Vector([0,0,4,0],orientation=row);
if (s7 = e)then
print("Synchronising")
elif (s7 = g) then
print ("Synchronising")
elif (s7 = h) then
print ("Synchronising")
elif (s7 = g) then
print("Synchronising")
else
print("Not Synchronising");
end if;

return s, s1, s2, s3, s4, s5, s6, s7, s8, s9, s10, s11, s12, s13, s14;
end proc:

 

Thanks 

Find a loop do define these matrices, even for large values of n. Record the Matrix Mn 

 

all i,j=1,2,13,..,n for example 

N := `<|>`(`<,>`(1, 2), `<,>`(2, 2))

P := `<|>`(`<,>`(1, 2, 3), `<,>`(2, 2, 3), `<,>`(3, 3, 3))

Q := `<|>`(`<,>`(1, 2, 3, 4), `<,>`(2, 2, 3, 4), `<,>`(3, 3, 3, 4), `<,>`(4, 4, 4, 4))

 

Please help!!!

hey friends I am attaching my worksheet. In it am working on wavelet method now problem is in the end of worksheet am having output of  four column vector matrixes and i want all these column vector matrixes in A matrix of order n x n but problem is am not getting my desire matirx form. I want all the entires from all the vector column matrix to be wriiten in the matrix A in column wise.

new.mw

Hi all,

restart:
a:=<1,2;3,4>;
whattype(a);
b:=Matrix(2,2,[1,2,3,4]);
whattype(b);
evalb(a=b);

Is there any real difference between these two methods of creathing matrices?

How do I check if two matrices are equal?

Thanks,

Casper

indexing a matrix...

February 25 2013 jmw193 20

I'm pretty new to Maple 16.

 

I need to calculate a 6x6 transformation matrix for 8 different angles (8 different plys...classical lamination theory problem in mechanics of composites class), and I don't want to cut and paste the matrix 8 times changing the angle.  I want to put the 8 angles in a column vector (I think that's right...maybe it needs to be an array) and then have the transformation matrix call the angles one at a time to calculate a transformation matrix for each angle ...

Hello,

I am trying to make a three dimensional cube in maple by using matrices. I also have to make it rotate about the x,y,z axes, project it with three different origins, and translate it (I'm assuming this means translate it across an axis?). I have already made the cube, but when I try to rotate it across one of the axes by multiplying my matrix by one of the rotation matrices (this is what I put un maple: P1:= matrix([[1,0,0,0],[0,cos(theta),-sin(theta),0],[0,sin(theta),cos(theta),0...

1 2 3 Page 1 of 3