主题:有人知道residue 这个函数吗??
liyt100
[专家分:0] 发布于 2007-11-21 20:20:00
这个residue函数不会、、、看了书没找到。。。help英语。、看的晕晕的。。题目是这样的
试用residue命令来计算下列表达式的部分分式展开:(3*s^2+5*s+3)/(((s+1)^3)*(s^2+1))
同时再利用residue命令,验算所得答案是否正确
什么意思。。怎么做?? 大哥大姐谢谢啦!!!
回复列表 (共3个回复)
沙发
abc886y365 [专家分:0] 发布于 2007-11-22 10:03:00
当:
时
在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
板凳
abc886y365 [专家分:0] 发布于 2007-11-22 10:11:00
一个分式为1/(m*(m+1)),用residue就是把上式化为1/m-1/(m+1)
3 楼
liyt100 [专家分:0] 发布于 2007-11-22 17:15:00
>> 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 =
[]
为何结论有虚数???
我来回复