График функции Вейерштрасса

@Фырфа · Регистрация: 28.02.2019

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

Добрый день.
Необходимо построить график функции Вейерштрасса в питоне, используя numpy и matplotlib. Возникла загвоздка с выражением самой функции: как ее описать в питоне. Получилось как-то так:

Python
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
x_data = np.arange(-2, 2.01, 0.001)
n = 100
b = 0.5
a = 3
y_data = np.array([])
y = np.array([])
for x in x_data:
    for i in range(n):
        y = np.append(y, np.cos(np.pi*x*a**i)*b**i)
    y_sum = np.sum(y)
    y_data = np.append(y_data, y_sum)
    y = np.array([])
 
plt.plot(x_data, y_data)
plt.axis('equal')
plt.grid = True
plt.show()

Мне такой вариант не очень нравится. Хотелось бы узнать, есть ли какой-то способ описать функцию суммы бесконечного ряда в питоне? И есть ли какой-то другой способ построить график этой функции?
Например, здесь при n=1000 вылетает ошибка, что переменная Int слишком большая, чтобы конвертировать во float.

@Gdez · 08.11.2020, 10:44

Фырфа,

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
import numpy as np
from math import pi, cos, log
 
def power(a, n):
    return (1 if n == 0
            else power(a * a, n // 2) if n % 2 == 0
            else a * power(a, n - 1))
 
x_data = np.arange(-2, 2.01, 0.01)
n = 100
b = 0.9
eps = 1e-4
print('b^n = {:5.3e}'.format(power(b,n)))
m = round(log(eps) / log(b)) + 1
print('достаточное количество членов ряда : ',m)
n = m
a = 3
y_data = np.array([])
for x in x_data:
    y = 0
    for i in range(n):
        ai = power(a,i) * x * pi
        bi = power(b,i)    
        y += cos(ai)*bi
    y_data = np.append(y_data, y)
#print(x_data,y_data)
print("end")

Новые блоги и статьи Все статьи Все блоги /
Знаешь почему 90% людей редко бывают счастливыми? kumehtar 14.04.2026 Потому что они ждут. Ждут выходных, ждут отпуска, ждут удачного момента. . . а удачный момент так и не приходит.	Фиксация колонок в отчете СКД Maks 14.04.2026 Фиксация колонок в СКД отчета типа Таблица. Задача: зафиксировать три левых колонки в отчете. Процедура ПриКомпоновкеРезультата(ДокументРезультат, ДанныеРасшифровки, СтандартнаяОбработка) / / . . .	Настройки VS Code Loafer 13.04.2026 { "cmake. configureOnOpen": false, "diffEditor. ignoreTrimWhitespace": true, "editor. guides. bracketPairs": "active", "extensions. ignoreRecommendations": true, . . .	Оптимизация кода на разграничение прав доступа к элементам формы Maks 13.04.2026 Алгоритм из решения ниже реализован на нетиповом документе, разработанного в конфигурации КА2. Задачи, как таковой, поставлено не было, проделанное ниже исключительно моя инициатива. Было так:. . .
Контроль заполнения и очистка дат в зависимости от значения перечислений Maks 12.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "ПланированиеПерсонала", разработанного в конфигурации КА2. Задача: реализовать контроль корректности заполнения дат назначения. . .	Архитектура слоя интернета для сервера-слоя. Hrethgir 11.04.2026 В продолжение https:/ / www. cyberforum. ru/ blogs/ 223907/ 10860. html Знаешь что я подумал? Раз мы все источники пишем в голове ветки, то ничего не мешает добавить в голову такой источник, который сам. . .	Подстановка значения реквизита справочника в табличную часть документа Maks 10.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "ПланированиеПерсонала", разработанного в конфигурации КА2. Задача: при выборе сотрудника (справочник Сотрудники) в ТЧ документа. . .	Очистка реквизитов документа при копировании Maks 09.04.2026 Алгоритм из решения ниже применим как для типовых, так и для нетиповых документов на самых различных конфигурациях. Задача: при копировании документа очищать определенные реквизиты и табличную. . .