暂态分析matlab计算,关于暂态分析的电流图分析（Prony算法程序调试）

搭建了输电线路的仿真图出来了故障电流的数据，把数据提出来用程序可以运行，但是我不知道如何设置这个算法的采样点，导致出来的频谱结果一直都是50HZ工频的。完全没有其他高频分量能帮忙看看程序那里改动可以从故障时间点开始采样吗？

%打开指定文件，并对信号进行Prony分析计算

function [M, Amp, Fre, damp, Pha, main_f, main_damp, x, xc, y, Amp_Response, er, all, N, dt]=x_svdprony(x_in, dt, fL, showfigure)

format long;

%     a=zeros(1,20);

load 76sanxiangduanlu.MAT iX0006aBega;

x_in=iX0006aBega;%xlsread('shiyantu1.xlsx');

%x_in=iX0006aBega';

x = x_in;

cpu=cputime;

%N=size(x,1);

N=500;

%取N/2的整数部分为初始的P

P=floor(N/2);

%去直流环节

x_Sum = 0;

for i=1:N

x_Sum = x_Sum + x(i);

end

x_av = x_Sum / N;

if x_av > 10E-10

for i=1:N

x(i) = x(i) - x_av;

end

end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%滤波环节,滤波环节

%fL=1;

fL=2;

if fL>1

for i=fL+1:N

x(i-fL)=0;

for j=1:fL

x(i-fL) = x(i-fL)+(1/fL)*x(i-j+1);

end

end

end

N=N-fL;

tt=0:1:N-1;

%P要求为偶数

if mod(P, 2) ~= 0

P = P - 1;

end

P1=P;

if mod(P, 2) == 0

% Generate R，生