【正文】 oising[J],Neural Networks and Learning Systems,2015,Volume(PP):11 .[13] Scanlon,Patricia。 Kennedy, Irwin O.。Liu,Yongheng. Feature extraction approaches to?RF fingerprinting for device identification in femtocells[J]. Bell Labs Technical Journal,2010, Volume(15):141151.附錄一close allclear allclcN=500。fs=500。n=1:N。dt=1/fs。a=[1 1 1 0 0 1 1 0 1 0 ]。g1=a。g2=~a。g11=(ones(1,500))39。*g1。g1a=g11(:)39。t=0:dt:10dt。t1=length(t)。ask1=g1a.*cos(2*pi*50*t)。%ask的頻率為50ask2=g1a.*cos(2*pi*45*t)。%ask的頻率為45ask3=g1a.*cos(2*pi*40*t)。%ask的頻率為40figure(1)。subplot(311)。plot(t,ask1)。title(39。ask1原始信號(hào)39。)。xlabel(39。時(shí)間 t39。)。ylabel(39。幅值 A39。)。grid on。subplot(312)。plot(t,ask2)。title(39。ask2原始信號(hào)39。)。xlabel(39。時(shí)間 t39。)。ylabel(39。幅值 A39。)。grid on。subplot(313)。plot(t,ask3)。title(39。ask3原始信號(hào)39。)。xlabel(39。時(shí)間 t39。)。ylabel(39。幅值 A39。)。grid on。figure(2)。subplot(311)。[tfr1,t,f1]=tfrwv(ask139。,1:N,N)。contour(t/fs,f1(1:length(f1)/2)*fs,abs(tfr1(1:length(f1)/2,:)))。title(39。ASK1的WignerVille分布時(shí)頻圖39。)。xlabel(39。時(shí)間 t39。)。ylabel(39。頻率 f39。)。grid on。subplot(312)。[tfr2,t,f2]=tfrwv(ask239。,1:N,N)。contour(t/fs,f2(1:length(f2)/2)*fs,abs(tfr2(1:length(f2)/2,:)))。title(39。ASK2的WignerVille分布時(shí)頻圖39。)。xlabel(39。時(shí)間 t39。)。ylabel(39。頻率 f39。)。grid on。subplot(313)。[tfr3,t,f3]=tfrwv(ask339。,1:N,N)。contour(t/fs,f3(1:length(f3)/2)*fs,abs(tfr3(1:length(f3)/2,:)))。title(39。ASK3的WignerVille分布時(shí)頻圖39。)。xlabel(39。時(shí)間 t39。)。ylabel(39。頻率 f39。)。grid on。figure(3)。mesh(t/fs,f1(1:length(f1)/2)*fs,abs(tfr1(1:length(f1)/2,:)))。title(39。ASK1的三維時(shí)頻圖39。)。xlabel(39。時(shí)間 t39。)。ylabel(39。頻率 f39。)。zlabel(39。幅值 A39。)。grid on。figure(4)。mesh(t/fs,f2(1:length(f2)/2)*fs,abs(tfr2(1:length(f2)/2,:)))。title(39。ASK2的三維時(shí)頻圖39。)。xlabel(39。時(shí)間 t39。)。ylabel(39。頻率 f39。)。zlabel(39。幅值 A39。)。grid on。figure(5)。mesh(t/fs,f3(1:length(f3)/2)*fs,abs(tfr3(1:length(f3)/2,:)))。title(39。ASK3的三維時(shí)頻圖39。)。xlabel(39。時(shí)間 t39。)。ylabel(39。頻率 f39。)。zlabel(39。幅值 A39。)。grid on。[COEFF1,SCORE1,latent1] = prinp(tfr1)。[COEFF2,SCORE2,latent2] = prinp(tfr2)。[COEFF3,SCORE3,latent3] = prinp(tfr3)。S1=sum(latent1)。S2=sum(latent2)。S3=sum(latent3)。C1=100*latent1/sum(latent1)。C2=100*latent2/sum(latent2)。C3=100*latent3/sum(latent3)。figure(6)。pareto(C1)。title(39。ask1的貢獻(xiàn)率39。)。figure(7)。pareto(C2)。 title(39。ask2的貢獻(xiàn)率39。)。figure(8)。pareto(C3)。 title(39。ask3的貢獻(xiàn)率39。)。SCORE1_test=SCORE139。SCORE2_test=SCORE239。SCORE3_test=SCORE339。pca1=SCORE1_test(1:500*2)。pca2=SCORE2_test(1:500*2)。pca3=SCORE3_test(1:500*2)。d12=pdist2(pca1,pca2,39。Euclidean39。)。d13=pdist2(pca1,pca3,39。Euclidean39。)。tfr11=tfr1(1:500*500)。tfr21=tfr2(1:500*500)。tfr31=tfr3(1:500*500)。D12=pdist2( tfr11,tfr21,39。Euclidean39。)。D13=pdist2( tfr11,tfr31,39。Euclidean39。)。disp(39。降維后ask1分別與askask3的距離:39。)。disp(d12)。disp(d13)。disp(39。降維前ask1分別與askask3的距離:39。)。disp(D12)。disp(D13)。附錄二function [tfr,t,f] = tfrwv(x,t,N,trace)%TFRWV WignerVille時(shí)頻分布.% [TFR,T,F]=TFRWV(X,T,N,TRACE) 計(jì)算維格納分布% 一個(gè)離散時(shí)間信號(hào)x, % 或表示兩信號(hào)之間互維格納. % % X : 信號(hào)如果自動(dòng)維納，或[X1，X2]如果交叉維納.% T : 瞬時(shí)時(shí)間 (默認(rèn) : 1:length(X)).% N : 矩陣頻率數(shù)量 (默認(rèn) : length(X)).% TRACE : 如果非零，算法的進(jìn)展顯示(默認(rèn) : 0).% TFR : 時(shí)頻表示。當(dāng)沒(méi)有輸出參數(shù)，tfrwv運(yùn)行tfrqview。% F : 歸一化頻率向量.% 例如:% sig=fmlin(128,)。 tfrwv(sig)。% % %%% if (nargin == 0), error(39。At least one parameter required39。)。end。[xrow,xcol] = size(x)。 if (nargin == 1), t=1:xrow。 N=xrow 。 trace=0。elseif (nargin == 2), N=xrow 。 trace=0。elseif (nargin == 3), trace = 0。end。 if (N0), error(39。N must be greater than zero39。)。end。[trow,tcol] = size(t)。if (xcol==0)||(xcol2), error(39。X must have one or two columns39。)。elseif (trow~=1), error(39。T must only have one row39。)。 elseif (2^nextpow2(N)~=N), fprintf(39。For a faster putation, N should be a power of two\n39。)。end。 tfr= zeros (N,tcol)。 if trace, disp(39。WignerVille distribution39。)。 end。for icol=1:tcol, ti= t(icol)。 taumax=min([ti1,xrowti,round(N/2)1])。 tau=taumax:taumax。 indices= rem(N+tau,N)+1。 tfr(indices,icol) = x(ti+tau,1) .* conj(x(titau,xcol))。 tau=round(N/2)。 if (ti=xrowtau)amp。amp。(ti=tau+1), tfr(tau+1,icol) = * (x(ti+tau,1) * conj(x(titau,xcol)) + ... x(titau,1) * conj(x(ti+tau,xcol))) 。 end。 if trace, disprog(icol,tcol,10)。 end。end。 tfr= fft(tfr)。 if (xcol==1), tfr=real(tfr)。 end 。 if (nargout==0), tfrqview(tfr,x,t,39。tfrwv39。)。elseif (nargout==3), f=(*(0:N1)/N)39。end。30