波尔 发表于 2014-4-9 19:41

怎么把采集的振动数据加入到Duffing方程

   怎么把采集的振动数据加入到Duffing方程从而进行故障诊断？

牛小贱 发表于 2014-4-9 20:41

由于之前没接触过这个duffing方程 ，故从网上下载的源程序，可以参考一下：从外部导入数据s4（没有数据s4！！）。
程序代码：1.定义duffing方程：
function dx=duffing(t,x)
load s4
w=1;
l=s4 (round(t/12000+1));
dx=;2.求微分方程
=ode45(@duffing,,);
x1=x(:,1);
x2=x(:,2);
x1=mod(x1,2*pi);
x1(x1>pi)=x1(x1>pi)-2*pi;
plot(t(1:50:end),x1(1:50:end))%频闪采样图形
figure
h=plot(x1,x2)希望对你有所帮助！！

波尔 发表于 2014-4-9 20:51

牛小贱 发表于 2014-4-9 20:41
由于之前没接触过这个duffing方程 ，故从网上下载的源程序，可以参考一下：从外部导入数据s4（没有数据s4！ ...

定义Duffing方程中的第四行代码第一个字符是一还是L?

牛小贱 发表于 2014-4-9 21:21

波尔 发表于 2014-4-9 20:51
定义Duffing方程中的第四行代码第一个字符是一还是L?

肯定是L啊……这个程序不能运行啊，因为没有数据，只是给你一个参考。但是你可以把你的数据导入，修改一下程序即可。

波尔 发表于 2014-4-10 15:30

还有人给点意见吗

liguangzhigong 发表于 2014-9-2 21:26

请问，这个问题你解决了吗？我看这个程序还是有些问题

liguangzhigong 发表于 2014-9-2 21:26

04.l=s4 (round(t/12000+1));这一句不对吧

liguangzhigong 发表于 2014-9-2 21:38

是不是应该是l=s4 (round(t*12000+1));

liguangzhigong 发表于 2014-9-2 21:39

04.x1=mod(x1,2*pi);

05.x1(x1>pi)=x1(x1>pi)-2*pi;

06.plot(t(1:50:end),x1(1:50:end))%频闪采样图形

07.figure
这几个也没没用吧

Vickyvictoria 发表于 2015-10-30 21:03

liguangzhigong 发表于 2014-9-2 21:26
04.l=s4 (round(t/12000+1));这一句不对吧

这个应该没什么问题，主要看s4的含义

Vickyvictoria 发表于 2015-10-30 21:04

liguangzhigong 发表于 2014-9-2 21:39
04.x1=mod(x1,2*pi);

05.x1(x1>pi)=x1(x1>pi)-2*pi;

这个对于程序估计是有用的，不过这么处理可以看出该程序太不严格了

ljq519781936 发表于 2019-1-19 21:09

x1=mod(x1,2*pi);
x1(x1>pi)=x1(x1>pi)-2*pi;这个有啥用请教下？
页: [1]
查看完整版本: 怎么把采集的振动数据加入到Duffing方程