如何解决这个积分,我有我想要的功能和整合这个函数t和无穷之间得到一个明确的解决方案,有人能帮助我吗

6视图(30天)
α= 0.2
β= 0.05
σ= 0.01
r0 = 0.05
ε= -0.46
NRepl = 60
NSteps = 60
T = 60
dt = T / NSteps
模型=“vasicek”
rng (75)
T = 60
% t = 0:1:60
r = pr_1 (r0,α,β,σ,T, NSteps, NRepl,模型)
=意味着(右)
% r1 =公关(r0,α,β,σ,T, NSteps, NRepl,模型)
% b =意味着(r1)
% r2 = pr_2 (r0,α,β,σ,T, NSteps, NRepl,模型)
% c =意味着(r2)
图;
情节(,“r”)
%抓住
%的阴谋(b,“k”)
%抓住
%的阴谋(c, b)
%推迟
% %
α1 = -0.095
φ= 0.001
sigma1 = 0.001
v0 = 0.00123
rng (1075)
v = mortalitysimulation(α1、φ,sigma1 v0, T, NSteps, NRepl,模型)
m =意味着(v)
v2 = mortalitysimulation1(α1,φ,sigma1 v0, T, NSteps, NRepl,模型)
z =意味着(v2)
v1 = mortalitysimulation2(α1,φ,sigma1 v0, T, NSteps, NRepl,模型)
l =意味着(v1)
图;
情节(m)
%抓住
%的阴谋(z,“—”)
%抓住
%情节(l '——')
%推迟
% %
B =αβ*
一个= - alpha
sigma1 =σ^ 2
h = 3
ρ= 0.05
H - H = / H
B1 =α1 *φ
A1 =α1
sigma11 = sigma1 ^ 2
k = 1 / h
% t = 1:1:60
d =(ε;0)
D = D * D '
% %
ZZ = k * (1 / A1) * (B1 - k (3/4) * * sigma11 / (A1 ^ 2))
XX = k * (B1 - k (1/2) * * sigma11 / (A1 ^ 2))
CC = k * (B1-k * sigma11 / (A1 ^ 2)) * (1 / A1)
NN = (1/4) * k ^ 2 * sigma11 / (A1 ^ 2) * (1 / A1)
毫米= H * (1 /) * (B - H (3/4) * * sigma1 / (^ 2))
你= H * (B - H (1/2) * * sigma1 / (^ 2))
KK = H * (b - H * sigma1 / (^ 2)) * (1 /)
JJ = 1/4 * ((H ^ 2) * sigma1 / (^ 2)) * (1 /)
信谊t
信谊年代
有趣= @ (s) exp (ZZ-XX * s-CC * exp (A1 * s)神经网络* exp (2 * A1 * (s) - k * (1 / A1) * (1-exp (A1 * s)) * m + MM-LL * (s - t) kk * exp (A * (s - t))剧照* exp (2 * * (s - t)) - h * (1 /) * (1-exp (A * (s - t)) * A *(1/2) *ε^ 2 * (s - t) - k *ρ* s))
int I = (t,有趣,年代,正无穷)

接受的答案

Devineni Aslesha
Devineni Aslesha 2020年8月25日
在给定的代码的方式 int 用于计算的积分 有趣的 关于变量s从正。然而, int 函数是不能解决给定的积分表达式。你可以参考 近似的不定积分 int 函数,它返回未解决的积分。被积函数可以近似函数( 有趣的 )或照顾括号而定义的乐趣,这样就可以得到解决。在计算积分与极限之前,试着计算积分w.r。t变量“年代”。
使用int I =(有趣的)而不是int I = (t,有趣,年代,正无穷))。
一旦解决了积分,然后从t正应用限制。

更多的答案(0)

类别

找到更多的在数值类型帮助中心文件交换

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!