频谱_功率谱_matlab功率谱_功率谱_信号功率谱_matlab功率谱_频谱功率谱

clc; clear all;close all; %% 加载数据原始 load DBN_no_noise.mat ; load DBN_noise_2.mat; load VCO_no_noise.mat; load VCO_noise_1; load VCO_noise_2; data1=DBN_no_noise';data2=DBN_noise_2'; data3=VCO_no_noise';data4=VCO_noise_1'; data5=VCO_noise_2'; fs=1; % %% 将信号放入同一个矩阵方便分析 data_all=[data1(2,:);data2(2,:);data3(2,:);data4(2,:);data5(2,:)]; Len=size(data_all,2); %% 计算频谱 f_x = (1:floor(Len/2))/Len*fs; %频谱范围，截取前半段（抽样频率高于最大频率的2倍) fft_data=zeros(1,floor(Len/2)); for i=1:size(data_all,1) temp_fre=abs((fft(data_all(i,:),Len))); temp_fre=10*log10(temp_fre(1:floor(Len/2))); %由于幅度响应是偶函数，所以截取前一半绘制 fft_data=[fft_data;temp_fre]; end fft_data=fft_data(2:end,:); %% Welch估计功率谱、‘hamming窗’ figure; windows=hamming(256); NFFT=512; noverlap=128; [pxx1,f1]=pwelch(data_all(1,:),windows,noverlap,NFFT,fs); [pxx2,f2]=pwelch(data_all(2,:),windows,noverlap,NFFT,fs); [pxx3,f3]=pwelch(data_all(3,:),windows,noverlap,NFFT,fs); [pxx4,f4]=pwelch(data_all(4,:),windows,noverlap,NFFT,fs); [pxx5,f5]=pwelch(data_all(5,:),windows,noverlap,NFFT,fs); pxx1=10*log10(pxx1); pxx2=10*log10(pxx2); pxx3=10*log10(pxx3); pxx4=10*log10(pxx4); pxx5=10*log10(pxx5); close all %% DBN无噪声加噪声2数据对比 figure; subplot 211;plot(f_x,fft_data(1,:)); title('DBN无噪声频谱');xlabel('频率');ylabel('幅值') subplot 212;plot(f_x,fft_data(2,:)); title('DBN加噪声2频谱');xlabel('频率');ylabel('幅值') figure; subplot 211;plot(f1,pxx1); title('DBN无噪声功率谱');xlabel('频率');ylabel('幅值') subplot 212;plot(f2,pxx2); title('DBN加噪声2功率谱');xlabel('频率');ylabel('幅值') %% VCO无噪声加噪声1数据对比 figure; subplot 211;plot(f_x,fft_data(3,:)); title('VCO无噪声频谱');xlabel('频率');ylabel('幅值') subplot 212;plot(f_x,fft_data(4,:)); title('VCO加噪声1频谱');xlabel('频率');ylabel('幅值') figure; subplot 211;plot(f3,pxx3); title('VCO无噪声功率谱');xlabel('频率');ylabel('幅值') subplot 212;plot(f4,pxx4); title('VCO加噪声1功率谱');xlabel('频率');ylabel('幅值') %% VCO 加噪声1加噪声2数据对比 figure; subplot 211;plot(f_x,fft_data(4,:)); title('DBN加噪声1频谱');xlabel('频率');ylabel('幅值') subplot 212;plot(f_x,fft_data(5,:)); title('DBN加噪声2频谱');xlabel('频率');ylabel('幅值') figure; subplot 211;plot(f4,pxx4); title('DBN加噪声1功率谱');xlabel('频率');ylabel('幅值') subplot 212;plot(f5,pxx5); title('DBN加噪声2功率谱');xlabel('频率');ylabel('幅值') %% DBN 加噪声2和VCO加噪声2数据对比 figure; subplot 211;plot(f_x,fft_data(2,:)); title('VCO加噪声2频谱');xlabel('频率');ylabel('幅值') subplot 212;plot(f_x,fft_data(5,:)); title('DBN加噪声2频谱');xlabel('频率');ylabel('幅值') figure; subplot 211;plot(f2,pxx2); title('VCO加噪声2功率谱');xlabel('频率');ylabel('幅值') subplot 212;plot(f5,pxx5); title('DBN加噪声2功率谱');xlabel('频率');ylabel('幅值')