1、Matlab程序:
s1=1/30:1/30:1;
s2=1:-1/40:-0、5;
s3=-0、5:0、05:-0、1;
s=[s1 s2 s3];
plot(s);
s_a=s;
for i=1:49
s_a=[s_a s];
end
figure();plot(s_a);
ss1=0:0、03/220:0、03;
ss2=0、03:1/2000:1;
ss3=1:-1/2500:0、03;
ss4=0、03:-0、03/411:0;
ss=[ss1 ss2 ss3 ss4];
figure();plot(ss);
ss_mu=s_a、*ss;
2、figure();plot(ss_mu);
p_press=200:-2/50:2/50;
figure();plot(p_press);
max=0;
k=0;
for i=1:5000
if(ss_mu(i)>max)
max=ss_mu(i);
k=i;
end
end
m=0;
min=10;
for i=11:k
f=0;
if min>abs(ss_mu(i)-0、55*ss_mu(k))
f=1;
for j=-10:10
i
3、f(ss_mu(i+j)>ss_mu(i))
f=0;
end
end
if f==1
min=abs(ss_mu(i)-0、55*ss_mu(k))
m=1;
end
end
end
n=0;
min=10;
for i=k:5000
f=0;
if min > abs(ss_mu(i)-0、82*ss_mu(k))
f=1;
for j=-10:10
if(ss_mu(i+j)>ss_mu(i))
f=0;
end
end
if f==1
min= abs(ss_mu(i)-0、82*ss_mu(k))
n=1;
end
end
end
OM=ss_mu(k)
OS=ss_mu(m)
OD=ss_mu(n)
PM=p_press(k)
PS=p_press(m)
PD=p_press(n)