function [f,p] = fft_plot(input_pre,Fs,SAM) L = FsSAM;%采样点数 input_p = input_pre(1:(floor(length(input_pre)/L)L)); P = zeros(L/2+1,floor(length(input_p)/L)); for i = 1:floor(length(input_p)/L) f = Fs(0:(L/2))/L; input = input_p((i-1)L+1:iL); primary_data = input-mean(input); primary_data = reshape(primary_data,[],1); primary_data = primary_data(1:L); %数据直接进行fft处理并plot pri_fft = fft(primary_data); P2 = abs(pri_fft/L); P1 = P2(1:L/2+1); P(2:end-1,i) = 2P1(2:end-1); end Pover = mean(P'); Pover = Pover'; plot(f,Pover); p = Pover; grid on hold on end
时间: 2024-03-13 20:44:29 浏览: 77
你修改后的代码存在两个错误:
1.第一行的 `L = FsSAM;%采样点数` 应该改为 `L = Fs*SAM;%采样点数`,因为 `Fs` 和 `SAM` 是两个变量,需要使用乘法运算符 `*`。
2.第 18 行的 `P(2:end-1,i) = 2P1(2:end-1);` 应该改为 `P(2:end-1,i) = 2*P1(2:end-1);`,因为 `2` 和 `P1(2:end-1)` 是两个变量,需要使用乘法运算符 `*`。
修改后的代码如下:
```
function [f,p] = fft_plot(input_pre,Fs,SAM)
L = Fs*SAM;%采样点数
input_p = input_pre(1:(floor(length(input_pre)/L)*L));
P = zeros(L/2+1,floor(length(input_p)/L));
f = Fs*(0:(L/2))/L; % 将 f 的定义放在循环外面
for i = 1:floor(length(input_p)/L)
input = input_p((i-1)*L+1:i*L);
primary_data = input-mean(input);
primary_data = reshape(primary_data,[],1);
primary_data = primary_data(1:L);
%数据直接进行fft处理并plot
pri_fft = fft(primary_data);
P2 = abs(pri_fft/L);
P1 = P2(1:L/2+1);
P(2:end-1,i) = 2*P1(2:end-1); % 修改这一行
end
Pover = mean(P');
Pover = Pover';
plot(f,Pover);
p = Pover;
grid on
hold on
end
```
阅读全文