跳转到内容

File:Dirichlet kernels.svg

页面内容不支持其他语言。
這個文件來自維基共享資源
维基百科,自由的百科全书

原始文件 (SVG文件,尺寸为720 × 495像素,文件大小:57 KB)


摘要

描述
English: Plots of the first few Dirichlet kernels made with matplotlib. This function has a period of 2π.
日期
来源 Dirichlet.png
作者
其他版本
这是一张修改过的图片,这意味着它已在原版本的基础上通过软件进行了编辑,改动内容:vectorized。其原始版本为:Dirichlet.png。修改者:Azaghal of Belegost

许可协议

我,本作品著作权人,特此采用以下许可协议发表本作品:
w:zh:知识共享
署名 相同方式共享
本文件采用知识共享署名-相同方式共享 4.0 国际许可协议授权。
您可以自由地:
  • 共享 – 复制、发行并传播本作品
  • 修改 – 改编作品
惟须遵守下列条件:
  • 署名 – 您必须对作品进行署名,提供授权条款的链接,并说明是否对原始内容进行了更改。您可以用任何合理的方式来署名,但不得以任何方式表明许可人认可您或您的使用。
  • 相同方式共享 – 如果您再混合、转换或者基于本作品进行创作,您必须以与原先许可协议相同或相兼容的许可协议分发您贡献的作品。

Source code

Source Code in python:
from __future__ import division # In case someone wants to use Python 2
import numpy as np
import matplotlib as mpl
import matplotlib.pyplot as plt
import seaborn as sns # Just for the colors

mpl.rcParams['axes.labelsize'] = 20
mpl.rcParams['axes.titlesize'] = 22
mpl.rcParams['legend.fontsize'] = 18
mpl.rcParams['xtick.labelsize'] = 14
mpl.rcParams['ytick.labelsize'] = 14
mpl.rcParams['lines.linewidth'] = 3

def dirichlet(n, x):
    return np.where(x!=0, np.sin((n + .5)*x)/np.sin(.5*x), 2*n + 1)

points = 2**10 + 1
kernels = 4
x_min = -np.pi
x_max = np.pi

x_vals = np.linspace(x_min, x_max, num=points)
n_vals = np.arange(1, kernels + 1)
D = dirichlet(np.resize(n_vals, (points, kernels)).T,
              np.resize(x_vals, (kernels, points)))

fig = plt.figure()
ax = fig.add_subplot(1, 1, 1)
for kernel, n in zip(D, n_vals):
    ax.plot(x_vals, kernel,
            label=r'$\frac{\sin((' + str(n) + r'+1/2)x)}{\sin(x/2)}$')

ax.xaxis.set_ticks(np.pi*np.arange(-4, 5)/4)
ax.xaxis.set_ticklabels([r'$-\pi$', r'$-\frac{3\pi}{4}$', r'$-\frac{\pi}{2}$',
                         r'$-\frac{\pi}{4}$', r'$0$', r'$\frac{\pi}{4}$',
                         r'$\frac{\pi}{2}$', r'$\frac{3\pi}{4}$', r'$\pi$'])
ax.legend()
ax.set_xlim([x_min, x_max])
ax.set_xlabel(r'$x$')
ax.set_ylabel(r'$D_n(x)$')
ax.xaxis.grid(color='gray')
ax.yaxis.grid(color='gray')
plt.tight_layout()
plt.savefig('dirichlet.svg')

说明

添加一行文字以描述该文件所表现的内容

此文件中描述的项目

描繪內容

文件历史

点击某个日期/时间查看对应时刻的文件。

日期/时间缩⁠略⁠图大小用户备注
当前2014年11月21日 (五) 21:382014年11月21日 (五) 21:38版本的缩略图720 × 495(57 KB)Azaghal of BelegostUser created page with UploadWizard

以下页面使用本文件:

全域文件用途

以下其他wiki使用此文件:

元数据