我试着做一个谐振子的简单例子,它将用Runge四阶法求解。待解的二阶常微分方程(ODE)及其初始条件如下:
Y‘+y=0
y(0) =0,y'(0) = 1/pi
范围在0到1之间,有100个步骤。我使用u作为辅助变量,将我的二阶ODE用两种一阶代码分隔开:
y‘=u
u‘= -y
解析解为正弦y(x) = (1/pi)^2 sin(pi*x)。
我的Python代码如下:
from math import pi
from numpy import arange
from matplotlib.pyplot import plot, show
# y' = u
# u'
在C/C++中,true标准化为1,false为0。虽然不是很好的做法
if variable:
#do something
在python中,决策似乎是可以的,如果变量是数字,则根据变量是零还是零来做出决定。以下内容是合理的。
a = 123123
if a:
print "if condition is true" # this prints
a = "a string"
if a:
print "if condition is true" # this prints
a = None
if a:
我试图使用python/scipy计算两个极限之间的积分。
我正在使用在线计算器双重检查我的结果(,),当我设定了某些限制时,我的结果不同意。
所用的代码是:
import scipy as sp
import numpy as np
def integrand(x):
return np.exp(-0.5*x**2)
def int_test(a,b):
# a and b are the lower and upper bounds of the integration
return sp.integrate.quad(integrand,a,b)
当设定(a,
我想用Python语言读取Fortran90二进制文件,因此有一个用Scipy开发的模块叫做"fortranFile",就像这里给出的
如何在python中包含此模块,以便可以从代码中调用它,如下所示
import fortranFile
f = fortranFile('file.bin')
x = f.readReals('f')
print x
现在,我得到了错误
Traceback (most recent call last):
File "readfortbin.py", line 5, in <m
我正在读R.G.Droomey的“如何用计算机解决它”一书,在第三章中,我发现了这个问题--“设计和实现iteratively compute the reciprocal of a number的算法”。我完全搞不懂如何做到这一点,因为他之前教的是如何to compute the square roots,然后突然提出了这个问题。有什么关系?
这方面的算法是什么?另外,当我们可以直接找到一个数字的倒数时,我们为什么需要这个呢?
我不知道如何让Python在执行重复正弦运算时像我的计算器一样工作。
我的Python代码:
number = float(.565)
pi = float(math.pi)
for j in range(25):
number = round(math.sin(number*pi), 16)
这将返回0.446445551265。
当我用我的计算器做同样的运算时,我得到0.4464454868。
无论我对round()使用多高的精度,我都不能得到与我的计算器相同的结果。
我也尝试过使用Decimal库,但没有得到更好的结果。有什么想法吗?
我想知道如何从零开始创建一个情感模型。我有我的数据,文本列表,没有关于情感的标签。
Author Quotes
Dan Brown “Everything is possible. The impossible just takes longer.”
Dan Brown “Great minds are always feared by lesser minds.”
Dan Brown “Men go to far greater lengths to avoid what they fear than to obtain what they desire.”
Dan
我已经设置了一个基本的C#应用程序,以便在我使用Google.Apis.Pagespeedonline.v2 nuget包指定的网站上运行PageSpeed测试。
设置非常简单,而且我有一个变量,可以指定url,该url随后进入服务
// Create the service.
var service = new PagespeedonlineService(new BaseClientService.Initializer
{
ApplicationName = "PageSpeed Sample",
我在windows 8上使用python3.864位Pycharm,在试图从列表中生成数据时出现了内存错误。
我要做的是使用.csv包将一个巨大的csv (25 do )读入列表中,使用pd.Dataframe生成一个数据文件,然后导出一个带有pd.to_stata函数的.dta文件。我的内存是64 My,比数据大得多。
下面是错误消息:
MemoryError: Unable to allocate 25.8 GiB for an array with shape (77058858, 45) and data type object
我发现了三个类似的问题,但没有一个对我有用。
:这个解决
我很难理解如何在Apache Commons Math中使用。
我有一个函数,我想要得到它的一阶导数。然后我想要得到多个不同值的导数的值。
Logit logit = new Logit(0.1, 10.0);
DerivativeStructure ds = // How to instanctiate?
DerivativeStructure dsRes = logit.value(ds);
// How to use dsRes to get the value of the derivative function applied on
// several values?
此外,如
我是Python的新手,在这方面我真的被难住了。我正在读一本书,代码运行得很好;我就是不明白!
T[i+1] = m*v[i+1]ˆ**/L
这段代码中的双星号部分是什么?它甚至后面跟着一个正斜杠。变量L的初始化值为1.0,但看起来像是有人摔倒在键盘上,但代码运行正常。这是一个数学表达式还是别的什么?如果能帮助我理解这一点,我将不胜感激。谢谢!
完整代码:
from pylab import *
g = 9.8 # m/sˆ2
dt = 0.01 # s
time = 10.0 # s
v0 = 2.0 # s
D = 0.05 #
L = 1.0 # m
m = 0.5 # kg
# Num
在找到最佳参数后,我试图绘制数据的Probability Density Function (PDF) plot,但图显示的是一条平线而不是一条曲线。
这是缩放的问题吗?
这是Continuous or Discrete data的问题吗?数据文件可用
这里的目的是得到最佳的分配配件,然后绘制PDF函数。
我的数据值非常小,比如:0.21, 1.117 .etc。数据统计和PDF地块如下:
我的剧本如下:
from time import time
from datetime import datetime
start_time = datetime.now()