1、内容简介
略
2、内容说明
略
3、仿真分析
clc
clear
close all
%企q号q:1762016542
load DATA.mat
fs=173.61;
t = 0:1/fs:177/fs;
x1 = X(1,1:178);
figure
plot(t,x1)
xlabel 时间/s
ylabel 幅值
title 时域
x_end = X(:,end);
disp('种类数目')
min(x_end)
max(x_end)
[f,y0]=change_fft(fs,x1);
figure
plot(f,y0)
xlabel 频率/Hz
ylabel 幅值
title 频域
num = length(x_end);
PSK = zeros(num,5);
RIR1 = zeros(num,5);
for i = 1:num
[f,y0]=change_fft(fs,X(i,1:178));
PSK(i,:) = PSI(f,y0);%调用函数
temp = PSI(f,y0);
RIR1(i,:) = RIR(temp);
end
figure
plot(PSK(:,1))
hold on
plot(PSK(:,2))
hold on
plot(PSK(:,3))
hold on
plot(PSK(:,4))
hold on
plot(PSK(:,5))
hold off
xlabel 行数
ylabel PSI
title 功率谱
figure
plot(RIR1(:,1))
hold on
plot(RIR1(:,2))
hold on
plot(RIR1(:,3))
hold on
plot(RIR1(:,4))
hold on
plot(RIR1(:,5))
hold off
xlabel 行数
ylabel RIR
title 功率谱
xlswrite('data.xls', [PSK RIR1], 'Sheet1')
4、参考论文
基于发作间期EEG的癫痫自动诊断系统的研究与设计_刘国权.caj