韦尔奇估计的For循环
5视图(30天)
显示旧的评论
我有一个问题在创建for循环计算韦尔奇估计不使用pwelch函数。谁能帮我做这个吗?
我们要计算韦尔奇估计使用一个100年代peridograms滑动时间间隔。为此,计算第一PSD估计peridograms 1到10,第二估计从2到11日第三估计3到12等。”
f = 1000;
M_welch = 10 * f + 1;%损害窗口长度多
w_welch =损害(M_welch);%的窗口函数直接应用于信号
w_welch = w_welch /规范(w_welch);%单位能源正常化
Lmax =地板(N / M_welch);% peridograms的最大数量
L2 = 10;%的peridograms数100年代间隔
y = x;
Y =重塑(Y (1: Lmax * M_welch), M_welch, Lmax) ';%信号分割成L片段长度M
%将每段信号放入coloumn并创建一个矩阵
Y = Y。* repmat (w_welch Lmax 1);%乘以窗函数——元素的乘法
% phi_sum = 0 (1、2 * m - 1);%初始化和对ACF
M_welch S_sum_welch = 0 (1);%初始化和PSD
f_Hz = (f / 2) *(0:地板(M_welch / 2)) / (M_welch / 2);%相应积极的频率向量(赫兹)
Rp = 1;
Lmax-9 deltaRatio_welch = 0 (1);
为我= 1:Lmax-9% 100年代第一次
为j =我:我+ 9%循环周期图1日至10日周期图
S_sum_welch = S_sum_welch + 2 * T * abs ((fft (Y (j,:)))) ^ 2;%更新PSD %为什么是2 * T
结束
结束
