ошибки в коде, как исправить

@Invisability · Регистрация: 20.12.2019

Студворк — интернет-сервис помощи студентам

Python
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
import numpy as np
import matplotlib.pyplot as plt
from sympy import *
import matplotlib as mpl
mpl.rcParams['font.family'] = 'fantasy'
mpl.rcParams['font.fantasy'] = 'Comic Sans MS, Arial'
T1=5;T2=4;T3=4#постооянные времени звеньев
K1= 1.5; K2= 1; K3= 1# коэффициенты усиления звеньев
dt= 0.01 #Шаг моделирования
N =10000 #Число точек кривой разгона
?st = 84.97 #Статическое значение температуры нагреваемой воды
Gst = 470.1#Статическое значение расхода греющей воды
Gmax =0.02* Gst #Максимальное отклонение расхода 
?dop =0.09* ?st #Допустимое отклонение температуры: 
x=[0 for w in  np.arange(0,N)]# начальное заполнение списка x
y1=[0 for w in  np.arange(0,N)]#начальное заполнение списка y1
y2=[0 for w in  np.arange(0,N)]#начальное заполнение списка y2
y3=[0 for w in  np.arange(0,N)]#начальное заполнение списка y3
mu=[0 for w in  np.arange(0,N)]#начальное заполнение списка mu
ep=[0 for w in  np.arange(0,N)]#начальное заполнение списка ep
ii=[0 for w in  np.arange(0,N)]#начальное заполнение списка ii
"""Переходные процессы в АСР при ступенчатом возмущении Gmax
по каналу регулирующего воздействия"""
def PP(Kp,Ki):
         def fa(K,T,dt,x,y):#Разностное уравнение  кривой разгона
                  return (1-dt/T)*y+K*x*(dt/T)
         def fi(K,dt,x,y):
                  return dt*K*x+y#Разностное уравнение И-звена
         def fp(K,x):                
                  return K*x#Разностное уравнение П-звена
         for i in np.arange(0,N):# численное решение разностных уравнений
                  if i+1>N-1:break
                  x[i+1]=mu[i]+Gmax
                  y1[i+1]=fa(K1,T1,dt,x[i],y1[i])
                  y2[i+1]=fa(K2,T2,dt,y1[i],y2[i])
                  y3[i+1]=fa(K3,T3,dt,y2[i],y3[i])
                  ep[i+1]=-y3[i]
                  ii[i+1]=fi(Ki,dt,ep[i],ii[i])
                  mu[i+1]= fp(Kp,ep[i])+ii[i+1]
         return y3
""" Оптимальные настройки регуляторов"""
KpП= 1.299; KiП= 0
KpИ= 0; KiИ= 0.053
KpПИ = 0.734; KiПИ = 0.105
Kpгр = 2.366; Kiгр = 0.352
""" Учёт ?st """
VПИ=[w+?st for w in PP(KpПИ,KiПИ)]
VП=[w+?st for w in PP(KpП,KiП)]
VИ=[w+?st for w in PP(KpИ,KiИ)]
Vгр=[w+?st for w in PP(Kpгр,Kiгр)]
V=[w+?st for w in PP(0,0)]
u=[?st+?dop for i in  np.arange(0,N)]
t=[dt*i for i in  np.arange(0,N)]
plt.title('Переходные процессы в АСР при ступенчатом возмущении ') 
plt.ylabel('VПИ,VП,VИ,Vгр,u') #Метка по оси y 
plt.xlabel('Время - t') #Метка по оси x 
plt.grid(True) #Сетка
plt.plot(t,VП,linewidth=2,label='АСР с П-регулятором')
plt.plot(t,VИ,linewidth=2,label='АСР с И-регулятором')
plt.plot(t,VПИ,linewidth=2,label='АСР с ПИ-регулятором') 
plt.plot(t,Vгр,linewidth=2,label='АСР на границе устойчивости') 
plt.plot(t,u,linewidth=2,label='Допустимое отклонение температуры') 
plt.plot(t,V,linewidth=2,label='Кривая разгона объекта регулирования ')
plt.legend(loc='best')
plt.show() #Показать график

Помогите исправить ошибки пожалуйста

Добавлено через 1 минуту

Python
1
2
3
4
5
6
7
VПИ=[w+?st for w in PP(KpПИ,KiПИ)]
VП=[w+?st for w in PP(KpП,KiП)]
VИ=[w+?st for w in PP(KpИ,KiИ)]
Vгр=[w+?st for w in PP(Kpгр,Kiгр)]
V=[w+?st for w in PP(0,0)]
u=[?st+?dop for i in  np.arange(0,N)]
t=[dt*i for i in  np.arange(0,N)]

не понятно что это за часть кода и что с ней делать

@Gdez · 20.12.2021, 14:52

Invisability,

Python
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
import numpy as np
 
T1=5;T2=4;T3=4#постоянные времени звеньев
K1= 1.5; K2= 1; K3= 1# коэффициенты усиления звеньев
dt= 0.01 #Шаг моделирования
N =10000 #Число точек кривой разгона
st = 84.97 #Статическое значение температуры нагреваемой воды
Gst = 470.1#Статическое значение расхода греющей воды
Gmax =0.02 * Gst #Максимальное отклонение расхода 
dop =0.09 * st #Допустимое отклонение температуры: 
 
"""
x=np.zeros(N)# начальное заполнение списка x
y1=np.zeros(N)#начальное заполнение списка y1
y2=np.zeros(N)#начальное заполнение списка y2
y3=np.zeros(N)#начальное заполнение списка y3
mu=np.zeros(N)#начальное заполнение списка mu
ep=np.zeros(N)#начальное заполнение списка ep
ii=np.zeros(N)#начальное заполнение списка ii
"""
x, y1, y2, y3, mu, ep, ii = np.atleast_1d(np.zeros((7, N)))
"""Переходные процессы в АСР при ступенчатом возмущении Gmax
по каналу регулирующего воздействия"""
def PP(Kp,Ki):
         def fa(K,T,dt,x,y):#Разностное уравнение  кривой разгона
                  return (1-dt/T)*y+K*x*(dt/T)
         def fi(K,dt,x,y):
                  return dt*K*x+y#Разностное уравнение И-звена
         def fp(K,x):                
                  return K*x#Разностное уравнение П-звена
         for i in range(N):# численное решение разностных уравнений
                  if i+1>N-1:break
                  x[i+1]=mu[i]+Gmax
                  y1[i+1]=fa(K1,T1,dt,x[i],y1[i])
                  y2[i+1]=fa(K2,T2,dt,y1[i],y2[i])
                  y3[i+1]=fa(K3,T3,dt,y2[i],y3[i])
                  ep[i+1]=-y3[i]
                  ii[i+1]=fi(Ki,dt,ep[i],ii[i])
                  mu[i+1]= fp(Kp,ep[i])+ii[i+1]
         return y3 + st
""" Оптимальные настройки регуляторов"""
KpП= 1.299; KiП= 0
KpИ= 0; KiИ= 0.053
KpПИ = 0.734; KiПИ = 0.105
Kpгр = 2.366; Kiгр = 0.352
""" Учёт ?st """
"""
VПИ=[w + st for w in PP(KpПИ,KiПИ)]
VП=[w + st for w in PP(KpП,KiП)]
VИ=[w + st for w in PP(KpИ,KiИ)]
Vгр=[w + st for w in PP(Kpгр,Kiгр)]
V=[w + st for w in PP(0,0)]
u=[st + dop for i in  np.arange(0,N)]
t=[dt*i for i in  np.arange(0,N)]
"""
VПИ=PP(KpПИ,KiПИ)
VП=PP(KpП,KiП)
VИ=PP(KpИ,KiИ)
Vгр=PP(Kpгр,Kiгр)
V=PP(0,0)
u=np.full(N, st + dop)
t=dt*np.arange(N)
...
...

Новые блоги и статьи Все статьи Все блоги /
[Owen Logic] Поддержание уровня воды в резервуаре количеством включённых насосов: моделирование и выбор регулятора ФедосеевПавел 14.03.2026 Поддержание уровня воды в резервуаре количеством включённых насосов: моделирование и выбор регулятора ВВЕДЕНИЕ Выполняя задание на управление насосной группой заполнения резервуара,. . .	делаю науч статью по влиянию грибов на сукцессию anaschu 13.03.2026 прикрепляю статью	SDL3 для Desktop (MinGW): Создаём пустое окно с нуля для 2D-графики на SDL3, Си и C++ 8Observer8 10.03.2026 Содержание блога Финальные проекты на Си и на C++: hello-sdl3-c. zip hello-sdl3-cpp. zip Результат:	Установка CMake и MinGW 13.1 для сборки С и C++ приложений из консоли и из Qt Creator в EXE 8Observer8 10.03.2026 Содержание блога MinGW - это коллекция инструментов для сборки приложений в EXE. CMake - это система сборки приложений. Здесь описаны базовые шаги для старта программирования с помощью CMake и. . .
Как дизайн сайта влияет на конверсию: 7 решений, которые реально повышают заявки Neotwalker 08.03.2026 Многие до сих пор воспринимают дизайн сайта как “красивую оболочку”. На практике всё иначе: дизайн напрямую влияет на то, оставит человек заявку или уйдёт через несколько секунд. Даже если у вас. . .	Модульная разработка через nuget packages DevAlt 07.03.2026 Сложившийся в . Net-среде способ разработки чаще всего предполагает монорепозиторий в котором находятся все исходники. При создании нового решения, мы просто добавляем нужные проекты и имеем. . .	Модульный подход на примере F# DevAlt 06.03.2026 В блоге дяди Боба наткнулся на такое определение: В этой книге («Подход, основанный на вариантах использования») Ивар утверждает, что архитектура программного обеспечения — это структуры,. . .	Управление камерой с помощью скрипта OrbitControls.js на Three.js: Вращение, зум и панорамирование 8Observer8 05.03.2026 Содержание блога Финальная демка в браузере работает на Desktop и мобильных браузерах. Итоговый код: orbit-controls-threejs-js. zip. Сканируйте QR-код на мобильном. Вращайте камеру одним пальцем,. . .