【Simulink系列】——动态系统仿真 之 离散系统线性离散系统

2024-02-04 18:36

本文主要是介绍【Simulink系列】——动态系统仿真 之 离散系统线性离散系统,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

一、离散系统定义

离散系统是指系统的输入与输出仅在离散的时间上取值,而且离散的时间具有相同的时间间隔。满足下列条件:

①系统(的输入输出)每隔固定时间间隔才更新一次。固定时间间隔称为采样时间。

②系统的输出依赖当前的输入、以往的输入和输出见数学描述

③系统具有离散的状态。状态是指系统前一时刻的输出量。

设系统的输入变量为u\left(n T_s\right), n=0,1,2, \cdots,其中n为采样时刻,Ts为系统的采样时间。由于Ts为固定值,因而系统输入u\left(n T_s\right)常被记为u(n)。设系统输出为y(nTs),也可记为y(n),则离散系统的数学描述为:

y(n)=f(u(n), u(n-1), \cdots, y(n-1), y(n-2), \cdots)

注:离散系统初始状态的确定非常重要。

二、离散Simulinlk仿真

有一个人口变化系统。设某一年的人口数目为p(n),其中n表示年份,它与上一年的人口数目p(n-1)、人口繁殖速率r以及新增资源所能满足的个体数目K之间的动力学模型如下差分方程:

p(n)=r p(n-1)\left[1-\frac{p(n-1)}{K}\right]

此差分方程可看出系统为非线性离散系统。设人口初始值p(0)=600000,人口繁殖速率r为1.08,新增资源满足个体数K=200000,建立Simulink模型,分析人口数目的变化趋势。

(1)建立系统模型

Unit Delay为单位延迟模块,主要功能是将输入信号延迟一个采样时间,这里将p(n)作为它的输入以得到p(n-1)。

(2)模块参数设置

人口繁殖速率r:取1.08,
1/K:取1/2000000,
Unit Delay的初始值为初始人口数目:取600000。

(3)仿真参数设置

仿真时间0~50s,求解器设置为discrete离散,即需要使用离散求解器;仿真步长选择Fixed-step定步长。

(4)仿真运行

注:若离散系统中所有模块的采样时间均相同,则此系统为单速率离散系统(选定步长求解器),否则为多速率离散系统(选变步长仿真更好)。

三、线性离散系统定义

在离散系统中,线性离散系统具有重要的地位。当离散系统同时满足齐次性和叠加性时,即:

T\left\{\alpha u_{1}(n)+\beta u_{2}\left(n\right)\right\}=\alpha T\left\{u_{1}\left(n\right)\right\}+\beta T\left\{u_{2}\left(n\right)\right\}

此离散系统称为线性离散系统。最一般的数学描述为:

y(n)=f\bigl(u(n),u(n-1),\cdots;y(n-1),y(n-2),\cdots\bigr)

也可用差分方程描述(上为状态方程,下为输出方程):

\begin{array}{rl}&x(n+1)=f\big(x(n),u(n),n\big)\\&y(n)=g\big(x(n),u(n),n\big)\end{array}

除了上述两种方式,针对线性离散系统本身的特点,常使用Z变换来描述线性离散系统。Z变换时对离散信号进行分析的一个强有力工具,尤其线性离散系统

对于一个离散信号u(n),其Z变换为U\bigl(z\bigr)=\sum_{k=-\infty}^{\infty}u\bigl(k\bigr)z^{-k},一般离散系统的时间大于0,则Z变换为U\bigl(z\bigr)=\sum_{k=0}^{\infty}u\bigl(k\bigr)z^{-k},简记为U(z)=Z\left\{u(n)\right\}

利用离散信号的Z变换结果得到原离散信号的过程,叫做Z变换的逆变换,简记为u\left(z\right)=Z^{-1}\left\{U(n)\right\}

Z变换有两个重要性质

①Z变换同时满足齐次性和叠加性。

Z\left\{\alpha u_{1}\left(n\right)+\beta u_{2}\left(n\right)\right\}=\alpha Z\left\{u_{1}\left(n\right)\right\}+\beta Z\left\{u_{2}\left(n\right)\right\}

②若u(n)的Z变换为U(z),则u(n-1)的Z变换为z^{-1}U(z)

四、线性离散系统Simulink仿真

低通滤波器模型。它可以滤除信号中的高频部分,以获取信号中有用的低频信号,一个低通数字滤波器的差分方程如下:

y(n)-1.6y(n-1)+0.7y(n-2)=0.04u(n)+0.08u(n-1)+0.04u(n-2)

其中u(n)为滤波器输入,y(n)为滤波器输出。可见其为线性离散系统。

(1)Z变换

对该滤波器系统的差分方程进行Z变换,可得下面描述:

\frac{Y(z)}{U(z)}=\frac{0.04+0.08z^{-1}+0.04z^{-2}}{1-1.6z^{-1}+0.7z^{-2}}

(2)建立系统模型

发送方首先使用高频正弦波(载波)对一低频锯齿波(原信号波)进行幅度调制,然后在无损信道中传递此幅度调制信号;接受方在接收到幅度调制信号后,首先进行解调,然后使用低通数字滤波器对解调后的信号进行滤波以获得低频锯齿波信号

Sine Wave用来产生高频载波信号,Sine Wave1用来产生解调信号;

Signal Generator产生低频锯齿波信号;

Discrete Filter表示数字滤波器。

(3)模块参数设置

高频载波(Sine Wave)频率设置为1000rad/sec,其余默认;Signal Generator波形选择锯齿波sawtooth,其余默认;解调正弦波(Sine Wave1)频率也为1000rad/sec,采样时间0.005s。

Discrete Filter设置如下:

注:数字滤波器的采样时间应该与解调信号的采样时间一致。

(4)仿真参数设置

仿真时间0~10s,选择求解器并设置为Variable-step变步长ode45(Dormand-Prince),最大步长设置为0.005

(5)仿真运行

数字滤波器的输出信号与原始锯齿信号并不完全一致,存在一定失真,这种失真不可避免,因为实际上并不存在理想的滤波器,调制时信号之间不可避免会相互干扰。

疑问:博友们,这里遇到各问题,锯齿波的scope显示周期与设定的信号周期不一致???

往期精彩

STM32专栏(付费9.9)icon-default.png?t=N7T8http://t.csdnimg.cn/E2F88

OpenCV-Python专栏(付费9.9)icon-default.png?t=N7T8http://t.csdnimg.cn/zK1jV

AI底层逻辑专栏(付费9.9)icon-default.png?t=N7T8http://t.csdnimg.cn/zic0f

机器学习专栏(免费)icon-default.png?t=N7T8http://t.csdnimg.cn/FaXzAFreeRTOS专栏(免费)icon-default.png?t=N7T8http://t.csdnimg.cn/SjIqU电机控制专栏(免费)icon-default.png?t=N7T8http://t.csdnimg.cn/FNWM7

这篇关于【Simulink系列】——动态系统仿真 之 离散系统线性离散系统的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/678442

相关文章

Spring Security 从入门到进阶系列教程

Spring Security 入门系列 《保护 Web 应用的安全》 《Spring-Security-入门(一):登录与退出》 《Spring-Security-入门(二):基于数据库验证》 《Spring-Security-入门(三):密码加密》 《Spring-Security-入门(四):自定义-Filter》 《Spring-Security-入门(五):在 Sprin

不懂推荐算法也能设计推荐系统

本文以商业化应用推荐为例,告诉我们不懂推荐算法的产品,也能从产品侧出发, 设计出一款不错的推荐系统。 相信很多新手产品,看到算法二字,多是懵圈的。 什么排序算法、最短路径等都是相对传统的算法(注:传统是指科班出身的产品都会接触过)。但对于推荐算法,多数产品对着网上搜到的资源,都会无从下手。特别当某些推荐算法 和 “AI”扯上关系后,更是加大了理解的难度。 但,不了解推荐算法,就无法做推荐系

基于人工智能的图像分类系统

目录 引言项目背景环境准备 硬件要求软件安装与配置系统设计 系统架构关键技术代码示例 数据预处理模型训练模型预测应用场景结论 1. 引言 图像分类是计算机视觉中的一个重要任务,目标是自动识别图像中的对象类别。通过卷积神经网络(CNN)等深度学习技术,我们可以构建高效的图像分类系统,广泛应用于自动驾驶、医疗影像诊断、监控分析等领域。本文将介绍如何构建一个基于人工智能的图像分类系统,包括环境

水位雨量在线监测系统概述及应用介绍

在当今社会,随着科技的飞速发展,各种智能监测系统已成为保障公共安全、促进资源管理和环境保护的重要工具。其中,水位雨量在线监测系统作为自然灾害预警、水资源管理及水利工程运行的关键技术,其重要性不言而喻。 一、水位雨量在线监测系统的基本原理 水位雨量在线监测系统主要由数据采集单元、数据传输网络、数据处理中心及用户终端四大部分构成,形成了一个完整的闭环系统。 数据采集单元:这是系统的“眼睛”,

第10章 中断和动态时钟显示

第10章 中断和动态时钟显示 从本章开始,按照书籍的划分,第10章开始就进入保护模式(Protected Mode)部分了,感觉从这里开始难度突然就增加了。 书中介绍了为什么有中断(Interrupt)的设计,中断的几种方式:外部硬件中断、内部中断和软中断。通过中断做了一个会走的时钟和屏幕上输入字符的程序。 我自己理解中断的一些作用: 为了更好的利用处理器的性能。协同快速和慢速设备一起工作

嵌入式QT开发:构建高效智能的嵌入式系统

摘要: 本文深入探讨了嵌入式 QT 相关的各个方面。从 QT 框架的基础架构和核心概念出发,详细阐述了其在嵌入式环境中的优势与特点。文中分析了嵌入式 QT 的开发环境搭建过程,包括交叉编译工具链的配置等关键步骤。进一步探讨了嵌入式 QT 的界面设计与开发,涵盖了从基本控件的使用到复杂界面布局的构建。同时也深入研究了信号与槽机制在嵌入式系统中的应用,以及嵌入式 QT 与硬件设备的交互,包括输入输出设

JAVA智听未来一站式有声阅读平台听书系统小程序源码

智听未来,一站式有声阅读平台听书系统 🌟 开篇:遇见未来,从“智听”开始 在这个快节奏的时代,你是否渴望在忙碌的间隙,找到一片属于自己的宁静角落?是否梦想着能随时随地,沉浸在知识的海洋,或是故事的奇幻世界里?今天,就让我带你一起探索“智听未来”——这一站式有声阅读平台听书系统,它正悄悄改变着我们的阅读方式,让未来触手可及! 📚 第一站:海量资源,应有尽有 走进“智听

动态规划---打家劫舍

题目: 你是一个专业的小偷,计划偷窃沿街的房屋。每间房内都藏有一定的现金,影响你偷窃的唯一制约因素就是相邻的房屋装有相互连通的防盗系统,如果两间相邻的房屋在同一晚上被小偷闯入,系统会自动报警。 给定一个代表每个房屋存放金额的非负整数数组,计算你 不触动警报装置的情况下 ,一夜之内能够偷窃到的最高金额。 思路: 动态规划五部曲: 1.确定dp数组及含义 dp数组是一维数组,dp[i]代表

科研绘图系列:R语言扩展物种堆积图(Extended Stacked Barplot)

介绍 R语言的扩展物种堆积图是一种数据可视化工具,它不仅展示了物种的堆积结果,还整合了不同样本分组之间的差异性分析结果。这种图形表示方法能够直观地比较不同物种在各个分组中的显著性差异,为研究者提供了一种有效的数据解读方式。 加载R包 knitr::opts_chunk$set(warning = F, message = F)library(tidyverse)library(phyl

【区块链 + 人才服务】可信教育区块链治理系统 | FISCO BCOS应用案例

伴随着区块链技术的不断完善,其在教育信息化中的应用也在持续发展。利用区块链数据共识、不可篡改的特性, 将与教育相关的数据要素在区块链上进行存证确权,在确保数据可信的前提下,促进教育的公平、透明、开放,为教育教学质量提升赋能,实现教育数据的安全共享、高等教育体系的智慧治理。 可信教育区块链治理系统的顶层治理架构由教育部、高校、企业、学生等多方角色共同参与建设、维护,支撑教育资源共享、教学质量评估、