关于正态分布绘图及其参数估计

lcczx

这个是代码，希望大家看看，不知道小弟还有什么改进的地方，小弟刚学MATLAB滴~
PS：这里有一个问题就是，Y坐标，我想在画线的地方标上40%60%这些概率，而不是等距的坐标值。
先谢谢各位大侠了，呵呵~~

1. clc;clf;clear           %清楚原始数据干扰
   
2. a=0;                    
   
3. b=15;                   %初始变量
   
4. for k=1:1:15
   
5. i=[0.0001 0.001 0.01 0.05 0.1 0.2 0.4 0.5 0.6 0.8 0.9 0.95 0.99 0.999 0.9999];
   
6.     y1=norminv(i(k));
   
7.     x1=[a,b];
   
8.     plot(x1,[y1,y1])
   
9.     hold on
   
10. end
    
11. title('正态坐标纸');         %绘制坐标图
    
12. load test1.dat              %数据导入
    
13. m=size(test1);n=m(2);       %求数据输入的数目
    
14. for k=1:1:n
    
15.     l(k)=(k-0.3)/(n+0.4);   %求失效概率中位秩
    
16. end
    
17. for i=1:1:n
    
18.     plot(test1(2,i),norminv(l(i)),'o')  %描点
    
19. hold on
    
20. end
    
21. x=test1(2,:);                  %赋值给x
    
22. y=norminv(l(:)');             %赋值给y
    
23. avgx=sum(x)/n;              %x平均数
    
24. avgy=sum(y)/n;              %y平均数
    
25. B=(sum(x.*y)-n*avgx*avgy)/(sum(x.*x)-n*avgx*avgx);A=avgy-B*avgx;    %最小二乘法求线性直线
    
26. for p=0:0.01:15
    
27.     plot(p,A+B*p)
    
28. end                         %绘制直线
    
29. rol=(sum(x.*y)-n*avgx*avgy)/(sqrt((sum(x.*x)-n*avgx*avgx)*(sum(y.*y)-n*avgy*avgy)));    %求相关系数
    
30. rola=2.576/sqrt(n+1);                                                                   %求标准相关系数
    
31. if rol>=rola                %判断是否线性相关
    
32.     disp('两者线性相关')
    
33.     cu='期望等于    ';cr='标准差等于     ';cum='Kh';
    
34.     cuu=(0.5-A)/B;crr=(-A+0.841)/B-(-A+0.5)/B;
    
35.     disp([cu,num2str(cuu),cum])
    
36.     disp([cr,num2str(crr)])     %参数估计
    
37. else
    
38.     disp('两者线性不相关')       %非线性相关，quit
    
39. end

复制代码

[ 本帖最后由 lcczx 于 2007-9-1 14:59 编辑 ]

lcczx

问题1：如果把坐标画在横线处，以40%60%等比较。。。。

问题2：现在想COMPLIER成EXE文件，MCC不可以了，因为是脚本文件；用MATCOM也没有成功。。。说是没有找到*.m文件或者文件没有起作用。。。

[ 本帖最后由 eight 于 2007-9-3 22:01 编辑 ]