solve residue by matlab

plz help me ....
how can I solve this problem by using residue by matlab :
25cos(πz)/((2z+1)^2 (z-2)^2 )
plz help me .. thank u .

Risposte (2)

Star Strider
Star Strider il 24 Mar 2015
If you have the Symbolic Math Toolbox, use the partfrac funciton:
syms z
f = 25*cos(pi*z)/((2*z+1)^2 * (z-2)^2 );
pf = partfrac(f, z)
produces:
pf =
cos(pi*z)/(z - 2)^2 - (4*cos(pi*z))/(5*(z - 2)) + (8*cos(pi*z))/(5*(2*z + 1)) + (4*cos(pi*z))/(2*z + 1)^2

4 Commenti

adrian zizo
adrian zizo il 24 Mar 2015
thank u ...but I want the value of residue at (z=2 and z=-0.5) how can I find it?
Star Strider
Star Strider il 24 Mar 2015
The ‘pf’ expression provides you with the partial fraction evaluated at exactly those values. If you substitute them in for ‘z’ in that expression, you will get an infinite result for each, as expected.
adrian zizo
adrian zizo il 24 Mar 2015
but when I enter : syms z
f = 25*cos(pi*z)/((2*z+1)^2 * (z-2)^2 );
pf = partfrac(f, z)
in matlab , matlab give me error ...my matlab version is R2010a
Star Strider
Star Strider il 24 Mar 2015
According to the online documentation, partfrac was introduced in R2015a, so you wouldn’t have it in R2010a.
I posted the result as ‘pf’ in my Answer, so you have the result of the partfrac operation without having to have the function in your version.

Accedi per commentare.

Konstantinos Sofos
Konstantinos Sofos il 24 Mar 2015

0 voti

I would suggest to use L'Hôpital's rule but you could use also the following Adaptive numerical limit (and residue) estimation

2 Commenti

adrian zizo
adrian zizo il 24 Mar 2015
I want residue of complex numbers .. plz help me
Konstantinos Sofos
Konstantinos Sofos il 24 Mar 2015
Modificato: Konstantinos Sofos il 24 Mar 2015
Have a look in the link that i attached you before...
for your given example
f = @(z) 25*cos(pi*z)./((2*z+1).^2 .* (z-2).^2 )

Accedi per commentare.

Categorie

Tag

Richiesto:

il 24 Mar 2015

Commentato:

il 24 Mar 2015

Community Treasure Hunt

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

Start Hunting!

Translated by