% 计算功率谱密度函数-----------
f0101=[]; P0101_Along_SC=[]; P0101_Acros_SC=[]; P0101_Verta_SC=[];
P0101_Along_VK=[]; P0101_Acros_VK=[]; P0101_Verta_VK=[];
NFFT=6000; window=boxcar(6000);
[SPA,n]=pwelch(F010101_Along_Vibr(:,2),window,[],NFFT,Fs); % (样本2)\1塔\10m高
semilogx(n,SPA);
grid on;
f=n.*F010101_Along_Jfcd(2)./F010101_Along_mean(2); f0101=[f0101 f];
P_Temp1=n.*SPA./(F010101_Along_Rmsd(2)).^2;
P_Temp2=4.*f./(1+70.8.*f.^2).^(5/6);
P0101_Along_SC=[P0101_Along_SC P_Temp1];
P0101_Along_VK=[P0101_Along_VK P_Temp2];
[SPA,n]=pwelch(F010101_Acros_Vibr(:,2),window,[],NFFT,Fs); % (样本2)\1塔\10m高
plot(n(1:600,1),SPA(1:600,1),'k')
P_Temp1=n.*SPA./(F010101_Acros_Rmsd(2)).^2;
P_Temp2=4.*f.*(1+755.2*f.^2)./(1+283.2.*f.^2).^(11/6);
P0101_Acros_SC=[P0101_Acros_SC P_Temp1];
P0101_Acros_VK=[P0101_Acros_VK P_Temp2];
[SPA,n]=pwelch(F010101_Verta_Vibr(:,2),window,[],NFFT,Fs); % (样本2)\1塔\10m高
plot(n,SPA,'k')
P_Temp1=n.*SPA./(F010101_Verta_Rmsd(2)).^2;
P_Temp2=4.*f.*(1+755.2*f.^2)./(1+283.2.*f.^2).^(11/6);
P0101_Verta_SC=[P0101_Verta_SC P_Temp1];
P0101_Verta_VK=[P0101_Verta_VK P_Temp2];
评论1