| 查看: 1358 | 回复: 10 | |||
[求助]
Matlab三重积分求助！！！！！！！！！！！！！！！！
|
|
对由x0 x1 y1 y0 z0 z1构成的立方体的整个区域积分，代码如下，为何算出来结果不对，谁知道怎么回事，该如何解决，谢谢了！！！ function F(m,n,q) clear syms x y z ; bohr=0.529177; Pi=3.1415; a=3/bohr; alat=5.43/bohr; x0=277.05286; x1=482.27720; y1=379.66503; y0=0.0000000; z0=0.0000000; z1=112.87339; X=0.5*(x0+x1); Y=0.5*(y0+y1); Z=0.5*(z0+z1); w=(x-X)^2+(y-Y)^2+(z-Z)^2; f=(2*Pi/a)^3*exp(-w/a^2); int1=int(f,z,z1,z0); int2=int(int1,y,y1,y0); int3=int(int2,x,x1,x0) |
回复此楼» 收录本帖的淘贴专辑推荐
 matlab典型案例 |
» 猜你喜欢
我标书代码变了
已经有45人回复
-
博士毕业高校和就业的相关问题
已经有10人回复
-
山东大学(青岛)“天然药物生物智造”课题组 招聘“博士后”(年薪20.4-55.6万元)
已经有11人回复
-
34岁读博士晚吗
已经有28人回复
-
硕博巨婴，也许才刚刚开始
已经有83人回复
-
2024国社科通讯评审
已经有18人回复
-
最新消息：2024国自然下载文件名变了
已经有13人回复
-
成都产品质量检测研究院
已经有10人回复
-
审稿人含糊拒稿，还需要回复吗？如何回复？
已经有18人回复
-
F口信息学部拿面上，大概需要什么样的成果
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助matlab问题
已经有5人回复
- 继续求助Matlab绘图！
已经有8人回复
- matlab 积分求助
已经有8人回复
- 求助四重积分的matlab求解程序
已经有5人回复
- 求助matlab编程
已经有13人回复
- 【求助】利用MATLAB求解的外围轮廓包络线为什么不与图形外围轮廓重合？怎样才可以精确
已经有8人回复
- 【求助】matlab求解一个简单的积分（急！！！）
已经有7人回复
- 求助！！关于MATLAB
已经有6人回复
- 【求助】三重积分求助
已经有5人回复
- 【求助】matlab 解矩阵方程 重金悬赏解决方案
已经有10人回复
- 【求助】积分微分方程matlab求解
已经有6人回复
- 【求助】初学matlab
已经有11人回复
- 【求助】matlab 积分的问题
已经有5人回复
- 【求助】三重积分
已经有9人回复
2楼2012-09-06 17:09:11
3楼2012-09-06 17:51:54
4楼2012-09-07 09:18:34
5楼2012-09-07 10:46:33
6楼2012-09-07 10:48:46
7楼2012-09-07 17:05:44
|
恩，这我知道，我把int用triplequad替换其他不变，如下 w=(x-X)^2+(y-Y)^2+(z-Z)^2; F=@(x,y,z)(2*Pi/a)^3*exp(-w/a^2); Q=triplequad(F,pssn_x0,pssn_x1,pssn_y0,pssn_y1,pssn_z0,pssn_z1) 但是运行后出现以下提示：该如何修改，我对这块不熟悉，希望你能帮下忙。先谢谢你了。 ??? Error using ==> innerintegral Inputs to innerintegral must be floats, namely single or double. Error in ==> triplequad>innerintegral at 66 Q = zeros(size(y), superiorfloat(fcl, xmax, y, z, varargin{:})); Error in ==> dblquad>innerintegral at 80 fcl = intfcn(xmin, y(1), varargin{:}); %evaluate only to get the class below Error in ==> quad at 77 y = f(x, varargin{:}); Error in ==> dblquad at 60 Q = quadf(@innerintegral, ymin, ymax, tol, trace, intfcn, ... Error in ==> triplequad at 53 Q = dblquad(@innerintegral, ymin, ymax, zmin, zmax, tol, quadf, intfcn, ... Error in ==> F at 32 Q=triplequad(F,pssn_x0,pssn_x1,pssn_y0,pssn_y1,pssn_z0,pssn_z1) |
8楼2012-09-07 21:07:17
walk1997
金虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 4676.2
- 红花: 22
- 帖子: 1066
- 在线: 797.8小时
- 虫号: 416039
- 注册: 2007-06-29
- 性别: GG
- 专业: 粒子物理学和场论
9楼2012-09-13 10:43:09
walk1997
金虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 4676.2
- 红花: 22
- 帖子: 1066
- 在线: 797.8小时
- 虫号: 416039
- 注册: 2007-06-29
- 性别: GG
- 专业: 粒子物理学和场论
|
In[1]:= X = (x0 + x1)/2; Y = (y0 + y1)/2; Z = (z0 + z1)/2; w = (x - X)^2 + (y - Y)^2 + (z - Z)^2; f = (2*Pi/a)^3*exp (-w/a^2); Integrate[f, {x, x0, x1}, {y, y0, y1}, {z, z0, z1}] 积分结果 我没copy错吧？ Out[6]= (1/(3 a^5))2 exp \[Pi]^3 (x0 - x1) (y0 - y1) ((x0 - x1)^2 + (y0 - y1)^2 + (z0 - z1)^2) (z0 - z1) |
10楼2012-09-13 10:44:54
10