log图绘图学习

%%数据载入

load starData

nObs = size(spectra,1)

lambdaStart = 630.02

lambdaDelta = 0.14

%%行向量

lambdaEnd = lambdaStart + (nObs-1)*lambdaDelta

lambda = (lambdaStart:lambdaDelta:lambdaEnd)' %%“ ' ”将向量、矩阵转置

%%提取数据

s = spectra(:,6)

log函数图

loglog(lambda,s,".-") %%x,y,点标记实现连接

xlabel("Wavelength") %%x轴标签

ylabel("Intensity")%%y轴标签

索引

[sHa,idx]=min(s)%%提取最小值sHa和对应位置idx

lambdaHa=lambda(idx)%%提取对应位置的数据

标注点

hold on%%原图上

loglog(lambdaHa,sHa,"rs","MarkerSize",8)%%x，y位置，红色方框，线条尺寸

hold off

循环

for c=(1:7)%%语句循环，从1到7，步进为1，

s=spectra(:,c)%%提取spectra的列向量

[sHa,idx]=min(s)%%y值的最小值及对应序号

lambdaHa=lambda(idx)%%对应序号对应的x值

z=lambdaHa/656.28-1

speed(c)=z*299792.458

if speed(c)<=0%%条件语句

loglog(lambda,s,"--")%%log图，虚线

else speed>0%%判断

loglog(lambda,s,"Linewidth",3)%%画log图,线宽3

end

hold on%%画在同一张图上

end

hold off%%解除

legend(starnames)%%添加图例

movaway=starnames(speed>0)%%从图中提取speed>0的星球名称