[MATLAB数学相关] 概率密度函数求解

[TTT_IOU]

不知下面程序是你所要的：
A=[15 3 17 18 8 30 23 30 28 29 25 26 27 5 8 4 27 13 15 20];
[a,b]=hist(A,3:6:27);
a1=a/length(A);
bar(b,a1);
axis([0,30,0,.45])
set(gca,'ytick',0:.05:.5,'yticklabel','0|10%|15%|20%|25%|30%|35%|40%|45%|50%')
希望有帮助！！！

[kfcrex]

作者: TTT_IOU

不知下面程序是你所要的：
A=[15 3 17 18 8 30 23 30 28 29 25 26 27 5 8 4 27 13 15 20];
[a,b]=hist(A,3:6:27);
a1=a/length(A);
bar(b,a1);
axis([0,30,0,.45])
set(gca,'ytick',0:.05:.5,'yticklabel','0|10%|15%|20%|25%|30%|35%|40%|45%|50%')
希望有帮助！！！

可是这个程序运行出来的结果不对啊，我对Matlab不大清楚，出来的横坐标为什么是3，9，15……27呢？应该是6，12，18……30相对应的概率也不对呀~
不过可以麻烦您解释下hist这个函数吗？最好能举点不同于hlelp文档里的例子，help文档里的hist函数说明看不懂，当初看别人做的时候也是用这个函数，一直没看明白，E文太差了+_+b

[zsy312]

%先编辑函数，a为已知数组，s为待求密度分布
function s=dens(a)
clc
N=10;
len=ceil(max(a)/N);
s=zeros(1,len);
for i=1:len
s(i)=length(find(a>N*(i-1)&a<=N*i));
end
s=s/length(a);
%%%%%%%%%%%%%%%%%%
主程序调用：
y=dens(a);
bar(100*y);%百分率
%区间我设成10了，可以改成别的。

[zsy312]

画图函数及自变量的区间自己改吧，plot可以做曲线，想光滑话的话，可以用interp1函数：
k=ceil(max(a)/10);%要和刚才的区间长度一致
y=interp1((1:k)*10-5,s*100,1:0.02:max(a),'spline');
plot(1:0.02:max(a),y)
axis([0,max(a),0,max(y)])

[kfcrex]

作者: zsy312

%先编辑函数，a为已知数组，s为待求密度分布
function s=dens(a)
clc
N=10;
len=ceil(max(a)/N);
s=zeros(1,len);
for i=1:len
s(i)=length(find(a>N*(i-1)&a<=N*i));
end
s=s/length(a);
%%%%%%%%%%%%%%%%%%
主程序调用：
y=dens(a);
bar(100*y);%百分率
%区间我设成10了，可以改成别的。

我运行了下，出现这个问题
??? Error: File: D:\Matlab7\work\wentishiyan.m Line: 12 Column: 1
Missing variable or function.
话说小弟其实是个初学者，您所写的对我来说就是天书一般的存在啊~
不过这个问题我已经得到答案了，就是用hist函数来写
A=[15 3 17 18 8 30 23 30 28 29 25 26 27 5 8 4 27 13 15 20];
x=3:6:27;
num=hist(A,x);
p=num/20;
plot(x,p)
再次感谢诸位的大力帮助~

[zsy312]

改了函数的文件名就运行不了了，因为你后边要同时该命令才行。
用hist能做出来再好不过了，我也有收获。