A problem about residues(is matlab wrong?)

  • Context: MATLAB 
  • Thread starter Thread starter oahsen
  • Start date Start date
  • Tags Tags
    Matlab
Click For Summary
SUMMARY

The MATLAB code provided calculates the residues of the function 1/(z-i/9)^3 using the residue function, yielding results of residues 0 and 1 at the pole z=i/9. The pole is confirmed to be of order three, yet the presence of two different residue values raises questions about MATLAB's algorithm. The residue at a pole is defined as the coefficient of 1/(z-a) in the Laurent expansion, which in this case should only yield a single residue of zero. The discrepancy suggests a misunderstanding of MATLAB's output or a potential issue with the implementation.

PREREQUISITES
  • Understanding of complex analysis, specifically residues and poles.
  • Familiarity with MATLAB programming and its residue function.
  • Knowledge of Laurent series and their application in complex functions.
  • Basic concepts of contour integration in complex analysis.
NEXT STEPS
  • Review MATLAB documentation on the residue function for detailed understanding.
  • Study the properties of residues in complex analysis, focusing on poles of higher order.
  • Learn about Laurent series expansions and their role in determining residues.
  • Explore contour integration techniques and their implications in complex function theory.
USEFUL FOR

Students and professionals in mathematics, particularly those studying complex analysis, as well as MATLAB users seeking to understand residue calculations and their implications in contour integrals.

oahsen
Messages
58
Reaction score
0
I have written this simple code fragment to the MATLAB for finding the residue(s) of the function 1/(z-i/9)^3;

b=[1];
a=[1 -i/3 -1/27 i/(729) ];
[r p k]=residue(b,a)

and get the following result;
r =

0
0
1


p =

0.0000 + 0.1111i
0.0000 + 0.1111i
0.0000 + 0.1111i


k =

[]

The poles are true. Function has a pole of order three at z=i/9 . However, there are two different values for the residue: 0 and 1. I could not get the meaning behind that? How can a function have different residue at the same point? Additionally if we assume MATLAB is right then the result of the contour integral should be 2pi*i at the unit circle. However, we can easily show that this contour integral is zero. So, is something wrong with the residue algorithm of MATLAB or there is something I do not consider?
 
Physics news on Phys.org
The residue of a pole in a function f(z) at z=a is the coefficient of 1/(z-a) in the Laurent expansion of f(z). In your case, a=i/9 and the coefficeint of 1/(z-i/9) is zero. So the residue (there is only one residue!) is zero.
 
Avodyne said:
The residue of a pole in a function f(z) at z=a is the coefficient of 1/(z-a) in the Laurent expansion of f(z). In your case, a=i/9 and the coefficeint of 1/(z-i/9) is zero. So the residue (there is only one residue!) is zero.

Yes, logically it should have only one residue which is zero. However, what I did not understand is why has MATLAB find also a residue which is 1 at the pol z=i/9 ?
 
I don't use MATLAB so I can't help.
 

Similar threads

  • · Replies 12 ·
Replies
12
Views
4K
  • · Replies 4 ·
Replies
4
Views
2K
  • · Replies 10 ·
Replies
10
Views
3K
  • · Replies 9 ·
Replies
9
Views
3K
  • · Replies 6 ·
Replies
6
Views
4K
  • · Replies 2 ·
Replies
2
Views
4K
  • · Replies 3 ·
Replies
3
Views
5K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 8 ·
Replies
8
Views
2K
  • · Replies 2 ·
Replies
2
Views
3K