求助：关于TRAPZ函数

jimin

Z = TRAPZ(X,Y,DIM) or TRAPZ(Y,DIM) integrates across dimension DIM
    of Y. The length of X must be the same as size(Y,DIM)).

    Example: If Y = [0 1 2
                     3 4 5]

    then trapz(Y,1) is [1.5 2.5 3.5] and trapz(Y,2) is [2
                                                        8];

xx=-10:0.5:10;
tt=0.01:0.1:1;
tau=0:0.01:1;
a=2;
[X,T,TAU]=meshgrid(xx,tt,tau);
F=1/2/2./sqrt(pi*T).*exp(-(X-TAU).^2/4/2^2./T);
js=trapz(F,3);
不大明白dim的含义，js=trapz(F,3);
这句中dim＝3是在做什么，麻烦指点一下，谢谢！！

fanghuikeer

trapz函数是采用复化梯形公式求积分，其应用格式为
I=trapz(x,y)
其中x分区间离散化数据点构成的向量，y与x维的向量,
返回值是积分的近似值．．．

jimin

先谢拉，这个我知道，我不明白的是
Z = TRAPZ(X,Y,DIM) or TRAPZ(Y,DIM) integrates across dimension DIM
    of Y. The length of X must be the same as size(Y,DIM)).
这里面dim的用法，如dim=3

[ 本帖最后由 ChaChing 于 2009-12-5 15:54 编辑 ]

eight

应该是对第三维积分，你生成的F是一个三维矩阵

[ 本帖最后由 ChaChing 于 2009-12-5 15:55 编辑 ]

jimin

Okey,I get it,thank you.

[ 本帖最后由 ChaChing 于 2009-12-5 15:55 编辑 ]

sdutnuaa

回复 1 # jimin 的帖子

clc
xx=-10:0.2:10;
tt=0.01:0.1:1;
tau=-0.5:0.001:0.5;
theata=2*(-tau./abs(tau)).*tau+1.1;
a=2;
[X,T,TAU]=meshgrid(xx, tt,tau);
F=theata.*1/2/2./sqrt(pi*T).*exp(-(X-TAU).^2/4/2^2./T);
js=0.2*trapz(tau,F);
waterfall(X(:,:,1),T(:,:,1),js)
运行后出现下列状况 是什么意思？该如何解决？谢谢啦
??? Error using ==> rdivide
Matrix dimensions must agree.

Error in ==> chuanre at 8
F=theata.*1/2/2./sqrt(pi*T).*exp(-(X-TAU).^2/4/2^2./T);

ChaChing

回复 6 # sdutnuaa 的帖子

建议一步步自行试看看!
sqrt(pi*T)
exp(-(X-TAU).^2/4/2^2./T)
sqrt(pi*T).*exp(-(X-TAU).^2/4/2^2./T)
...

ChaChing

回复 6 # sdutnuaa 的帖子

Ref 5F, 常见的程序出错问题整理 http://forum.vibunion.com/thread-46001-1-1.html
(from http://forum.vibunion.com/home-s ... -blog-id-18250.html)

sdutnuaa

回复 7 # ChaChing 的帖子

我是刚学matlab的菜鸟，  我那个程序是编写的一维热传导方程在初始条件为一三角形的下的传导过程，能否详细指点下？

ChaChing

回复 9 # sdutnuaa 的帖子

aa=sqrt(pi*T).*exp(-(X-TAU).^2/4/2^2./T); 为10*101*1001
theata为1*1001
请问两者如何点乘?:@)

sdutnuaa

回复 10 # ChaChing 的帖子

恩 我也在考虑这个问题 谢谢  我再自己好好研究下