MATLAB高阶系统参数辨识程序

%高阶惯性辨识  调用point函数
clear all;
close all;
clc;
[FileName] = uigetfile（‘*.mat‘‘Select the Mat-file‘）%弹出选择窗口       新技能
load（FileName）;
plot（ty‘y‘‘linewidth‘4）
hold on
[lpm]=size（y）;
if m>lp lp=m;end
dt=t（3）-t（2）;
sum=0;num=0;
for i=lp-20:lp;
    sum=sum+y（i）;
    num=num+1;
end
sum=sum/num;
ys=sum;
k=ys        %定K
[e1post1]=point（0.4yslpy）;
[e2post2]=point（0.8yslpy）;
n=round（（1.075*post1/（post2-post1）+0.5）^2）%四舍五入
T=（post1+post2）*dt/2.16/n
%辨识完成，画图比较
a=exp（-dt/T）;
b=1-a;
y1=0;
x（1:n）=0;
for i=2:lp
    x（1）=a*x（1）+k*b*1;
    x（2:n）=a*x（2:n）+b*x（1:n-1）;
    y1（i）=x（n）;
end
n
plot（ty1‘k‘‘linewidth‘1）

    
        tle>鏈皢瀵硅薄寮曠敤璁剧疆鍒板璞＄殑瀹炰緥銆?/title>
        ta name=“viewport“ content=“width=device-width“ />