Цикл Python

@legenda32 · Регистрация: 09.09.2014

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

Не получается сделать цикл. Идея такая, нужно посчитать q[N] и pogr в точках x0=от 0 до 0.9. Например в точке 0.1, 0.2 и тд.
Я попробовал сделать, но у меня получаются не те результаты. Если задать x0=0.9 и сделать все без цикла, то результат q[N] должен быть 1.219...

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
66
67
68
69
70
71
72
73
74
75
76
77
78
#!/usr/bin/python
import numpy as np 
import pylab as py
import sympy as sp
myeps=0.001
xn=sp.Symbol('xn')
xm=sp.Symbol('xm')
yn=sp.Symbol('yn')
ym=sp.Symbol('ym')
Phi=(1/(2*sp.pi))*sp.log(1/(sp.sqrt((xn-xm)**2+(yn-ym)**2)))
#pprint(Phi)
N=100;Rp=1; Rc=0.001; x0=0; y0=0; xc=x0+Rc; yc=y0; C1=1;
theta=np.linspace(2*np.pi, 0, N+1, endpoint=True)
x=np.array([Rp*np.sin(theta),Rp*np.cos(theta)])
 
 
 
def withoutZero(sFunct=None,minDist=myeps):
    def vectorFunct(A,B,C,D):
        gd=((A-B)**2+(C-D)**2)>minDist**2
        AA=A.copy(); AA[np.logical_not(gd)]+=1
        result=np.zeros_like(A)
        result[gd]=sp.lambdify((xn,xm,yn,ym),sFunct,modules='numpy')(AA,B,C,D)[gd]
        return result
    return vectorFunct
  
def Normal(x1):
    n=x1[0].shape[0]-1
    return np.array([-x1[1,1:]+x1[1,:n],x1[0,1:]-x1[0,:n]])
def last(sFunct=None):
    def vectorFunct(A,B,C,D):
       return sp.lambdify((xn,xm,yn,ym),sFunct,modules='numpy')(A,B,C,D)
    return vectorFunct
with open ('skv.txt', 'w') as skv:
    with open ('pogr.txt', 'w') as pog:
        for x0 in np.linspace(0,1,10, endpoint=False):
      h1=np.array([x[1,:N]/2.+x[1,1:]/2.,x[0,:N]/2.+x[0,1:]/2.])
      Xn,Xm=np.meshgrid(h1[1],h1[1])
      Yn,Ym=np.meshgrid(h1[0],h1[0])
      KF0=withoutZero(sp.diff(Phi,xn),myeps)(Xn,Xm,Yn,Ym)
      KF1=withoutZero(sp.diff(Phi,yn),myeps)(Xn,Xm,Yn,Ym)
      norm=Normal(x)
      Matrix=KF0*norm[0]+KF1*norm[1]
      Matrix=-2*Matrix+np.eye(N)
      XN=h1[1]
      YN=h1[0]
      Kf0=last(sp.diff(Phi,xn))(XN,xc,YN,yc)
      Kf1=last(sp.diff(Phi,yn))(XN,xc,YN,yc)
      stroka=Kf0*norm[0]+Kf1*norm[1]
      Phi0=1/(2*np.pi)*np.log(1/Rc)
      Phi2=2*(1/(2*np.pi)*np.log(1/(np.sqrt((h1[1]-x0)**2+(h1[0]-y0)**2))))
      Phi2=np.append(Phi2,Phi0)
      c = np.row_stack((Matrix,stroka))
      c = np.column_stack((c, Phi2))
      b = np.zeros(N)
      b = np.append(b,C1)
      q = np.linalg.solve(c, b)
      q0=(2*np.pi*C1)/np.log(Rp/Rc)
      pogr=(q[N]/q0-1)
      print q[N]
      #print q0
      #print pogr
      pog.write(str(pogr)+'\n')
      skv.write(str(x0)+'\n')
    
 
skv.close()
pog.close()
data_x = []
with open("skv.txt") as f:
   for line in f:
     data_x.append([float(x) for x in line.split()])
data_y = []
with open("pogr.txt") as f:
   for line in f:
     data_y.append([float(y) for y in line.split()])
#py.plot(data_x, data_y, 'r')
#py.show()

@passant · 29.05.2018, 15:42

Строка 36 и последующие. Что с отступами?? Что в теле цикла for ?

@legenda32 · 29.05.2018, 15:45 **[ТС]**

Сбились при копировании, скорее всего

@passant · 29.05.2018, 16:07

Не знаю, но так как приведено в посте - цикл for не выполняется ни разу!

@legenda32 · 29.05.2018, 16:18 **[ТС]**

Исправил отступы, но этот цикл все равно не работает не так как надо

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
66
67
68
69
70
71
72
73
74
75
76
77
78
#!/usr/bin/python
import numpy as np 
import pylab as py
import sympy as sp
myeps=0.001
xn=sp.Symbol('xn')
xm=sp.Symbol('xm')
yn=sp.Symbol('yn')
ym=sp.Symbol('ym')
Phi=(1/(2*sp.pi))*sp.log(1/(sp.sqrt((xn-xm)**2+(yn-ym)**2)))
#pprint(Phi)
N=100;Rp=1; Rc=0.001; x0=0; y0=0; xc=x0+Rc; yc=y0; C1=1;
theta=np.linspace(2*np.pi, 0, N+1, endpoint=True)
x=np.array([Rp*np.sin(theta),Rp*np.cos(theta)])
 
 
 
def withoutZero(sFunct=None,minDist=myeps):
    def vectorFunct(A,B,C,D):
        gd=((A-B)**2+(C-D)**2)>minDist**2
        AA=A.copy(); AA[np.logical_not(gd)]+=1
        result=np.zeros_like(A)
        result[gd]=sp.lambdify((xn,xm,yn,ym),sFunct,modules='numpy')(AA,B,C,D)[gd]
        return result
    return vectorFunct
  
def Normal(x1):
    n=x1[0].shape[0]-1
    return np.array([-x1[1,1:]+x1[1,:n],x1[0,1:]-x1[0,:n]])
def last(sFunct=None):
    def vectorFunct(A,B,C,D):
       return sp.lambdify((xn,xm,yn,ym),sFunct,modules='numpy')(A,B,C,D)
    return vectorFunct
with open ('skv.txt', 'w') as skv:
    with open ('pogr.txt', 'w') as pog:
        for x0 in np.linspace(0,1,10, endpoint=False):
           h1=np.array([x[1,:N]/2.+x[1,1:]/2.,x[0,:N]/2.+x[0,1:]/2.])
           Xn,Xm=np.meshgrid(h1[1],h1[1])
           Yn,Ym=np.meshgrid(h1[0],h1[0])
           KF0=withoutZero(sp.diff(Phi,xn),myeps)(Xn,Xm,Yn,Ym)
           KF1=withoutZero(sp.diff(Phi,yn),myeps)(Xn,Xm,Yn,Ym)
           norm=Normal(x)
           Matrix=KF0*norm[0]+KF1*norm[1]
           Matrix=-2*Matrix+np.eye(N)
           XN=h1[1]
           YN=h1[0]
           Kf0=last(sp.diff(Phi,xn))(XN,xc,YN,yc)
           Kf1=last(sp.diff(Phi,yn))(XN,xc,YN,yc)
           stroka=Kf0*norm[0]+Kf1*norm[1]
           Phi0=1/(2*np.pi)*np.log(1/Rc)
           Phi2=2*(1/(2*np.pi)*np.log(1/(np.sqrt((h1[1]-x0)**2+(h1[0]-y0)**2))))
           Phi2=np.append(Phi2,Phi0)
           c = np.row_stack((Matrix,stroka))
           c = np.column_stack((c, Phi2))
           b = np.zeros(N)
           b = np.append(b,C1)
           q = np.linalg.solve(c, b)
           q0=(2*np.pi*C1)/np.log(Rp/Rc)
           pogr=(q[N]/q0-1)
           print q[N]
           #print q0
           #print pogr
           pog.write(str(pogr)+'\n')
           skv.write(str(x0)+'\n')
    
 
skv.close()
pog.close()
data_x = []
with open("skv.txt") as f:
   for line in f:
     data_x.append([float(x) for x in line.split()])
data_y = []
with open("pogr.txt") as f:
   for line in f:
     data_y.append([float(y) for y in line.split()])
#py.plot(data_x, data_y, 'r')
#py.show()

@passant · 30.05.2018, 14:40

А "как надо"? Сие известно только тому, кто эту программу писал. С точки зрения остальных программа работает в точности так, как это описано в коде и никак иначе.

@legenda32 · 30.05.2018, 14:52 **[ТС]**

Как она должна работать по моей задумке, я описал выше, если в ручную перебирать значения x0, то все работает норм, но когда я делаю это в цикле, то прога выдает совсем другое

@passant · 30.05.2018, 14:55

Дебагером пользовались? Изменение переменных по шагам проследили?
Что показала отладка? Где конкретно возникло расхождение ожидаемых и полученных результатов? На каком шаге цикла, в какой операции?

но когда я делаю это в цикле, то прога выдает совсем другое

- подсказка: внутри цикла у вас ошибка в формулах

@legenda32 · 30.05.2018, 17:54 **[ТС]**

Не знаком со всеми этими вещами, я пробовал через np.linspace разбить x0 от 0 до 0.9, потом задавал цикл for x0 in x0.flat, но все равно результат не верен, хотя х0 вроде бы принимает те значения которые я ему даю

@passant · 31.05.2018, 09:09

Повторю - ошибка в ваших формулах.
Если решили заняться программированием - дебегингом придется заниматься, на каком языке вы бы не программировали. Благо на Python это просто и даже приятно.
По-русски это называется "процесс отладки". Отладка - это 80% работы программиста. Делать это за вас врят-ли кто захочет - мало того, что надо вникать в ваш не шибко красивый код, искать где вместо плюса вы минус поставили, где отсупы неправильно вывели, где лишнюю запятую втулили. Так еще и в вашей задаче разбираться. В общем - самое время начать осваивать.
Удачи.

Новые блоги и статьи Все статьи Все блоги /
Программный контроль заполнения реквизита табличной части документа Maks 02.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "СписаниеМатериалов", разработанного в конфигурации КА2. Задача: реализовать контроль заполнения реквизита табличной части. . .	wmic не является внутренней или внешней командой Maks 02.04.2026 Решение: DISM / Online / Add-Capability / CapabilityName:WMIC~~~~ Отсюда: https:/ / winitpro. ru/ index. php/ 2025/ 02/ 14/ komanda-wmic-ne-naydena/	Программная установка даты и запрет ее изменения Maks 02.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "СписаниеМатериалов", разработанного в конфигурации КА2. Задача: при создании документов установить период списания автоматически. . .	Вывод данных в справочнике через динамический список Maks 01.04.2026 Реализация из решения ниже выполнена на примере нетипового справочника "Спецтехника" разработанного в конфигурации КА2. Задача: вывести данные из ТЧ нетипового документа. . .
Функция заполнения текстового поля в реквизите формы документа Maks 01.04.2026 Алгоритм из решения ниже реализован на нетиповом документе "ВыдачаОборудованияНаСпецтехнику" разработанного в конфигурации КА2, в дополнении к предыдущему решению. На форме документа создается. . .	К слову об оптимизации kumehtar 01.04.2026 Вспоминаю начало 2000-х, университет, когда я писал на Delphi. Тогда среди программистов на форумах активно обсуждали аккуратную работу с памятью: нужно было следить за переменными, вовремя. . .	Идея фильтра интернета (сервер = слой+фильтр). Hrethgir 31.03.2026 Суть идеи заключается в том, чтобы запустить свой сервер, о чём я если честно мечтал давно и давно приобрёл книгу как это сделать. Но не было причин его запускать. Очумелые учёные напечатали на. . .	Модель здравосоХранения 6. ESG-повестка и устойчивое развитие; углублённый анализ кадрового бренда anaschu 31.03.2026 В прикрепленном документе раздумья о том, как можно поменять модель в будущем

@passant 578 / 411 / 69 Регистрация: 09.01.2018 Сообщений: 1,363
	29.05.2018, 15:42
	Строка 36 и последующие. Что с отступами?? Что в теле цикла for ? 0

@legenda32 1 / 1 / 1 Регистрация: 09.09.2014 Сообщений: 39
	29.05.2018, 15:45 [ТС]
	Сбились при копировании, скорее всего 0

@passant 578 / 411 / 69 Регистрация: 09.01.2018 Сообщений: 1,363
	29.05.2018, 16:07
	Не знаю, но так как приведено в посте - цикл for не выполняется ни разу! 0

@passant 578 / 411 / 69 Регистрация: 09.01.2018 Сообщений: 1,363
	30.05.2018, 14:40
	А "как надо"? Сие известно только тому, кто эту программу писал. С точки зрения остальных программа работает в точности так, как это описано в коде и никак иначе. 0

@legenda32 1 / 1 / 1 Регистрация: 09.09.2014 Сообщений: 39
	30.05.2018, 14:52 [ТС]
	Как она должна работать по моей задумке, я описал выше, если в ручную перебирать значения x0, то все работает норм, но когда я делаю это в цикле, то прога выдает совсем другое 0

@legenda32 1 / 1 / 1 Регистрация: 09.09.2014 Сообщений: 39
	30.05.2018, 17:54 [ТС]
	Не знаком со всеми этими вещами, я пробовал через np.linspace разбить x0 от 0 до 0.9, потом задавал цикл for x0 in x0.flat, но все равно результат не верен, хотя х0 вроде бы принимает те значения которые я ему даю 0

@passant 578 / 411 / 69 Регистрация: 09.01.2018 Сообщений: 1,363
	31.05.2018, 09:09
	Повторю - ошибка в ваших формулах. Если решили заняться программированием - дебегингом придется заниматься, на каком языке вы бы не программировали. Благо на Python это просто и даже приятно. По-русски это называется "процесс отладки". Отладка - это 80% работы программиста. Делать это за вас врят-ли кто захочет - мало того, что надо вникать в ваш не шибко красивый код, искать где вместо плюса вы минус поставили, где отсупы неправильно вывели, где лишнюю запятую втулили. Так еще и в вашей задаче разбираться. В общем - самое время начать осваивать. Удачи. 0