西莫电机圈

 找回密码
 立即注册

QQ登录

只需一步,快速开始

手机号码,快捷登录

手机号码,快捷登录

查看: 1593|回复: 1

[求助] 求助从一个叠加了扰动量的正弦信号中提取出扰动量的MATLAB程序

[复制链接]

该用户从未签到

发表于 2012-5-17 12:55 | 显示全部楼层 |阅读模式

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

您需要 登录 才可以下载或查看,没有账号?立即注册

x
本帖最后由 血色黄昏 于 2012-5-17 13:07 编辑

调速 扰动时间20-21s 相电压 相电流瞬时值.png     部分程序如下,我看的不是很懂。具体情况是一个电力系统模型稳定运行后,20s在发电机侧添加一个振荡扰动,时长6s,扰动消失后系统重新稳定。在workspace中保存了线路中一个节点采样得到的电压值,如u1,采样点为60000个,现在要提取u1中在添加扰动后叠加的扰动量,即通过u1电压值分离得到扰动量的波形。右图中上部分为u1,稳定时为一个正弦信号。有没有哪位高手能帮忙看看这个问题,优化下程序,万分感激,先谢谢了。
















plot (u1(1:60000,1), 'DisplayName', 'u1(1:60000,1)', 'YDataSource', 'u1(1:60000,1)'); figure(gcf)
a1=hilbert(u1)
hu1=abs(a1)
plot (hu1(1:60000,1), 'DisplayName', 'hu1(1:60000,1)', 'YDataSource', 'hu1(1:60000,1)'); figure(gcf)
a2=hilbert(u2)
hu2=abs(a2)
plot (hu2(1:60000,1), 'DisplayName', 'hu2(1:60000,1)', 'YDataSource', 'hu2(1:60000,1)'); figure(gcf)
a3=hilbert(u3)
hu3=abs(a3)
plot (hu3(1:60000,1), 'DisplayName', 'hu3(1:60000,1)', 'YDataSource', 'hu3(1:60000,1)'); figure(gcf)
plot (hu2(1:60000,1), 'DisplayName', 'hu2(1:60000,1)', 'YDataSource', 'hu2(1:60000,1)'); figure(gcf)
plot (hu1(1:60000,1), 'DisplayName', 'hu1(1:60000,1)', 'YDataSource', 'hu1(1:60000,1)'); figure(gcf)
hju1=hu1(30000:50000)
mm1=max(hju1)
mi1=min(hju1)
hjmu1=hju1-(mm1-mi1)/2-mi1
plot (hjmu1(1:20001,1), 'DisplayName', 'hjmu1(1:20001,1)', 'YDataSource', 'hjmu1(1:20001,1)'); figure(gcf)
hju3=hu3(30000:50000)
mm3=max(hju3)
mi3=min(hju3)
hjmu3=hju3-(mm3-mi3)/2-mi3
plot (hjmu3(1:20001,1), 'DisplayName', 'hjmu3(1:20001,1)', 'YDataSource', 'hjmu3(1:20001,1)'); figure(gcf)
plot (hjmu1(1:20001,1), 'DisplayName', 'hjmu1(1:20001,1)', 'YDataSource', 'hjmu1(1:20001,1)'); figure(gcf)
a=xcorr(hjmu3,hjmu1)
[m13,t13]=max(a)
plot (hjmu1(1:20001,1), 'DisplayName', 'hjmu1(1:20001,1)', 'YDataSource', 'hjmu1(1:20001,1)'); figure(gcf)
hjmju1=hjmu1(2500:12500)
plot (hjmju1(1:10001,1), 'DisplayName', 'hjmju1(1:10001,1)', 'YDataSource', 'hjmju1(1:10001,1)'); figure(gcf)
hjmju3=hjmu3(2500:12500)
plot (hjmju3(1:10001,1), 'DisplayName', 'hjmju3(1:10001,1)', 'YDataSource', 'hjmju3(1:10001,1)'); figure(gcf)
plot (hjmju1(1:10001,1), 'DisplayName', 'hjmju1(1:10001,1)', 'YDataSource', 'hjmju1(1:10001,1)'); figure(gcf)
a=xcorr(hjmju1,hjmju3)
[m13,t13]=max(a)
plot (a(1:20001,1), 'DisplayName', 'a(1:20001,1)', 'YDataSource', 'a(1:20001,1)'); figure(gcf)
r13=a
plot (r13(1:20001,1), 'DisplayName', 'r13(1:20001,1)', 'YDataSource', 'r13(1:20001,1)'); figure(gcf)
rr13=r13/100000
plot (rr13(1:20001,1), 'DisplayName', 'rr13(1:20001,1)', 'YDataSource', 'rr13(1:20001,1)'); figure(gcf)
西莫电机论坛微信公众平台正式上线!★详情请点击★ 西莫电机论坛会员交流专用群欢迎您西莫电机论坛加群请注明论坛用户名及所从事专业,否则不予通过

签到天数: 7 天

连续签到: 1 天

[LV.3]偶尔看看II

发表于 2012-6-2 22:16 | 显示全部楼层
没看懂,不过我看思路好像是有傅立叶分析。
西莫电机论坛微信公众平台正式上线!★详情请点击★ 西莫电机论坛会员交流专用群欢迎您西莫电机论坛加群请注明论坛用户名及所从事专业,否则不予通过
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

QQ|小黑屋|手机版|Archiver|西莫电机圈 ( 浙ICP备10025899号-3 浙公网安备:33028202000436号

GMT+8, 2024-11-5 17:23 , Processed in 0.088727 second(s), 32 queries .

Powered by Discuz! X3.4

© 2001-2023 Discuz! Team.

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