声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1180|回复: 3

[编程技巧] 求极限

[复制链接]
发表于 2009-5-10 16:22 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
function    q=j_itae(x,y,z)
Kp=x; Ki=y; Kd=z;
Ti=Kp/Ki; Td=Kd/Kp; T=0.5;
numpid=[Ti*Td,Kp*Ti,1]; denpid=[0,Ti,0];
numd=[0,0,1]; dend=[1,2,1]; numr=[0,1]; denr=[1,0];
G_pid=tf(numpid,denpid);
G_d=tf(numd,dend); Rs=tf(numr,denr);
G=G_pid*G_d;
e2=(Rs/(1+G))^2;
syms  s
q=limt(e2,s,0)

大家好,这是一个传递函数,e2是误差的平方,我想通过上面的程序求出e2在s趋近于0的时候的极限,用到的是limt函数,不太会用,请求指点,谢谢!

[ 本帖最后由 ChaChing 于 2009-5-11 21:26 编辑 ]
回复
分享到:

使用道具 举报

发表于 2009-5-10 17:42 | 显示全部楼层
help limit
 楼主| 发表于 2009-5-11 15:43 | 显示全部楼层

回复 沙发 ChaChing 的帖子

我用了,但是遇到了很多问题,可以帮忙看一下吗?非常感谢!
发表于 2009-5-11 19:09 | 显示全部楼层

syms  s
q=limt(e2,s,0)
改为
[num,den]=tfdata(e2);
len=length(num{1});
syms x
x1=sum(num{1}.*(x.^(len-1:-1:0)));
x2=sum(den{1}.*(x.^(len-1:-1:0)));
limit(x1/x2,x,0)
PS:limit不能直接对e2求极限。

评分

1

查看全部评分

您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-19 05:05 , Processed in 0.070720 second(s), 22 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表