雨流计数法matlab程源代码高等教育大学_高等教育-大学课件.pdf

《雨流计数法matlab程源代码高等教育大学_高等教育-大学课件.pdf》由会员分享，可在线阅读，更多相关《雨流计数法matlab程源代码高等教育大学_高等教育-大学课件.pdf（3页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、%RAIN FLOW cycle counting.%RAINFLOW counting function allows you to extract%cycle from random loading%SYNTAX%rf=RAINFLOW(ext)%rf=RAINFLOW(ext,dt)%rf=RAINFLOW(ext,extt)%OUTPUT%rf-rainflow cycles:matrix 3xn or 5xn dependend on input,%rf(l/)Cycles amplitude,%rf(2/)Cycles mean value”%rf:)Number of cycle

2、s(0.5 or 1.0),%rf(4/)Begining time(when input includes dt or extt data),%rf(5/)Cycle period(whe n in put in eludes dt or extt data),%INPUT%ext-signal points,vector nxl,ONLY TURNING POINTS!,%dt-sampling time,positive number,when the turning points%spaced equally,%extt-signal time,vector nxl,exact tim

3、e of occurrenee of turning points.%See also SIG2EX1；RFHIST,RFMATRIX,RFPDF3D.%RAINFLOW%Copyright(c)1999-2002 by Adam Nieslony,%MEXfu nction.function rfdemol(ext)%function rfdemol(ext)%RFDEMO1 shows cycles extracted from signal%using rainflow algoritm%INPUT:ext-option,number or vectors with turning%po

4、ints or time history.Default ext=16.%OUTPUT:no enable.%SYNTAX:%rfdemol%rfdemol(10)%rfdemol(2 3 2 42 5 1 6)%By Adam Nies 程 ny%Revised,10-Nov-2009%Visit the MATLAB Central File Exchange for latest version.error(nargchk(0,2znargi n)if n argi n=0,%turning points from 16 random numbers ext=sig2ext(rand n

5、(4);elseif length(ext(:)=l/%turning points from n random numbers ext=sig2ext(rand n(l,ext);else%turning points from vector ext ext=sig2ext(ext);end a=rainflow(ext,l);m n=size(a);%if n100,%button=questdlg(Rainflow found 1 num2str(sum(a(3/:)1 cycles!Do you want to con tinue?,.%Continue Operation/Yes/N

6、o/No1)；%if strcmp(button/No)%error(*Function aborted by user.1)%end%end coymcrgb;plot(0:length(ext)-l,ext,k：)hold on wyk=0:0.05:l;for c=l:n”coin r=rem(c-l,6)+l;nrl=round(a(4,c)+l);nr2=round(a(4,c)+l+a(5/c)*a(3/c);if a(3/c)=1.0/if ext(nrl)ext(nr2)z plot(wyk.*a(5/c)*0.5+a(4,c)/cos(wyk.*pi)*a(l,c)+a(2/

7、c)/col(colnr)text(a(4,c)/a(2/c)+a(l/c),int2str(c)：Half-cycle,down1,.Colorcollcolnrl/VerticalAlig nment/bottom1)else plot(wyk.*a(5/c)*0.5+a(4/c)/cos(pi+wyk.*pi)*a(l/c)+a(2/c),col(colnr)text(a(4,c)/a(2/c)-a(l/c),int2str(c)Half-cycle,up1,.Colorcollcolnrl/VerticalAlig nment/top)end end end xlabepeaks,co

8、unted from O)ylabel(value)title(Rainflow cycles extracted from sign al)legend(peaks from signal：。)hold off dispCRow 1:amplitude)dispCRow 2:mearT)dispCRow 3:number of cycles(cycle or half cycle)1)disp(*Row 4:begin time of extracted cycle or half cycle)dispCRow 5:period of a cycle)disp(a)最新文件-仅供参考-已改成 wod 文本-改