clear?all
close?all
clc

fid=fopen（‘1.txt‘‘r‘）;
[xcount]=fscanf（fid‘%f‘[2inf]）;???%讀入應力時程
[MN]=size（x）;
label=2;????%label表示分析哪一列應力
t1=x（11:N）;
subplot（311）
plot（t1x（label1:N））;
grid?on
title（‘應力時程曲線‘）;
xlabel（‘time（s）‘）;
ylabel（‘LU‘）;

%剔除非極值點
S（1）=x（21）;n=1;t（1）=t1（1）;
for?i=2:（N-1）%雨滴
????if?（x（labeli）>x（labeli-1））&&（x（labeli）>=x（labeli+1））%求極大值
????????n=n+1;S（n）=x（labeli）;t（n）=t1（i）;
????elseif?（x（labeli）????????n=n+1;S（n）=x（labeli）;t（n）=t1（i）;
????end
end
n=n+1;S（n）=x（labelN）;t（n）=t1（N）;
subplot（312）;
plot（tS）;
title（‘應力時程極值標記‘）;
xlabel（‘time（s）‘）;
ylabel（‘LU‘）;

n1=n;%儲存極值點個數(shù)
Se=0;%極值平均值
for?j=1:n1
????Se=Se+S（j）;
end
Se=Se/n1;
%n1記錄應力峰值點數(shù)?????
%n記錄剩余點的個數(shù)，當n<4時，循環(huán)結束
%i記錄四個點