matlab 求已知概率密度函数的随机数生成

N=10000; %需要随机数的个数
a=zeros(N,1); %存放随机数的数列
n=0;
f1=@(t) 1./(1.2*pi*(1+5*(t-7.3).^2));
f2=@(t) 1./(1.05*pi*(1+6*(t-8.2).^2));
tt=linspace(0,24,1000);
ff=f1(tt).*(tt<=7.5)+f2(tt).*(tt>7.5);%根据公式计算概率密度
s=trapz(tt,ff);  %计算整个区间概率密度的积分
ff=ff/s;         %归一化概率密度while n<Nt=rand(1)*24;%生成[0,24]均匀分布随机数if t<=7.5f=f1(t)/s;elsef=f2(t)/s;end         %计算对应密度函数值f(t)r=rand(1);  %生成[0,1]均匀分布随机数if r<=f     %如果随机数r小于f(t)，接纳该t并加入序列a中n=n+1;a(n)=t;end
end%以上为生成随机数列a的过程，以下为统计检验随机数列是否符合分布
num=100;         %分100个区间统计
[x,c]=hist(a,num);    %统计不同区间出现的个数
dc=24/num;        %区间大小
x=x/N/dc;         %根据统计结果计算概率密度bar(c,x,1); hold on;  %根据统计结果画概率密度直方图
plot(tt,ff,'r'); hold off; %根据公式画概率密度曲线

matlab 求已知概率密度函数的随机数生成相关推荐

matlab 求函数绝对值,概率密度函数估计matlab命令 A a abs 绝对值、(4)
isnumeric 若是数值数组则为真 isobject 若是对象则为真 isprime 若是质数则为真 isreal 若是实数则为真 isspace 若是空格则为真 issparse 若是稀疏矩阵则 ...
用matlab求已知分布的概率密度函数,已知某分布地概率密度函数，怎样得到服从该分部的随机数...
我想用谢中华老师的<MATLAB统计分析与应用39个案例分析>中的crnd模块来求. 该模块程序为: function y = crnd(pdffun, pdfdef, m, n) %生成 ...
matlab状态空间法算反馈阵,matlab中已知系统的状态方程怎样绘制系统阶跃响应曲线...
matlab中,怎样设置坐标系统的原点在左上角? axisijx=-pi:0.1:pi;y=sin(x);plot(x,y,'g:'); 谁能帮我做Matlab的作业:分析系统能否通过状态反馈实现镇定 ...
matlab已知函数表达式画函数图像,怎么用matlab画已知函数表达式的一个函数图像？函数比较复杂的……...
点击查看怎么用matlab画已知函数表达式的一个函数图像?函数比较复杂的--具体信息答:要具体哦,如果一个自变量,函数图像就是曲线,用plot 如果是2个自变量,函数图像就是曲面了,如f=3*exp ...
[MATLAB]从已知矩阵中取出子阵
首先要知道冒号表示法,如下: EX: 如上有4×6矩阵,假设我们要取出第一行,第二行,第四行,则可以有B=A([1 2 4],:) 假设要取出第一行,第二行,第四行的1-2列和4-6列,可以有B=A( ...
matlab 椭圆参数传递,用MATLAB绘图已知椭圆参数方程
用MATLAB绘图已知椭圆参数方程为 cos ,[0,2] sin x a t t y b t π= ? ∈ ? = ? 一.实验名称 MATLAB绘图功能二.实验目的及要求 1.理解MATLAB绘 ...
matlab 参数方程绘图,用matlab绘图已知椭圆参数方程
用matlab绘图已知椭圆参数方程 1 用MATLAB绘图已知椭圆参数方程为 cos , [0,2 ] sin x a t t y b t        一.实验名称 MATLAB绘图功 ...
用C语言求已知集合的笛卡儿乘积（离散数学）
用C语言求已知集合的笛卡儿乘积(离散数学) 实验要求: 通过编程实现求给定集合A和B的笛卡儿乘积C(C=A×B)D(D=A×B×A)的运算. 实验内容已知所给集合A和B,求A与B的笛卡儿乘积C(C= ...
Matlab利用已知方阵创建斜对角矩阵
Matlab利用已知方阵创建斜对角矩阵 A=[ 1 2 3 4 ;7 8 9 10 ;13 14 15 16 ;19 20 21 22 ]; 对于小矩阵,直接利用for循环解决 n = size(A, ...

matlab 求已知概率密度函数的随机数生成

matlab 求已知概率密度函数的随机数生成相关推荐

最新文章

热门文章