MATLAB求助,用FFT直接计算x=cos(2*pi*40*t)+2sin(2*pi*120*t)+w(t)的功率谱.用FFT直接计算x=cos(2*pi*40*t)+2sin(2*pi*120*t)+w(t)的功率谱.最好带上注释,
来源：学生作业帮助网 编辑：六六作业网 时间：2024/12/19 23:49:27
MATLAB求助,用FFT直接计算x=cos(2*pi*40*t)+2sin(2*pi*120*t)+w(t)的功率谱.用FFT直接计算x=cos(2*pi*40*t)+2sin(2*pi*120*t)+w(t)的功率谱.最好带上注释,
MATLAB求助,用FFT直接计算x=cos(2*pi*40*t)+2sin(2*pi*120*t)+w(t)的功率谱.
用FFT直接计算x=cos(2*pi*40*t)+2sin(2*pi*120*t)+w(t)的功率谱.
最好带上注释,
MATLAB求助,用FFT直接计算x=cos(2*pi*40*t)+2sin(2*pi*120*t)+w(t)的功率谱.用FFT直接计算x=cos(2*pi*40*t)+2sin(2*pi*120*t)+w(t)的功率谱.最好带上注释,
clc
clear
Fs=500;%采样频率
t=0:1/Fs:1;
xn=cos(2*pi*40*t)+2*sin(2*pi*120*t);%自己加入噪声
nfft=512;
cxn=xcorr(xn,'unbiased');%计算自相关函数
CXk=fft(cxn,nfft);
Pxx=abs(CXk);
index=0:round(nfft/2-1);
k=index*Fs/nfft;
%plot_Pxx=10*log(Pxx(index+1));
plot_Pxx=Pxx(index+1);
plot(k,plot_Pxx)
没有加入噪声,就是你函数中的w(t) 部分,自己加吧.
Fs=1000; %采样频率
n=0:1/Fs:1;
%产生含有噪声的序列
xn=cos(2*pi*40*n)+3*cos(2*pi*100*n)+randn(size(n));
window=boxcar(length(xn)); %矩形窗
nfft=1024;
[Pxx,f]=periodogram(xn,window,nfft,Fs); %直接...
全部展开
Fs=1000; %采样频率
n=0:1/Fs:1;
%产生含有噪声的序列
xn=cos(2*pi*40*n)+3*cos(2*pi*100*n)+randn(size(n));
window=boxcar(length(xn)); %矩形窗
nfft=1024;
[Pxx,f]=periodogram(xn,window,nfft,Fs); %直接法
plot(f,10*log10(Pxx));
%间接法
cxn=xcorr(symbol,'unbiased'); %计算序列的自相关函数
P2=fft(cxn,128);
subplot(1,3,2);
plot(f,10*log10(P2(1:65)));
xlabel('frequency')
ylabel('power')
title('间接法')
%加窗法
window2=blackman(100); %blackman窗
noverlap=20; %数据无重叠
range='onesided'; %频率间隔为[0 1000/2]，只计算一半的频率
[P3(1:65),f]=pwelch(symbol,window2,noverlap,128,1000,range);
plot_P3=10*log10(P3(1:65));
subplot(1,3,3);
plot(f,plot_P3(1:65));
xlabel('frequency')
ylabel('power')
title('加窗法')
收起