楼主: 164544298
3878 1

[程序分享] 变上限二重积分 MATLAB [推广有奖]

  • 0关注
  • 4粉丝

已卖:1份资源

本科生

39%

还不是VIP/贵宾

-

威望
0
论坛币
15 个
通用积分
1.6392
学术水平
0 点
热心指数
0 点
信用等级
0 点
经验
152 点
帖子
27
精华
0
在线时间
154 小时
注册时间
2017-4-5
最后登录
2025-7-10

楼主
164544298 在职认证  发表于 2019-1-27 13:45:19 |AI写论文
30论坛币
变上限二重积分 假定这个二重积分的值,要求解积分上限,公式中的v=5.2823,ρ=[1 0.9774;0.9774 1], ej=0.5840, X=[xi;xj]
我利用interal2和fzero函数会出现报错,大家可以帮我看一错在哪了吗?

v=5.2823

p=[1 0.9774;0.9774 1]

X=[xi;xj]

syms xi xj ei

f=@(xj,xi) ((det(p)^(-0.5))*gamma((v+2)/2)*(1+(1/v)*X’*inv(p)*X)^(-(v+2)/2) )/(gamma(v/2)*v*pi)

[ei,fval,exitflag] = fzero(@(ei)integral2(f,-Inf, 0.5840, -Inf, ei) - 0.95,-1e3)

会出现报错:

错误使用 fzero (line 306)

FZERO 无法继续,因为用户提供的 function_handle ==> @(ei)integral2(f,-Inf,0.5840,-Inf,ei)-0.95 失败,出现下面的错误。


输入函数必须返回 'double' 或 'single' 值。找到 'sym'。




沙发
164544298 在职认证  发表于 2019-1-27 18:05:31
琢磨出来了,把矩阵形式转化为数字的形式式就可以了

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群
GMT+8, 2025-12-26 06:49