Search All of the Math Forum:

Views expressed in these public forums are not endorsed by NCTM or The Math Forum.

Notice: We are no longer accepting new posts, but the forums will continue to be readable.

Topic: need help
Replies: 3   Last Post: Jan 28, 2013 7:47 AM

 Messages: [ Previous | Next ]
 Nasser Abbasi Posts: 6,677 Registered: 2/7/05
Re: need help
Posted: Jan 28, 2013 7:47 AM

On 1/27/2013 9:56 PM, Muhammad Ramzan wrote:
>

>> syms p
>> A = [1 p;0 1];
>> B = [1 p;0 1-p];
>>
>> mat = npermutek(['A' 'B'],4);
>>
>> for i=1:length(mat)
>> a=eval(sym(mat(i,1)));
>> b=eval(sym(mat(i,2)));
>> c=eval(sym(mat(i,3)));
>> d=eval(sym(mat(i,3)));
>> kron(kron(kron(a,b),c),d)
>> end
>> ------------------------------------
>>

>
>
> Kindly clarify and help further
>
> 1. d=eval(sym(mat(i,3))); whether it should be d=eval(sym(mat(i,4))); or not?

yes ofocurse. It was just a typo. It is clear it should have been mat(i,4).

>
> 2. I have to further use these 16 matrices, generated by the code
>you have supplied through kron(kron(kron(a,b),c),d), for multiplication with
>another 16x16 matrix, say the first one you have given in your reply and to sum up the result.
> i.e. tensor(A,A,A,A).(16x16 matrix).ctranspose(tensor(A,A,A,A))+.......+
>tensor(B,B,B,B).(16x16 matrix).ctranspose(tensor(B,B,B,B))
>
> where tensor(A,A,A,A) means the first 16x16 matrix generated by your code and given by you above.
>

I do not have time now. But why do not you try to work it out.
I suggest doing one step at a time. You are working in symbolics,
so you need to stick to Matlab functions that supports symbolics.

Take the result from one row, and do the operations above on it.
Add the code in the loop shown above.

Sorry, I have a HW to finish myself and no time.

good luck
--Nasser

Date Subject Author
1/27/13 Nasser Abbasi
1/28/13 Nasser Abbasi