使用近似技术在一个或多个维度上集成函数的算法,而不是使用符号代数和微积分的精确闭合形式解决方案。包括自适应求积法,蒙特卡罗方法,有限元分析,马尔可夫链等概念。
我使用不同的数值方法来理解不同类型的积分器在不同时间步长产生的结果。我通过
我尝试编写一个代码,但它没有给我所期望的结果,因为我期望一个更接近准确值的值。 将 numpy 导入为 np 从 scipy.special 导入 root_hermitenorm def Gauss_hermite(f...
我尝试编写一个代码,但它没有给我所期望的结果,因为我期望一个更接近准确值的值。 将 numpy 导入为 np 从 scipy.special 导入 root_hermitenorm def Gauss_hermite(f...
我正在致力于将模型从 MATLAB 转换为 Python。该模型的关键在于 MATLAB 的 ode15。在 MATLAB 执行中,ode15s 有标准选项: 选项 = odeset() [t P] = ode15s(@
将 numpy 导入为 np 从 scipy import 集成 将 matplotlib.pyplot 导入为 plt 导入日期时间 开始 = datetime.datetime.now() plt.rcParams['axes.grid'] = True XX=[无、11.3、14.8、7.6、10.5、12.7...
假设: z = [0.4 0.5 0.75]' 函数 y = myfunct(x) y = 四边形(@sin, 0, x) 我想计算 sin(x) 从 0 到 0.4、0.5 和 0.75 的定积分,使用: 我的函数(z) 然而,Matlab
正方形内两条线之间的面积 [-1,+1] x [-1,+1]
我正在Matlab中进行一个项目,需要找到正方形[-1,+1]x[-1,+1]内相交于一点(xIntersection,yIntersection)的两条线之间的面积。所以我们的想法是减去...
我正在生成此图(简化代码版本): k = 1000 r = [100 220 470 1*k 2200 4700 10*k 22*k 47*k 100*k 220*k 470*k 1000*k ] unModNB = [0.72 0.746 0.801 0.92 1.16 1.69 2.78 4.6 6.45 9.1 11...
我正在尝试定义一个看起来有点像这样的函数: 函数 u = 预期效用(x) u = int(效用(x,I)*密度(I),I,-10,10) 函数 t = 效用(x,I) m=x+I 如果 m>=0 u = m^(1/2) 否则...
Matlab 集成拟合的“linearinterp”返回错误“第一个输入参数必须是函数句柄”
尝试对一些实验收集的数据进行积分。 使用信封和abs函数后,我使用拟合函数来获得我想要积分的方程(不幸的是“poly”......
有人可以帮助我解决这个问题。我必须使用 Matlab 的函数 trapz 来计算积分,增加网格间隔 N。我需要这样做,直到达到容差。 N = 1; %初始
我有一个用户定义的函数 函数[Rout2] = SSP2(a1,a2,b1,b2,pB) 如果a1/b1>a2/b2 如果pB<=a1-b1*a2/b2 Rout2 = (a1-pB)/b1; else if pB<=a1+a2 Rout2 = (a1+a...
我有向量 top、bottom、mid,我试图将它们与 i 积分,将原点保持在索引 11。 i=0:0.01e-06:1.3e-06; 底部=cumtrapz(i(11:结束),b(11:结束)); 顶部=cumtrapz(i(11:结束),t(11:结束));...
我正在使用integral2 matlab命令对二重积分进行数值计算。我的代码如下: l1=0; m1=0; funn = @(theta,phi)(Y_hsph_new(l1,m1,theta,phi)); q = 积分2(funn,0,pi/2,0,2*...
我无法对 x= 0 到 pi/2 和 y = 0 到 pi/2 的 heaviside(y-f) 进行积分(仅 syms 积分“int”)。运行代码后,输出如下所示,但我无法获取
通过在 matlab 中将变量限制与未知参数集成进行数据拟合
我想用方程拟合实验数据: y = C(1) * 74.826 * (x/T(1))^3 * 积分 (t^4*exp(t)/(exp(t)-1)^2), 0, T(1)/ x) + C(2) * 24.942 * (T(2)/x)^2 * exp(T(2)/x)/(exp(T(2)/x)-1)^2 + C...
我在matlab中有以下代码: 伽玛=1.01; [t,phi] = ode45(@(t,x)(伽玛-F(x,pi/6, 0.5)), [0,100], 0); 绘图(t,phi) 推迟; 标题('\伽玛 = 1.01') 函数输出 = F(p, a, b) phi = mod(p,2*...
我有一个概率分布函数,如下: f = @(x) h0 * (1/sqrt(2*pi*sig2)) .* exp(-(x-mh0).^2 / (2*sig2))+ (1-h0) * (1/ sqrt(2*pi*sig2)) .* exp(-(x-mh1).^2 / (2*sig2)); 我...
我正在 Matlab 中对螺线管建模。我正在评估柱塞运动每一步的螺线管电感 L 的值。 使用力-行程特征数据,我...
这个 DOFOR 伪代码在 MATLAB 中的等价物是什么?
我按照教科书中的以下伪代码进行操作,并且能够在 MATLAB 中成功实现左侧的函数。现在我正在研究右边的函数,但我没有