经典PSD法定子电流分析的matlab程序

以下是一个使用Matlab实现经典PSD法进行定子电流分析的简单示例程序：

% 假设已经获取到了电流波形数据，存储在变量current中

fs = 1000; % 采样频率
t = (0:length(current)-1)/fs; % 时间向量
N = length(current); % 数据长度

% 计算比例特征（Proportional）
proportional = current - mean(current);

% 计算积分特征（Integral）
integral = cumsum(current - mean(current))/fs;

% 计算微分特征（Derivative）
derivative = diff(current - mean(current))*fs;

% 绘制原始信号和特征信号的波形图
figure;
subplot(4,1,1);
plot(t, current);
title('原始信号');
xlabel('时间 (s)');
ylabel('电流 (A)');

subplot(4,1,2);
plot(t, proportional);
title('比例特征');
xlabel('时间 (s)');
ylabel('电流 (A)');

subplot(4,1,3);
plot(t(1:N-1), integral);
title('积分特征');
xlabel('时间 (s)');
ylabel('电流积分 (A*s)');

subplot(4,1,4);
plot(t(1:N-1), derivative);
title('微分特征');
xlabel('时间 (s)');
ylabel('电流变化率 (A/s)');

这个程序假设你已经将电流波形数据存储在名为`current`的变量中，然后分别计算了比例特征、积分特征和微分特征。最后，使用subplot函数将原始信号和特征信号的波形图绘制在一张图上。

请注意，这只是一个简单的示例程序，你可以根据实际需要进行修改和扩展。另外，根据你的数据和分析目的，可能需要进行更多的数据预处理和特征提取步骤。