FFT及周期图法算功率谱密度

你是否正在寻找关于功率谱密度的内容？让我把最内涵的东西奉献给你：

FFT及周期图法算功率谱密度

t=0:0.0001:0.1; %时间间隔为0.0001，说明采样频率为10000Hz x=square(2*pi*1000*t); %产生基频为1000Hz的方波信号 n=randn(size(t)); %白噪声 f=x; %在信号中加入白噪声 figure(1); subplot(2,1,1); plot(t,f); %画出原始信号的波形图 ylabel('幅值(V)'); xlabel('时间(s)'); title('原始信号'); y=fft(f,1024); %对原始信号进行离散傅里叶变换，参加DFT采样点的个数为1024 subplot(2,1,2); m=abs(y); f1=(0:length(y)/2-1)'*10000/length(y);%计算变换后不同点对应的幅值 plot(f1,m(1:length(y)/2)*2/length(y)); ylabel('幅值的模');xlabel('时间(s)');title('原始信号傅里叶变换'); %用周期图法估计功率谱密度 p=y.*conj(y)/1024; %计算功率谱密度 ff=10000*(0:499)/1024; %计算变换后不同点对应的频率值 figure(2);plot(ff,p(1:500));ylabel('幅值');xlabel('频率(Hz)'); title('功率谱密度（周期图法）');

求FFT变换时，纵坐标和横坐标的问题

clear all;clc fs=500; t=(0:1:100)/fs; s=cos(2*pi*50*t); n=length(t); figure(1),subplot(211);plot(t,s) f=fftshift(fft(s)); subplot(212);plot((-n/2:n/2-1)*fs/n,abs(f)*2/n);

posted on

以上就是关于功率谱密度的全部内容，相信你一定会非常满意，。

本文来自电脑杂谈，转载请注明本文网址：
http://www.pc-fly.com/a/jisuanjixue/article-6623-1.html