27694021
2007-07-18, 15:35
load sunspot.dat %打开sunspot文件
year=sunspot(:,1): %第一列为year
wolfer=sunspot(:,2); %第二列为wolfer
plot(year,wolfer); %以year为横轴,以wolfer为纵轴画图
title(‘sunspot Data’) %标题为sunspot Data
Y=fft(wolfer);傅立叶变换
N=length(Y); 计算Y的长度
Y(1)=[];
power=abs(Y(1:N/2).^2);
nyquist=1/2;
freq=(1:N/2)/(N/2)*nyquist;
plot(freq,power),grid on 画图 产生网格线
xlable(‘年’); X轴为年
title(‘能量图象’) 标题为能量图象
priod=1./freq ;
plot(period,power); 以period,power为 XY轴画图
axis([0 40 0 2e7])
grid on 产生网格线
ylabel(‘能量’) Y轴标为能量
xlabel(‘年’) X轴标为年
还没翻译完全，好心人补下！
year=sunspot(:,1): %第一列为year
wolfer=sunspot(:,2); %第二列为wolfer
plot(year,wolfer); %以year为横轴,以wolfer为纵轴画图
title(‘sunspot Data’) %标题为sunspot Data
Y=fft(wolfer);傅立叶变换
N=length(Y); 计算Y的长度
Y(1)=[];
power=abs(Y(1:N/2).^2);
nyquist=1/2;
freq=(1:N/2)/(N/2)*nyquist;
plot(freq,power),grid on 画图 产生网格线
xlable(‘年’); X轴为年
title(‘能量图象’) 标题为能量图象
priod=1./freq ;
plot(period,power); 以period,power为 XY轴画图
axis([0 40 0 2e7])
grid on 产生网格线
ylabel(‘能量’) Y轴标为能量
xlabel(‘年’) X轴标为年
还没翻译完全，好心人补下！