求助:关于TRAPZ函数
Z = TRAPZ(X,Y,DIM) or TRAPZ(Y,DIM) integrates across dimension DIMof 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 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;
=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是在做什么,麻烦指点一下,谢谢!! trapz函数是采用复化梯形公式求积分,其应用格式为
I=trapz(x,y)
其中x分区间离散化数据点构成的向量,y与x维的向量,
返回值是积分的近似值... 先谢拉,这个我知道,我不明白的是
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 编辑 ] 应该是对第三维积分,你生成的F是一个三维矩阵
[ 本帖最后由 ChaChing 于 2009-12-5 15:55 编辑 ] Okey,I get it,thank you.
[ 本帖最后由 ChaChing 于 2009-12-5 15:55 编辑 ] 回复 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;
=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); 回复 6 # sdutnuaa 的帖子
建议一步步自行试看看!
sqrt(pi*T)
exp(-(X-TAU).^2/4/2^2./T)
sqrt(pi*T).*exp(-(X-TAU).^2/4/2^2./T)
... 回复 6 # sdutnuaa 的帖子
Ref 5F, 常见的程序出错问题整理 http://forum.vibunion.com/thread-46001-1-1.html
(from http://forum.vibunion.com/home-space-uid-63979-do-blog-id-18250.html) 回复 7 # ChaChing 的帖子
我是刚学matlab的菜鸟,我那个程序是编写的一维热传导方程在初始条件为一三角形的下的传导过程,能否详细指点下? 回复 9 # sdutnuaa 的帖子
aa=sqrt(pi*T).*exp(-(X-TAU).^2/4/2^2./T); 为10*101*1001
theata为1*1001
请问两者如何点乘?:@) 回复 10 # ChaChing 的帖子
恩 我也在考虑这个问题 谢谢我再自己好好研究下
页:
[1]