%高階慣性辨識??調用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“?/>