| 查看: 403 | 回复: 0 | |||
[交流]
求解二重积分，总出现各种各样问题！！！！跪求大神指导！！！！！！！
|
|
求二重积分没办法求出来，总是显示，??? Inputs to dblquad must be floats, namely single or double. Error in ==> dblquad>innerintegral at 85 Q = zeros(size(y), superiorfloat(fcl, xmax, y)); Error in ==> quad at 62 y = f(x, varargin{:}); Error in ==> dblquad at 64 Q = feval(quadf, @innerintegral, ymin, ymax, tol, trace, intfcn, ... 代码如下，烦请各位大神给看看到底是出了什么问题呀？急求！！！！ format long g data0=[2556 2592 2688 2832 3087 3408 3552 4124 4248 4296 4368 4872 5064 5376 5568 6744 7096 7248 7392 7584 7800 7920 8520 8880 8904 9210 9264 9480 10032 10056 10115 11496 11520 11808 12000 12168 12186 12264 12408 12528 12624 14544 14928 15144 15336 15768 16176 16176 16440 16656 16800 17160 17232 17568 17592 17760 17832 18120 18576 19008 20112 20232 20256 20544 20664 21168 21432 21576 22872 22992]; [parmhat0, parmci0] = wblfit(data0); mbar = 2.10815930057601; sigmam = 0.0273906928410007; nbar =13604.6417116356; sigman =503437.078749078; c=70; d=prod(data0); e=sum(data0); syms t int(int((x^c)*d^(x-1))/(2*pi*sigmam*sigman*y^(c*x))*exp(-(t/y)^x-e/(y^x)-(y-nbar)^2/2*sigman-(x-mbar)^2/2*sigmam), 0, inf),0,inf) dblquad(@(x,y)((x^c)*d^(x-1))/(2*pi*sigmam*sigman*y^(c*x))*exp(-(t/y)^x-e/(y^x)-(y-nbar)^2/2*sigman-(x-mbar)^2/2*sigmam), 0, inf,0,inf) |
回复此楼» 猜你喜欢
津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生一名
已经有0人回复
-
【原创讨论】从电子约束到物质编辑：一套可迭代的环形磁场科技树
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有121人回复
-
【方案分享】单环磁场+轴心控制+偏转导出电子束约束系统（可行性实验）
已经有6人回复
-
【修正版】单环用磁约束低速电子实验方案（简化版）
已经有0人回复
-
桂林理工大学物理学专业招收调剂，还有三个名额！！！
已经有22人回复
-
考博自荐
已经有4人回复
-
山东大学第二批博士研究生招生
已经有0人回复
-
中国科学院东莞材料科学与技术研究所-2026年博士招生-吴昊研究员-磁学与自旋电子学
已经有0人回复
-
《电磁学》教材推荐
已经有1人回复
-
【急招】合肥工大核聚变材料计算方向2026级工程博士生
已经有4人回复