Расчет статических характеристик объекта с помощью модели динамики

N.

G1.

G2.

P2.

P3.

K.

H.

2.6.

2.4.

1.4.

3.1.

7.71 068.

2.2.

Файл-сценарий:

clc.

clear.

close all.

global G1 G2 k P2 P3.

G1=2.6; G2=2.4; P2=1.4; P3=3.1; k=7.71 068;

nach=[0 2.2]; H₀=0;

n=10;del=0.1*k;k1=k/2;

for i=1:n.

if i==1.

k=k1;

[t, H]=ode45(@fun3,nach, H₀);

h (1)=real (H (end));K (1)=k;

figure (1).

plot (t, H).

title ('H=f (t)').

xlabel ('t');ylabel ('H');

grid on;

else.

k=k1+del*i;

[t, H]=ode45(@fun3,nach, H (end));

h (i)=real (H (end));

K (i)=k;

% figure (i).

% plot (t, H).

end.

end.

figure (2).

plot (t, H).

title ('H=f (t)').

xlabel ('t');ylabel ('H');

grid on.

figure (3).

plot (K, h).

title ('H=f (k)').

xlabel ('k');ylabel ('H');

grid on.

k=7.71 068;

G₁=G½;h=[];

for i=1:n.

if i==1.

G1=G₁;

[t, H]=ode45(@fun3,nach, H₀);

h (1)=real (H (end));G3(1)=G1;

else.

G1=G₁+del*i;

[t, H]=ode45(@fun3,nach, H (end));

h (i)=real (H (end));

G3(i)=G1;

end.

end.

figure (4).

plot (G3,h).

title ('H=f (G1)').

xlabel ('G1');ylabel ('H');

grid on.

hold on.

G1=2.6;

G₂=G2/2;G=[]; h=[];

for i=1:n.

if i==1.

G2=G₂;

[t, H]=ode45(@fun3,nach, H₀);

h (1)=real (H (end));G (1)=G2;

else.

G2=G₂+del*i;

[t, H]=ode45(@fun3,nach, H (end));

h (i)=real (H (end));

G (i)=G2;

end.

end.

figure (5).

plot (G, h).

title ('H=f (G2)').

xlabel ('G2');ylabel ('H');

grid on.

Файл-функция:

function dydt = fun3(t, H).

global k.

global G2.

global G1.

global P2.

global P3.

G3=k*sqrt (H+P2-P3);

dydt=[G1+G2-G3];

end.