当前位置: 首页 > news >正文

Simulink中Matlab function使用全局变量

目录

  • 一. 引言
  • 二. 普通Matlab function使用全局变量
  • 三. Simulink中的Matlab function使用全局变量
  • 四. 如何利用Matlab function的全局变量施加随机噪声

一. 引言

最近发现了之前仿真中的一个问题,记录一下备忘。

Matlab function中有时候需要用到全局变量,其用法较为简单,而Simulink中的Matlab function全局变量用法稍复杂一点,按照步骤来即可,但第一次用还是需要一个详细教程,以免出错。

二. 普通Matlab function使用全局变量

普通Matlab function全局变量定义:

function W = weight(m)
global GG;
W = m * GG;

将上述代码保存成文件weight.m.

注意调用的时候也要声明全局变量:

global GG;
GG = 10;
a = weight(4)

三. Simulink中的Matlab function使用全局变量

详细的多种用法参考这篇文章: Simulink高级应用—— Matlab Function中共享或使用全局变量或全局数据

这里只说其中一种用法,使用data sotre memory来实现全局变量共享。

比如我这里要给控制器输出端加一个噪音信号,用Matlab function,就可以定义为时间的函数。

在这里插入图片描述
1. MATLAB Function的代码中 声明全局变量

考虑最简单的情况,t=5s时加一个幅值为noise_amp的噪音,定义:

function y = fcn(t)
global noise_amp;
if t > 5y = noise_amp;
else y = 0;
end

2. MATLAB Function 注册全局变量

在MATLAB Function的代码编辑界面,点击Edit Data

点击

点击Add Data,成功添加一个变量,命名为noise_amp(与代码中的必须一致)。变量类型由input改为data store memory

在这里插入图片描述
3. Simulink中添加 Data Store Memory模块

从library中拖入或直接在空白页面输入Data Store Memory,添加一个Data Store Memory模块。命名与Matlab function中一致,必须设置好明确的数据类型,不能为auto(比如我设置成double)。设置信号的类型(real)和初始值。我这里设置的初始值0.3,即噪音幅值0.3

在这里插入图片描述

现在运行一下,结果如图所示。可以看到t=5s时直接在输出端叠加了一个0.3的噪音,之后在PID的调节下回到设置的期望值1。这样就仿真了一个脉冲形式的抗扰控制。

在这里插入图片描述

四. 如何利用Matlab function的全局变量施加随机噪声

上面是在某一时刻施加了一个噪音,现在我想以10Hz频率给输出端施加随机噪音(控制周期为0.001s, 1000Hz),应该怎么施加。

这个问题乍一看很简单,其实仔细想一想不是那么简单。比如下面这这么写就是错的:

function y = fcn(t,u)if mod(1000*t,100) == 0   % 噪声频率为0.1s给一次, 仿真周期是0.001sy = 0.1 * (-1+2*rand);
elsey = 0;
end

它仿真出来的是下图形式,每一次的噪音只在那一时刻起了作用,表现在输出信号上就是一个个的“尖刺”,实际上每一次施加的噪音都应该是持续性的,所以应该用累加形式。

在这里插入图片描述

正确的写法如下(新定义了一个全局变量noise,步骤参照前面):

function y = fcn(t, data1)global noise;if mod(1000*t,100) == 0   % 噪声频率为0.1s给一次, 仿真周期是0.001sy = noise + 0.1 * (-1+2*rand);
elsey = noise;
end
noise = y;

这样子输出就是下面这样的,可以看到是连续信号(只演示原理,控制效果忽略~)。

在这里插入图片描述

http://www.lryc.cn/news/483626.html

相关文章:

  • WPF-控件的属性值的类型转化
  • 海思Hi3516DV300上播放G711U音频文件
  • Linux源码阅读笔记-V4L2框架基础介绍
  • 列表(list)
  • 使用Python抓取数据的实战指南
  • GIC寄存器介绍
  • c++实现B树(下)
  • 外星人入侵
  • 【数据仓库】hbase的安装与简单操作
  • 为什么RNN(循环神经网络)存在梯度消失和梯度爆炸?
  • 【数据库】数据库迁移的注意事项有哪些?
  • MQTT协议解析 : 物联网领域的最佳选择
  • pycharm中from[本地包]import文件/模块出现问题(最最最全方法!)
  • MongoDB在现代Web开发中的应用
  • Python Bokeh 数据可视化教程
  • (一)<江科大STM32>——软件环境搭建+新建工程步骤
  • 内存大小的单位转换
  • 如何在 Spring MVC 中使用 `@PostMapping`? 如何在 Spring MVC 中使用 `@PutMapping`?
  • AIGC Agent(智能体)应用开发高级工程师实战培训 —— 线上8周系统教学课程学习路线图
  • GDSC、CTRP数据库学习
  • 【嵌入式】ESP32开发(一)ESP-IDF概述
  • 最新6.7分非肿瘤纯生信,使用机器学习筛选慢阻肺中的关键基因。机器学习在非肿瘤生信文章中正火,可重复!
  • vue 提交表单抹除字段为空的数据
  • web实验3:虚拟主机基于不同端口、目录、IP、域名访问不同页面
  • 英伟达Isaac Manipulator产品体验
  • 网安加·百家讲坛 | 仝辉:金融机构鸿蒙应用安全合规建设方案
  • PHP Session
  • 泷羽sec学习打卡-Linux基础2
  • # 【STM32F1】——无线收发模块RF200与串口通信
  • 计算机网络:运输层 —— TCP 协议概述与 TCP 报文段首部格式