风振分析命令.docVIP

下载本文档

2
0
约3.09千字
约 5页
2018-05-26 发布于河南
举报
版权申诉

风振分析命令.doc

1、本文档共5页，可阅读全部内容。
2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
5、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
6、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
7、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
8、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

风振分析命令

！画对比曲线 n=0.001:0.001:0.5 v10=34 x=1200*n/v10 s1=0.667*x.^2./n./(1+x.^2).^(4/3) %规格化功率谱 plot(n,s1) a=(x.^2-0.7729^2).^2+(2*1.477*0.7729*x).^2 b=800*2.182*x.^2/v10 s2=b./a %近似谱 plot(n,s2) plot(n,s1,n,s2) plot(n,s1,k-,n,s2,b:) clear Pij=[ 8.6719 8.9066 11.2417 8.3669 7.1518] h=[16.7 49.6 85 112.8 146] for i=1:5 for j=1:5 Spij(i,j)=exp(-abs(h(i)-h(j))/60).*Pij(i).*Pij(j) end end goal=chol(Spij) goal=goal N=4001 afa=2.182 Wf=0.1376 kexi=1.477 dw=0.01256 V=34 x=sqrt(-1) fai=2*pi*rand(1,N) for time=1:4001 t=(time-1)*0.05 for j=1:5 sum1=0 for m=1:j sum2=0 for k=1:N w=(k-1)*dw F1=sqrt(afa*pi*V/900)*w*x F2=(Wf^2-w^2)+2*kexi*Wf*w*x GS=F1/F2 GD(j,m)=sqrt(2*dw)*goal(j,m)*abs(GS) if k=1 theta=0.0 end if k1 theta=atan2(imag(GS),real(GS)) end Result=GD(j,m)*cos(w*t+theta+fai(k)) sum2=sum2+Result end xie(j,m)=sum2 sum1=sum1+xie(j,m) end ft(time,j)=sum1 end ok(time)=t end for i=1:4001 tu1(i)=ft(i,1) tu2(i)=ft(i,2) tu3(i)=ft(i,3) tu4(i)=ft(i,4) tu5(i)=ft(i,5) end n=0.001:0.001:10 v10=34 x=1200*n/v10 a=(x.^2-0.7729^2).^2+(2*1.477*0.7729*x).^2 b=800*2.182*x.^2/v10 s2=b./a %近似谱 c=2.182*v10*pi*pi*pi*n.^2/225 d=(4*pi*pi*n.^2-0.1376^2).^2+ 6.5226*n.^2 s3=c./d [pxx,f]=psd(tu1,1002,20,hanning(1002),501) f=f*0.2 loglog(f,pxx,n,s2) clear Pij=[ 8.6719 8.9066] h=[16.7 49.6 ] for i=1:2 for j=1:2 Spij(i,j)=exp(-abs(h(i)-h(j))/60).*Pij(i).*Pij(j) end end goal=chol(Spij) goal=goal N=501 afa=2.182 Wf=0.1376 kexi=1.477 dw=0.01256 V=34 x=sqrt(-1) fai=2*pi*rand(1,N) for time=1:501 t=(time-1)*0.05 for j=1:2 sum1=0 for m=1:j sum2=0 for k=1:N w=(k-1)*dw F1=sqrt(afa*pi*V/900)*w*x F2=(Wf^2-w^2)+2*kexi*Wf*w*x GS=F1/F2 GD(j,m)=sqrt(2*dw)*goal(j,m)*abs(GS) if k=1 theta=0.0 end if k1 theta=atan2(imag(GS),real(GS)) end Result=GD(j,m)*cos(w*t+theta+fai