Question

matlab中积分的问题

p=linspace(0,0.05,100);
b=（很长的关于p的公式）
t=b.^2*pi^2*10^(-3)*r^4*0.25.*p/(r^2+p.^2);
f=trapz(p,t);
这样的t关于p的积分会显示有错误，但是当我将t等式中的p提出来单独乘一遍（如下）就可以积分了。求大神指导。
p=linspace(0,0.05,100);
b=……；
t=b.^2*pi^2*10^(-3)*r^4*0.25/(r^2+p.^2);
t=t.*p;
f=trapz(p,t);
上式可以进行积分

Answer 1

你的这个问题，错在这一句：t=b.^2*pi^2*10^(-3)*r^4*0.25.*p/(r^2+p.^2);
“b.^2*pi^2*10^(-3)*r^4*0.25.*p” 是个向量吧？“r^2+p.^2”也是向量吧？
那应该用“./”呀，你前面都用了，怎么到这儿粗心了呢。
改一下就可以了：t=b.^2*pi^2*10^(-3)*r^4*0.25.*p./(r^2+p.^2);

你后面的运算没有报错，是因为用“/”也能算出t的值来的
相当于2个向量相除，解出的t值只有一个值，并不是个向量
而你后面再用t=t.*p;运算，相当于人为又把t变成向量了
但不是t本来的值。虽然没报错，但结果应该不是很准确
也有可能当r取一定值时，结果比较接近。
反正，我试了一下，2种积分的结果相差倒是不大，主要是不知道你b的定义。

Answer 2

% 向量乘除的表示问题 t=b.^2*pi^2*10^(-3)*r^4*0.25.*p/(r^2+p.^2); 改为
t=b.^2.*pi^2.*10^(-3).*r^4.*0.25.*p./(r^2+p.^2);