ngth (x_pl1) (% g , x_pl1 (index)) ( n )=- A ; (mn, T/2); _ min2=pc;=- A; (mn, T/2);=k + 1;
x_pl2=pc; pf == 1
for index=1: 1: length (x_pl2) ( %graquo;,x_pl2(index))(laquo;nraquo;)=clock;c1==c2=0;=time(c1,c2);=c(1)*31104000+c(2)*2592000+c(3)*86400+c(4)*3600+c(5)*60+c(6)*1;((tgt;=max_rtd)amp;(kgt;=min_pn)amp;(k*Tgt;=min_mtd))=1;((sum(equality(x_pl2,x_pl1,m,n))~=length(equality(x_pl2,x_pl1,m,n)))|...
(sum(equality(x_min2,x_min1,m,n))~=length(equality(x_min2,x_min1,m,n))))amp;(F==0)_min1=x_min2;_pl1=x_pl2;=-A;(mn,T/2);_min2=pc;=-A;(mn,T/2);=k+1;
x_pl2=pc; pf == 1
for index=1: 1: length (x_pl2) ( %graquo;,x_pl2(index))(laquo;nraquo;)=clock;c1==c2=0;=time(c1,c2);=c(1)*31104000+c(2)*2592000+c(3)*86400+c(4)*3600+c(5)*60+c(6)*1;
if ((t gt;=max_rtd) amp; (k gt;=min_pn) amp; (k * T gt;=min_mtd))=1;=- A; F == 0pf == 1 ( n ) ( Phase hodograph point is defined successfully n ) ( Passed period number is% g , k) ( n ) ( Passed model time is% g , T * k) ( s n ) ( Passed real time is% g , t) ( s n )=(x_pl2 + x_min2) ./ 2;
x_pl2=x_pl2-d; _pl2 (length (x_pl2) +1)=T; _pl2 (length (x_pl2) +1)=1;
if pf == 1 ( n ) ( Warning: n ) ( Function PH_point=f (mn, A, T, sv, pc, m, n, min_pn, min_mtd, max_rtd, df, pf) n ) ( failed to find R-characteristic in period T =% g , T) ( s n ) ( You should increase max_rtd and parameters Relative tolerance raquo ; n ) ( and Absolute tolerance in Simulation parametres to get correct n ) ( result here. ) ( Also don t set m and n too much. n ) ( n ) ( Passed period number is% g , k) ( n ) ( Passed model time is% g , T * k) ( s n ) ( Passed real time is% g , t) ( s n )
df == 1_pl2=(x_pl1 + x_pl2) ./ 2; _min2=(x_min1 + x_min2) ./ 2;=(x_pl2 + x_min2) ./ 2; _pl2=x_pl2-d; _pl2(length(x_pl2)+1)=T;_pl2(length(x_pl2)+1)=2;_pl2=x_pl2.*NaN;_pl2(length(x_pl2)+1)=T;_pl2(length(x_pl2)+1)=3;
_point=x_pl2;
else ( n )
fprintf ( Error: n ) ( Incorrect parameters in function n ) ( PH_point=f (mn, A, T, sv, pc, m, n, min_pn, min_mtd, max_rtd , df, pf) n ) _ point=NaN;
-функція Script
% Script.m m-файл для визначення всього фазового годографа
global A sv pc ( off raquo ;, laquo;allraquo;);(laquo;ph_point.txtraquo;)on;=laquo;hydraulic_driveraquo;;=27;=5;=2.99;=0.01;=Ts;=[0 0 0 0 0 0 0 0 0 0];=0;=[1e - 6 1e - 2 1e - 12 1e - 8 1e - 6 1e - 1 вересня 1e - 2 1e - 3 1e - 9 1e - 5 1e - 15 ];=[5 5 5 5 5 5 5 5 5 5 5 5];_pn=5;_mtd=6;_rtd=10;=1;=1;=1;Tgt;=Tf=PH_point(mn,A,T,sv,pc,m,n,min_pn,min_mtd,max_rtd,df,pf);j=1:1:length(a)(i,j)=a(j);=T-dT;=i+1;
