Not getting expected results

Het guys, i attached a code here. The problem is that the very first 'R' matrix i have to get is a Identity Matrix. Unfortunately, i don't know what went wrong there is some mistake. Could anyone shed some light, what i did wrong?. Thanking you in advance.

Answer 1

[James Tursa]

I get the identity matrix with your code:

>> phi = 0;
>> psi = 0;
>> theta = 0;
>>         R=[cosd(phi).*cosd(psi)-sind(phi).*cosd(theta).*sind(psi)   sind(phi).*cosd(psi)+cosd(phi).*cosd(theta).*sind(psi)   sind(theta).*sind(psi)
            -cosd(phi).*sind(psi)-sind(phi).*cosd(theta).*cosd(psi) -sind(phi).*sind(psi)+cosd(phi).*cosd(theta).*cosd(psi) sind(theta).*cosd(psi)
            sind(phi).*sind(theta)  -cosd(phi).*sind(theta)  cosd(theta)]
R =
     1     0     0
     0     1     0
     0     0     1

[satish kumar]

Thanks James for the response, yes for me also if i take independently my psi and theta values. But, i want to get that results by exactly passing through the for loop. i have to check for lots of variation in theta and psi values there I expect the first value of the 'R' as Identitymatrix which is not happening at now.

When,thetaV = (0:10:90); % euler angles for a rotation

psiV = (0:1:90); % euler angles for a rotation

and these going through for loop should get me Identity Matrix which i am not getting. Could you just see through this again. Thanks again James.

[James Tursa]

I copied your entire code and ran it in the debugger, pausing it at the first iteration. R is the identitiy matrix.

[satish kumar]

Really, james. I need one more help from you redarding the next part of this code. should i write here?

[James Tursa]

Sure.

[satish kumar]

Basically, i am expecting plots like screenshot(26).png which i am attaching here by running the code practice3.m . Could you just help me through.