| 查看: 1949 | 回复: 4 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
pyypyy37金虫 (正式写手)
|
[求助]
滑动T检验 Matlab程序 检验一个降雨序列的变异。求助！谢谢
|
||
|
求滑动T检验程序和资料。 我想绘图，类似于下图这种，不知道怎么编程。应该不难，但是没有接触过，还望有人能指点一二。谢谢各位高手。  @QCL($DK(8)B_5G{RL4)IR1.jpg |
回复此楼» 猜你喜欢
职称评审没过，求安慰
已经有56人回复
-
最近几年招的学生写论文不引自己组发的文章
已经有5人回复
-
26申博自荐
已经有3人回复
-
A期刊撤稿
已经有4人回复
pyypyy37
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1055.6
- 散金: 153
- 红花: 2
- 沙发: 13
- 帖子: 959
- 在线: 40.6小时
- 虫号: 1102381
- 注册: 2010-09-18
- 专业: 环境工程
送红花一朵 |
我想问一下，这只是绘图的命令吧？我把我的程序发给你。麻烦你帮我看看，如果要画出我贴出的那个图，还要加什么命令。谢谢！！你太好了！检测变异是不是看h值。0的地方就是变异点？？再次感谢。   clear Xn=load ('E:\pyy\shuju\rain.txt'); % X=X./max(X); X=Xn'; N=length(X); m=0; for i=2:N-2 left=X(1:i); right=X(i+1:N); N_left=length(left); N_right=length(right); mu_left=mean(left); mu_right=mean(right); std_left=std(left); std_right=std(right); Sw=sqrt(((N_left-1)*std_left^2+(N_right-1)*std_right^2)/N_left+N_right+1); T(i)=(mu_left-mu_right)/(Sw*sqrt(1/N_left+1/N_right)); if abs(T(i))>0.025 m=m+1; a(m)=X(i); b(m)=T(i); c(m)=i; end [h(i),significance(i),ci]=ttest2(left,right,0.05); |
3楼2013-04-27 09:26:39
libralibra
至尊木虫 (著名写手)
骠骑将军
- 程序强帖: 40
- 应助: 817 (博后)
- 金币: 12914.1
- 红花: 64
- 帖子: 2238
- 在线: 287.3小时
- 虫号: 696514
- 注册: 2009-02-05
- 专业: 计算机软件
2楼2013-04-26 23:43:33