Question

matlab求数值解

求大侠帮忙求下sigma数值解，顺便给下程序学习一下！

Answer 1

matlab为了处理正态分布的概率函数，专门定义一个函数erf，

用这个函数可以减少很多工作量。

回到你的问题，这是一个非线性方程，先要画个图看看0点的大致位置

你从图中可以看出，方程有2个零点，分别在1附近和8附近，这个这就是非线性方程迭代的初值，分别求得两个解为：

x1 =   1.28476428662986
x2 =   7.52563996056696

具体给你matlab代码。

function my_solve
clc;
format long g;
close all;
figure;
x = 0:0.1:20;
y = Fs(x);
plot(x,y);
grid on;
x1 = fsolve( @Fs , 1 ) ;
x2 = fsolve( @Fs , 8 ) ;
x1
x2
end
function y = Fs( sigma )
y=(erf(3/sqrt(2)./sigma)-erf(2/sqrt(2)./sigma))/2-0.05;
end

Answer 2

norm=@(x,sigma)1/2/pi/sigma*exp(-x.^2./2./sigma.^2);
qn=@(sigma) quad(norm,2,3,1e-6,[],sigma)-0.02; %0.05无解
sigma=fzero(qn,1)

追问

你说无解是不是因为积分变量弄错了，应该是qn=@(x) quad(norm,2,3,1e-6,[],x)-0.02;
这才是表示对x积分吧？

追答

1、你的式子和我的式子是等价的。
2、必然是sigma为未知数

quad是对第一个参数积分的，具体请help