完善资料让更多小伙伴认识你,还能领取20积分哦, 立即完善>
|
|
相关推荐
1个回答
|
|
数学过程推导
具'j体的数i学推导过程a就不做详绍了,在现代控制理论教材中都有,这儿贴两张图片 根据上面公式推导,假设存在系统,其中, 其离散系统表达式为 求系统的状态转移矩阵,式中为拉氏反变换 matlab求解过程为 A = [0 1;0 -2]; B1 = [0 ;1]; syms s t Ls; % 求状态转移矩阵 利用拉氏变换,syms为符号函数用来定义数学函数 I = eye(size(A)); Ls = inv(s*I - A); STM = ilaplace(Ls,s,t) %状态转移矩阵,ilaplace为拉氏反变换函数 同理可以求H syms s t Ls; % 求状态转移矩阵 利用拉氏变换 I = eye(size(A)); Ls = inv(s*I - A);% collect 函数为合并同类项 STM = ilaplace(Ls,s,t) %状态转移矩阵 syms T HLs = int(STM,t,0,T); H = HLs*B1 运行结果如下 在simulink仿真过程中,我们需要G和H为已知的矩阵,当T固定时就可以求出G和H的具体数值了,这是需要用到符号函数求解函数subs,设步长T为0.001;matlab命令为 dt = 0.001;A = [0 1;0 -2];B1 = [0 ;1];syms s t Ls; % 求状态转移矩阵 利用拉氏变换 I = eye(size(A)); Ls = inv(s*I - A);% collect 函数为合并同类项 STM = ilaplace(Ls,s,t); %状态转移矩阵 G = double(subs(STM,t,dt)) % 符号函数求解syms T HLs = int(STM,t,0,T); H = HLs*B1; H = real(double(subs(H,T,dt))) % 符号函数求解 结果如下: 到这matlab部分就完成了剩下的就是在simulink中的具体仿真了,这个有需要再写吧,贴下自己做的一个简单仿真和连续系统同时进行的以便对数据对比,结果是连续系统和离散系统仿真数据是一样的。 |
|
|
|
只有小组成员才能发言,加入小组>>
调试STM32H750的FMC总线读写PSRAM遇到的问题求解?
1771 浏览 1 评论
X-NUCLEO-IHM08M1板文档中输出电流为15Arms,15Arms是怎么得出来的呢?
1619 浏览 1 评论
1070 浏览 2 评论
STM32F030F4 HSI时钟温度测试过不去是怎么回事?
724 浏览 2 评论
ST25R3916能否对ISO15693的标签芯片进行分区域写密码?
1673 浏览 2 评论
1936浏览 9评论
STM32仿真器是选择ST-LINK还是选择J-LINK?各有什么优势啊?
729浏览 4评论
STM32F0_TIM2输出pwm2后OLED变暗或者系统重启是怎么回事?
569浏览 3评论
594浏览 3评论
stm32cubemx生成mdk-arm v4项目文件无法打开是什么原因导致的?
552浏览 3评论
小黑屋| 手机版| Archiver| 电子发烧友 ( 湘ICP备2023018690号 )
GMT+8, 2024-12-22 21:36 , Processed in 0.788793 second(s), Total 79, Slave 62 queries .
Powered by 电子发烧友网
© 2015 bbs.elecfans.com
关注我们的微信
下载发烧友APP
电子发烧友观察
版权所有 © 湖南华秋数字科技有限公司
电子发烧友 (电路图) 湘公网安备 43011202000918 号 电信与信息服务业务经营许可证:合字B2-20210191 工商网监 湘ICP备2023018690号