关于本站
人大经济论坛-经管之家：分享大学、考研、论文、会计、留学、数据、经济学、金融学、管理学、统计学、博弈论、统计年鉴、行业分析包括等相关资源。
经管之家是国内活跃的在线教育咨询平台!
获取电子版《CDA一级教材》
完整电子版已上线CDA网校，累计已有10万+在读~ 教材严格按考试大纲编写，适合CDA考生备考，也适合业务及数据分析岗位的从业者提升自我。
毕业论文
- 开题报告 | 【独家发布】论文 ...
- 开题报告 | 周五双学位论文开 ...
- 开题报告 | 还是找开题报告的 ...
- 开题报告 | 求浙江大学MBA论文 ...
- 开题报告 | 交开题报告
- 开题报告 | 本科毕业论文开题 ...
- 开题报告 | 开题报告、文献检 ...
- 开题报告 | 写开题报告中嘤嘤 ...
TOP热门关键词
连玉君stata课程：http://www.peixun.net/author/3.html连玉君优酷视频：http://i.youku.com/arlionStata现场培训-初级班【New,四天，北京】http://www.peixun.net/view/307_detail.htmlStata现场培训-高级班【New,四 ...
免费学术公开课,扫码加入 |
连玉君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&fromuid=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
「经管之家」APP：经管人学习、答疑、交友，就上经管之家！
免流量费下载资料----在经管之家app可以下载论坛上的所有资源，并且不额外收取下载高峰期的论坛币。
涵盖所有经管领域的优秀内容----覆盖经济、管理、金融投资、计量统计、数据分析、国贸、财会等专业的学习宝库，各类资料应有尽有。
来自五湖四海的经管达人----已经有上千万的经管人来到这里，你可以找到任何学科方向、有共同话题的朋友。
经管之家（原人大经济论坛），跨越高校的围墙，带你走进经管知识的新世界。
扫描下方二维码下载并注册APP
免流量费下载资料----在经管之家app可以下载论坛上的所有资源，并且不额外收取下载高峰期的论坛币。
涵盖所有经管领域的优秀内容----覆盖经济、管理、金融投资、计量统计、数据分析、国贸、财会等专业的学习宝库，各类资料应有尽有。
来自五湖四海的经管达人----已经有上千万的经管人来到这里，你可以找到任何学科方向、有共同话题的朋友。
经管之家（原人大经济论坛），跨越高校的围墙，带你走进经管知识的新世界。
扫描下方二维码下载并注册APP
您可能感兴趣的文章
人气文章
1.凡人大经济论坛-经管之家转载的文章,均出自其它媒体或其他官网介绍,目的在于传递更多的信息,并不代表本站赞同其观点和其真实性负责；
2.转载的文章仅代表原创作者观点,与本站无关。其原创性以及文中陈述文字和内容未经本站证实,本站对该文以及其中全部或者部分内容、文字的真实性、完整性、及时性，不作出任何保证或承若；
3.如本站转载稿涉及版权等问题,请作者及时联系本站,我们会及时处理。
2.转载的文章仅代表原创作者观点,与本站无关。其原创性以及文中陈述文字和内容未经本站证实,本站对该文以及其中全部或者部分内容、文字的真实性、完整性、及时性，不作出任何保证或承若；
3.如本站转载稿涉及版权等问题,请作者及时联系本站,我们会及时处理。
