请求计算机绘星系视星等、绝对星等图

huxn

依据下面截图条件绘图

happysxyf

前辈理论深奥， 我只能用自己发明的酷语言绘图如下：采用静波兰解析算式，每秒千万行级解析能力，绘制精度在0.000001左右。确实不是双曲线！


酷语言代码

1. 
   
2. /*
   
3. @ECHO OFF&CLS&COLOR 0F&MODE CON COLS=100 LINES=30&CLGO "%~F0"&PAUSE>NUL&EXIT /B
   
4. COPYRIGHT@2017~2019 BY HAPPY
   
5. */
   
6. 
   
7. TO FUNCTIONS_IMG :M :X :Y
   
8. 
   
9.         $PENS=2 $PENW=3 $PENC=RAND(15)+1
   
10.         FOR ::Z 0.1 3
    
11.         [
    
12.                 :X=(:M-(5*LN(1.22624*(10^10)*SINH(::Z/2))-40)*EXP(::Z))*EXP(0.016444)
    
13.                 :Y=3*:M*EXP(-::Z)
    
14.                 LINE  :X       :Y
    
15.                 $PENX=:X $PENY=:Y
    
16.                 EVAL  :X
    
17.                 EVAL  :Y
    
18.         ]0.01
    
19. 
    
20. TO MAIN
    
21. 
    
22.         FUNCTIONS_IMG 17.323 0 0
    
23.         FUNCTIONS_IMG 19.426 0 0
    
24.         FUNCTIONS_IMG 23.403 0 0
    
25.         FUNCTIONS_IMG 25.677 0 0
    

复制代码

——————————————————————————————————————
——————————————————————————————————————
——————————————————————————————————————
下面演示绘制其他曲线，绝对没有问题,绘制精度在0.000001左右，用的gdiplus库。
蝶形线：



酷语言代码

1. 
   
2. /*
   
3. @ECHO OFF&CLS&COLOR 0F&MODE CON COLS=100 LINES=30&CLGO "%~F0"&PAUSE>NUL&EXIT /B
   
4. COPYRIGHT@2017~2019 BY HAPPY
   
5. */
   
6. 
   
7. TO FUNCTIONS_IMG :A
   
8. 
   
9.         $PENS=2 $PENW=13 $PENC=4
   
10.         FOR ::T 0 3600
    
11.         [       
    
12.                 IF (::T)%360=0
    
13.                 [
    
14.                         $PENC=rand(15)+1
    
15.                 ]
    
16.                 $TP0=(exp(cos(deg(::T)))-2*cos(deg(4*::T))+(sin(deg(::T/12)))^5)*70
    
17.                 LINE  $TP0*cos(deg(::T))       $TP0*sin(deg(::T))
    
18.                 $PENX=$TP0*cos(deg(::T)) $PENY=$TP0*sin(deg(::T))
    
19.         ]
    
20. 
    
21. TO MAIN
    
22.        
    
23.         FUNCTIONS_IMG 0.0003
    
24.         PAUSE

复制代码

huxn

说明:绘天文物理图，一般两坐标轴不可能同一原点0，此图纵轴原点在起始点0，横轴起始点为30，纵轴可以放在右边，4条类似双曲线簇在第二象限內，其余象限不存在，是宇宙红移Z的定义域限制.m0(0为下标)标在纵轴图外边,M标在横轴下面.这个图m约取13至27时,有无穷多的曲线簇,现在仅绘4条曲线,已很好反映图的分布规律.单元格请柵去.

huxn

补充：请下载星系表II/326/zcatev ,表中仅点取红移redshifts zpho，视星等gmag，rmag。即方程中m=gmag或rmag。表中16<gmag，rmag<26，
0.1<zpho<5。表左边有小方框max，点击出现数字，取999,下载后表自有行数据，只下载约3千即可保存计算机内。以 zpho和m为方程(1)(2)中的对应红移、视星等。计算机绘出m0―M图,此图点集,是上面绘出曲线簇。是上面绘图高一级应用。网趾：http://vizier.u-strasbg.fr/viz-b ... e&-oc.form=sexa

huxn

广大读者:
我发出4个帖子已近20天了,不知什么原因,我国已毕业应用数学专业的青年才俊很多.我出的天体物浬学的前沿基础,是应用数学专业的青年大展身手的内容,给出的定义域符合天文观测数据,且都是光滑函数,还看不到解决问题的好囬复.我的4个问题本来是天体物浬学教材中最基础的问题.
谁人绘得好图,我写论文要用你们绘的图及绘图作者名及网站.我可以绘图,可惜年岁限制,绘图不好,无法达到发表论文标准,才发帖求助,请青年才俊发挥你们才华!