24小时热门版块排行榜

返回列表

pyypyy37

金虫 (正式写手)

应助: 0 (幼儿园)
金币: 1055.6
散金: 153
红花: 2
沙发: 13
帖子: 959
在线: 40.6小时
虫号: 1102381
注册: 2010-09-18
专业: 环境工程

[求助] 滑动T检验 Matlab程序检验一个降雨序列的变异。求助！谢谢

求滑动T检验程序和资料。
我想绘图，类似于下图这种，不知道怎么编程。应该不难，但是没有接触过，还望有人能指点一二。谢谢各位高手。

@QCL($DK(8)B_5G{RL4)IR1.jpg

回复此楼

» 猜你喜欢

英文综述是否需要润色及查重已经有5人回复
救命帖已经有5人回复
限项规定已经有5人回复
为什么nbs上溴没有产物点出现呢已经有9人回复
招博士已经有3人回复
存款400万可以在学校里躺平吗已经有35人回复
最失望的一年已经有18人回复
求推荐博导已经有4人回复
求推荐英文EI期刊已经有5人回复

1楼 2013-04-26 22:22:11

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

» 本帖已获得的红花（最新10朵）

pyypyy37

matlab/VB/python/c++/Java写程序请发QQ邮件:790404545@qq.com

2楼2013-04-26 23:43:33

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

pyypyy37

金虫 (正式写手)

应助: 0 (幼儿园)
金币: 1055.6
散金: 153
红花: 2
沙发: 13
帖子: 959
在线: 40.6小时
虫号: 1102381
注册: 2010-09-18
专业: 环境工程

送红花一朵

引用回帖:

2楼: Originally posted by libralibra at 2013-04-26 23:43:33
没有数据,只能大概给个例子

T = 2.31;
t = 1:20;
x = randi(,1,length(t));
y1 = T*ones(size(x));
y2 = -y1;
plot(t,x,'k*-',t,y1,'r',t,y2,'b');
ylim();

...

我想问一下，这只是绘图的命令吧？我把我的程序发给你。麻烦你帮我看看，如果要画出我贴出的那个图，还要加什么命令。谢谢！！你太好了！检测变异是不是看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

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖