声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2494|回复: 10

[控制系统类] 一个奇怪的问题,同一个系统的传函表示和差分表示仿真结果不一样??

[复制链接]
发表于 2008-10-29 23:29 | 显示全部楼层 |阅读模式

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

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

x
我发现一个奇怪的问题,请大家帮忙指点下,谢谢各位:

被控系统是一个稳定的离散系统,输入就是简单的恒值信号(1),无反馈。

(1)、当直接用simulink仿真时(系统用Z传函表示),系统如下:

Simulink.JPG
仿真结果如下:
simulink_result.JPG

(2)、当使用Matlab直接写程序仿真时:
      被控系统用其差分方程表示:y(i) = 3.7564 * y(i-1) - 5.2922 * y(i-2)  + 3.3146 * y(i-3)  - .7788 * y(i-4) + .0000598 * u(i-1)  + .0001722 * u(i-2)  - .0001589 * u(i-3) - .000051 * u(i-4)
      输入为恒值:u(i) = 1
      可是仿真结果却是:
Matlab_result.JPG

怎么会差这么多?用matlab直接仿真的应该没多少问题啊,这个程序应该是很简单的。。。。
各位帮忙看看吧,谢谢:)
回复
分享到:

使用道具 举报

发表于 2008-10-30 08:42 | 显示全部楼层
差分方程写错了
发表于 2008-10-30 10:38 | 显示全部楼层
建议楼主列出程序
发表于 2008-11-2 22:34 | 显示全部楼层
我试了下,方程没问题,原因在于参数0.7788,末尾少了一个2,导致系统不稳定了。

评分

1

查看全部评分

发表于 2008-11-3 08:52 | 显示全部楼层
那个2 应该不影响,他的差分方程写错了
发表于 2008-11-3 10:33 | 显示全部楼层

回复 地板 xianfa110 的帖子

乖乖0.7788与0.77882误差不到0.01%, 就会导致系统不稳定了!
只不过直觉怪怪的, 除非系统本身邻近不稳定区域, 不然差那麽少会有那麽大影响吗? 感觉或许尚有其他设定有问题!
对不起, 控制这东东, 个人是门外汉, 且隔行如隔山
各位大牛, 麻烦帮帮释疑!
发表于 2008-11-4 00:00 | 显示全部楼层
呵呵,不信你试试,不用他那个差分方程,直接把那个离散传函的参数少个2,看看发散不。
看看那个传函的参数都那么小,就知道这个2的影响有多大了。

发表于 2008-11-4 00:06 | 显示全部楼层
看看系统的特征根
>> c=[1 -3.7564 5.2922 -3.3146 0.77882];
>> roots(c)

ans =

   0.9902 + 0.0324i
   0.9902 - 0.0324i
   0.8880 + 0.0701i
   0.8880 - 0.0701i

>> c=[1 -3.7564 5.2922 -3.3146 0.7788];
>> roots(c)

ans =

   1.0000         
   0.9593         
   0.8986 + 0.0667i
   0.8986 - 0.0667i
发表于 2008-11-4 09:12 | 显示全部楼层
确实发散 , 但是比较奇怪的是他的差分方程确实写错了,可是和我的正确模型的结果图像一样,怪!!!
发表于 2008-11-5 00:01 | 显示全部楼层
差分方程没写错啊!z^(-1)就表示延迟一拍。没问题的。

原帖由 yufeng 于 2008-11-4 09:12 发表
确实发散 , 但是比较奇怪的是他的差分方程确实写错了,可是和我的正确模型的结果图像一样,怪!!!
发表于 2008-11-6 10:57 | 显示全部楼层
我曾经遇到过这样的问题,在simulink里用示波器看到的波形,和输出到workspace再用plot画出来的,两个一点都不一样。就像你你这样大的差距。
我这样解决:在smulink仿真开始前按ctrl+E,在workspace I/O下,把limit data points to last 的值1000后面添三个0,就解决了。
所以我想没准是plot画图采样或者精度的问题?
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-10 22:48 , Processed in 0.078803 second(s), 23 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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