用 Stata 画蝴蝶和花瓣

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

---

    连玉君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)
*-----------------------------------------------------------------------------


*--------------------用 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)      
*------------------------------------------------------------------------------------------

---

    连玉君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

---

扫码加我 拉你入群

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

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

分享2 收藏16 回帖

关键词：Stata tata forvalues forvalue replace 蝴蝶 暑期培训 fromuid 花蝴蝶 local

蝴蝶Piao<<liang！小花差点。

收藏了

学习学习

牛逼啊

太美了    ！

精致作品！

精品，太漂亮了~

漂亮

连老师威武