matlab概率⽣成函数求概率，已知某概率密度函数，如何产⽣
⼀服从该分布的随机数...
谢⽼师：
我想⽤您说的这个程序来求解问题。您编写的程序的M⽂件就是下⾯的这个crnd函数。我将我⾃⼰编的概率密度函数放⼊后出现了好多问题，不知道该如何解决，不知道您能不能帮忙看⼀下，⾮常感谢！
function y = crnd(pdffun, pdfdef, m, n)
%⽣成任意⼀元连续分布随机数
%  y = crnd(pdffun, pdfdef, m, n)，产⽣指定⼀元连续分布的随机数，m⾏n列。pdffun为密度
%  函数表达式，pdfdef为密度函数定义域。注意：pdfdef只能是有限区间，若密度函数定义域为⽆限区
%  间，应设成⽐较⼤的有限区间，例如[-10000,10000]
%
%  example:
%  pdffun = 'x*(x>=0 & x<1)+(2-x)*(x>=1 & x<2)';
%  x = crnd(pdffun, [0 2], 1000, 1);
%  [fp,xp] = ecdf(x);
%  ecdfhist(fp, xp, 20);
%  hold on
matlab生成随机数%  fplot(pdffun, [0 2], 'r')
%
%  Copyright 2009 - 2010 xiezhh.
%  $Revision: 1.0.0.0 $  $Date: 2009/12/17 12:10:00 $
fun = vectorize(['(' pdffun ')' '*x']);    % x*f(x)运算向量化
try              %try模块
xm = quadl(fun, min(pdfdef), max(pdfdef));    % 计算x的数学期望xm，对fun进⾏积分，后⾯两个参数分别为上下限
catch
xm = mean(pdfdef);    % 计算定义区间的平均值xm
end
pdffun = vectorize(['(' pdffun ')' '*x/x']);    % x*f(x)/x运算向量化
cdfrnd = rand(m*n, 1);    % 产⽣[0,1]上均匀分布随机数
y = zeros(m*n, 1);        % 产⽣0向量作为变量y的初值
options = optimset;      % 产⽣⼀个控制迭代过程的结构体变量
options.Display = 'off';  % 不显⽰中间迭代过程
% 通过循环计算指定⼀元连续分布的随机数
for i = 1:m*n
funcdf = @(x)[quadl(pdffun, min(pdfdef), x) - cdfrnd(i)];
y(i) = fsolve(funcdf, xm, options);
end
y=reshape(y,[m,n]);    % 把向量y变为矩阵
在commend window中，有如下程序：
>> x=-10:0.1:10;
u1=0;
u2=0;
q=5;
b=q/sqrt(2);
pdffun='0.4*(1/(q*sqrt(2*pi)))*exp(-(x-u1).^2/2*q^2)+0.6*(1/(2*b))*exp(-abs(x-u2)/b)'; % 调⽤crnd函数⽣成100个服从指定⼀元连续分布的随机数
x = crnd(pdffun, [-10 10], 100, 1);
[fp,xp] = ecdf(x);    % 计算经验累积概率分布函数值
ecdfhist(fp,xp,20);  % 绘制频率直⽅图
hold on
fplot(pdffun, [-10 10], 'r')    % 绘制真实密度函数曲线
xlabel('x');      % 为X轴加标签
ylabel('f(x)');    % 为Y轴加标签
legend('频率直⽅图', '密度函数曲线')    % 为图形加标注框
运⾏后出现了很多问题：
Error using ==> inline.feval
Not enough inputs to inline function.
Error in ==> quadl at 64
y = feval(f,x,varargin{:}); y = y(:).';
Error in ==> crnd>@(x)[quadl(pdffun, min(pdfdef), x) - cdfrnd(i)] at 32
funcdf = @(x)[quadl(pdffun, min(pdfdef), x) - cdfrnd(i)];
Error in ==> fsolve at 180
fuser = feval(funfcn{3},x,varargin{:});
Error in ==> crnd at 33
y(i) = fsolve(funcdf, xm, options);

版权声明：本站内容均来自互联网，仅供演示用，请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系QQ:729038198，我们将在24小时内删除。