本文作者在MATLAB中遇到使用cpsd和pwelch函数进行频谱分析的问题，主要疑惑在于函数参数的理解及是否正确应用。代码中计算自频谱密度和互频谱密度，并探讨了(cpsd替代csd, pwelch替代psd)的说法。作者提问了空参数[]的含义，以及数据量N和transform length L的选择是否合适。"
123538346,9125126,Python图像处理：直方图均衡化与高斯滤波详解,"['Python', '图像处理', '计算机视觉', '直方图', '滤波']
本帖最后由 多多会迈特来布 于 2016-2-12 05:24 编辑
这是我遇到问题的代码：
[code]N=2^14; %Number of RecordsL=2^14;%Transform Length
SmoothingK=2^0;
fd=zeros(L,1); %Frequency division matrix
%Sampling Frequency
Fs=(N-1)/(Times(N,1)-Times(1,1));
%Find the resolution from the samplingfrequency
f=Fs/L;
%Auto spectal density for each probe
SP1=pwelch(P1,(L/SmoothingK),[],L)/(Fs/12.566);
SP2=pwelch(P2,(L/SmoothingK),[],L)/(Fs/12.566);
SP3=pwelch(P3,(L/SmoothingK),[],L)/(Fs/12.566);
%Generate cross spectral densities
CSD12=cpsd(P1,P2,(L/SmoothingK),[],L);
CSD13=cpsd(P1,P3,(L/SmoothingK),[],L);
PH11=zeros(L,1);
PH12=zeros(L,1);
PH13=zeros(L,1);
for i = 1:(L/2)
PH12(i,1)=atan2(imag(CSD12(i,1)),real(CSD12(i,1)));
PH13(i,1)=atan2(
最低0.47元/天 解锁文章
7451
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
