- 积分
- 1828
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-11-19
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
clc
clear
s1=1:1:18000;
s2=1:1:300;
ss2=7680:1:7979;
ss1=460800:1:478799;
data=load('g:\03013\1002\ud.txt');
data1=load('g:\03013\1002\oneminute\0201u.txt');
u=data(:,7);
uu=data1(:,1);
%subplot(211)
%a1=mean(uu(ss2));
%a2=max(uu(ss2));
%a3=max(u(ss1));
%a4=min(uu(ss2));
%scatter(s1,u(ss1),'.k')
%set(gca,'xtick',[1 3600 7200 10800 14400 18000],'xticklabel',{'08h','09h','10h','11h','12h','13h'})
%title('2.6秒量级短期爬坡样本')
%xlabel('时间(h)')
%ylabel('风速(m/s)')
%subplot(212)
%scatter(s2,uu(ss2),'.k')
v=polyfit(s2,uu(ss2),1);
t=1:0.5:300;
w=polyval(t,w);
plot(t,w,s2,uu(ss2),'.k')
set(gca,'xtick',[1 60 120 180 240 300],'xticklabel',{'08h','09h','10h','11h','12h','13h'})
title('2.6分钟量级短期爬坡样本')
xlabel('时间(h)')
ylabel('风速(m/s)') |
|