《2022年随机载荷模拟及雨流计数方法在matlab中的实现 .pdf》由会员分享，可在线阅读，更多相关《2022年随机载荷模拟及雨流计数方法在matlab中的实现 .pdf（6页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、近日来的工作一、仿真随机过程sw=2; A=sqrt(2*sw*detaw)（1）中心频率为10pi，detaw=0.01，带宽为 2pi 的低频窄带随机过程for n=9:0.01:11; xn=A*sin(n*pi*t+rand(1,1)*2*pi); x=x+xn; end （2）中心频率为 100pi，detaw=0.01，带宽为 2pi的高频窄带随机过程for n=99:0.01:101; yn=A*sin(n*pi*t+rand(1,1)*2*pi); y=y+yn;end（3）低频与高频组合成的宽带随机过程z=x+y 0102030405060708090100-10-50510
2、低 频 窄 带 随 机 过 程0102030405060708090100-10-50510高 频 窄 带 随 机 过 程0102030405060708090100-10-50510组 合 成 的 宽 带 随 机 过 程问题一：关于带宽。带宽为2pi 似乎太宽了，若带宽改为0.2pi，detaw 依然不变，得到图如下。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 1 页，共 6 页 - - - - - - - - - 0102030405060708090100-1-0.500.
3、51低 频 窄 带 随 机 过 程0102030405060708090100-0.4-0.200.20.4高 频 窄 带 随 机 过 程0102030405060708090100-1-0.500.51组 合 成 的 宽 带 随 机 过 程二、用雨流计数法截取循环（1）得到随机过程x的极点值（储存在矩阵S中） 。n=length(x);s=x(1);for i=2:n-1;if (x(i)x(i-1)&(x(i)x(i+1)|(x(i)x(i-1)&(x(i)x(i+1); s=s,x(i);endendS=s,x(n); （2）对极值点间变程进行比较分析，截取循环while i+3leng
4、th(S) i=i+1; detaS1=abs(S(i+1)-S(i); %获得四个极值点间的三段变程 detaS2=abs(S(i+2)-S(i+1); detaS3=abs(S(i+3)-S(i+2);if (detaS2=detaS1)&(detaS2高频窄带 低频窄带（自己这么想的，不知是否行的通）三、对剩余载荷的处理re=ZSr,ZSr;名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 4 页，共 6 页 - - - - - - - - - reside,Ba,Bm,Q=Ra
5、inflow(re);050100150200250300350400-1001020residue0100200300400500600700800-1001020residue+residue050100150200250300-1001020residue+cycles问题三：我的理解是对于剩余载荷residue，在复制一个一样的residue，让它们首尾相接作为一个新的过程再进行雨流处理，得到residue2+cycles 。我对文章的理解是residue2 要和原本的residue 一致，但这又如何保证呢？总之，对这里理解的不是很清晰。四、载荷的重构（将循环插入剩余载荷中）（1）筛选
6、出要插入循环Q中高于门槛值threshold的循环存于 S中c,r=size(Q); S=;Abs=;for i=1:c if abs(Q(i,2)-Q(i,1)threshold S=S;Q(i,1),Q(i,2); Abs=Abs,abs(Q(i,2)-Q(i,1);endend （2）对 S中循环按幅值大小降序排列Abs,index=sort(Abs);index=fliplr(index);k1=S(:,1);S(:,1)=k1(index);k2=S(:,2);S(:,2)=k2(index);（3）按幅值由大到小插入载荷序列中（这部分编出来程序运行没有问题，但结果显示似乎并没有插进
7、去，还在寻找问题中 ）名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 5 页，共 6 页 - - - - - - - - - c,r=size(S); no1=;no2=;while m=cif S(m,1)=S(m,2)&(reside(i+1)=S(m,1) %筛选出第 m个循环能够插入的序列rde=rde;reside(i),reside(i+1); indexrde=indexrde;i,i+1;endend k,g=size(rde);if k=0 m=m+1;else c
8、hu=randi(k,1,1); %在所有能够插入的序列中随机选取一个插入 p=indexrde(chu,1); reside=reside(1:p),S(m,1),S(m,2),reside(p+1:length(reside); m=m+1;endelse n=length(reside); rin=;indexrin=;for i=1:n-1if (reside(i)=S(m,1) rin=rin;reside(i),reside(i+1); indexrin=indexrin;i,i+1;endend k,g=size(rin);if k=0 m=m+1;else chu=randi(k,1,1); p=indexrin(chu,1); reside=reside(1:p),S(m,1),S(m,2),reside(p+1:length(reside); m=m+1;endendend问题四：这种以逆雨流法重构载荷的意义就在于去除一些不必要的载荷吗？名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 6 页，共 6 页 - - - - - - - - -