作业帮matlab 做积分函数曲线
来源：学生作业帮 编辑：拍题作业网作业帮 分类：数学作业 时间：2024/04/25 23:26:00
matlab 做积分函数曲线
函数如图：
想画出此函数曲线 r为横轴
自己写了一个 不对 求改正
x=0:0.01:15;
mu=0;
b=1;
s=1;
b0=b.^2;
d0=s.^2;
f0=@(z)(1./z).*exp(-((log(z) - mu).^2./(2*d0))-(((x.^2)+(z.^2))./ (2*b0)));
f1=quad(f0,0,inf);
f=(f1.*x)./(b0.*sqrt(2*pi*d0));
plot(x,f);
函数如图：
想画出此函数曲线 r为横轴
自己写了一个 不对 求改正
x=0:0.01:15;
mu=0;
b=1;
s=1;
b0=b.^2;
d0=s.^2;
f0=@(z)(1./z).*exp(-((log(z) - mu).^2./(2*d0))-(((x.^2)+(z.^2))./ (2*b0)));
f1=quad(f0,0,inf);
f=(f1.*x)./(b0.*sqrt(2*pi*d0));
plot(x,f);
1、quad的积分表达式中
2、quad的积分限不能为无穷大,换用integral函数（2012a以上版本）.
把
f0=@(z)(1./z).*exp(-((log(z) - mu).^2./(2*d0))-(((x.^2)+(z.^2))./ (2*b0)));
f1=quad(f0,0,inf);
改成
f0=@(z,x)(1./z).*exp(-((log(z) - mu).^2./(2*d0))-(((x.^2)+(z.^2))./ (2*b0)));
f1=arrayfun(@(x)integral(@(z)f0(z,x),0,inf),x);
2、quad的积分限不能为无穷大,换用integral函数（2012a以上版本）.
把
f0=@(z)(1./z).*exp(-((log(z) - mu).^2./(2*d0))-(((x.^2)+(z.^2))./ (2*b0)));
f1=quad(f0,0,inf);
改成
f0=@(z,x)(1./z).*exp(-((log(z) - mu).^2./(2*d0))-(((x.^2)+(z.^2))./ (2*b0)));
f1=arrayfun(@(x)integral(@(z)f0(z,x),0,inf),x);