gonglvpu.rar_welch_力窗函数

clc clear fs=3200; N=1600; % t=0:0.0003125:0.3; % x0=sin(2*pi*50*t); % % x01=0.4*exp(-250*t).*sin(2*pi*1000*t); % r=zeros(1,length(t)); % for z=1:400 % r(z+400)=x01(z); % end % x=x0+r; [x t]=PQgenerator1(5,fs,N,0,0.1162,0.1286,0.45,8,5); y=awgn(x,10,'measured'); nfft=1024; %cxn=xcorr(y,'unbiased'); cxk=fft(y,nfft); pxx=abs(cxk); index=0:round(nfft/2-1); k=index*fs/nfft; plotpxx=pxx(index+1); figure(2); % subplot(2,1,1); % plot(t,x); % axis([0,0.5,-2,2]); % xlabel('t/s'); % ylabel('电压/mv'); % subplot(2,1,2); plot(k,plotpxx); %axis([0,800,0,300]); xlabel('f/Hz'); ylabel('功率'); h2=axes('Position',[0.5 0.5 0.3 0.3]) plot(k(75:175),plotpxx(75:175))