import numpy as np
import matplotlib.pyplot as plt
import math
from scipy.interpolate import make_interp_spline #平滑曲线用的函数
plt.axis([0, 130,-16, 16])
xss=np.array([0,5,10,15,20,25,30,35,40,45,50,55,60,65,70,75,80,85,90,95,100,105,110,115,120,125,130])
ys=np.array([-1.56699, -1.55284, -1.51299, -1.44526, -1.34541, -1.20726, -1.02088, -0.775333, -0.463055,
-0.0576436, 0.475682, 1.15845, 2.00164, 2.99478, 4.13066, 5.38871, 6.73956, 8.10639, 9.43044,
10.6792, 11.8363, 12.8766, 13.7775, 14.5165, 15.0723, 15.4162, 15.5427])
plt.ion()
for k in range(20):
for i in range(101):
A=math.sin(i*math.pi/50)
yss=A*ys
x_smooth = np.linspace(xss.min(), xss.max(), 300)
y_smooth = make_interp_spline(xss, yss)(x_smooth)
plt.cla()
plt.plot(x_smooth, y_smooth,linewidth = '2', color= '#FF0000')
plt.axis([0, 130,-16, 16])
plt.show()
plt.pause(0.005)
if i==0 and k==0:input()

清明上河图》画卷简介
《清明上河图》画卷，北宋风俗画作品。作者张择端，字正道，东武(今山东诸城)人。宋徽宗时为宫廷画家。少年时到京城汴梁(今河南开封)游学，后习绘画，尤喜画舟车、市桥、郭径，自成一家。《清明上河图》是他的代表作，曾经为宣和内府所收藏。绢本、淡设色，纵24．8厘米，横528．7厘米，现藏北京故宫博物院。该图描绘了清明时节，北宋京城汴梁以及汴河两岸的繁华景象和自然风光。作品以长卷形式，采用散点透视的构图法，将繁杂的景物纳入统一而富于变化的画面中，画中人物500多，衣着不同，神情各异，其间穿插各种活动，注重戏剧性，构图疏密有致，注重节奏感和韵律的变化，笔墨章法都很巧妙。全图分为三个段落。
首段，汴京郊野的春光：在疏林薄雾中，掩映着几家茅舍、草桥、流水、老树、扁舟。两个脚夫赶着五匹驮炭的毛驴，向城市走来。一片柳林，枝头刚刚泛出嫩绿，使人感到虽是春寒料峭，却已大地回春。路上一顶轿子，内坐一位妇人。轿顶装饰着杨柳杂花，轿后跟随着骑马的、挑担的，从京郊踏青扫墓归来。环境和人物的描写，点出了清明时节的特定时间和风俗，为全画展开了序幕。
中段，繁忙的汴河码头：汴河是北宋国家漕运枢纽，商业交通要道，从画面上可以看到人烟稠密，粮船云集，人们有在茶馆休息的，有在看相算命的，有在饭铺进餐的。还有“王家纸马店”，是扫墓卖祭品的，河里船只往来，首尾相接，或纤夫牵拉，或船夫摇橹，有的满载货物，逆流而上，有的靠岸停泊，正紧张地卸货。横跨汴河上的是一座规模宏大的木质拱桥，它结构精巧，形式优美。宛如飞虹，故名虹桥。有一只大船正待过桥。船夫们有用竹竿撑的；有用长竿钩住桥梁的；有用麻绳挽住船的；还有几人忙着放下桅杆，以便船只通过。邻船的人也在指指点点地象在大声吆喝着什么。船里船外都在为此船过桥而忙碌着。桥上的人，也伸头探脑地在为过船的紧张情景捏了一把汗。这里是名闻遐迩的虹桥码头区，车水马龙，熙熙攘攘，名副其实地是一个水陆交通的会合点。
后段，热闹的市区街道：以高大的城楼为中心，两边的屋宇鳞次栉比，有茶坊、酒肆、脚店、肉铺、庙宇、公廨等等。商店中有绫罗绸缎、珠宝香料、香火纸马等的专门经营，此外尚有医药门诊，大车修理、看相算命、修面整容，各行各业，应有尽有，大的商店门首还扎着“彩楼欢门”，悬挂市招旗帜，招揽生意，街市行人，摩肩接踵，川流不息，有做生意的商贾，有看街景的士绅，有骑马的官吏，有叫卖的小贩，有乘座轿子的大家眷属，有身负背篓的行脚僧人，有问路的外乡游客，有听说书的街巷小儿，有酒楼中狂饮的豪门子弟，有城边行乞的残疾老人，男女老幼，士农工商，三教九流，无所不备。交通运载工具：有轿子、骆驼、牛马车、人力车，有太平车、平头车，形形色色，样样俱全。绘色绘形地展现在人们的眼前。总计在五米多长的画卷里，共绘了五百五十多个各色人物，牛、马、骡、驴等牲畜五、六十匹，车、桥二十多辆，大小船只二十多艘。房屋、桥梁、城楼等也各有特色，体现了宋代建筑的特征。张择端的《清明上河图》是一幅描写北宋汴京城一角的现实主义的风俗画，具有很高的历史价值和艺术水平。

package com.zuidaima.swing;
import java.awt.Color;
import java.awt.Graphics;
import java.awt.Image;
import java.awt.Toolkit;
import javax.swing.JFrame;
@SuppressWarnings("serial")
class Cardioid extends JFrame {
// 定义窗口大小
private static final int WIDTH = 480;
private static final int HEIGHT = 600;
// 获取屏幕大小
private static int WINDOW_WIDTH = Toolkit.getDefaultToolkit()
.getScreenSize().width;
private static int WINDOW_HEIGHT = Toolkit.getDefaultToolkit()
.getScreenSize().height;
// 构造函数
public Cardioid() {
// 设置窗口标题
super("❤形线");
// 设置背景色
this.setBackground(Color.BLACK);
// 设置窗口位置
this.setLocation((WINDOW_WIDTH - WIDTH) / 2,
(WINDOW_HEIGHT - HEIGHT) / 2);
// 设置窗口大小
this.setSize(WIDTH, HEIGHT);
// 设置窗口布局
this.setLayout(getLayout());
// 设置窗口可见
this.setVisible(true);
// 设置窗口默认关闭方式
this.setDefaultCloseOperation(DISPOSE_ON_CLOSE);
}
public void paint(Graphics g) {
double x, y, r; // 横纵坐标及半径
Image image = this.createImage(WIDTH, HEIGHT);
Graphics pic = image.getGraphics();
// 绘制图形
for (int i = 0; i < 100; i++) {
for (int j = 0; j < 100; j++) {
r = Math.PI / 45 + Math.PI / 45 * i
* (1 - Math.sin(Math.PI / 45 * j)) * 18;
x = r * Math.cos(Math.PI / 45 * j) * Math.sin(Math.PI / 45 * i)
+ WIDTH / 2;
y = -r * Math.sin(Math.PI / 45 * j) + HEIGHT / 2;
pic.setColor(Color.MAGENTA);
pic.fillOval((int) x, (int) y, 2, 2);
}
// 生成图片
g.drawImage(image, 0, 0, this);
}
}
public static void main(String[] args) {
new Cardioid();
}
}
转自:http://www.oschina.net/code/snippet_1591393_36195

原始方程源于此贴一楼：直通车
整理修改后：
被窝修改成这样：x
yx[p,i]
xx[,i]
par(bg="pink")
plot(main="中国(标题可更改)", c(rev(xx),xx),c(rev(y),yx),type='l',xlim = c(-1,1),xlab = "LOVE YOU",col.lab="red",font.lab=9,cex.lab=3,ylab = "",axes=F,col=q)
Sys.sleep(1)} #Sys.sleep(1)动态效果
据说笛卡尔死前寄出的最后一封信，里面只有短短的一行：r=a(1−sinθ)'>r=a(1−sinθ) r=a(1−sinθ) ，这就是有名的心形函数。情人节将至，我用R语言的grid包画了几幅图片，希望借此平台赠与我相恋五年的男友，也希望与各位统计爱好者分享快乐。
首先，我利用grid.lines()把转化为直角坐标系的(x, y)两两相连围成心形，构建了一个heart function作为基本图形。x=16(sint)3y=13cost–5cos2t–2cos3t–cos4t'>xy =16(sint) 3 =13cost–5cos2t–2cos3t–cos4t x=16(sint)3y=13cost–5cos2t–2cos3t–cos4t
为了得到嵌套心形图案，我使用了grid包创建了多个viewport。viewport是grid包的一个重要特色，此概念类似于photoshop的图层。
创建一个viewport，我们需要设置它的位置、长度和宽度，下图虚线实际上并不出现在R的output里面，但这个矩形区域图层会成为接下来画图的区域。构建了新的viewport以后，我们可以用pushViewport()命令锁定该图层，使之成为目标区域。我们也可以构建多个viewport，几个viewport之间可以通过命令相互切换。
例如，在第一个图层的基础上在新建一个图层，调整新图层的长度与宽度使之稍微小于第一个图层， 用pushViewport()锁定新图层，再调用一次心形函数，以此类推，循环创建多个嵌套图层，并依次在各个图层上画心形函数，于是我们可以得到一系列嵌套的心形。
此外，grid包允许我们对图形进行复制、旋转、放缩等修改。要旋转心形函数，我们并不需要修改函数本身，而是可以通过旋转viewport的方式旋转我们所需要绘制的图形。设置新viewport，调整angle函数，那么在此图层下绘制的任何图形将会被旋转。
利用viewport对图形进行修改，我们可以绘制各种有趣的图形pattern，本人只是grid包的初学者，如有偏颇之处望多多包涵。最后，祝愿大家情人节快乐！
以下是“情人节礼物”的代码：直通车