基于混沌集成决策树的电能质量复合扰动识别(matlab代码)

news2024/11/16 3:32:04

目录

1 主要内容

2 部分代码

3 程序结果

4 程序链接


主要内容

该程序参考《基于混沌集成决策树的电能质量复合扰动识别》,主要做的是S变换电能质量扰动识别,通过S变换对电能质量扰动(谐波,闪变,暂升等单一扰动和复合扰动)进行变换得到时频图,并对其进行特征提取。

部分代码

%% 电能质量扰动函数
V1=sin(w*x)+a_3*sin(3*w*x+b_3)+a_5*sin(5*w*x+b_5)+a_7*sin(7*w*x+b_7);
% plot(t,V1);    title('1:谐波'); 谐波
V2=(1-a2.*(stepfun(x,0.08)-stepfun(x,0.16))).*sin(w*x);
%   plot(t,V2)  电压暂降   0.1-0.9
V3=(1+a2.*(stepfun(x,0.08)-stepfun(x,0.16))).*sin(w*x);
  % 电压暂升    
V4=(1-a4.*(stepfun(x,0.08)-stepfun(x,0.16))).*sin(w*x);
%   plot(t,V4)      %0.9-1  电压中断
V5=(1+a_f*sin(b*w*x)).*sin(w*x);
%subplot(2,1,2);  plot(t,V5);   title('5:闪变');
V6=sin(w*x)+a6*exp(-(x-0.02)/tao).*sin(Wn*w*x).*(stepfun(x,0.02)-stepfun(x,0.03));
%plot(t,V6)   暂态振荡
V7=sin(w*x)+a7*exp(-(x-0.02)/tao).*(stepfun(x,0.02)-stepfun(x,0.06)); 
%plot(t,V7)   暂态脉冲
​
V8=(1-a2.*(stepfun(x,0.18)-stepfun(x,0.26))).*sin(w*x)+a_3*sin(3*w*x+b_3)+a_5*sin(5*w*x+b_5)+a_7*sin(7*w*x+b_7);
%电压暂降+谐波
V9=(1+a2.*(stepfun(x,0.18)-stepfun(x,0.26))).*sin(w*x)+a_3*sin(3*w*x+b_3)+a_5*sin(5*w*x+b_5)+a_7*sin(7*w*x+b_7);
%电压暂升+谐波
V10=(1-a4.*(stepfun(x,0.08)-stepfun(x,0.16))).*sin(w*x)+a_3*sin(3*w*x+b_3)+a_5*sin(5*w*x+b_5)+a_7*sin(7*w*x+b_7);
%电压中断+谐波
V11=(1+a_f*sin(b*w*x)).*sin(w*x)+a_3*sin(3*w*x+b_3)+a_5*sin(5*w*x+b_5)+a_7*sin(7*w*x+b_7);
%电压闪变+谐波
​
V12=(1-a2.*(stepfun(x,0.08)-stepfun(x,0.16))).*sin(w*x)+a6*exp(-(x-0.02)/tao).*sin(Wn*w*x).*(stepfun(x,0.02)-stepfun(x,0.03));
%电压暂降+暂态振荡
V13=(1+a2.*(stepfun(x,0.08)-stepfun(x,0.16))).*sin(w*x)+a6*exp(-(x-0.02)/tao).*sin(Wn*w*x).*(stepfun(x,0.02)-stepfun(x,0.03));
%电压暂升+暂态振荡
V14=(1+a_f*sin(b*w*x)).*sin(w*x)+a6*exp(-(x-0.02)/tao).*sin(Wn*w*x).*(stepfun(x,0.02)-stepfun(x,0.03));
%电压闪变+暂态振荡
V15=sin(w*x)+a_3*sin(3*w*x+b_3)+a_5*sin(5*w*x+b_5)+a_7*sin(7*w*x+b_7)+a6*exp(-(x-0.02)/tao).*sin(Wn*w*x).*(stepfun(x,0.02)-stepfun(x,0.03));
%谐波+暂态振荡
​
V16=(1-a2.*(stepfun(x,0.08)-stepfun(x,0.16))).*sin(w*x)+a7*exp(-(x-0.02)/tao).*(stepfun(x,0.02)-stepfun(x,0.06)); 
%电压暂降+暂态脉冲
V17=(1+a2.*(stepfun(x,0.08)-stepfun(x,0.16))).*sin(w*x)+a7*exp(-(x-0.02)/tao).*(stepfun(x,0.02)-stepfun(x,0.06)); 
%电压暂升+暂态脉冲
V18=(1+a_f*sin(b*w*x)).*sin(w*x)+a7*exp(-(x-0.02)/tao).*(stepfun(x,0.02)-stepfun(x,0.06)); 
%闪变+暂态脉冲
V19=sin(w*x)+a_3*sin(3*w*x+b_3)+a_5*sin(5*w*x+b_5)+a_7*sin(7*w*x+b_7)+a7*exp(-(x-0.02)/tao).*(stepfun(x,0.02)-stepfun(x,0.06)); 
%谐波+暂态脉冲
​
V20=(1-a2.*(stepfun(x,0.08)-stepfun(x,0.16))).*sin(w*x)+a_3*sin(3*w*x+b_3)+a_5*sin(5*w*x+b_5)+a_7*sin(7*w*x+b_7)+a6*exp(-(x-0.02)/tao).*sin(Wn*w*x).*(stepfun(x,0.02)-stepfun(x,0.03));
%暂降+谐波+暂态振荡
V21=(1+a2.*(stepfun(x,0.08)-stepfun(x,0.16))).*sin(w*x)+a_3*sin(3*w*x+b_3)+a_5*sin(5*w*x+b_5)+a_7*sin(7*w*x+b_7)+a6*exp(-(x-0.02)/tao).*sin(Wn*w*x).*(stepfun(x,0.02)-stepfun(x,0.03));
%暂升+谐波+暂态振荡
V22=(1+a_f*sin(b*w*x)).*sin(w*x)+a_3*sin(3*w*x+b_3)+a_5*sin(5*w*x+b_5)+a_7*sin(7*w*x+b_7)+a7*exp(-(x-0.02)/tao).*(stepfun(x,0.02)-stepfun(x,0.06));
%闪变+谐波+暂态脉冲
V23=(1-a2.*(stepfun(x,0.08)-stepfun(x,0.16))).*sin(w*x)+a_3*sin(3*w*x+b_3)+a_5*sin(5*w*x+b_5)+a_7*sin(7*w*x+b_7)+a6*exp(-(x-0.02)/tao).*sin(Wn*w*x).*(stepfun(x,0.02)-stepfun(x,0.03))+a7*exp(-(x-0.02)/tao).*(stepfun(x,0.02)-stepfun(x,0.06));
%暂降+谐波+暂态振荡+暂态脉冲
%plot(x,V20)
​
%%  V1 谐波  S变换  做时频图

程序结果

4 程序链接

 点击直达!

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/728943.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

15 | 边界:微服务的各种边界在架构演进中的作用?

目录 演进式架构 微服务还是小单体? 微服务边界的作用 逻辑边界 物理边界 代码边界 正确理解微服务的边界 总结 那重点落到边界的时候,总结一下就是,微服务的设计要涉及到逻辑边界、物理边界和代码边界等等。 那么这些边界在微服务架…

vue清除地址栏参数(可以单个,可以多个)

需求 vue跳转新页面后,清除url里面的参数,就是上一个页面带过来的参数,只用一次 mounted () {this.$nextTick(() > {let url this.getnewurl()window.history.replaceState(null, null, url);})},methods: { //根据参数名去清除&#xf…

深度学习训练营之调用Gensim来训练Word2Vec模型

深度学习训练营之调用Gensim来训练Word2Vec模型 原文链接环境介绍前置工作下载Gensim库对于原始语料进行分词添加停用词 模型训练模型介绍模型正式训练 计算词频 原文链接 🍨 本文为🔗365天深度学习训练营 中的学习记录博客🍦 参考文章&#…

C#:了解LINQ,简化数据查询和操作的强大工具

文章目录 linq关键字fromwhereselectorderbyjoingroupletinto linq方法筛选方法WhereOfType 排序方法:OrderByOrderByDescendingThenByThenByDescending 投影方法:SelectSelectMany 分组方法:GroupBy 连接方法:JoinGroupJoin 聚合…

高效方案:30万条数据插入 MySQL 仅需13秒

本文主要讲述通过MyBatis、JDBC等做大数据量数据插入的案例和结果。 30万条数据插入插入数据库验证 实体类、mapper和配置文件定义 User实体 mapper接口 mapper.xml文件 jdbc.properties sqlMapConfig.xml 不分批次直接梭哈 循环逐条插入 MyBatis实现插入30万条数据 J…

SpringBoot 集成 MybatisPlus 一——介绍

MybatisPlus 是 Mybatis 的升级版本,是对 Mybatis 的简化,因为他们的口号就是“为简化开发而生”。 1、创建数据表 CREATE TABLE ​​User​​ ( ​​id​​ INT NOT NULL, ​​username​​ VARCHAR(50) NULL DEFAULT NULL, ​​gendar​​ CHAR(2) NU…

了解微服务架构

微服务架构 软件架构历史 Software As A Service,这不仅仅是⼀个理念,它更多释放的是企业在新⼀轮的市场竞争中,如何使⽤轻量级的组织架构和新的软件架构来更好的服务企业⾯向未来的战略调整和市场定位,从⽽赢得未来的市场空间。…

【计算机视觉 | 图像分割】arxiv 计算机视觉关于图像分割的学术速递(7 月 6 日论文合集)

文章目录 一、分割|语义相关(15篇)1.1 Prompting Diffusion Representations for Cross-Domain Semantic Segmentation1.2 ZJU ReLER Submission for EPIC-KITCHEN Challenge 2023: Semi-Supervised Video Object Segmentation1.3 Multi-Modal Prototypes for Open-Set Semanti…

基于YOLO的3D人脸关键点检测方案

目录 前言一、任务列表二、3D人脸关键点数据H3WB2.下载方法3.任务4.评估5.使用许可 3DFAWAFLW2000-3D 三、3D关键点的Z维度信息1.基于3DMM模型的方法2.H3WB 四、当前SOTA的方法1.方法1 五、我们的解决方法1.数据转为YOLO格式2.修改YOLO8Pose的入口出口3.开始训练,并…

aardio 的addHeaders请求

aardio群 625494397 废话不多说 直接开干! import web.json; import console; import inet.whttp; web.json.parse( json );//转化json格式 h.addHeaders {Accept-Encoding gzip, deflate, br;Accept-Language zh-CN,zh;q0.9;User-Agent Mozilla/5.0 (Windows …

浅谈智能照明控制管理系统的功能介绍

安科瑞电气股份有限公司 上海嘉定 201801 摘要:智能照明控制系统较好地实现了智能控制、人性化照明和节能降耗的功能,使其在楼宇控制领域变得越来越重要,越来越受到人们的重视。本文介绍了智能照明控制系统的概念、特点、优势、发展方向等内容,并着重对智能照明控制…

多个项目的进度管理:掌握这3点,轻松应对

在现代的企业竞争中,许多企业组织都会同时进行多个项目。而有效的进度管理对于确保项目按时交付和优化资源利用至关重要。那么有哪些方法可以帮助项目管理者掌握多个项目的进度管理呢? 一:明确目标和优先级 在管理多个项目时,明确…

windows系统使用nvm实现多版本node切换

介绍nvm 是 node version manager(node 版本管理工具)的缩写,是一个命令行工具,用于管理和切换到不同版本的 node.js。 不同的项目可能需要不同版本的 node.js 和 npm(node 包管理器),例如&…

增强型视觉系统 (EVS)

增强型视觉系统 EVS 1、增强型视觉系统概览2、车载相机 HAL2.1 EVS 应用2.2 EVS 管理器2.3 EVS HIDL 接口2.4 内核驱动程序 《增强型视觉系统 (EVS) 1.1 集成指南》 车载相机 HAL 1、增强型视觉系统概览 为了增强视频串流管理和错误处理,Android 11 更新了车载相机…

OD查找偏移表达式

1、在CE中先找到血地址2DF94C28 2、在od数据段dd 血地址 dd 2DF94C28查找到数据 3、在血地址位置右键-断点-硬件访问-DWORD,再操作游戏会再改变这个地址代码段下断点,查找到 4、验证查找的是否正确,在8199FA处下断点,dd edi288处…

github port 22 connection refused 不通过科学上网连接不上问题

原因是用密钥连接github走的是ssh,ssh默认端口是22没有加密,github不让走22 解决办法是在 ~/.ssh/config里添加 Host github.com Hostname ssh.github.com Port 443 让ssh连github的时候走443端口,TLS加密的 ssh -T gitgithub.com 说明已…

Docker深入

一、简介 Docker是一个用于构建、运行、传送应用程序的平台。 2、为什么使用Docker 在部署服务器过程中,配置各种环境变量、第三方依赖等等,耗费时间太长,使用Docker可以将他们打包成一个集装箱,只要在开发环境中运行成功&…

Hive 复习重点

文章目录 特点SQL查询转换成MR作业的过程内部表 & 外部表 & 分区表 & 分桶表内部表外部表分区表 分桶表DQL语法性能优化SQL语句优化数据格式化 ORC, Parquet列式存储小文件过多优化并行执行优化JVM 优化推测执行行列过滤limit 限制本地模式(开启MR,单机…

stm32mp157——通过按键中断实现LED灯的亮灭

EXTI章节初始化 void hal_exit_init(EXTI编号,GPIOF组编号,中断触发方式) 调用:hal_exit_init(EXTI9,0x05,下降沿触发) 调用: hal_exit_init(EXTI7,0x05,下降沿触发) 调用: hal_exit_init(EXTI8,0x05,下降沿触发) GIC章节初始化 void hal_gic_init(中断号,中断优先级) 调用: ha…

Redis - 一篇讲解根据 Key 前缀统计分析内存占用

问题描述 今天遇到一个 Redis 内存打挂了的问题,想看看哪个前缀 Key 占用内存比较大?! 原因分析 我们都知道如果直接用 Keys 参数去做统计很危险,而且也只能统计数量,当然也可以排序去前几名的占用内存 Key 对应的大…