声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2295|回复: 1

[综合讨论] 关于正态分布绘图及其参数估计

[复制链接]
发表于 2007-9-1 11:59 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
这个是代码,希望大家看看,不知道小弟还有什么改进的地方,小弟刚学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 编辑 ]

test1.dat

65 Bytes, 下载次数: 27

这个是TEST1.DAT

回复
分享到:

使用道具 举报

 楼主| 发表于 2007-9-3 21:23 | 显示全部楼层
问题1:如果把坐标画在横线处,以40%60%等比较。。。。

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

[ 本帖最后由 eight 于 2007-9-3 22:01 编辑 ]
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-29 01:20 , Processed in 0.072387 second(s), 21 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表