1

我有一组数据x(1000,2)

2   18.1199158925616
3   9.56472607516328
7   34.3629828698699
4   30.9865760010538
7   54.4869549747429
8   54.7696742935546
2   1.45832407805944
3   15.4231700658276
8   57.0479651194063
1   5.34909555648445
2   11.5593766410445
2   17.4392978010484
6   25.4072607558367
4   15.1248062077068
2   13.2035176450795
7   33.2507069131440
5   21.6859046133920
3   18.4458236149053
1   0.824962982951128
8   46.6322862041701
8   43.6889420657037
2   17.0930689315346
3   13.9754117410459
2   1.99260423380962
11  86.3172411757665
8   50.6244869641866
1   1.16602336459361
2   12.6604166373800
2   17.4519078987572
5   27.4799515521804

plot(x(:,1),x(:,2),'*r');以散点图的形式绘制: 在此处输入图像描述

但希望normpdf对 x(:,1) 的选定值使用如下所示的内容: 在此处输入图像描述

x(:,2)要在对应于的列中查找值x(:,1)==3

[x3]=x(find(x(:,1)==3),2);
mu=mean(x3);
sig=std(x3);
Y=normpdf(x3,mu,sig);

figure('color',[1,1,1]);
hhh=plot3(3*ones(length(xi)),x(find(x(:,1)==3),2),Y,'-r')
p = polyfit(x(:,1),x(:,2),1);
plin=plot(x(:,1),polyval(p,x(:,1)),'-b');

我的版本看起来像裤子 但是我如何让红色虚线 pdf 成为平滑曲线,因为 '-r' 给了我之字形线。我必须以某种方式订购数据吗?另外,如何使法线曲线中心的红点下降到蓝线?也许最重要的是,作为一个见多识广的读者,您希望看到这种情节还是其他不那么花哨的情节?

4

1 回答 1

1

是的,对x-values 的排序应该会让你看到红色曲线。你也不需要find()在这里。我还没有测试过代码,但这些更改应该可以工作:

[x3]=sort(x(x(:,1)==3,2));
mu=mean(x3);
sig=std(x3);
Y=normpdf(x3,mu,sig);

figure('color',[1,1,1]);
hhh=plot3(3*ones(length(xi)),x3,Y,'-r')
p = polyfit(x(:,1),x(:,2),1);
plin=plot(x(:,1),polyval(p,x(:,1)),'-b');

至于这样一个情节的价值,我猜它会很混乱。但我想它可以用来给出数据分布的一般视图。

于 2013-03-26T15:32:39.283 回答