小编dus*_*tin的帖子

如何以非常高的质量保存python图？

也就是说,当我继续放大保存在pdf文件中的对象时,没有模糊？

另外,保存它的最佳模式是什么？

png,eps？还是其他一些？我不能这样做,pdf因为有一个隐藏的数字发生混乱Latexmk编译.

我有两个数据向量,我已将它们放入matplotlib.scatter().现在,我想过度拟合这些数据的线性拟合.我该怎么做？我尝试过使用scikitlearn和np.scatter.

在MySQL命令行界面中,当您执行查询时,它将告诉您在打印结果后执行查询所需的时间.

在Postgres命令行界面(psql)中它没有告诉你.我知道如何配置日志记录,以便我可以从日志中获取信息,但将它打印到标准输出会更方便MySQL.

可以这样做吗？

我想multiprocessing在大型数据集上使用以找到两个gps点之间的距离.我构建了一个测试集,但是我无法multiprocessing在这个集合上工作.

import pandas as pd
from geopy.distance import vincenty
from itertools import combinations
import multiprocessing as mp

df = pd.DataFrame({'ser_no': [1, 2, 3, 4, 5, 6, 7, 8, 9, 0],
                'co_nm': ['aa', 'aa', 'aa', 'bb', 'bb', 'bb', 'bb', 'cc', 'cc', 'cc'],
                'lat': [1, 2, 3, 4, 5, 6, 7, 8, 9, 10],
                'lon': [21, 22, 23, 24, 25, 26, 27, 28, 29, 30]})



def calc_dist(x):
    return pd.DataFrame(
               [ [grp,
                  df.loc[c[0]].ser_no,
                  df.loc[c[1]].ser_no,
                  vincenty(df.loc[c[0], x], 
                           df.loc[c[1], x])
                 ] …

我如何在Python中用数字解决ODE？

考虑

\ddot{u}(\phi) = -u + \sqrt{u}

具备以下条件

u(0) = 1.49907

和

\dot{u}(0) = 0

有约束

0 <= \phi <= 7\pi.

最后,我想生成一个参数图,其中x和y坐标是作为u的函数生成的.

问题是,我需要运行odeint两次,因为这是二阶微分方程.我试过让它在第一次运行后再次运行,但它回来时出现雅可比错误.必须有一种方法可以同时运行它两次.

这是错误:

odepack.error:该函数及其雅可比行列式必须是可调用函数

以下代码生成的代码.有问题的行是sol = odeint.

import numpy as np
from scipy.integrate import odeint
import matplotlib.pyplot as plt
from numpy import linspace


def f(u, t):
    return -u + np.sqrt(u)


times = linspace(0.0001, 7 * np.pi, 1000)
y0 = 1.49907
yprime0 = 0
yvals = odeint(f, yprime0, times)

sol = odeint(yvals, y0, times)

x = 1 / sol * …

当我尝试使用该命令在ubuntu 13.04上安装Scrapy时sudo easy_install -U Scrapy,出现以下错误:

Searching for Scrapy
Reading pypi.python.org/simple/Scrapy/
Best match: Scrapy 0.16.5
Processing Scrapy-0.16.5-py2.7.egg
Scrapy 0.16.5 is already the active version in easy-install.pth
Installing scrapy script to /usr/local/bin

Using /usr/local/lib/python2.7/dist-packages/Scrapy-0.16.5-py2.7.egg
Processing dependencies for Scrapy
Searching for Twisted>=8.0
Reading pypi.python.org/simple/Twisted/
Best match: Twisted 13.1.0
Downloading pypi.python.org/packages/source/T/Twisted/Twisted-13.1.0.tar.bz2#md5=5609c91ed465f5a7da48d30a0e7b6960
Processing Twisted-13.1.0.tar.bz2
Writing /tmp/easy_install-xYafJF/Twisted-13.1.0/setup.cfg
Running Twisted-13.1.0/setup.py -q bdist_egg --dist-dir /tmp/easy_install-xYafJF/Twisted-13.1.0/egg-dist-tmp-ybetxt
twisted/runner/portmap.c:10:20: fatal error: Python.h: No such file or directory
compilation terminated.
error: Setup script exited with error: command 'x86_64-linux-gnu-gcc' failed with …

我试图用柯蒂斯复制轨道力学中的一个情节,但我不能完全理解它.但是,我已经np.arctan2从切换到np.arctan.

也许我的实施arctan2不正确？

import pylab
import numpy as np

e = np.arange(0.0, 1.0, 0.15).reshape(-1, 1)

nu = np.linspace(0.001, 2 * np.pi - 0.001, 50000)
M2evals = (2 * np.arctan2(1, 1 / (((1 - e) / (1 + e)) ** 0.5 * np.tan(nu / 2) -
           e * (1 - e ** 2) ** 0.5 * np.sin(nu) / (1 + e * np.cos(nu)))))

fig2 = pylab.figure()
ax2 = fig2.add_subplot(111)

for Me2, _e in zip(M2evals, e.ravel()):
    ax2.plot(nu.ravel(), …

我是Python的新手,我从未学过任何其他编程语言.我似乎得到了这个错误,我已经阅读了其他帖子,但他们说在[dollar = 0]之前放置全局,这会产生语法错误,因为它不允许[= 0].我正在使用[美元]作为计数器,所以我可以跟踪我添加到它的内容并在需要时将其显示回来.有人能帮助我吗？谢谢.

<>代码<>

    dollars = 0

    def sol():
        print('Search or Leave?')
        sol = input()
        if sol == 'Search':
            search()
        if sol == 'Leave':
            leave()

    def search():
        print('You gain 5 bucks')
        dollars = dollars + 5
        shop()

    def leave():
        shop()

    def shop():
        shop = input()
        if shop == 'Shortsword':
            if money < 4:
                print('I\'m sorry, but you don\'t have enough dollars to buy that item.')
                shop1()
            if money > 4:
                print('Item purchased!')
                print('You now have ' + dollars + …

import numpy as np
from scipy.optimize import fsolve

musun = 132712000000
T = 365.25 * 86400 * 2 / 3
e = 581.2392124070273


def f(x):
    return ((T * musun ** 2 / (2 * np.pi)) ** (1 / 3) * np.sqrt(1 - x ** 2)
        - np.sqrt(.5 * musun ** 2 / e * (1 - x ** 2)))


x = fsolve(f, 0.01)
f(x)

print x

这段代码有什么问题？它似乎行不通.

我查看了scipy.integrate.ode但我无法找到如何实际使用这些集成方法,dorpi5以及dop853.

我想尝试将ode集成python与mathematica我的python代码集成到这两种方法中,以了解它如何影响结果但不知道如何.