回 帖 发 新 帖 刷新版面

主题:有人知道residue 这个函数吗??

这个residue函数不会、、、看了书没找到。。。help英语。、看的晕晕的。。题目是这样的 

试用residue命令来计算下列表达式的部分分式展开:(3*s^2+5*s+3)/(((s+1)^3)*(s^2+1)) 
同时再利用residue命令,验算所得答案是否正确 

什么意思。。怎么做?? 大哥大姐谢谢啦!!!

回复列表 (共3个回复)

沙发

当:

 时 

在MATLAB中直接用分子/分母的系数表示时有
    num=[b0,b1,…,bm];
    den = [a0,a1,…,an];
则命令
    [r,p,k] = residue(num,den)
将求出两个多项式Y(s)和X(s)之比的部分分式展开的留数、极点和直接项。Y(s)/X(s)的部分分式展开由下式给出: 

 

例2-A2 考虑下列传递函数:

  

命令 [r,p,k] = residue(num,den)
将给出下列结果: 

[r,p,k]=residue(num,den)
r=
-6.000
-4.000
3.000
p=
-3.000
-2.000
-1.000
k=
2

留数为列向量r,极点位置为列向量p,直接项是行向量k。以下是Y(s)/X(s)的部分分式展开的MATLAB表达形式:

  


命令
[num,den] = residue(r,p,k)


执行后得到如下结果:

[num,den]=residue(r,p,k)
num=
2.0000 5.0000 3.0000 6.0000
den=
1.0000 6.0000 11.0000 6.0000

板凳

一个分式为1/(m*(m+1)),用residue就是把上式化为1/m-1/(m+1)

3 楼


>> b=[0 0 0 3 5 3] 
a=[1 3 4 4 3 1] 

[r, p, k] = residue(b,a) 

b =

     0     0     0     3     5     3


a =

     1     3     4     4     3     1


r =

   1.2500          
  -0.0000          
   0.5000          
  -0.6250 - 0.6250i
  -0.6250 + 0.6250i


p =

  -1.0000          
  -1.0000          
  -1.0000          
  -0.0000 + 1.0000i
  -0.0000 - 1.0000i


k =

     []


为何结论有虚数???

我来回复

您尚未登录,请登录后再回复。点此登录或注册