我需要编写一个python代码来计算1/(1+x^2)的积分(-5,5)的确切值。我知道答案是2弧坦(5),大致相当于2.746801.
下面是我所写的代码,但是我得到了一个稍微不同的答案,我想知道我是否可以做些什么来使这段代码更准确?谢谢你的帮助!
## The function to be integrated
def func(x):
return 1/(1 + x**2)
## Defining variables
a = -5.0
b = 5.0
dx = 1.0
Area = 0
## Number of trapezoids
n = int((b-a)/dx)
我试图使用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,
我想在下面的代码中找到输出功率Po的积分: Vo = 54.6
# defining a function for duty cycle, output current and output power
def duty_cycle(output_voltage, array_voltage):
duty_cycle = np.divide(output_voltage, array_voltage)
return duty_cycle
def output_current(array_current, duty_cycle):
output_current =
我有一个问题,因为这个工作适用于许多函数,但当我试图绘制正弦的积分时(我正在使用matlab 2010),我遇到了麻烦:
clear all
close all
clc
x = linspace(-10, 10, 100);
f = @(x) sin(x);
I = arrayfun(@(x) quad(f, 0, x), x);
plot(x, f(x),'r', x, I, 'b')
我希望有一个-cos(x),但是我得到了偏移量为1的东西,为什么会发生这种情况?该如何解决这个问题呢?
当我想通过RCaller运行一些代码时,Java会抛出以下异常:
Exception in thread "JavaFX Application Thread" com.github.rcaller.exception.ExecutionException: Can not send the source code to R file due to: java.io.IOException: The pipe is being closed Maximum number of retries exceeded.
这是我的密码:
protected void initialize
我已经创建了一个python脚本,它可以找到christoffel符号并求解测地方程。今天我想,在弯曲的2D空间中添加计算弧长的功能会很有用。我已经创建了使用渐近的代码,一切似乎都是正确的,但由于某种原因,它只是没有计算积分!它只是简单地显示了积分。谁能帮帮我。我已经尝试了所有的方法,甚至使用dsolve,但都不起作用。这是我的代码的当前状态: import sympy as sp
from sympy.interactive import printing
printing.init_printing(use_latex=True)
g = [[[11], [12]], [[21],
我将查找$I_m=\int_0^{2\pi} \prod_{k=1}^m cos(kx){}dx$,其中$m=1,2,3\ldots$
简单的SAGE代码:
x=var('x')
f = lambda m,x : prod([cos(k*x) for k in range(1,m+1)])
for m in range(1,15+1):
print m, numerical_integral(f(m,x), 0, 2*pi)[0],integrate(f(m,x),x,0,2*pi).n()
输出:
1 -1.47676658757e-16 0.000000000000
我试图运行一个SymPy集成,在这里我创建了自己的变量,积分是
import numpy as np
import seaborn
import math
import pandas as pd
import matplotlib as plt
import scipy.integrate as integrate
import scipy.special as special
from sympy import init_session
init_session()
%matplotlib inline
dP, dT = symbols('dP dT', real = T
我试图构建一些简单的运算,比如一个导数和积分函数来操作lambda函数,因为渐近和was都在努力地集成一些我传递给它们的东西。
导数函数没有给出任何问题,而是在绘制时返回输入函数的导数,但积分函数不返回,也没有绘制输入的正确积分。
import matplotlib.pyplot as plt
import numpy as np
from phys_func import func
sr = [-10,10]
x = np.linspace(sr[0],sr[1], 100)
F = lambda x: x**2
f = func.I(F,x)
plt.plot(x,F(x), lab
用Monte积分估计下列积分:
我试图对下面的问题进行蒙特卡罗积分,其中p(x)是高斯分布,均值为1,方差为2 (见图)。
我被告知,一旦我们从正态分布中提取样本,pdf就会在积分中消失。请解释这个概念,以及如何用Python解决这个问题。下面是我的尝试。
def func(x):
return (math.exp(x))*x
mu = 1
sigma = sqrt(2)
N = 1000
areas = []
for i in range(N):
xrand = np.zeros(N)
for i in range (len(xrand)):
x
我试图计算一些积分,例如:
a = 1/sqrt(2);
b = -5;
c = 62;
d = 1;
f = exp(-x^2-y^2)*(erfc((sym(a) + 1/(x^2+y^2)*(sym(b)*x+sym(d)*y))*sqrt((x^2+y^2)*sym(10.^(c/10))))...
+ erfc((sym(a) - 1/(x^2+y^2)*(sym(b)*x+sym(d)*y))*sqrt((x^2+y^2)*sym(10.^(c/10)))));
h = int(int(f,x,-Inf,Inf),y,-Inf,Inf);
我正在计算可见光谱中太阳辐射的百分比。
我的代码是:
h=6.626e-34; %planck constant
c=3e8; %speed of light
k=1.38066e-23; %boltzman constant
T1=5777; %temperature in surface of sun
syms f
J1 =(2*pi*h/(c^2))*((f^3)/(exp((h*f)/(k*T1))-1)); %planck radiation law
hold on;
ezplot(J1,[0,1.5e15])
a=int(J1,f,0,inf)
我想解决这样的问题:
dy/dt = 0.01*y*(1-y), find t when y = 0.8 (0<t<3000)
我尝试过Python中的ode函数,但它只能在给定t时才能计算y。
那么,在Python中有什么简单的方法来解决这个问题呢?
PS:这个函数只是一个简单的例子。我真正的问题太复杂了,无法用分析的方法解决。所以我想知道如何用数字来解决这个问题。我认为这个问题更像是一个优化问题:
Objective function y(t) = 0.8, Subject to dy/dt = 0.01*y*(1-y), and 0<t<3000
我真正的问
我想要积分一个三维函数。让我们将该函数定义为:
def pi(x, y, z): return x + y ** 2 + z ** 3
作为一个简单的例子。让我们选择作为[0,1]x[0,2]x[0,3]的域。是我尝试的第一件事。我创建了一个pi(x,y,z)评估的3D张量,然后进行3个1D积分:
from scipy.integrate import trapz
import numpy as np
x = np.linspace(0, 1)
y = np.linspace(0, 2)
z = np.linspace(0, 3)
print(trapz(trapz(trapz(pi(x[:,