|
楼主 |
发表于 2010-4-22 11:35
|
显示全部楼层
总算是写出来了
我也发上来 希望有空的人 帮我修改修改
%导入数据
x1=1:1:100;
y1=sqrt(x1);
%相关系数
n=length(x1);
c=cell(n-1,1);
x=c;
y=x;
for ii=1:n-1;
x{ii} = x1(1:ii+1);
y{ii} = y1(1:ii+1);
c{ii} = corrcoef(x{ii},y{ii});
end
for ii=1:n-2;
b(ii)=c{ii+1}(1,2);
end
bb=[1,1,b];
%斜率
d=y1./x1;
%绘制原图 相关系数 斜率 对比
figure
subplot(3,1,1)
plot(x1,y1)
title('原图')
grid on
subplot(3,1,2)
plot(x1,bb)
title('相关系数')
grid on
subplot(3,1,3)
plot(x1,d)
title('斜率')
grid on
%在相关系数中找开始小于a的点
a=0.99;
for iii=1:n;
flag=(bb(iii)<=a);
if flag;
break;
end
end
num=iii-1;
x2=x1(1:num); %整理后的x
y2=y1(1:num); %整理后的y |
|