如何在python中给出一个任意的初始条件
在Python中给出一个任意的初始条件，可以通过以下步骤实现：
- 定义变量：首先，你需要定义一个变量来存储初始条件的值。例如，你可以使用以下代码创建一个名为"initial_condition"的变量，并将其值设置为任意的初始条件。
initial_condition = "任意的初始条件"- 根据需要选择数据类型：根据你的需求，选择适当的数据类型来存储初始条件的值。Python提供了多种数据类型，包括字符串、整数、浮点数、布尔值等。根据具体情况选择合适的数据类型。
- 使用初始条件：一旦你定义了初始条件的变量，你可以在代码中使用它。根据你的需求，可以将初始条件用于计算、条件判断、函数调用等各种操作。
下面是一个示例，展示了如何在Python中给出一个任意的初始条件：
# 定义初始条件变量
initial_condition = "任意的初始条件"
# 使用初始条件
print("初始条件是：" + initial_condition)
# 进行计算
result = len(initial_condition)
print("初始条件的长度是：" + str(result))这个示例中，我们首先定义了一个名为"initial_condition"的变量，并将其值设置为"任意的初始条件"。然后，我们使用该初始条件进行了一些操作，例如打印初始条件的值和计算初始条件的长度。你可以根据具体需求修改和扩展这个示例。
相关·内容
我正在用pdepe解决一个二阶偏微分方程。我的问题的初始条件是一个二阶常微分方程的解，我用ode45解决了它。如何在pdepe中传递ode45解决方案作为我的初始条件？
浏览 1提问于2014-12-21得票数 0
1回答
我希望能够在数字上区分和集成Python中的数组。我意识到，这是有作用的，在裸体和参与。我注意到一个抵消，然而，当集成。
例如，我从一个初始函数y=cos(x)开始。
然后我用numpy.gradient取导数。它按预期工作(图为-sin(x))：
当我将导数与scipy.cumtrapz集成时，我希望得到初始函数。然而，也有一些抵消。我知道-sin(x)的积分是cos(X)+常数，所以这个常数不是用累积积分计算的吗？
我关心的是，如果你有一些任意的信号，并且不知道初始条件/边界条件，那么+常量项会不会与累积项一起消失呢？有解决这个问题的方法吗？
我使用的代码如下：
import num
浏览 0提问于2020-12-07得票数 1
回答已采纳
当我给出alpha>1的值时，例如，alpha=1.05，它给出错误：
??? Error using ==> fde12 at 114
A not sufficient number of assigned initial conditions.
Order ALPHA = 1.050000 requires 2 initial conditions. See FDE12
当我给出y0作为2个初始条件的向量时，它给出了错误的输出。我认为这是因为我给出了同一组参数的ODE45返回的两个值中的初始条件。我知道第二个初始条件只能迭代地来自fde12，但不确定怎么做。
浏览 3提问于2013-04-28得票数 0
这个问题是关于statsmodels.tsa.stattools python库中的增广Dickey-Fuller测试实现的- adfuller()。
原则上，AIC和BIC应该为一组可用的模型计算信息准则，并选择最优(信息损失最小的模型)。
但是，它们是如何在增强的Dickey-Fuller的背景下运作的呢？
我不明白的是:我已经设置了maxlag=30，BIC选择了带有一些信息标准的lags=5。我已经设置了maxlag=40 - BIC仍然选择lags=5，但是信息标准已经改变了！为什么世界上相同数量的滞后的信息标准会随着最大滞后的变化而不同？
有时，这会导致模型的选择发生变化，当最大时
浏览 0提问于2015-11-01得票数 3
回答已采纳
我试着写一个算法来解决非线性常微分方程。
dr/dt = r(t)+r²(t)
它有(一种可能的)解决方案
r(t) = r²(t)/2+r³(t)/3
为此，我在python中实现了euler方法和RK4方法。为了检查错误，我在rosettacode上使用了这个例子：
dT/dt = -k(T(t)-T0)
使用该解决方案
T0 + (Ts − T0)*exp(−kt)
因此，我的代码现在看起来像这样
import numpy as np
from matplotlib import pyplot as plt
def test_func_1(t, x):
return x
浏览 7提问于2017-02-16得票数 2
回答已采纳
1回答
我有两个测试用例，但其中只有一个给出了正确的输出。第二个是给出数组的第一个元素。请告诉我哪里出了问题。谢谢。
#include <iostream>
using namespace std;
void thirdLargest(int arr[], int arr_size)
{
int m1 = arr[0], m2 = arr[0], m3 = arr[0], i;
for(i = 0; i < arr_size; i++)
{
if(m1 < arr[i])
{
m3 = m2;
m2 = m1;
浏览 0提问于2019-08-21得票数 0
回答已采纳
下面的模拟是一个弹跳球模型。我没有将状态端口连接到确定的初始条件，而是连接积分器输出。给出了第一步后的结果，初始条件块的输出为零。为什么？
浏览 1提问于2017-04-10得票数 0
回答已采纳
1回答
用FiPy求解柱面几何中扩散方程的稳态解与用其它软件(如Mathematica )得到的解有很大不同。 方程式是： $0 = \frac{1}{r}\frac{d}{dr}\left(\frac{r}{T^{1/2}}\frac{dT}{dr}\right) + cte*T^{3/2}$ 这意味着，通过使用CylindricalGrid1D网格，我们可以将方程写成： mesh = CylindricalGrid1D(nr=100, dr=0.01, origin=0.0)
T = CellVariable(name='temperature', mesh=mesh, hasO
浏览 25提问于2019-02-13得票数 0
回答已采纳
1回答
我相信我在试图导入和测试在Simulink中创建的系统识别工具箱(SIT)中的传递函数时，做了一些根本错误的事情。
举一个简单的例子来说明我在做什么。我有一个输入，它是一个从12秒到25秒的偏移正弦波，振幅为1，频率为1.5rad/s，给出了测量的输出。
我使用SIT创建了一个简单的2极1零传递函数，该函数给出了以下协议：
然后，我尝试将此传递函数导入Simulink，以便在以下配置中进行研究，该结构具有频率1.5rad/s的正弦输入和启动t=12。
当我运行这个模拟13秒时，块的输入与预期相同，但后传递函数信号显示的结果与预期的基本一致，并且是一个数量级的输出。
pre：
浏览 2提问于2014-02-25得票数 1
回答已采纳
我试着在SimMechanics中旋转机器的旋转关节。执行器来自SimElectronics，所以我使用旋转-旋转接口块将电机和关节连接在一起。整个块模型如下(也是可用的)
但是，当我试图运行模拟时，我得到了以下错误：
没有为一个或多个Simulink转换器块提供足够的输入导数，与突出显示的解决程序配置块相关联，用于所选的求解器。隐式求解器(ode23t、ode15s和ode14x)通常比显式求解器需要更少的输入导数，而本地求解器不需要任何输入导数。您可以通过以下链接的Simulink转换器块的输入处理选项卡上选择不同的选项(包括打开输入过滤或增加用户提供的输入导数)来提供更多的输
浏览 2提问于2013-04-28得票数 2
回答已采纳
我是这个网站的新手，我不是英国人，所以如果这篇文章不像你预期的那样好，请原谅。让我们进入正题:为了给你一个我所说的“运行按钮”的例子，以python one为例。在visual studio代码中安装python扩展后，您可以按下屏幕右上方的按钮，然后vsc会自动在对活动.py文件运行python的外壳中执行命令。给出这个例子，我想知道如何为vsc制作一个按钮，它可以自动创建一个新的终端，执行一个我想在按钮的制作过程中插入的命令，当然终端不应该在给出期望的输出后消失。只是给你更多的信息，我在ubuntu上，如果这有任何帮助的话。请随时告诉我，我是否应该编辑此问题。
浏览 96提问于2021-04-06得票数 2
回答已采纳
1回答
我有兴趣解决下面的系统
x′=2x+3y+z
y′=x−3y+2z
z′=x+2y−4z
在不同的时刻，甚至在无限大的时间里，有条件。
x(0)=1
y(0)=1
z(+inf)=0
我试图谷歌，但它似乎是没有方法的PYTHON，MATLAB适合这个系统。有什么想法吗？谢谢!
浏览 1提问于2021-03-08得票数 0
我正在将一些废水代谢模型转移到Python，特别是关于聚磷酸盐积累生物体(PAO)的一些模型。我已经成功地将模型移到Python中，并使用Gekko进行求解，这是很棒的！(在待实现的方程中仍有一些调整)。
我可以得到O2的变化到wokr和简单的氧气供应，但是，我不太知道如何让它通过多个循环重复来观察储存聚合物在许多循环中的变化。
这个氧气步骤(O2)的代码在初始循环中工作：
t = np.linspace(0, 19, 19) m.time = t
O2_step = np.zeros(19)
O2_step[17:] = 0.05
O2 = m.Param(value=O2_step)
我
浏览 8提问于2022-08-29得票数 2
我正在尝试将Simulink模型(为求解二阶ODE而构建并如下图所示)“转换”为python。
inW和inT是维数1601x2的矢量，增益MFW和MFT具有维数2x1，增益MC和MK具有维数2x2。所有初始条件都设置为0。
out_acc、out_rate和out_disp的维数均为1601x2。
在python中设置这个问题有帮助吗？
谢谢!
西蒙妮
浏览 2提问于2015-05-19得票数 0
回答已采纳
在上面的模拟中，我使用的是标准的SHM微分方程，a= -w^2*x，我使用的是带odeint的Python。尽管编辑了几次，我还是把输出作为一条直线，而不是正弦曲线。守则是：
from scipy.integrate import odeint
from pylab import *
k = 80 #Spring Constant
m = 8 #mass of block
omega = sqrt(k/m) #angular velocity
def deriv(x,t):
return array([x[1],(-1)*(k/m)*x[0]])
t = linspace(0,
浏览 0提问于2014-12-12得票数 1
回答已采纳
1回答
本程序用于求解方程使用二分法给出错误“函数显示返回值”。
在这种方法中，我们给出了一个函数f(x)，并且我们近似该函数的2个根a和b，使得f(a).f(b)<0。
然后我们找到另一个点
c=(a+b)/2
if f(c)==0
then root=c;
else
if f(a).f(c)<0
b=c;
if f(b).f(c)<0
a=c;
对于给定的迭代次数，我们重复这些步骤
#include<stdio.h>
#include<math.h>
#define e 0.000001/*e is the prescribed accura
浏览 2提问于2015-08-27得票数 1
我有一个用python编写的应用程序，用于计算函数的最小返回值。Am采用scipy.optimize.mminimize和SLSQP作为优化方法。它在循环中运行，为了节省时间，从查找局部极小值开始，我需要它使用我提供的x0。问题似乎是，它不在乎我给它什么x0。它只是在随机值的情况下开始优化。我做错什么了？
我编写了一个小型测试应用程序来测试最小化程序上的x0：
import numpy as np
import matplotlib.pyplot as plt
from scipy.optimize import minimize
global log
log = []
counter = 0
浏览 1提问于2020-02-21得票数 0
回答已采纳
1回答
我试图利用方法来解决这个问题，这个问题的表述方式类似于at。
最小化群内平方和(WCSS)：
但在我的公式中，必须最小化模块集群内的和。
给出了整数集X和簇k的个数。需要选择k个簇整数的值，这样就可以最小化差异模块的聚类和。
我是以交互的方式进行的，随机获取初始的mu值，然后将其调整为分配给集群的元素的平均值。
然而，这种方法只给出简单测试用例的正确答案。
浏览 0提问于2013-05-03得票数 1
1回答
我正在尝试求解一元一阶线性微分方程，目前使用的是scipy.integrate中的odient模块。但是，它在$y_0$中采用的初始条件是在域$x_0$的初始边界上求值的，而我得到的是某个随机点$x$处的$y$的值。 关于类似问题的建议是使用solve_bvp，这也不能完全解决我的问题。 我该怎么做呢？
浏览 21提问于2021-07-16得票数 0
1回答
我对python完全陌生，并尝试集成以下ode：
$\dot{x} = -2x-y^2$
$\dot{y} = -y-x^2
这会导致数组中的所有内容都为0，但我做错了什么呢？它主要是复制代码，与另一个，而不是耦合的ode，它工作得很好。
import numpy as np
import matplotlib.pyplot as plt
from scipy.integrate import ode
def fun(t, z):
"""
Right hand side of the differential equations
dx/d
浏览 1提问于2016-10-27得票数 0
回答已采纳
我是一个经验丰富的Matlab程序员，但我不知道如何解决这个(似乎)简单的问题。我有一个由5个方程和5个未知数组成的系统：
当不涉及代数方程时，我知道如何使用ode45求解ODE。在上述系统中，V (速度)和C (加速度)都是已知的常数，C是航天器的横向加速度。
解决这一问题的办法如下：
在t=0，我们知道Theta(0)，x(0)和y(0)。请记住，V和C都是常数和已知的。
给定Theta(0)和C/V，我们得到了Theta(t1)积分第4方程。有了Theta的这个新值，我们应该能够计算出新的Vx(t1)和Vy(t1)，这将给出x(t1)和y(t2)的新值。
重复
浏览 0提问于2019-02-09得票数 0
回答已采纳
1回答
我想使用Python的odeINT来集成由循环生成的多组方程。这些方程都是耦合的，因此必须通过对odeINT的一次调用同时进行积分。问题是初始条件("y0")必须是列表列表或矩阵(而不是简单的列表)。odeINT给出这个错误：“初始条件y0必须是一维的”。我想知道如何解决这个问题。以下是代码的一个示例；非常感谢您的任何想法。
class network:
def __init__(self):
self.i_range = 3
## INITIAL CONDITIONS WILL BE A LIST OF LISTS.
浏览 0提问于2016-12-27得票数 0
2回答
我正在做一个项目，我试图检测特定表面(竞技场)上的绿色和红色圆圈。当我尝试使用该竞技场的数字版本(PNG图像)时，我可以成功地检测到这两个彩色圆圈。 这是表面的数字图像： ? 现在，我把这个竞技场打印在一个柔盘上(没有那两个彩色圆圈)，并手动在上面放置彩色圆形硬币。但在通过一个130万像素的网络摄像头捕捉到它的图像后，颜色检测不起作用，并给出了错误的结果。 以下是通过网络摄像头捕捉到的打印的竞技场： ? 为什么没有检测到颜色？是否需要对摄像头图像进行后处理？我尝试过通过cv2.filter2D锐化图像，但也不起作用。 下面是从我的Python代码中检测红圈的代码片段： ip_im
浏览 33提问于2019-12-26得票数 1
回答已采纳
1回答
我问了一个关于matlabFunction如何工作()的问题，这引发了一个与ode45函数相关的问题。使用我在matlabFunction上的文章中给出的例子，当我通过ode45传递这个函数时，带一些初始条件，ode45是否将导数-s.*(X)读取为逼近未知函数x；-y+r.*x.*z和y，以及-b.*z+x.*y和z？更确切地说，我们有
matlabFunction([s*(y-x);-x*z+r*x-y; x*y-b*z],
'vars',{t,[x;y;z],[s;r;b]},'file','Example2');
浏览 0提问于2014-03-12得票数 0
回答已采纳
我在想象蟒蛇的振子系统。只要在其中一个变量中没有构建(使用dirac-delta函数)，系统就能工作。然而，当我试图包括它时，系统似乎没有变化。在下面，你可以找到我和一些问题一起经历的过程。
我用python编写了一些代码，能够正确地可视化规则振子的行为。但是，当在其中一个变量中包含一个逐步增量时，系统不会对此作出反应。下面，您可以找到包含增量的版本。代码运行，但是给出输出，就好像增量不存在一样。在代码中，您可以找到我获得正确解决方案的问题。
import numpy as np
from scipy.integrate import odeint
import matplotlib.pypl
浏览 1提问于2019-07-11得票数 0
1回答
人们可以使用Sympy来求解常微分方程。我在符号计算主题上的问题是:人们可以使用SymPy来求解一阶变系数方程吗？注意，这是ODE的一个特例。例如，如果我有一个像下面这样的等式 ? 如果可能的话，我该如何设置使用SymPy来求解这样的方程呢？
浏览 11提问于2020-11-16得票数 0
1回答
我写了一个代码来获得斐波纳契数列，使用斐波那契矩阵的概念，灵感来自于这篇Instagram帖子
def f(n):
import numpy as np
m=np.array([[1,1],[1,0]],dtype=np.int64)
m1=m
for i in range(0,n):
m1=np.matmul(m1,m)
return m1[0,1]
但在n=93之后，它开始给出负数。如果我使用np.int32作为数据类型，那么在n=47之后，它开始给出否定的和错误的结果。我使用的是Python3.9，我希望我的结果是整型的(而不是浮点型的)，怎么做才能得到n=1000或更
浏览 31提问于2021-05-29得票数 0
1回答
我使用Pipenv来管理Django项目的虚拟环境。我希望在Pipenv中设置一个别名，以执行一个冗长(并且经常使用)的python命令。
我尝试通过在我的.env文件中添加以下行来设置这个别名：
alias launch='python manage.py collectstatic --noinput; heroku local'
但是，当我在Pipenv中运行新的launch命令时，无法识别它：
$ pipenv shell
$ launch
bash: launch: command not found
有没有办法在我的Pipenv .env文件中添加别名，以便每次执行
浏览 1提问于2018-11-19得票数 0
回答已采纳
1回答
所以，我正在学习Python，我试着想出一段代码，用于在闰年创建一个图表，以便了解我的选择。我想可能是出了什么问题，我就是搞不清楚是什么。我的代码如下：
year = int(input("Which year do you want to check?: "))
if year % 4 == 0:
if year % 100 != 0:
if year % 400 == 0:
print("Leap year")
else:
print("Not leap")
else:
print("Not
浏览 5提问于2021-01-06得票数 0
我模拟了一个带有滞后控制器和锯子的控制回路。我确认了锯子是PT4系统。所以这个PT4街区是我的工厂。现在控制器给出了锯子的相对位置值。缩回0,3毫米或延伸0，3毫米。工厂想要有绝对的地位。因此，我从锯子的实际位置向具有相对位置的add块进行了反馈。到目前一切尚好。控制回路起作用了。唯一的问题是开始，因为模拟值从零开始。因此，需要一段时间才能将锯条延伸到所需位置。随后，在控制开始之前，将锯条延伸到8.6 to。因此，我也想从8.6 my开始模拟。
我该怎么做呢？
我尝试使用初始条件块，但这并没有改变任何事情。
图1显示了控制循环的输出，其中问题在开头。
图2显示了反馈
图3显示了IC块设置为8.
浏览 1提问于2016-02-12得票数 0
1回答
给出了一个函数，如何证明它在性质上是σ型的。到目前为止，我的方法是验证乙状结肠函数的属性是否成立：
1)它是单调的
2)受一对水平渐近线的约束。
3)它有一个“钟”形的一阶导数。
4)对于小于0的值是凸的，对于大于0的值是凹的。
5) Sigmoid函数及其仿射成分具有多重最优性。
例如，这是我试图验证的函数之一：
📷
试图证明所有属性似乎很费时，特别是在考试场景中。我认为有一个更好的方法存在，我应该遵循。请建议我可以使用什么其他方法。
谢谢。
浏览 0提问于2020-03-25得票数 2
1回答
我用FiPy来解决一个受生物学启发的问题。
本质上，我想代表一个二维平面，在不同的点，我有源和汇。源以固定的速率发射基板(不同的源可以有不同的速率)，汇以固定的速率消耗基板(不同的汇可以有不同的速率)。我的代码：
import numpy.matlib
from fipy import CellVariable, Grid2D, Viewer, TransientTerm, DiffusionTerm, ImplicitSourceTerm, ExplicitDiffusionTerm
from fipy.tools import numerix
from time import *
nx
浏览 2提问于2016-04-03得票数 2
回答已采纳
1回答
模拟卫星接近火星的运动
、、、、
我对Python非常陌生，并且已经生成了以下代码，但它不能工作。
这段代码试图绘制一颗卫星以可变的y和x距离(即2维空间)接近火星的运动。目前，火星被认为是静止的。
import scipy as sp
import numpy as np
import pylab as pl
import scipy.integrate as spi
G=6.67*(10**-11)
mm=6.4*(10**23)
#^mm is the mass of mars and G is the gravitational constant
def f(b,t):
xx=b[0]
浏览 1提问于2014-02-05得票数 0
1回答
我使用matlabFunction定义了一个函数。以下是代码：
matlabFunction([a16;-((1+x16^2)/(2*x16))*a16],'vars',{x16,[a16]},'file','DE_19')
似乎起作用了。然而，当我试图使用ode45来求解由matlabFunction定义的微分方程时，我会得到一个错误。以下是代码：
[x,y] = ode45(@(x16,Y) DE_19(x16,Y),[1,11],[2,7,5]);
我得到的错误是
Error using odearguments (line 93)
浏览 2提问于2014-03-03得票数 1
回答已采纳
我面临的问题如下:我想表示一个解决方案，其中给出了(t0,η)。比如说，(t0,η) = (0,1)。但是，我希望看到在整个(-1,1)域中表示的图形。
由于tspan与我的t0有着内在的联系，我如何使求解器解决给定t0的问题，而在tspan中的间隔中的初始值不同
因此，在给定的示例中，我想观察(-0.88, 0.88)中的ODE，但是要解决条件(t0, y0) = (0,1)的问题。
Ex8(c)
using DifferentialEquations
using Plots
function example2!(dy, y, t, p)
c1 = p
dy[1] = c
浏览 4提问于2020-12-20得票数 2
回答已采纳
1回答
我写了一些代码来解决二阶微分方程，但它给出了一个完全错误的结果。我认为问题在于以正确的方式表达Euler方法。我也尝试了另一个二阶常微分方程，但在近似y(x)时也失败了。你能指出错误在哪里吗？请看图表和代码：
解算ODE：
y"(x)=-1001y'(x)-1000y(x)
y(0)=1, y'(0)=0
解析解决方案：
y(x)=(1000exp(-x)-exp(-1000x))/999
重写为2个ODE的系统；
y'(x)=v(x)
v'(x)=-1001v(x)-1000y(x)=f(x,y,v)
y(0)=1,v(0)=0
Jupyter
浏览 1提问于2017-04-25得票数 0
1回答
我基本上需要在我的模拟中的每个时间步做'x =x+ a*b‘操作。B的值在整个模拟过程中都在变化。
请建议我应该使用什么块；如果可能的话，也要有一个图表。我想我需要一个区块，它可以作为一个来源，但也有一个选项，以更新。或者执行上述任务的任何其他方式就足够了。
我已经在Matlab脚本中做到了这一点，这是非常容易的，但我需要一个Simulink模型解决方案。
浏览 0提问于2020-03-31得票数 1
3回答
我有一份if/else声明。
if (something) {
//run if logic
} else {
//run else logic
}
在我的if logic phase期间，我可能会发现一些东西让我想要运行if/ the语句的其他部分。逻辑阶段太复杂，不适合于初始条件something。
是否有类似break;的内容可用于跳转到if/else语句的else部分。
我正在做这件事，但是我的团队不喜欢使用goto
if (something) {
//run if logic
if (somethingComplex) {
goto
浏览 3提问于2016-06-27得票数 2
回答已采纳
在命令行上，
python -c "z = [i for i in range(5)]; print (z)"给出
[0, 1, 2, 3, 4]
;的意思是应该有一个换行符。
但是换行符后面跟着标签呢？还是两个标签？
下面的代码是非功能性的，但它说明了我想要做的事情：
python -c "z = [i for i in range(5)]; for i in z: print (i)"
作为输出：
1
2
3
4
5
因此，基本上问题是:如何使用python -c执行以下代码
z = [i for i in range(5)]
for i in z:
浏览 2提问于2020-05-06得票数 1
回答已采纳
这就是我到目前为止得到的问题和结果，但我不知道如何才能获得下一个C1。Solve the IVP (1/x + 2y^2x)dx + (2yx^2 - cos(y))dy = 0, y(1) = pi.给出了一个隐式解决方案 from sympy import *
x = symbols('x')
y = Function('y')
deq = diff(y(x),x) + (1/x + 2*y(x)**2*x)/(2*y(x)*x**2 - cos(y(x)))
ysoln = dsolve(deq, y(x))
浏览 14提问于2020-09-14得票数 0
回答已采纳
在大多数解释语言中(例如在psuedo PHP中给出的例子)，我可以这样做
function f($name, &$errors = null) {
if(something_wrong) {
if(!is_null($errors) {
$errors[] = 'some error';
}
return false;
}
}
if(!f('test', $errs = array()) {
print_r($errs);
}
得到这样的结果
Array
浏览 6提问于2015-11-09得票数 0
回答已采纳
当我从shell脚本调用Python脚本时，它运行得很好。 python script.py 但在我从Gerrit中拉出相同的脚本，然后添加用于调用Python脚本的代码后，它给出了以下错误： ./script.sh: line 126: python: command not found
浏览 206提问于2021-04-21得票数 0
回答已采纳
1回答
数值求解微分方程
、、、、
我试着求解一个非常简单的方程f= t**2。我编写了一个for -循环，以便第一次使用f，然后使用每个循环的解作为下一个循环的输入函数。
我不确定我的数值求解方法是否正确，出于某种原因，我的循环只工作了两次(一次是通过if-然后是else-语句)，然后只给出零。
任何帮助都非常感谢。谢谢！
## IMPORT PACKAGES
import numpy as np
import math
import sympy as sym
import matplotlib.pyplot as plt
## Loop to solve numerically
for i in range(1,4,1)
浏览 2提问于2020-06-05得票数 1
回答已采纳
2回答
根据中给出的示例，我们可以从python脚本中调用MATLAB文件triarea.m，如下所示：
import matlab.engine
eng = matlab.engine.start_matlab()
eng.triarea(nargout=0)
如果我们在Python中有一个由用户file_name = 'triarea'输入MATLAB文件名的变量，
eng.file_name(nargout=0)
返回错误
MatlabExecutionError:用于输入'file_name‘类型的参数的未定义函数’file_name‘。
如何在Python中运
浏览 8提问于2019-10-02得票数 1
回答已采纳
1回答
我有一个在OSX10.8Xcode中生成c++11 clang服务器的qt
它对名称(单个挖掘线程)进行加密证明的工作挖掘。
当我单击.app时，过程需要4. 1/2小时
当我在.app文件夹中运行确切的exe时，从终端运行进程需要30分钟
问题，我如何调试这个？
谢谢
====================================
更糟的是：
在终端中运行的挖掘服务器。
如果我启动了连接到服务器的GUI程序，并且只发送(ipc)它就可以使用“it”命令:4小时
如果我启动了一个连接到服务器并只发送(ipc)的CL命令，那么“it”命令: 30分钟。
这两种情况下，服务器都是在一个紧
浏览 2提问于2014-04-16得票数 0
回答已采纳
我试图演示如何使用系统传递函数的定义和python-control模块来“求解”微分方程初值问题的“解”(模拟解)。事实上，我真的是个新手。
我有一个简单的微分作为例子：y'' - 4y' + 13y = 0，具有这些初始条件：y(0) = 1和y'(0) = 0。
我用手实现了这个传递函数：Y(s) = (s - 4)/(s^2 - 4*s + 13)。
所以，在Python中，我正在编写这段代码(请注意，y_ans是这个微分IVP作为的答案。
import numpy as np
import control as ctl
import matplotlib
浏览 14提问于2020-10-29得票数 1
回答已采纳
我正在开发一个模型预测控制器(MPC)与移动地平线估计(MHE)插件的动态仿真程序。
我的问题是，模拟程序在每个时间步骤中执行Python脚本。因此，每一个时间步骤，一个新的模型，在GEKKO被产生。是否有可能重新加载模型和数据文件？例如，给出数据到GEKKO的路径？
诚挚的问候,
莫里茨
浏览 15提问于2022-01-12得票数 2
我正在尝试在页面上显示XML提要中的项目
这是我的代码。
$feed = 'http://awebsite.co.uk/directory/rssfeed.php?thecondition=(title%20LIKE%20'%british%'%20OR%20description%20LIKE%20'%british%')';
$xml = simplexml_load_file("$feed");
if ( $xml->channel->item !='' ) {
foreach ($xml->
浏览 0提问于2013-12-13得票数 0
1回答
我想使用python来计算一些数据的校验和，问题是在sha1中，hashlib输入是以字符串形式给出的。
有没有可能在python中计算sha1，但以某种方式给出原始字节作为输入？
我之所以这样问，是因为如果我想要计算文件的哈希，在C中我会使用openssl库，只传递普通的字节，但在Python语言中，我需要传递字符串，所以如果我要计算某个特定文件的哈希，在两种语言中我会得到不同的结果。
浏览 4提问于2013-03-11得票数 0
回答已采纳
2回答
我有一个非常简单的问题:生成y_1，.，y_n
y_t = x_t +a* y_{t-1} +b* y_{t-2}
给出了a，b，x_1，.，x_n，以及y_0和y_{-1}。我可以用MATLAB以两种方式(在使用filter之前使用一个循环和一些手动调整)来完成这个任务。如何使用filter的初始条件内建备抵来完成此操作？
我的代码如下：
clc; clear; tic;
n = 20;
x = (1:1:n)'; %//' this is here so code looks good in SO.com
yzero = -3;
ymi
浏览 0提问于2016-03-10得票数 0
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
实时音视频
即时通信 IM活动推荐
腾讯云开发者
