pinpu.zip_桥梁_桥梁matlab_频谱

%%%%%康妮采样频率128Hz，东华5920 clear clc close all hidden data=load('F:\毕业终极资料\新建文件夹 (2)\温州\2013\龙湾基地台风菲特20131005\处理后数据\feite201310061745_1925\feite201310061745g.TXT');%文件路径 xx1=data(:,8); %10x xx2=data(:,9); %10y % xx3=data(:,7); %18x % xx4=data(:,8); %18y % xx5=data(:,9); %25x % xx6=data(:,10); %25y % xx7=data(:,11); %32x % xx8=data(:,12); %32y % xx9=data(:,13); %41x % xx10=data(:,14); %41y % xx1=xx1-mean(xx1) % xx2=xx2-mean(xx2) xx1=xx1-mean(xx1) xx2=xx2-mean(xx2) fs=25.6; nfft=2^18; window=hamming(2048); noverlap=25.6; [Pxx1,f]=pwelch(xx1,window,noverlap,nfft,fs); [Pxx2,f]=pwelch(xx2,window,noverlap,nfft,fs); % [Pxx3,f]=pwelch(xx3,window,noverlap,nfft,fs); % [Pxx4,f]=pwelch(xx4,window,noverlap,nfft,fs); % [Pxx5,f]=pwelch(xx5,window,noverlap,nfft,fs); % [Pxx6,f]=pwelch(xx6,window,noverlap,nfft,fs); % [Pxx7,f]=pwelch(xx7,window,noverlap,nfft,fs); % [Pxx8,f]=pwelch(xx8,window,noverlap,nfft,fs); % [Pxx9,f]=pwelch(xx9,window,noverlap,nfft,fs); % [Pxx10,f]=pwelch(xx10,window,noverlap,nfft,fs); % p=0.95; % % x1=xx1-xx9;%%%以41层为准，相同方向的相减求互谱密度 % x2=xx3-xx9; % x3=xx5-xx9; % x4=xx7-xx9; % % x5=xx2-xx10; %%10 % x6=xx4-xx10; %%18 % x7=xx6-xx10; %%25 % x8=xx8-xx10; %%32 % % % % [Px1,f1]=csd(xx9,x1,nfft,fs,hanning(1000),0,p);%%互谱密度 % [Px2,f1]=csd(xx9,x2,nfft,fs,hanning(1000),0,p);%%互谱密度 % [Px3,f1]=csd(xx9,x3,nfft,fs,hanning(1000),0,p);%%互谱密度 % [Px4,f1]=csd(xx9,x4,nfft,fs,hanning(1000),0,p);%%互谱密度 % % [Px5,f1]=csd(xx10,x5,nfft,fs,hanning(1000),0,p);%%互谱密度 % [Px6,f1]=csd(xx10,x6,nfft,fs,hanning(1000),0,p);%%互谱密度 % [Px7,f1]=csd(xx10,x7,nfft,fs,hanning(1000),0,p);%%互谱密度 % [Px8,f1]=csd(xx10,x8,nfft,fs,hanning(1000),0,p);%%互谱密度 ls = 12; %level %dbs = 5:20; dbs = 15; wn = strcat( 'db', num2str( dbs )); [ctl ltl] = wavedec( Pxx2, ls, wn ); tla = wrcoef( 'a', ctl, ltl, wn, ls ); tlds=[]; j=1; while j<=ls tld = wrcoef( 'd', ctl, ltl, wn, j ); if j>1 tlds = tlds + tld; else tlds=tld; end j=j+1; end i=i+1; Pxx=tla; % figure(1); % plot(f,Pxx,'linewidth',2); % grid on % axis square % set(gca, 'XLim', [0 14],'FontSize',18,'FontName','Times New Roman'); % set(gca, 'XTick', [0 0.4 0.8 1.2 1.6 2 2.4 2.8 3.2 3.6 4]); % xlabel('f/Hz','FontSize',18,'FontName','Times New Roman'); % ylabel('10层楼板X自谱密度/(w/Hz)','FontSize',18,'FontName','Times New Roman'); figure(2); plot(f,Pxx,'linewidth',2); grid on axis square % set(gca, 'XLim', [0 4],'FontSize',18,'FontName','Times New Roman'); % set(gca, 'XTick', [0 0.4 0.8 1.2 1.6 2 2.4 2.8 3.2 3.6 4]); xlabel('f/Hz','FontSize',18,'FontName','Times New Roman'); ylabel('纵向自谱密度/(w/Hz)','FontSize',18,'FontName','Times New Roman'); figure(3); plot(f,10*log10(Pxx3),'linewidth',2); grid on axis square set(gca, 'XLim', [0 4],'FontSize',18,'FontName','Times New Roman'); set(gca, 'XTick', [0 0.4 0.8 1.2 1.6 2 2.4 2.8 3.2 3.6 4]); xlabel('f/Hz','FontSize',18,'FontName','Times New Roman'); ylabel('18层楼板X自谱密度/(w/Hz)','FontSize',18,'FontName','Times New Roman'); figure(4); plot(f,10*log10(Pxx4),'linewidth',2); grid on axis square set(gca, 'XLim', [0 4],'FontSize',18,'FontName','Times New Roman'); set(gca, 'XTick', [0 0.4 0.8 1.2 1.6 2 2.4 2.8 3.2 3.6 4]); xlabel('f/Hz','FontSize',18,'FontName','Times New Roman'); ylabel('18层楼板Y自谱密度/(w/Hz)','FontSize',18,'FontName','Times New Roman'); figure(5); plot(f,10*log10(Pxx5),'linewidth',2); grid on axis square set(gca, 'XLim', [0 4],'FontSize',18,'FontName','Times New Roman'); set(gca, 'XTick', [0 0.4 0.8 1.2 1.6 2 2.4 2.8 3.2 3.6 4]); xlabel('f/Hz','FontSize',18,'FontName','Times New Roman'); ylabel('25层楼板X自谱密度/(w/Hz)','FontSize',18,'FontName','Times New Roman'); figure(6); plot(f,10*log10(Pxx6),'linewidth',2); grid on axis square set(gca, 'XLim', [0 4],'FontSize',18,'FontName','Times New Roman'); set(gca, 'XTick', [0 0.4 0.8 1.2 1.6 2 2.4 2.8 3.2 3.6 4]); xlabel('f/Hz','FontSize',18,'FontName','Times New Roman'); ylabel('25层楼板Y自谱密度/(w/Hz)','FontSize',18,'FontName','Times New Roman'); figure(7); plot(f,10*log10(Pxx7),'linewidth',2); grid on axis square set(gca, 'XLim', [0 4],'FontSize',18,'FontName','Times New Roman'); set(gca, 'XTick', [0 0.4 0.8 1.2 1.6 2 2.4 2.8 3.2 3.6 4]); xlabel('f/Hz','FontSize',18,'FontName','Times New Roman'); ylabel('32层楼板X自谱密度/(w/Hz)','FontSize',18,'FontName','Times New Roman'); figure(8); plot(f,10*log10(Pxx8),'linewidth',2); grid on axis square set(gca, 'XLim', [0 4],'FontSize',18,'FontName','Times New Roman'); set(gca, 'XTick', [0 0.4 0.8 1.2 1.6 2 2.4 2.8 3.2 3.6 4]); xlabel('f/Hz','FontSize',18,'FontName','Times New Roman'); ylabel('32层楼板Y自谱密度/(w/Hz)','FontSize',18,'FontName','Times New Roman'); figure(9); plot(f,10*log10(Pxx9),'linewidth',2); grid on axis square set(gca, 'XLim', [0 4],'FontSize',18,'FontName','Times New Roman'); set(gca, 'XTick', [0 0.4 0.8 1.2 1.6 2 2.4 2.8 3.2 3.6 4]); xlabel('f/Hz','FontSize',18,'FontName','Times New Roman'); ylabel('41层楼板X自谱密度/(w/Hz)','FontSize',18,'FontName','Times New Roman'); figure(10); plot(f,10*log10(Pxx10),'linewidth',2); grid on axis square set(gca, 'XLim', [0 4],'FontSize',18,'FontName','Times New Roman'); set(gca, 'XTick', [0 0.4 0.8 1.2 1.6 2 2.4 2.8 3.2 3.6 4]); xlabel('f/Hz','FontSize',18,'FontName','Times New Roman'); ylabel('41层楼板Y自谱密度/(w/Hz)','FontSize',18,'FontName','Times New Roman'); %下面为求互谱和相位的图形，自己再改下 %y向互谱 figure(11); plot(f1,10*log10(Px5)); grid on set(gca, 'XLim', [0 4],'FontSize',18,'FontName','Times New Roman'); set(gca, 'XTick', [0 0.4 0.8 1.2 1.6 2 2.4 2.8 3.2 3.6 4]); xlabel('f/Hz','FontSize',18,'FontName','Times New Roman'); ylabel('10层楼板Y互谱密度/(w/Hz)','FontSize',18,'FontName','Times New Roman'); figure(12); plot(f1,10*log10(Px6)); grid on set(gca, 'XLim', [0 4],'FontSize',18,'FontName','Times New Roman'); set(gca, 'XTick', [0 0.4 0.8 1.2 1.6 2 2.4 2.8 3.2 3.6 4]); xlabel('f/Hz','FontSize',18,'FontName','Times New Roman'); ylabel('18层楼板Y互谱密度/(w/Hz)','FontSize',18,'FontName','Times New Roman'); figure(13); plot(f1,10*log10(Px7)); grid on set(gca, 'XLim', [0 4],'FontSize',18,'FontName','Times New Roman'); set(gca, 'XTick', [0 0.4 0.8 1.2 1.6 2 2.4 2.8 3.2 3.6 4]); xlabel('f/Hz','FontSize',18,'FontName','Times New Roman'); ylabel('25层楼板Y互谱密度/(w/Hz)','FontSize',18,'FontName','Times New Roman'); figure(14); plot(f1,10*log10(Px8)); grid on set(gca, 'XLim', [0 4],'FontSize',18,'FontName','Times New Roman'); set(gca, 'XTick', [0 0.4 0.8 1.2 1.6 2 2.4 2.8 3.2 3.6 4]); xlabel('f/Hz','FontSize',18,'FontName','Times New Roman'); ylabel('32层楼板Y互谱密度/(w/Hz)','FontSize',18,'FontName','Times New Roman'); %x向互谱 figure(15); plot(f1,10*log10(Px1)); grid on set(gca, 'XLim', [0 4],'FontSize',18,'FontName','Times New Roman'); set(gca, 'XTick', [0 0.4 0.8 1.2 1.6 2 2.4 2.8 3.2 3.6 4]); xlabel('f/Hz','FontSize',18,'FontName','Times New Roman'); ylabel('10层楼板X互谱密度/(w/Hz)','FontSize',18,'FontName','Times New Roman'); figure(16); plot(f1,10*log10(Px2)); grid on set(gca, 'XLim', [0 4],'FontSize',18,'FontName','Times New Roman'); set(gca, 'XTick', [0 0.4 0.8 1.2 1.6 2 2.4 2.8 3.2 3.6 4]); xlabel('f/Hz','FontSize',18,'FontName','Times New Roman'); ylabel('18层楼板X互谱密度/(w/Hz)','FontSize',18,'FontName','Times New Roman'); figure(17); plot(f1,10*log10(Px3)); grid on set(gca, 'XLim', [0 4],'FontSize',18,'FontName','Times New Roman'); set(gca, 'XTick', [0 0.4 0.8 1.2 1.6 2 2.4 2.8 3.2 3.6 4]); xlabel('f/Hz','FontSize',18,'FontName','Times New Roman'); ylabel('25层楼板X互谱密度/(w/Hz)','FontS