本篇文章，主要是为上篇文章作补充，给出完整代码即效果。
clear all clc csi_trace = read_bf_file(sample_data/lie1.dat); pac_num=size(csi_trace,1); subcarrier=zeros(1,990); Hall=zeros(990,30); Z=zeros(990,30); Zk=zeros(990,1); for j=1:30 for i=1:990; csi_entry=csi_trace{i}; csi=get_scaled_csi(csi_entry); csi1=squeeze(csi(1,:,:)).;% 30*3 complex csiabs=db(abs(csi1)); csiabs=csiabs(:,2); csi1=csi1(:,2); subcarrier(i)=csiabs(j); if(subcarrier(i)>=35) subcarrier(i)=35; else if(subcarrier(i)<=1) subcarrier(i)=1; end end end yd=wden(subcarrier,heursure,s,one,10,sym3); vs=movvar(yd,50) figure(1) subplot(2,1,1); plot(subcarrier); title(beforeDWT); subplot(2,1,2); plot(yd); title(afterDWT); figure(9) plot(vs); xlabel(time(s)); ylabel(variance); title(movvar_Origin); hold on Hall(:,j)=yd.; end %uncentralize for k=1:30 for i=1:990 Zk(i)=Hall(i,k)-mean(Hall(:,k)); end Z(:,k)=Zk; end %Z covariance C=cov(Z); % Eigenvector e, eigenvalue R [E,R]=eig(C); R=ones(1,30)*R;%The eigenvalues generated by the eig function are transformed into row vectors TZ_ER=[R;E]; TZ_ER=TZ_ER; %The eigenvectors are arranged in descending order according to the size of R TZ_ER=sortrows(TZ_ER,1,descend); R=TZ_ER(:,1);%Separating eigenvalues E=TZ_ER(:,2:end);%Separating eigenvectors Pca=zeros(30,1); E=E Pca=E(:,1);%chose the 1st pinciple component Y=Hall*Pca; res=sum(R(1:2))/sum(R); apv=movvar(Y,50) figure(2) plot(apv,k,linewidth,3); xlabel(time(s)); ylabel(variance); title(movvar_PCA); figure(15) plot (Y,linewidth,3); hold on title(afterPCA); fprintf(consists the percent\n,res*100); for i=1:30 plot(Hall(:,i)); end Butterworth hfc = 300; lfc = 0; fs = 1000; order = 10; [b,a] = butter(order, hfc/(fs/2)); figure(5) freqz(b,a) noisy_sig = D{7}; dataIn = randn(1000,1); dataOut = filter(b,a,Y); figure(6) subplot(2,1,1) plot(Y) hold on title(original); subplot(2,1,2) plot(dataOut); hold on title(afterbutterworth); Y=dataOut; FFT F=fft(Y,1024); figure(9) plot(abs(F)); axis([0 1800 0 1500]); STFT fs = 1000;% window = 512; noverlap = window/2; nfft=1024; f_len = window/2 + 1; f = linspace(0, 150e3, f_len); % s= spectrogram(Y, window, noverlap); % figure(1); % imagesc(20*log10((abs(s))));xlabel(Samples); ylabel(Freqency); % colorbar; [s, f, t, p] = spectrogram(Y, window, nfft, f, fs); figure(8); imagesc(t, f, p);xlabel(Samples); ylabel(Freqency); colorbar; [s, f, t] = spectrogram(Y, 512,256,f, fs); figure; imagesc(t, f, 20*log10((abs(s))));xlabel(Samples); ylabel(Freqency); colorbar; figure(16) subplot(2,1,1) spectrogram(Y,128,127,128,fs,yaxis); title(STFT); hold on [s,f,t,ps]=spectrogram(Y,64,63,64,fs,yaxis); m=4:33 n=1:927 pmax=max(ps(m,n),[],1); subplot(2,1,2) % % figure(17) plot(pmax,linewidth,3,MarkerEdgeColor,k,MarkerFaceColor,g,MarkerSize,10); hold on title(power curve);<
整段代码效果如下图：
