作业帮matlab中积分的问题
来源：学生作业帮 编辑：作业帮 分类：数学作业 时间：2024/07/18 11:01:43
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);
上式可以进行积分
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);
上式可以进行积分
你的这个问题,错在这一句：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的定义.
“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的定义.