prony prony算法在电力系统的低频振荡中的应用

prony算法在电力系统的低频振荡中的应用，分析，判断，给出频率相位等

%prony算法模态参数识别
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear
clc
format long
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
fni=input（‘prony算法模态识别-输入数据文件名:‘‘s‘）;
fid=fopen （fni‘r‘）;               %打开只读文件
mn=fscanf （fid‘%d‘1）;            %模态阶数
%定义输入实测数据类型：
%ig=1时时域数据如冲击响应，自由振动，相互关函数，随机减量法处理结果；
%ig=2频域数据如频响函数实部虚部数据
ig=fscanf （fid‘%d‘1）;
%ig=1时f为采样频率sf，ig=2时f为采样间隔df
f=fscanf （fid‘%f‘1）;
fno=fscanf （fid‘%s‘1）;           %输出数据文件名
b=fscanf （fid‘%f‘[iginf]）;      %实测时域或频域数据
status=fclose （fid）;              
nm=2*mn;                           %建立复指数函数的阶数（为模态阶数的两倍）

%组织识别计算所用的时域数据及参数
if ig==1                           %实测时域数据
    sf=f;                          %取采样频率
    n=fix（length（b）/2）;            %取时域数据二分之一的长度
    h=b（

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       3995  2011-04-13 22:43  prony算法\prony.m

     文件     242176  2011-04-18 19:10  prony算法\扩展Prony方法_清华大学《现代信号处理》讲义_-张贤达.ppt

     目录          0  2011-04-19 20:34  prony算法

----------- ---------  ---------- -----  ----

               246171                    3