资源描述：
《风振分析命令》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、！画对比曲线n=0.001:0.001:0.5v10=34x=1200*n/v10s1=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).^2b=800*2.182*x.^2/v10s2=b./a%近似谱plot(n,s2)plot(n,s1,n,s2)plot(n,s1,'k-',n,s2,'b:')clearPij=[8.67198.906611.24178.36697.1518]h=[16.749.
2、685112.8146]fori=1:5forj=1:5Spij(i,j)=exp(-abs(h(i)-h(j))/60).*Pij(i).*Pij(j)endendgoal=chol(Spij)goal=goal'N=4001afa=2.182Wf=0.1376kexi=1.477dw=0.01256V=34x=sqrt(-1)fai=2*pi*rand(1,N)fortime=1:4001t=(time-1)*0.05forj=1:5sum1=0form=1:jsum2=0fork=1:Nw=(k-1)*dwF1=sqrt
3、(afa*pi*V/900)*w*xF2=(Wf^2-w^2)+2*kexi*Wf*w*xGS=F1/F2GD(j,m)=sqrt(2*dw)*goal(j,m)*abs(GS)ifk<=1theta=0.0endifk>1theta=atan2(imag(GS),real(GS))endResult=GD(j,m)*cos(w*t+theta+fai(k))sum2=sum2+Resultendxie(j,m)=sum2sum1=sum1+xie(j,m)endft(time,j)=sum1endok(time)=tendf
4、ori=1:4001tu1(i)=ft(i,1)tu2(i)=ft(i,2)tu3(i)=ft(i,3)tu4(i)=ft(i,4)tu5(i)=ft(i,5)endn=0.001:0.001:10v10=34x=1200*n/v10a=(x.^2-0.7729^2).^2+(2*1.477*0.7729*x).^2b=800*2.182*x.^2/v10s2=b./a%近似谱c=2.182*v10*pi*pi*pi*n.^2/225d=(4*pi*pi*n.^2-0.1376^2).^2+6.5226*n.^2s3=c./d
5、[pxx,f]=psd(tu1,1002,20,hanning(1002),501)f=f*0.2loglog(f,pxx,n,s2)clearPij=[8.67198.9066]h=[16.749.6]fori=1:2forj=1:2Spij(i,j)=exp(-abs(h(i)-h(j))/60).*Pij(i).*Pij(j)endendgoal=chol(Spij)goal=goal'N=501afa=2.182Wf=0.1376kexi=1.477dw=0.01256V=34x=sqrt(-1)fai=2*pi*ra
6、nd(1,N)fortime=1:501t=(time-1)*0.05forj=1:2sum1=0form=1:jsum2=0fork=1:Nw=(k-1)*dwF1=sqrt(afa*pi*V/900)*w*xF2=(Wf^2-w^2)+2*kexi*Wf*w*xGS=F1/F2GD(j,m)=sqrt(2*dw)*goal(j,m)*abs(GS)ifk<=1theta=0.0endifk>1theta=atan2(imag(GS),real(GS))endResult=GD(j,m)*cos(w*t+theta+fai(
7、k))sum2=sum2+Resultendxie(j,m)=sum2sum1=sum1+xie(j,m)endft(time,j)=sum1endok(time)=tendfori=1:501tu1(i)=ft(i,1)tu2(i)=ft(i,2)tu3(i)=ft(i,3)tu4(i)=ft(i,4)tu5(i)=ft(i,5)endclearPij=[8.67198.9066]h=[16.749.6]fori=1:2forj=1:2Spij(i,j)=exp(-abs(h(i)-h(j))/60).*Pij(i).*Pi
8、j(j)endendgoal=chol(Spij)goal=goal'N=1001afa=2.182Wf=0.1376kexi=1.477dw=0.01256V=34x=sqrt(-1)fai=2*pi*rand(1,N)fortime=1:1001t=(time-1)*0.
