Info

Cette question est clôturée. Rouvrir pour modifier ou répondre.

I am using the following matrix to find out its determinant....but during this process iam getting the following expression and getting the error "Subscript indices must either be real positive integers or logicals."..kindly help me to remove error.

1 vue (au cours des 30 derniers jours)
Ritesh Chaurasia
Ritesh Chaurasia le 30 Août 2014
Clôturé : MATLAB Answer Bot le 20 Août 2021
m2 =[besselj(n,t), 0, -(besselj(n,x)*((pi*y*besselk(n,z)*(-bessely(n+1,y)+(n/y)*bessely(n,y)))/2 - (pi*y*bessely(n,y)*(-besselk(n+1,z)+(n/z)*besselk(n,z))*e3*eta^2*omega^2*p2)/(2*e2*q)) - bessely(n,x)*((pi*y*besselk(n,z)*(-besselj(n+1,y)+(n/y)*besselj(n,y)))/2 - (pi*y*besselj(n,y)*(-besselk(n+1,z)+(n/z)*besselk(n,z))*e3*eta^2*omega^2*p2)/(2*e2*q))), -(- bessely(n,x)*((pi*besselj(n,y)*besselk(n,z)*beta1*n*1i)/(2*e2*omega) + (pi*besselj(n,y)*besselk(n,z)*beta1*eta^2*n*omega*p2^2*1i)/(2*e2*q^2)) + besselj(n,x)*((pi*besselk(n,z)*bessely(n,y)*beta1*n*1i)/(2*e2*omega) + (pi*besselk(n,z)*bessely(n,y)*beta1*eta^2*n*omega*p2^2*1i)/(2*e2*q^2)));...
0, eta*besselj(n,t), -(bessely(n,x)*eta*((pi*besselj(n,y)*besselk(n,z)*beta1*n*1i)/(2*omega) + (pi*besselj(n,y)*besselk(n,z)*beta1*n*p2^2*1i)/(2*omega*q^2)) - besselj(n,x)*eta*((pi*besselk(n,z)*bessely(n,y)*beta1*n*1i)/(2*omega) + (pi*besselk(n,z)*bessely(n,y)*beta1*n*p2^2*1i)/(2*omega*q^2))), -(besselj(n,x)*eta*((pi*y*besselk(n,z)*(-bessely(n+1,y)+(n/y)*bessely(n,y))/2 + (pi*bessely(n,y)*(-besselk(n+1,z)+(n/z)*besselk(n,z))*p2^2*r2)/(2*q)) - bessely(n,x)*eta*((pi*y*besselk(n,z)*(-besselj(n+1,y)+(n/y)*besselj(n,y))/2 + (pi*besselj(n,y)*(-besselk(n+1,z)+(n/z)*besselk(n,z))*p2^2*r2)/(2*q)))));...
((beta1*n)/p1^2*r1)*besselj(n,t), ((1i*omega)/p1)*(-besselj(n+1,t)+(n/t)*besselj(n,t)), -((besselj(n,x)*beta1*n*((pi*y*besselk(n,z)*(-bessely(n+1,y)+(n/y)*bessely(n,y))/2 - (pi*y*bessely(n,y)*(-besselk(n+1,z)+(n/z)*besselk(n,z))*e3*eta^2*omega^2*p2)/(2*e2*q)))/(p2^2*r1) - (bessely(n,x)*beta1*n*((pi*y*besselk(n,z)*(-besselj(n+1,y)+(n/y)*besselj(n,y)))/2 - (pi*y*besselj(n,y)*(-besselk(n+1,z)+(n/z)*besselk(n,z))*e3*eta^2*omega^2*p2)/(2*e2*q)))/(p2^2*r1) - ((-besselj(n+1,x)+(n/x)*besselj(n,x))*omega*((pi*besselk(n,z)*bessely(n,y)*beta1*n*1i)/(2*omega) + (pi*besselk(n,z)*bessely(n,y)*beta1*n*p2^2*1i)/(2*omega*q^2))*1i)/p2 + (-bessely(n+1,x)+(n/x)*bessely(n,x))*omega*((pi*besselj(n,y)*besselk(n,z)*beta1*n*1i)/(2*omega) + (pi*besselj(n,y)*besselk(n,z)*beta1*n*p2^2*1i)/(2*omega*q^2))*1i)/p2), -(-((-bessely(n+1,x)+(n/x)*bessely(n,x))*omega*((pi*y*besselk(n,z)*(-besselj(n+1,y)+(n/y)*besselj(n,y)/2 + (pi*y*besselj(n,y)*(-besselk(n+1,z)+(n/z)*besselk(n,z))*p2)/(2*q))*1i)/p2 + ((-besselj(n+1,x)+(n/x)*besselj(n,x))*omega*((pi*y*besselk(n,z)*(-bessely(n+1,y)+(n/y)*bessely(n,y)))/2 + (pi*y*bessely(n,y)*(-besselk(n+1,z)+(n/z)*besselk(n,z))*p2)/(2*q))*1i)/p2 - (bessely(n,x)*beta1*n*((pi*besselj(n,y)*besselk(n,z)*beta1*n*1i)/(2*e2*omega) + (pi*besselj(n,y)*besselk(n,z)*beta1*eta^2*n*omega*p2^2*1i)/(2*e2*q^2)))/(p2^2*r1) + (besselj(n,x)*beta1*n*((pi*besselk(n,z)*bessely(n,y)*beta1*n*1i)/(2*e2*omega) + (pi*besselk(n,z)*bessely(n,y)*beta1*eta^2*n*omega*p2^2*1i)/(2*e2*q^2)))/(p2^2*r1))));...
((-1i*omega*e1*eta)/p1)*(-besselj(n+1,t)+(n/t)*besselj(n,t)), ((eta*beta1*n)/p1^2*r1)*besselj(n,t), -(((-bessely(n+1,x)+(n/x)*bessely(n,x))*e1*eta*omega*((pi*y*besselk(n,z)*(-besselj(n+1,y)+(n/y)*besselj(n,y))/2 - (pi*y*besselj(n,y)*(-besselk(n+1,z)+(n/z)*besselk(n,z))*e3*eta^2*omega^2*p2)/(2*e2*q))*1i)/p2 - ((-besselj(n+1,x)+(n/x)*besselj(n,x))*e1*omega*eta((pi*y*besselk(n,z)*(-bessely(n+1,y)+(n/y)*bessely(n,y)))/2 - (pi*y*bessely(n,y)*(-besselk(n+1,z)+(n/z)*besselk(n,z))*e3*eta^2*omega^2*p2)/(2*e2*q))*1i)/p2 + (bessely(n,x)*beta1*eta*n*((pi*besselj(n,y)*besselk(n,z)*beta1*n*1i)/(2*omega) + (pi*besselj(n,y)*besselk(n,z)*beta1*n*p2^2*1i)/(2*omega*q^2)))/(p2^2*r1) - (besselj(n,x)*beta1*eta*n*((pi*besselk(n,z)*bessely(n,y)*beta1*n*1i)/(2*omega) + (pi*besselk(n,z)*bessely(n,y)*beta1*n*p2^2*1i)/(2*omega*q^2)))/(p2^2*r1))), -((besselj(n,x)*beta1*eta*n*((pi*y*besselk(n,z)*(-bessely(n+1,y)+(n/y)*bessely(n,y))/2 + (pi*y*bessely(n,y)*(-besselk(n+1,z)+(n/z)*besselk(n,z))*p2)/(2*q)))/(p2^2*r1) - (bessely(n,x)*beta1*eta*n*((pi*y*besselk(n,z)*(-besselj(n+1,y)+(n/y)*besselj(n,y)))/2 + (pi*y*besselj(n,y)*(-besselk(n+1,z)+(n/z)*besselk(n,z))*p2)/(2*q)))/(p2^2*r1) - ((-besselj(n+1,x)+(n/x)*besselj(n,x)))*e1*eta*omega*((pi*besselk(n,z)*bessely(n,y)*beta1*n*1i)/(2*e2*omega) + (pi*besselk(n,z)*bessely(n,y)*beta1*eta^2*n*omega*p2^2*1i)/(2*e2*q^2))*1i)/p2 + ((-bessely(n+1,x)+(n/x)*bessely(n,x))*e1*omega*eta((pi*besselj(n,y)*besselk(n,z)*beta1*n*1i)/(2*e2*omega) + (pi*besselj(n,y)*besselk(n,z)*beta1*eta^2*n*omega*p2^2*1i)/(2*e2*q^2))*1i)/p2)];
  5 commentaires
Geoff Hayes
Geoff Hayes le 30 Août 2014
Ritesh - what line in the code does the error message correspond to? What are the dimensions of x, y, z, and t? What is n, eta, omega? I'm guessing that
m2=[m11 m12 m13 m14;m21 m22 m23 m24;m31 m32 m33 m34;m41 m42 m43 m44];
is the last line of code. What are the values for each of these mZY? Are they scalars or something else?
dpb
dpb le 30 Août 2014
Continue to break it down into smaller and smaller pieces until you find the offending term.
Which of the expressions specifically is(are) the offending one(s) and the exact error text...I forget whether it will echo some of the subscript it thinks is the bad one or not...hmmm, a simple test here just echos the numeric value of a bum index, not the expression.
So, to proceed
a) start with the term that errors b) take each term inside and check it
Oh, and first try
whos bessel*
to see if you've got a previous typo which created an array of one of the function neames hanging around then
clear bessel*
before trying again and then recheck for the existance of an array.

Réponses (1)

per isakson
per isakson le 31 Août 2014
Modifié(e) : per isakson le 31 Août 2014
It looks like missing "*" after eta in
...*omega*eta((pi*y*besselk(n,z)*...
...*omega*eta((pi*besselj(n,y)...
Search for "eta(("
&nbsp
Take away
  • Keep the line length shorter than the width of the editor window
  • [Mark to matching ...] in Notepad++ was useful.
New on my MATLAB wishlist
  • [Mark to matching ... and Evaluate]

Community Treasure Hunt

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

Start Hunting!

Translated by