功率谱计算

clc
clear all
close all
s = load(‘C:\Users\lenovo\Desktop\f4_zy\f4_2_djg\EEG_channel(6).txt’);
s=s/1000/220;
len =length(s);
t = s(2:100001);
len =length(t);
max = max(t);
min = min(t);
fs = 1000;
figure(1)
plot(t);
[pxx,f] = pwelch(t,len,len-1,len,fs);
figure(2)
plot(f,pxx);
m=max-min;