全部版块 我的主页
› 论坛 › 计量经济学与统计论坛 五区 › 计量经济学与统计软件 › Stata专版
用 Stata 画蝴蝶和花瓣
arlionn
62331 951
2013-06-15

    连玉君stata课程：http://www.peixun.net/author/3.html

     连玉君优酷视频： http://i.youku.com/arlion

     Stata现场培训-初级班【New, 四天，北京】http://www.peixun.net/view/307_detail.html

     Stata现场培训-高级班【New, 四天，北京】http://www.peixun.net/view/308_detail.html


guo.bailing 用 SAS 画了一只漂亮的蝴蝶
（https://bbs.pinggu.org/forum.php?mod=viewthread&tid=2137037&page=1&from^^uid=40938）。
按其思路，我画了一只 Stata 蝴蝶，程序略改后，又画了两朵花。

*--------------------------------- Stata 花蝴蝶-------------------------------
*-蝴蝶身体部分
  clear
  local b0 = -400
  local b1 = 4
  local d  = 0.01
  local N = int((`b1'-`b0')/`d')
  set obs `N'
  gen t = -400 + _n*`d'
  gen a = cos(t)
  gen b = cos(4*t)
  gen x = sin(t)*(exp(a)-2*b-(sin(t/12)*sin(t/12)*sin(t/12)*sin(t/12)*sin(t/12)))
  gen y = cos(t)*(exp(a)-2*b-(sin(t/12)*sin(t/12)*sin(t/12)*sin(t/12)*sin(t/12)))
           
*-胡须
  dropvars x2 y2 x3
  gen y2 =.
  gen x2 =.        
  local j = 1
  local end = (_pi/4+0.20)
  forvalues i = 0.20(0.01)`end'{
    qui replace x2 = `i' in `j++'
  }
  replace y2 = 3*sin(2*(x2-0.20))+2.5
  gen x3 = -x2  //左半只

*-绘图
  twoway (line y x if y>0, lc(pink*0.55) lw(*2.5))    ///
         (line y x if y<0, lc(pink*0.70) lw(*2.5))    ///
         (line y2 x2     , lc(yellow*1.2) lw(*3.5))   ///
         (line y2 x3     , lc(yellow*1.2) lw(*3.5)),  ///
         yscale(off) xscale(off) legend(off)
*-----------------------------------------------------------------------------
butterfly.JPG

*--------------------用 Stata 绘制花儿: 忽悠小朋友用------------------------------------
  clear
  local b0 = -400
  local b1 = 4
  local d  = 0.005  //0.01
  local N = int((`b1'-`b0')/`d')
  set obs `N'
  gen t = -400 + _n*`d'
  gen a = cos(t)
  
*-主要参数  
  local i = 30  // 4   花瓣的数量
  local j = 24  // 12  花瓣叠入程度
  local s = 20  // 2   内瓣和外瓣比例
  gen b = cos(`i'*t)
  gen x = sin(t)*(exp(a)-`s'*b-(sin(t/`j')*sin(t/`j')*sin(t/`j')*sin(t/`j')*sin(t/`j')))
  gen y = cos(t)*(exp(a)-`s'*b-(sin(t/`j')*sin(t/`j')*sin(t/`j')*sin(t/`j')*sin(t/`j')))
   
*-绘图
  set scheme s2color
  local c "pink"
  local c "red*0.3"
  local z "0.8"     //线宽
  twoway (line y x if y>0, lc(`c') lw(*`z')) ///
         (line y x if y<0, lc(`c') lw(*`z')) ///
         , yscale(off) xscale(off) legend(off)      
*------------------------------------------------------------------------------------------
flower2.JPG


    连玉君stata课程：http://www.peixun.net/author/3.html

     连玉君优酷视频： http://i.youku.com/arlion

     Stata现场培训-初级班【New, 四天，北京】http://www.peixun.net/view/307_detail.html

     Stata现场培训-高级班【New, 四天，北京】http://www.peixun.net/view/308_detail.html
二维码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

全部回复
voodoo
2013-6-15 22:10:35
蝴蝶Piao<<liang！小花差点。












二维码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

zzl19910906
2013-6-15 23:30:43
收藏了
二维码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

改革同步
2013-6-16 18:09:11
学习学习
二维码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

weiha
2013-6-17 00:50:12
牛逼啊
二维码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

angelatao
2013-6-17 10:40:32
太美了    ！
二维码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

微信QQ空间QQ微博
扫码加好友，拉您进群
各岗位、行业、专业交流群