公告:1)网站程序升级:Q2A升级到1.8.6,Wordpress升级到5.7.2
2)修复了头像加载慢与提交问题反应慢等问题
2021-06-16

欢迎来到 Mathematica 问答社区

提问时请贴上文本代码

语法高亮:在编辑器中点击

被禁止的话题:广告破解

请阅读:《提问的智慧》

备用域名:mma.ooo

支持LaTex数学公式
行内公式标识符:\$ 或“$\backslash ($”+“$\backslash )$”,
行间公式标识符:\$\$ 或 “$\backslash [$”+“$\backslash ]$”

社区建议QQ群:365716997

分类

0 投票
338 浏览
j1 = 0.9

j3 = 0.9
j4 = 0.1
jz = 0.5
d1 = 0.1
d2 = 0.1

z = 4


r1 = 1/2 (Cos[x] + Cos[y])
r2 = Cos[x]*Cos[y]
a = 2*(j1 - j2 + j2*r2 - j1*r1)*z*sa + d1*(2*sa - 1) + jz*sb
b = 2*(j3 - j4 + j4*r2 - j3*r1)*z*sb + d2*(2*sb - 1) + jz*sa
c = jz*Sqrt[sa*sb]
ecl = (j2 - j1)*n*z*sa^2 - d1*n*sa^2 + (j4 - j3)*n*z*sb^2 - d2*n*sb^2 - jz*n*sa*sb

(*A+B>0*)
w1 = ((a - b)*((a + b)^2 - 4*c^2) + ((a + b)^2 + 4*c^2)*
     Sqrt[(a + b)^2 - 4*c^2])/(2*((a + b)^2 - 4*c^2))
w2 = ((-a + b)*((a + b)^2 - 4*c^2) + ((a + b)^2 + 4*c^2)*
     Sqrt[(a + b)^2 - 4*c^2])/(2*((a + b)^2 - 4*c^2))
q1 = w1/(Exp[w1/t] - 1)
q2 = w2/(Exp[w2/t] - 1)

sa = 0.5
sb = 0.5
j2 = 0.3

(*A+B>0*)
D[(1/(4*Pi^2))*NIntegrate[q1 + q2, {y, -Pi, Pi}, {x, -Pi, Pi}], t]

求助:1、最后一句求导中如果用NIntegrate,会导致出现以下错误并返回0值。N

Integrate::inumr: The integrand ((0. +4. (0.8 +Times[<<3>>]+Times[<<2>>])-4. (0.6 +Times[<<3>>]+Times[<<2>>])) ..... Plus[<<3>>]+4. Plus[<<3>>])^2)) has evaluated to non-numerical values for all sampling points in the region with boundaries {{0,3.14159},{0,3.14159}}.

怀疑应该是变量t没赋值导致的,但是我得对t求导,所以不能赋值啊,是不是就无法使用NIntegrate?

2.最后我用Integrate,但是一直无法运算出来,有什么便捷的方法么?

(PS:我最后只需要一个数值,不需要式子,所以才尝试用NIntegrate)

纠结了我好久了。crying

分类:函数 | 用户: 535 (36 分)
NIntegrate 里面不可以有未知符号的。
但由于 q1 + q2 太复杂,所以Integrate运行时间很长(很可能无法计算出结果)
所以还是想办法化简要求解的问题吧
化简也很难啊,我想到头都大了,都不知道咋化简。崩溃!~

登录 或者 注册 后回答这个问题。

...