matlab画积分曲线问题clearfigure;% b=-6;% s=3.5; % u=-17;u=10^(-17/10); %参数u,b,s不确定是上面一组还是这一组b=10^(-6/10);s=10^(3.5/10);x=0:0.01:30; %x也不一定要到30b0=b.^2;d0=s.^2;f0=@(z,x)(1./z).*exp(-((log(z) - u).^2./(2*

来源:学生作业帮助网 编辑:作业帮 时间:2024/05/04 15:35:22
matlab画积分曲线问题clearfigure;% b=-6;% s=3.5; % u=-17;u=10^(-17/10); %参数u,b,s不确定是上面一组还是这一组b=10^(-6/10);s=10^(3.5/10);x=0:0.01:30; %x也不一定要到30b0=b.^2;d0=s.^2;f0=@(z,x)(1./z).*exp(-((log(z) - u).^2./(2*

matlab画积分曲线问题clearfigure;% b=-6;% s=3.5; % u=-17;u=10^(-17/10); %参数u,b,s不确定是上面一组还是这一组b=10^(-6/10);s=10^(3.5/10);x=0:0.01:30; %x也不一定要到30b0=b.^2;d0=s.^2;f0=@(z,x)(1./z).*exp(-((log(z) - u).^2./(2*
matlab画积分曲线问题
clear
figure;
% b=-6;
% s=3.5;
% u=-17;
u=10^(-17/10); %参数u,b,s不确定是上面一组还是这一组
b=10^(-6/10);
s=10^(3.5/10);
x=0:0.01:30; %x也不一定要到30
b0=b.^2;
d0=s.^2;
f0=@(z,x)(1./z).*exp(-((log(z) - u).^2./(2*d0))-(((x.^2)+(z.^2))./ (2*b0))).* besseli(0,x.*z./b0);
f1=arrayfun(@(x)integral(@(z)f0(z,x),0,inf),x);
f=(f1.*x)./(b0.*sqrt(2*pi*d0));
plot(x,f,'r');
大概是因为贝塞尔或者指数函数的关系,不过还是出不来 ,

matlab画积分曲线问题clearfigure;% b=-6;% s=3.5; % u=-17;u=10^(-17/10); %参数u,b,s不确定是上面一组还是这一组b=10^(-6/10);s=10^(3.5/10);x=0:0.01:30; %x也不一定要到30b0=b.^2;d0=s.^2;f0=@(z,x)(1./z).*exp(-((log(z) - u).^2./(2*
问题出在Bessel函数上,具体我再分析一下.
请追问.