| 查看: 2150 | 回复: 4 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
pyypyy37金虫 (正式写手)
|
[求助]
滑动T检验 Matlab程序 检验一个降雨序列的变异。求助！谢谢
|
||
|
求滑动T检验程序和资料。 我想绘图，类似于下图这种，不知道怎么编程。应该不难，但是没有接触过，还望有人能指点一二。谢谢各位高手。  @QCL($DK(8)B_5G{RL4)IR1.jpg |
回复此楼» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有65人回复
-
时间戳又变了8-15
已经有23人回复
-
filecode=后面第一个是大写字母
已经有7人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有10人回复
-
filecode
已经有12人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有13人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有7人回复
-
关于Filecode分析方法
已经有13人回复
-
各位道友，我要去昆明玩几天，回来见。
已经有8人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有22人回复
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
