update

C/001/main.cpp

-//例20：C语言实现打印出心形，初学者的表白神器
-
-#include<stdio.h>//头文件 
-int main()//主函数入口 
-{
-	printf("      ****        ****\n");//打印第一行 
-	printf("   *********    *********\n");//打印第二行 
-	printf("*************  *************\n");//打印第三行 
-	int i, j;//定义变量 
-	for (i = 0; i < 3; i++)//打印4-6行，一共3行，因此i小于3 
-	{
-		for (j = 0; j < 29; j++)//限制每行输出*的个数 
-		{
-			printf("*");//这三行只打印*号，无空格输出 
-		}
-		printf("\n");//打印完一行需要进行换行 
-	}
-	for (i = 0; i < 7; i++) //打印7-13行，一共7行，因此i小于7 
-	{
-		for (j = 0; j < 2 * (i + 1) - 1; j++)//这个for循环和下面的for是并列的 
-		{
-			printf(" ");//打印空格 
-		}
-		for (j = 0; j < 27 - i * 4; j++)//读者可以带入几个数找出条件 
-		{
-			printf("*");//打印*
-		}
-		printf("\n");
-	}
-	for (i = 0; i < 14; i++)//打印最后一行的*
-	{
-		printf(" ");//打印空格 
-	}
-	printf("*\n");//打印* 
-	return 0;
-}
-

C/002/main.cpp

-#include<stdio.h>
-#include<math.h>
-int main()
-{
-	float y, x, z;
-
-	printf("那一天\n");
-	printf("第一次遇见你\n");
-	printf("忘不了\n");
-	printf("你的容颜\n");
-	printf("若轻云之蔽月,如流风之回雪\n");
-	printf("\n\n\n");
-	printf("其实\n");
-	printf("有一句话\n");
-	printf("我一直想对你说:\n");
-
-	for (double y = 2.5; y >= -1.6; y = y - 0.2)
-	{
-		for (double x = -3; x <= 4.8; x = x + 0.1)
-		{
-			(pow((x * x + y * y - 1), 3) <= 3.6 * x * x * y * y * y
-				|| (x > -2.4 && x < -2.1 && y<1.5 && y>-1)
-				|| (((x < 2.5 && x>2.2) || (x > 3.4 && x < 3.7)) && y > -1 && y < 1.5)
-				|| (y > -1 && y < -0.6 && x < 3.7 && x>2.2)) ? printf("*") : printf(" ");
-		}
-
-		printf("\n");
-	}
-
-	getchar();
-}

C/003/main.cpp

-#include <stdio.h>
-#include <math.h>
-#include <windows.h>
-#include <tchar.h>
-
-float f(float x, float y, float z) {
-    float a = x * x + 9.0f / 4.0f * y * y + z * z - 1;
-    return a * a * a - x * x * z * z * z - 9.0f / 80.0f * y * y * z * z * z;
-}
-
-float h(float x, float z) {
-    for (float y = 1.0f; y >= 0.0f; y -= 0.001f)
-        if (f(x, y, z) <= 0.0f)
-            return y;
-    return 0.0f;
-}
-
-int main() {
-    HANDLE o = GetStdHandle(STD_OUTPUT_HANDLE);
-    _TCHAR buffer[25][80] = { _T(' ') };
-    _TCHAR ramp[] = _T(".:-=+*#%@");
-
-    for (float t = 0.0f;; t += 0.1f) {
-        int sy = 0;
-        float s = sinf(t);
-        float a = s * s * s * s * 0.2f;
-        for (float z = 1.3f; z > -1.2f; z -= 0.1f) {
-            _TCHAR* p = &buffer[sy++][0];
-            float tz = z * (1.2f - a);
-            for (float x = -1.5f; x < 1.5f; x += 0.05f) {
-                float tx = x * (1.2f + a);
-                float v = f(tx, 0.0f, tz);
-                if (v <= 0.0f) {
-                    float y0 = h(tx, tz);
-                    float ny = 0.01f;
-                    float nx = h(tx + ny, tz) - y0;
-                    float nz = h(tx, tz + ny) - y0;
-                    float nd = 1.0f / sqrtf(nx * nx + ny * ny + nz * nz);
-                    float d = (nx + ny - nz) * nd * 0.5f + 0.5f;
-                    *p++ = ramp[(int)(d * 5.0f)];
-                }
-                else
-                    *p++ = ' ';
-            }
-        }
-
-        for (sy = 0; sy < 25; sy++) {
-            COORD coord = { 0, sy };
-            SetConsoleCursorPosition(o, coord);
-            WriteConsole(o, buffer[sy], 79, NULL, 0);
-        }
-        Sleep(33);
-    }
-}
\ No newline at end of file

C/004/Readme.md

-## 考虑到仓库大小，此程序单独存放，请前往
-# https://github.com/sun0225SUN/C-Love-Code
\ No newline at end of file

C/005/Readme.md

-## 考虑到仓库大小，此程序单独存放，请前往
-# https://github.com/sun0225SUN/C-Meteor-Shower
\ No newline at end of file

C/006/Readme.md

-## 考虑到仓库大小，此程序单独存放，请前往
-# https://github.com/sun0225SUN/C-Fireworks-Code
\ No newline at end of file

Csharp/001/Readme.md

-## 考虑到仓库大小，此程序单独存放，请前往
-# https://github.com/sun0225SUN/Be-My-Girlfriend
\ No newline at end of file

MatLab/001/just_heart.m

-t=[0:0.1:2*pi]
-x=16 * sin(t).^3
-y=13 * cos(t) - 5 * cos(2 * t) - 2 * cos(3 * t) - cos(4 * t)
-a=(x - min(x))/(max(x) - min(x))
-b=(y - min(y))/(max(y) - min(y))
-plot(x,y,'r')
\ No newline at end of file

MatLab/001/readme.md

-<img src="https://cdn.jsdelivr.net/gh/sun0225SUN/Awesome-Love-Code/assets/img/matlab/001.jpg"/>

MatLab/002/heart.m

-% matlab动态心形曲线（保存GIF格式）
-clc;clear all;
-a=10;
-x=-2:0.01:2;
-
-figure(1);
-set(gcf,'position',[0,0,800,600],'color','w');
-for i=1:100
-    str_title = strcat('\color{red}a=',num2str(a));
-    y = abs(x).^(2/3)+(0.9*sqrt((3.3-x.^2))).*sin(a*pi*x);
-    a=i/10;
-    figure(i);
-    
-    set(gcf,'position',[0,0,800,600],'color','w');
-    p1 = plot(x,y,'r','LineWidth',3);
-    hold on;
-    title('\color{red}f(x)=x^2^/^3+e/3*(π-x^2)^1^/^2*sin(a*π*x)','fontsize',15);
-    text(-0.2,2.3,str_title,'FontName','Times New Roman','FontSize',20);
-    xlim([-2 2]);
-    ylim([-1.5 2.5]);
-    frame = getframe(gcf);
-    im = frame2im(frame);
-    [I,map] = rgb2ind(im,256);
-    if i==1
-        imwrite(I,map,'心形图.gif','gif','Loopcount',Inf,'DelayTime',0.05);
-    else
-        imwrite(I,map,'心形图.gif','gif','WriteMode','append','DelayTime',0.05);
-    end
-    close(figure(i));
-end
-

MatLab/002/readme.md

-<img src="https://cdn.jsdelivr.net/gh/sun0225SUN/Awesome-Love-Code/assets/img/matlab/002.gif"/>

MatLab/003/heart_3d.m

-f=@(x,y,z)(x.^2+ (9./4).*y.^2 + z.^2 - 1).^3 - x.^2.*z.^3 - (9./80).*y.^2.*z.^3;
-[x,y,z]=meshgrid(linspace(-3,3));
-val=f(x,y,z);
-[p,v]=isosurface(x,y,z,val,0);
-patch('faces',p,'vertices',v,'facevertexcdata',jet(size(v,1)),'facecolor','w','edgecolor','flat');
-view(3);grid on;axis equal;
\ No newline at end of file

MatLab/003/readme.md

-<img src="https://cdn.jsdelivr.net/gh/sun0225SUN/Awesome-Love-Code/assets/img/matlab/003.jpg"/>

MatLab/004/my_heart.m

-[X,Y,Z] = meshgrid(-1.2:.02:1.2,-3:.03:3,-1:.02:1.3);
-f = smooth3((X.^2+9*Y.^2/9+Z.^2-1).^3-X.^2.*Z.^3-Y.^2.*Z.^3/20,'gauss');
-v = uniquetol(getfield(isosurface(X,Y,Z,f,0),'vertices'),.15,'byrows',1);
-n = isonormals(X,Y,Z,f,v);
-ar = bsxfun(@rdivide,n*[0 -1 0;1 0 0;0 0 0],sqrt(sum(n.^2,2)));
-
-[x,t] = meshgrid(0:.125:1,(-2:.2:10)*pi);
-[p,q] = deal(pi/2*exp(-t/(8*pi)),1-(1-mod(3.6*t,2*pi)/pi).^4/2);
-y = 2*(x.^2-x).^2.*sin(p);
-[p,q] = deal(q.*(x.*sin(p)+y.*cos(p)),q.*(x.*cos(p)-y.*sin(p)));
-
-figure color w, axis image vis3d off
-h = surface(p.*cos(t),p.*sin(t),q,'EdgeColor','n','FaceColor','r');
-arrayfun(@(t,i)copyobj(h,hgtransform('Mat',makehgtform('translate',5*v(i,:),...
-   'axisrotate',ar(i,:),t))),asin(sqrt(sum(ar.^2,2)))',1:size(v,1))
-view(32,12), set(camlight('head'),'color',[1 0.84 0.6]), lighting gouraud
\ No newline at end of file

MatLab/004/readme.md

-<img src="https://cdn.jsdelivr.net/gh/sun0225SUN/Awesome-Love-Code/assets/img/matlab/004.jpg"/>

MatLab/005/readme.md

-<img src="https://cdn.jsdelivr.net/gh/sun0225SUN/Awesome-Love-Code/assets/img/matlab/005.jpg"/>

MatLab/005/red_heart.m

-f=@(x,y,z)(x.^2+ (9./4).*y.^2 + z.^2 - 1).^3 - x.^2.*z.^3 - (9./80).*y.^2.*z.^3;
-[x,y,z]=meshgrid(linspace(-1.5,1.5));
-val=f(x,y,z);
-isosurface(x,y,z,val,0); 
-axis equal;view(3);colormap([1 0.2 0.2])

MatLab/006/love_heart.m

-%跳动爱心
-clear; clc; close all; 
-% NOTICE: Your MATLAB version should be at least R2019a !
-
-% --------------------------- functions -------------------------------
-f = @(x, y, z)(x.^2 + 2.25*y.^2 + z.^2 - 1).^3 -  ...
-    x.^2.* z.^3 - 0.1125*y.^2.*z.^3;
-g = @(x, y, z)(sqrt(x.^2+y.^2)-2.5).^2 + z.^2 - 0.4^2;
-
-% -------------------------- generate data ----------------------------
-t = linspace(-5, 5);
-
-[x1, y1, z1] = meshgrid(t);
-[x2, y2, z2] = meshgrid(t);
-
-val1 = f(x1, y1, z1);
-val2 = g(x2, y2, z2);
-
-[p1, v1] = isosurface(x1, y1, z1, val1, 0);
-[p2, v2] = isosurface(x2, y2, z2, val2, 0);
-
-% --------------------------- basic plot ------------------------------
-figure()
-subplot(1, 1, 1)
-
-h = patch('faces',p1,'vertices',v1,'facevertexcdata',jet(size(v1,1)),...
-    'facecolor','w','edgecolor','flat'); hold on;
-patch('faces',p2,'vertices',v2,'facevertexcdata',jet(size(v2,1)),...
-    'facecolor','w','edgecolor','flat');
-grid on; axis equal; axis([-3,3,-3,3,-1.5,1.5]); view(3)
-title(["$(x^2+\frac{9}{4}y^2+z^2-1)^3-x^2z^3-\frac{9}{80}y^2z^3=0$",...
-    "$(\sqrt{x^2+y^2}-R)^2 +z^2 = r^2$"],'Interpreter','latex','position',[3.3,4])
-warning('off');
-% 请在此处进行您的演讲！% 请在此处进行您的演讲！
-T = suptitle("$I\ Love\ U\ !$"); 
-% 请在此处进行您的演讲！% 请在此处进行您的演讲！
-set(T,'Interpreter','latex','FontSize',24)
-
-% -------------------------- generate gif -----------------------------
-pic_num = 1;
-for i = 1:20
-    v1 = 0.98 * v1;
-    set(h, 'vertices', v1); drawnow;
-    F = getframe(gcf);
-    I = frame2im(F);
-    [I,map]=rgb2ind(I,256);
-    if pic_num == 1
-        imwrite(I,map,'BeatingHeart.gif','gif','Loopcount',inf,'DelayTime',0.05);
-    else
-        imwrite(I,map,'BeatingHeart.gif','gif','WriteMode','append','DelayTime',0.05);
-    end
-    pic_num = pic_num + 1;
-end
-for i = 1:20
-    v1 = v1 / 0.98;
-    set(h, 'vertices', v1); drawnow;
-    F = getframe(gcf);
-    I = frame2im(F);
-    [I,map] = rgb2ind(I,256);
-    imwrite(I,map,'BeatingHeart.gif','gif','WriteMode','append','DelayTime',0.05);
-    pic_num = pic_num + 1;
-end
-
-% ----------------------------- dynamic -------------------------------
-while true
-    for i = 1:20
-        v1 = 0.98 * v1;
-        set(h, 'vertices', v1); drawnow;
-    end
-    for i = 1:20
-        v1 = v1 / 0.98;
-        set(h, 'vertices', v1); drawnow;
-    end
-end
-% ------------------------------- end ---------------------------------

MatLab/006/readme.md

-<img src="https://cdn.jsdelivr.net/gh/sun0225SUN/Awesome-Love-Code/assets/img/matlab/006.jpg"/>

MatLab/007/readme.md

-<img src="https://cdn.jsdelivr.net/gh/sun0225SUN/Awesome-Love-Code/assets/img/matlab/007.jpg"/>

MatLab/007/rose.m

-function rose
-grid on
-[x,t]=meshgrid((0:24)./24,(0:0.5:575)./575.*20.*pi-4*pi);
-p=(pi/2)*exp(-t./(8*pi));
-change=sin(20*t)/150;
-u=1-(1-mod(3.3*t,2*pi)./pi).^4./2+change;
-y=2*(x.^2-x).^2.*sin(p);
-
-r=u.*(x.*sin(p)+y.*cos(p)).*1.5;
-h=u.*(x.*cos(p)-y.*sin(p));
-
-map=[1.0000    0.6471    0.8275
-    0.9984    0.6353    0.8130
-    0.9969    0.6236    0.7985
-    0.9953    0.6118    0.7840
-    0.9937    0.6000    0.7695
-    0.9921    0.5882    0.7550
-    0.9906    0.5765    0.7404
-    0.9890    0.5647    0.7259
-    0.9874    0.5529    0.7114
-    0.9859    0.5412    0.6969
-    0.9843    0.5294    0.6824
-    0.9757    0.5149    0.6730
-    0.9670    0.5004    0.6636
-    0.9584    0.4859    0.6541
-    0.9498    0.4714    0.6447
-    0.9411    0.4568    0.6353
-    0.9325    0.4423    0.6259
-    0.9239    0.4278    0.6165
-    0.9153    0.4133    0.6070
-    0.9066    0.3988    0.5976
-    0.8980    0.3843    0.5882
-    0.8937    0.3780    0.5756
-    0.8894    0.3718    0.5631
-    0.8851    0.3655    0.5505
-    0.8808    0.3592    0.5380
-    0.8764    0.3529    0.5254
-    0.8721    0.3467    0.5129
-    0.8678    0.3404    0.5003
-    0.8635    0.3341    0.4878
-    0.8592    0.3279    0.4752
-    0.8549    0.3216    0.4627
-    0.8561    0.3165    0.4596
-    0.8573    0.3114    0.4564
-    0.8584    0.3063    0.4533
-    0.8596    0.3012    0.4502
-    0.8608    0.2961    0.4471
-    0.8620    0.2910    0.4439
-    0.8632    0.2859    0.4408
-    0.8643    0.2808    0.4377
-    0.8655    0.2757    0.4345
-    0.8667    0.2706    0.4314
-    0.8549    0.2620    0.4165
-    0.8432    0.2533    0.4016
-    0.8314    0.2447    0.3867
-    0.8196    0.2361    0.3718
-    0.8078    0.2274    0.3569
-    0.7961    0.2188    0.3420
-    0.7843    0.2102    0.3271
-    0.7725    0.2016    0.3122
-    0.7608    0.1929    0.2973
-    0.7490    0.1843    0.2824
-    0.7553    0.1827    0.2855
-    0.7616    0.1812    0.2887
-    0.7678    0.1796    0.2918
-    0.7741    0.1780    0.2949
-    0.7804    0.1764    0.2980
-    0.7867    0.1749    0.3012
-    0.7930    0.1733    0.3043
-    0.7992    0.1717    0.3074
-    0.8055    0.1702    0.3106
-    0.8118    0.1686    0.3137
-    0.7977    0.1631    0.3023
-    0.7836    0.1576    0.2910
-    0.7694    0.1521    0.2796
-    0.7553    0.1466    0.2682
-    0.7412    0.1411    0.2569
-    0.7271    0.1357    0.2455
-    0.7130    0.1302    0.2341
-    0.6988    0.1247    0.2227
-    0.6847    0.1192    0.2114
-    0.6706    0.1137    0.2000
-    0.6686    0.1141    0.1996
-    0.6667    0.1145    0.1992
-    0.6647    0.1149    0.1988
-    0.6628    0.1153    0.1984
-    0.6608    0.1157    0.1981
-    0.6588    0.1160    0.1977
-    0.6569    0.1164    0.1973
-    0.6549    0.1168    0.1969
-    0.6530    0.1172    0.1965
-    0.6510    0.1176    0.1961];
-set(gca,'CameraPosition',[2 2 2])
-hold on
-Xset=r.*cos(t);Yset=r.*sin(t);
-sf=surface(Xset,Yset,h,'EdgeAlpha',0.1,...
-    'EdgeColor',[0.5 0.5 0.5],'FaceColor','interp');
-colormap(map)
-
-
-theta=0;
-while 1
-    theta=theta+0.02;
-    set(sf,'XData',Xset.*cos(theta)-Yset.*sin(theta),...
-        'YData',Xset.*sin(theta)+Yset.*cos(theta))
-    pause(0.01)
-end
-
-
-end

MatLab/008/readme.md

-<img src="https://cdn.jsdelivr.net/gh/sun0225SUN/Awesome-Love-Code/assets/img/matlab/008.jpg"/>

MatLab/008/roseBall.m

-function roseBall
-clear;clc
-%曲面数据计算
-%==========================================================================
-[x,t]=meshgrid((0:24)./24,(0:0.5:575)./575.*20.*pi+4*pi);
-p=(pi/2)*exp(-t./(8*pi));
-change=sin(15*t)/150;
-u=1-(1-mod(3.6*t,2*pi)./pi).^4./2+change;
-y=2*(x.^2-x).^2.*sin(p);
-
-r=u.*(x.*sin(p)+y.*cos(p));
-h=u.*(x.*cos(p)-y.*sin(p));
-
-%颜色映射表
-%==========================================================================
-hMap=(h-min(min(h)))./(max(max(h))-min(min(h)));
-col=size(hMap,2);
-colorList=[0.0200    0.0400    0.3900
-         0    0.0900    0.5800
-         0    0.1300    0.6400
-    0.0200    0.0600    0.6900
-         0    0.0800    0.7900
-    0.0100    0.1800    0.8500
-         0    0.1300    0.9600
-    0.0100    0.2600    0.9900
-         0    0.3500    0.9900
-    0.0700    0.6200    1.0000
-    0.1700    0.6900    1.0000];
-
-% colorList=[0.2100    0.0900    0.3800
-%     0.2900    0.0700    0.4700
-%     0.4000    0.1100    0.4900
-%     0.5500    0.1600    0.5100
-%     0.7500    0.2400    0.4700
-%     0.8900    0.3200    0.4100
-%     0.9700    0.4900    0.3700
-%     1.0000    0.5600    0.4100
-%     1.0000    0.6900    0.4900
-%     1.0000    0.8200    0.5900
-%     0.9900    0.9200    0.6700
-%     0.9800    0.9500    0.7100];
-
-% colorList1=[0.2000    0.0800    0.4300
-%     0.2000    0.1300    0.4600
-%     0.2000    0.2100    0.5000
-%     0.2000    0.2800    0.5300
-%     0.2000    0.3700    0.5800
-%     0.1900    0.4500    0.6200
-%     0.2000    0.4800    0.6400
-%     0.1900    0.5400    0.6700
-%     0.1900    0.5700    0.6900
-%     0.1900    0.7500    0.7800
-%     0.1900    0.8000    0.8100
-% ];
-
-% colorList2=[0.1300    0.1000    0.1600
-%     0.2000    0.0900    0.2000
-%     0.2800    0.0800    0.2300
-%     0.4200    0.0800    0.3000
-%     0.5100    0.0700    0.3400
-%     0.6600    0.1200    0.3500
-%     0.7900    0.2200    0.4000
-%     0.8800    0.3500    0.4700
-%     0.9000    0.4500    0.5400
-%     0.8900    0.7800    0.7900
-% ];
-
-% colorList3=[0.3200    0.3100    0.7600
-%     0.3800    0.3400    0.7600
-%     0.5300    0.4200    0.7500
-%     0.6400    0.4900    0.7300
-%     0.7200    0.5500    0.7200
-%     0.7900    0.6100    0.7100
-%     0.9100    0.7100    0.6800
-%     0.9800    0.7600    0.6700
-% ];
-
-% colorList4=[0.9500    0.2300    0.6600
-%     0.7500    0.2100    0.6000
-%     0.6200    0.2000    0.5700
-%     0.4500    0.1800    0.5200
-%     0.3200    0.2100    0.5200
-%     0.2700    0.3100    0.6000
-%     0.2500    0.3600    0.6400
-%     0.1900    0.4800    0.7400
-% ];
-
-
-colorFunc=colorFuncFactory(colorList);
-dataMap=colorFunc(hMap');
-colorMap(:,:,1)=dataMap(:,1:col);
-colorMap(:,:,2)=dataMap(:,col+1:2*col);
-colorMap(:,:,3)=dataMap(:,2*col+1:3*col);
-
-    function colorFunc=colorFuncFactory(colorList)
-        xx=(0:size(colorList,1)-1)./(size(colorList,1)-1);
-        y1=colorList(:,1);y2=colorList(:,2);y3=colorList(:,3);
-        colorFunc=@(X)[interp1(xx,y1,X,'linear')',interp1(xx,y2,X,'linear')',interp1(xx,y3,X,'linear')'];
-    end
-
-
-%曲面旋转及绘制
-%==========================================================================
-surface(r.*cos(t),r.*sin(t),h+0.35,'EdgeAlpha',0.05,...
-    'EdgeColor',[0 0 0],'FaceColor','interp','CData',colorMap)
-
-hold on
-
-surface(r.*cos(t),r.*sin(t),-h-0.35,'EdgeAlpha',0.05,...
-    'EdgeColor',[0 0 0],'FaceColor','interp','CData',colorMap)
-Xset=r.*cos(t);
-Yset=r.*sin(t);
-Zset=h+0.35;
-
-yaw_z=pi*72/180;
-roll_x=pi-acos(-1/sqrt(5));
-R_z_2=[cos(yaw_z),-sin(yaw_z),0;
-    sin(yaw_z),cos(yaw_z),0;
-    0,0,1];
-R_z_1=[cos(yaw_z/2),-sin(yaw_z/2),0;
-    sin(yaw_z/2),cos(yaw_z/2),0;
-    0,0,1];
-R_x_2=[1,0,0;
-     0,cos(roll_x),-sin(roll_x);
-     0,sin(roll_x),cos(roll_x)];
- 
-[nX,nY,nZ]=rotateXYZ(Xset,Yset,Zset,R_x_2);
-surface(nX,nY,nZ,'EdgeAlpha',0.05,...
-'EdgeColor',[0 0 0],'FaceColor','interp','CData',colorMap)
-
-
-for k=1:4
-    [nX,nY,nZ]=rotateXYZ(nX,nY,nZ,R_z_2);
-    surface(nX,nY,nZ,'EdgeAlpha',0.05,...
-    'EdgeColor',[0 0 0],'FaceColor','interp','CData',colorMap)
-end   
-
-[nX,nY,nZ]=rotateXYZ(nX,nY,nZ,R_z_1);
-
-for k=1:5
-    [nX,nY,nZ]=rotateXYZ(nX,nY,nZ,R_z_2);
-    surface(nX,nY,-nZ,'EdgeAlpha',0.05,...
-    'EdgeColor',[0 0 0],'FaceColor','interp','CData',colorMap)
-end   
- 
-%--------------------------------------------------------------------------
-    function [nX,nY,nZ]=rotateXYZ(X,Y,Z,R)
-        nX=zeros(size(X));
-        nY=zeros(size(Y));
-        nZ=zeros(size(Z));
-        for i=1:size(X,1)
-            for j=1:size(X,2)
-                v=[X(i,j);Y(i,j);Z(i,j)];
-                nv=R*v;
-                nX(i,j)=nv(1);
-                nY(i,j)=nv(2);
-                nZ(i,j)=nv(3);
-            end
-        end
-    end
-%axes属性调整
-%==========================================================================
-ax=gca;
-grid on
-ax.GridLineStyle='--';
-ax.LineWidth=1.2;
-ax.XColor=[1,1,1].*0.4;
-ax.YColor=[1,1,1].*0.4;
-ax.ZColor=[1,1,1].*0.4;
-ax.DataAspectRatio=[1,1,1];
-ax.DataAspectRatioMode='manual';
-ax.CameraPosition=[-6.5914  -24.1625   -0.0384];
-
-
-
-end

MatLab/009/LoveFunc.m

-function LoveFunc
-LoveFunchdl=@(x,a)(x.^2).^(1/3)+0.9.*((3.3-x.^2).^(1/2)).*sin(a.*pi.*x);
-hold on
-grid on
-axis([-3 3,-2 4])
-x=-1.8:0.005:1.8;
-text(0,3.3,'$f(x)=x^{\frac{2}{3}}+0.9(3.3-x^2)^{\frac{1}{2}}\sin(\alpha\pi x)$',...
-    'FontSize',13,'HorizontalAlignment','center','Interpreter','latex');
-txt2=text(-0.35,2.9,'','FontSize',13,'HorizontalAlignment','left','Interpreter','latex','tag','alphadata');
-for a=1:0.01:20
-    delete(findobj('type','line'))
-    AlphaString=['$\alpha=',num2str(a),'$'];
-    Color=([1.0000 0.4902 0.6627]-[0.2118 0.4667 0.9961]).*(a/20)+[0.2118 0.4667 0.9961];
-    set(txt2,'string',AlphaString)
-    plot(x,LoveFunchdl(x,a),'color',Color,'LineWidth',1.2);
-    pause(0.003)
-end
-end

MatLab/009/readme.md

-<img src="https://cdn.jsdelivr.net/gh/sun0225SUN/Awesome-Love-Code/assets/img/matlab/009.jpg"/>

Python/001/love.py

-# coding:utf-8
-import sys
-import os
-import random
-import pygame
-from pygame.locals import *
-
-
-WIDTH, HEIGHT = 640, 480
-BACKGROUND = (255, 255, 255)
-
-def button(text, x, y, w, h, color, screen, size):
-	pygame.draw.rect(screen, color, (x, y, w, h))
-	font = pygame.font.Font('carton.ttf', size)
-	textRender = font.render(text, True, (0, 0, 0))
-	textRect = textRender.get_rect()
-	textRect.center = ((x+w/2), (y+h/2))
-	screen.blit(textRender, textRect)
-
-
-def title(text, screen, scale, color=(0, 0, 0)):
-	font = pygame.font.Font('carton.ttf', WIDTH//(len(text)*2))
-	textRender = font.render(text, True, color)
-	textRect = textRender.get_rect()
-	textRect.midtop = (WIDTH/scale[0], HEIGHT/scale[1])
-	screen.blit(textRender, textRect)
-
-
-def get_random_pos():
-	x, y = random.randint(20, 620), random.randint(20, 460)
-	return x, y
-
-
-def show_like_interface(text, screen, color=(255, 0, 0)):
-	screen.fill(BACKGROUND)
-	font = pygame.font.Font('carton.ttf', WIDTH//(len(text)))
-	textRender = font.render(text, True, color)
-	textRect = textRender.get_rect()
-	textRect.midtop = (WIDTH/2, HEIGHT/2)
-	screen.blit(textRender, textRect)
-	pygame.display.update()
-	while True:
-		for event in pygame.event.get():
-			if event.type == QUIT:
-				pygame.quit()
-				sys.exit()
-
-def show_like_interface1(screen):
-	screen.fill(BACKGROUND)
-	img3 = pygame.image.load("3.png")
-	imgRect = img3.get_rect()
-	imgRect.midtop = WIDTH // 2, HEIGHT // 4
-	screen.blit(img3, imgRect)
-	pygame.display.update()
-	while True:
-		for event in pygame.event.get():
-			if event.type == QUIT:
-				pygame.quit()
-				sys.exit()
-
-def show_like_interface2(text, screen, color=(255, 0, 0)):
-	screen.fill(BACKGROUND)
-	font = pygame.font.Font('carton.ttf', WIDTH//(len(text)))
-	textRender = font.render(text, True, color)
-	textRect = textRender.get_rect()
-	textRect.midtop = (WIDTH/2, HEIGHT/2)
-	screen.blit(textRender, textRect)
-	pygame.display.update()
-	while True:
-		for event in pygame.event.get():
-			if event.type == QUIT:
-					pygame.quit()
-					sys.exit()
-
-
-def main():
-	pygame.init()
-	screen = pygame.display.set_mode((WIDTH, HEIGHT), 0, 32)
-	pygame.display.set_caption('来自一个喜欢你很久的小哥哥')
-	clock = pygame.time.Clock()
-    
-	unlike_x_pos,unlike_y_pos = 370,380
-	unlike_pos_width, unlike_pos_height = 100,50
-	like_x_pos,like_y_pos = 180,370
-	like_pos_width,like_pos_height = 100,50
-
-	running = True
-	button_color = (192, 192, 192)
-	while running:
-		screen.fill(BACKGROUND)
-		img = pygame.image.load("1.png")
-		imgRect = img.get_rect()
-		imgRect.midtop = WIDTH//2, HEIGHT//4
-		screen.blit(img, imgRect)
-		for event in pygame.event.get():
-			if event.type == pygame.MOUSEBUTTONDOWN:
-				mouse_pos = pygame.mouse.get_pos()
-				if mouse_pos[0] < like_x_pos+like_pos_width+5 and mouse_pos[0] > like_x_pos-5 and\
-					mouse_pos[1] < like_y_pos+like_pos_height+5 and mouse_pos[1] > like_y_pos-5:
-					button_color = BACKGROUND
-					running = False
-		mouse_pos = pygame.mouse.get_pos()
-		if mouse_pos[0] < unlike_x_pos+unlike_pos_width+5 and mouse_pos[0] > unlike_x_pos-5 and\
-			mouse_pos[1] < unlike_y_pos+unlike_pos_height+5 and mouse_pos[1] > unlike_y_pos-5:
-			while True:
-				unlike_x_pos, unlike_y_pos = get_random_pos()
-				if mouse_pos[0] < unlike_x_pos+unlike_pos_width+5 and mouse_pos[0] > unlike_x_pos-5 and\
-					mouse_pos[1] < unlike_y_pos+unlike_pos_height+5 and mouse_pos[1] > unlike_y_pos-5:
-					continue
-				break
-		title('小姐姐，我观察你很久了', screen, scale=[2, 10])
-		title('做我女朋友好不好呀? *^_^*', screen, scale=[2, 6])
-		button('好呀', like_x_pos, like_y_pos, like_pos_width, like_pos_height, button_color, screen, 20)
-		button('算了吧', unlike_x_pos, unlike_y_pos, unlike_pos_width/2, unlike_pos_height/2, button_color, screen, 10)
-		pygame.display.flip()
-		pygame.display.update()
-		clock.tick(60)
-
-	if not os.path.exists("3.png"):
-		show_like_interface2('我就知道小姐姐你也喜欢我 *^_^*', screen, color=(0, 0, 0))
-	else:
-		show_like_interface1(screen=screen)
-
-if __name__ == '__main__':
-	main()
\ No newline at end of file

Python/001/love.spec

-# -*- mode: python -*-
-
-block_cipher = None
-
-
-a = Analysis(['love.py'],
-             pathex=['/Users/MING/.virtualenvs/PythonCodingTime-8iFaVP-J/', '/Users/MING/Github/love-with-python/be_my_girlfriend_mac'],
-             binaries=[],
-             datas=[('01.ming','.'),('PingFang.ttc','PingFang.ttc')],
-             hiddenimports=[],
-             hookspath=[],
-             runtime_hooks=[],
-             excludes=[],
-             win_no_prefer_redirects=False,
-             win_private_assemblies=False,
-             cipher=block_cipher)
-             
-pyz = PYZ(a.pure, a.zipped_data,
-             cipher=block_cipher)
-exe = EXE(pyz,
-          a.scripts,
-          a.binaries,
-          a.zipfiles,
-          a.datas,
-          name='love',
-          debug=False,
-          strip=False,
-          upx=True,
-          runtime_tmpdir=None,
-          console=False )

Python/001/readme.md

-# 执行
-
-```shell
-$ python3 -m pip install pygame
-$ cd be_my_girlfriend
-$ python3 love.py
-```
-
-# 效果
-
-<img src="https://cdn.jsdelivr.net/gh/sun0225SUN/Awesome-Love-Code/assets/img/python/001.jpg"/>

Python/002/index.py

-print('\n'.join([line for line in [''.join([('Love'[(x-y) % len('Love')] if ((x*0.05)**2+(y*0.1)**2-1)**3-(x*0.05)**2*(y*0.1)**3 <= 0 else ' ') for x in range(-30, 30)]) for y in range(30, -30, -1)] if line != '                                                            ']))
\ No newline at end of file

Python/002/readme.md

-# 执行
-
-```shell
-python index.py
-```
-
-# 效果
-
-```txt
-                veLoveLov           veLoveLov
-            eLoveLoveLoveLove   eLoveLoveLoveLove
-          veLoveLoveLoveLoveLoveLoveLoveLoveLoveLov
-         veLoveLoveLoveLoveLoveLoveLoveLoveLoveLoveL
-        veLoveLoveLoveLoveLoveLoveLoveLoveLoveLoveLov
-        eLoveLoveLoveLoveLoveLoveLoveLoveLoveLoveLove
-        LoveLoveLoveLoveLoveLoveLoveLoveLoveLoveLoveL
-        oveLoveLoveLoveLoveLoveLoveLoveLoveLoveLoveLo
-        veLoveLoveLoveLoveLoveLoveLoveLoveLoveLoveLov
-        eLoveLoveLoveLoveLoveLoveLoveLoveLoveLoveLove
-         oveLoveLoveLoveLoveLoveLoveLoveLoveLoveLove
-          eLoveLoveLoveLoveLoveLoveLoveLoveLoveLove
-          LoveLoveLoveLoveLoveLoveLoveLoveLoveLoveL
-            eLoveLoveLoveLoveLoveLoveLoveLoveLove
-             oveLoveLoveLoveLoveLoveLoveLoveLove
-              eLoveLoveLoveLoveLoveLoveLoveLove
-                veLoveLoveLoveLoveLoveLoveLov
-                  oveLoveLoveLoveLoveLoveLo
-                    LoveLoveLoveLoveLoveL
-                       LoveLoveLoveLov
-                          LoveLoveL
-                             Lov
-                              v
-```
-
-<img src="https://cdn.jsdelivr.net/gh/sun0225SUN/Awesome-Love-Code/assets/img/python/002.jpg"/>

Python/003/index.py

-import turtle
-import time
- 
- 
-# 清屏函数
-def clear_all():
-    turtle.penup()
-    turtle.goto(0, 0)
-    turtle.color('white')
-    turtle.pensize(800)
-    turtle.pendown()
-    turtle.setheading(0)
-    turtle.fd(300)
-    turtle.bk(600)
- 
- 
-# 重定位海龟的位置
-def go_to(x, y, state):
-    turtle.pendown() if state else turtle.penup()
-    turtle.goto(x, y)
- 
- 
-# 画线
-# state为真时海龟回到原点，为假时不回到原来的出发点
-def draw_line(length, angle, state):
-    turtle.pensize(1)
-    turtle.pendown()
-    turtle.setheading(angle)
-    turtle.fd(length)
-    turtle.bk(length) if state else turtle.penup()
-    turtle.penup()
- 
- 
-# 画箭羽
-def draw_feather(size):
-    angle = 30                          # 箭的倾角
-    feather_num = size//6               # 羽毛的数量
-    feather_length = size // 3          # 羽毛的长度
-    feather_gap = size//10              # 羽毛的间隔
-    for i in range(feather_num):
-        draw_line(feather_gap, angle+180, False)            # 箭柄，不折返
-        draw_line(feather_length, angle + 145, True)        # 羽翼，要折返
-    draw_line(feather_length, angle + 145, False)
-    draw_line(feather_num*feather_gap, angle, False)
-    draw_line(feather_length, angle + 145 + 180, False)
-    for i in range(feather_num):
-        draw_line(feather_gap, angle+180, False)            # 箭柄，不折返
-        draw_line(feather_length, angle - 145, True)        # 羽翼，要折返
-    draw_line(feather_length, angle - 145, False)
-    draw_line(feather_num*feather_gap, angle, False)
-    draw_line(feather_length, angle - 145 + 180, False)
- 
- 
-# 画爱心
-def draw_heart(size):
-    turtle.color('red', 'pink')
-    turtle.pensize(2)
-    turtle.pendown()
-    turtle.setheading(150)
-    turtle.begin_fill()
-    turtle.fd(size)
-    turtle.circle(size * -3.745, 45)
-    turtle.circle(size * -1.431, 165)
-    turtle.left(120)
-    turtle.circle(size * -1.431, 165)
-    turtle.circle(size * -3.745, 45)
-    turtle.fd(size)
-    turtle.end_fill()
-
-
-def hart_arc():
-    for i in range(200):
-        turtle.right(1)
-        turtle.forward(2)
- 
-# 画箭
-def draw_arrow(size):
-    angle = 30
-    turtle.color('black')
-    draw_feather(size)
-    turtle.pensize(4)
-    turtle.setheading(angle)
-    turtle.pendown()
-    turtle.fd(size*2)
- 
- 
-# 一箭穿心
-# 箭的头没有画出来，而是用海龟来代替
-def arrow_heart(x, y, size):
-    go_to(x, y, False)
-    draw_heart(size*1.15)
-    turtle.setheading(-150)
-    turtle.penup()
-    turtle.fd(size*2.2)
-    draw_heart(size)
-    turtle.penup()
-    turtle.setheading(150)
-    turtle.fd(size * 2.2)
-    draw_arrow(size)
- 
- 
-# 画出发射爱心的小人
-def draw_people(x, y):
-    turtle.penup()
-    turtle.goto(x, y)
-    turtle.pendown()
-    turtle.pensize(2)
-    turtle.color('black')
-    turtle.setheading(0)
-    turtle.circle(60, 360)
-    turtle.penup()
-    turtle.setheading(90)
-    turtle.fd(75)
-    turtle.setheading(180)
-    turtle.fd(20)
-    turtle.pensize(4)
-    turtle.pendown()
-    turtle.circle(2, 360)
-    turtle.setheading(0)
-    turtle.penup()
-    turtle.fd(40)
-    turtle.pensize(4)
-    turtle.pendown()
-    turtle.circle(-2, 360)
-    turtle.penup()
-    turtle.goto(x, y)
-    turtle.setheading(-90)
-    turtle.pendown()
-    turtle.fd(20)
-    turtle.setheading(0)
-    turtle.fd(35)
-    turtle.setheading(60)
-    turtle.fd(10)
-    turtle.penup()
-    turtle.goto(x, y)
-    turtle.setheading(-90)
-    turtle.pendown()
-    turtle.fd(40)
-    turtle.setheading(0)
-    turtle.fd(35)
-    turtle.setheading(-60)
-    turtle.fd(10)
-    turtle.penup()
-    turtle.goto(x, y)
-    turtle.setheading(-90)
-    turtle.pendown()
-    turtle.fd(60)
-    turtle.setheading(-135)
-    turtle.fd(60)
-    turtle.bk(60)
-    turtle.setheading(-45)
-    turtle.fd(30)
-    turtle.setheading(-135)
-    turtle.fd(35)
-    turtle.penup()
- 
- 
-# 第一个画面，显示文字
-def page0():
-    turtle.penup()
-    turtle.goto(-350, 0)
-    turtle.color('black')
-    turtle.write('专属于我们的情人节', font=('宋体', 60, 'normal'))
-    time.sleep(3)
- 
- 
-# 第二个画面，显示发射爱心的小人
-def page1():
-    turtle.speed(10)
-    draw_people(-250, 20)
-    turtle.penup()
-    turtle.goto(-150, -30)
-    draw_heart(14)
-    turtle.penup()
-    turtle.goto(-20, -60)
-    draw_heart(25)
-    turtle.penup()
-    turtle.goto(250, -100)
-    draw_heart(45)
-    turtle.hideturtle()
-    time.sleep(3)
- 
- 
-# 最后一个画面，一箭穿心
-def page2():
-    turtle.speed(1)
-    turtle.penup()
-    turtle.goto(-200, -200)
-    turtle.color('blue')
-    turtle.pendown()
-    turtle.write('大笨蛋         小笨蛋', font=('wisdom', 50, 'normal'))
-    turtle.penup()
-    turtle.goto(0, -190)
-    draw_heart(10)
-    arrow_heart(20, -60, 51)
-    turtle.showturtle()
- 
- 
-def main():
-    turtle.setup(900, 500)
-    page0()
-    clear_all()
-    page1()
-    clear_all()
-    page2()
-    turtle.done()
- 
- 
-main()

Python/003/readme.md

-# 执行
-
-```shell
-python index.py
-```
-
-# 效果
-
-<img src="https://cdn.jsdelivr.net/gh/sun0225SUN/Awesome-Love-Code/assets/img/python/003.jpg"/>

Python/004/code/flower_1.py

-"""
--*- coding: utf-8 -*-
-@Time : 2022/2/13 下午 1:34
-@Author : SunGuoqi
-@Website : https://sunguoqi.com
-@Github: https://github.com/sun0225SUN
-"""
-
-import numpy as np
-import matplotlib.pyplot as plt
-
-fig = plt.figure()
-ax = fig.gca(projection='3d')
-[x, t] = np.meshgrid(
-    np.array(range(25)) / 24.0,
-    np.arange(0, 575.5, 0.5) / 575 * 30 * np.pi - 4 * np.pi)
-p = (np.pi / 2) * np.exp(-t / (8 * np.pi))
-change = np.sin(20 * t) / 50
-u = 1 - (1 - np.mod(3.3 * t, 2 * np.pi) / np.pi)**4 / 2 + change
-y = 2 * (x**2 - x)**2 * np.sin(p)
-r = u * (x * np.sin(p) + y * np.cos(p)) * 1.5
-h = u * (x * np.cos(p) - y * np.sin(p))
-c = plt.get_cmap('magma')
-surf = ax.plot_surface(r * np.cos(t),
-                       r * np.sin(t),
-                       h,
-                       rstride=1,
-                       cstride=1,
-                       cmap=c,
-                       linewidth=0,
-                       antialiased=True)
-plt.show()

Python/004/code/flower_2.py

-import numpy as np
-import matplotlib.pyplot as plt
-from matplotlib import cm
-
-fig = plt.figure()
-ax = fig.gca(projection='3d')
-[x, t] = np.meshgrid(
-    np.array(range(25)) / 24.0,
-    np.arange(0, 575.5, 0.5) / 575 * 17 * np.pi - 2 * np.pi)
-p = (np.pi / 2) * np.exp(-t / (8 * np.pi))
-u = 1 - (1 - np.mod(3.6 * t, 2 * np.pi) / np.pi)**4 / 2
-y = 2 * (x**2 - x)**2 * np.sin(p)
-r = u * (x * np.sin(p) + y * np.cos(p))
-h = u * (x * np.cos(p) - y * np.sin(p))
-c = cm.gist_rainbow_r
-surf = ax.plot_surface(r * np.cos(t),
-                       r * np.sin(t),
-                       h,
-                       rstride=1,
-                       cstride=1,
-                       cmap=c,
-                       linewidth=0,
-                       antialiased=True)
-plt.show()

Python/004/code/flower_3.py

-import numpy as np
-import matplotlib.pyplot as plt
-from matplotlib import cm
-
-fig = plt.figure()
-ax = fig.gca(projection='3d')
-[x, t] = np.meshgrid(
-    np.array(range(25)) / 24.0,
-    np.arange(0, 575.5, 0.5) / 575 * 17 * np.pi - 2 * np.pi)
-p = (np.pi / 2) * np.exp(-t / (8 * np.pi))
-u = 1 - (1 - np.mod(3.6 * t, 2 * np.pi) / np.pi)**4 / 2
-y = 2 * (x**2 - x)**2 * np.sin(p)
-r = u * (x * np.sin(p) + y * np.cos(p))
-h = u * (x * np.cos(p) - y * np.sin(p))
-c = cm.get_cmap('spring_r')
-surf = ax.plot_surface(r * np.cos(t),
-                       r * np.sin(t),
-                       h,
-                       rstride=1,
-                       cstride=1,
-                       cmap=c,
-                       linewidth=0,
-                       antialiased=True)
-plt.show()

Python/004/code/flower_4.py

-import numpy as np
-import matplotlib.pyplot as plt
-
-fig = plt.figure()
-ax = fig.gca(projection='3d')
-# 将相位向后移动了6*pi
-[x, t] = np.meshgrid(
-    np.array(range(25)) / 24.0,
-    np.arange(0, 575.5, 0.5) / 575 * 20 * np.pi + 4 * np.pi)
-p = (np.pi / 2) * np.exp(-t / (8 * np.pi))
-# 添加边缘扰动
-change = np.sin(15 * t) / 150
-# 将t的参数减少，使花瓣的角度变大
-u = 1 - (1 - np.mod(3.3 * t, 2 * np.pi) / np.pi)**4 / 2 + change
-y = 2 * (x**2 - x)**2 * np.sin(p)
-r = u * (x * np.sin(p) + y * np.cos(p))
-h = u * (x * np.cos(p) - y * np.sin(p))
-c = plt.get_cmap('Reds')
-surf = ax.plot_surface(r * np.cos(t),
-                       r * np.sin(t),
-                       h,
-                       rstride=1,
-                       cstride=1,
-                       cmap=c,
-                       linewidth=0,
-                       antialiased=True)
-plt.show()

Python/004/code/flower_5.py

-import numpy as np
-import matplotlib.pyplot as plt
-
-fig = plt.figure()
-ax = fig.gca(projection='3d')
-[x, t] = np.meshgrid(
-    np.array(range(25)) / 24.0,
-    np.arange(0, 575.5, 0.5) / 575 * 6 * np.pi - 4 * np.pi)
-p = (np.pi / 2) * np.exp(-t / (8 * np.pi))
-change = np.sin(10 * t) / 20
-u = 1 - (1 - np.mod(5.2 * t, 2 * np.pi) / np.pi)**4 / 2 + change
-y = 2 * (x**2 - x)**2 * np.sin(p)
-r = u * (x * np.sin(p) + y * np.cos(p)) * 1.5
-h = u * (x * np.cos(p) - y * np.sin(p))
-c = plt.get_cmap('spring_r')
-surf = ax.plot_surface(r * np.cos(t),
-                       r * np.sin(t),
-                       h,
-                       rstride=1,
-                       cstride=1,
-                       cmap=c,
-                       linewidth=0,
-                       antialiased=True)
-plt.show()

Python/004/readme.md

-# 安装所需库
-
-```
-pip install pillow
-pip install matplotlib
-pip install numpy
-```
-
-# 运行
-
-```
-python 文件名.py
-```
-
-![](./img/eg.png)
-
-# 效果
-
-![01](./img/01.png)
-
-![02](./img/02.png)
-
-![03](./img/03.png)
-
-![04](./img/04.png)
-
-![05](./img/05.png)

Python/005/love.py

-import matplotlib.pyplot as plt
-import numpy as np
-
-
-def heart_3d(x, y, z):
-    return (x**2 + (9 / 4) * y**2 + z**2 -
-            1)**3 - x**2 * z**3 - (9 / 80) * y**2 * z**3
-
-
-def plot_implicit(fn, bbox=(-1.5, 1.5)):
-    ''' create a plot of an implicit function
-    fn ...implicit function (plot where fn==0)
-    bbox ..the x,y,and z limits of plotted interval'''
-    xmin, xmax, ymin, ymax, zmin, zmax = bbox * 3
-    fig = plt.figure()
-    ax = fig.add_subplot(111, projection='3d')
-    A = np.linspace(xmin, xmax, 100)  # resolution of the contour
-    B = np.linspace(xmin, xmax, 40)  # number of slices
-    A1, A2 = np.meshgrid(A, A)  # grid on which the contour is plotted
-
-    for z in B:  # plot contours in the XY plane
-        X, Y = A1, A2
-        Z = fn(X, Y, z)
-        cset = ax.contour(X, Y, Z + z, [z], zdir='z', colors=('red', ))
-        # [z] defines the only level to plot
-        # for this contour for this value of z
-
-    for y in B:  # plot contours in the XZ plane
-        X, Z = A1, A2
-        Y = fn(X, y, Z)
-        cset = ax.contour(X, Y + y, Z, [y], zdir='y', colors=('red', ))
-
-    for x in B:  # plot contours in the YZ plane
-        Y, Z = A1, A2
-        X = fn(x, Y, Z)
-        cset = ax.contour(X + x, Y, Z, [x], zdir='x', colors=('red', ))
-
-    # must set plot limits because the contour will likely extend
-    # way beyond the displayed level. Otherwise matplotlib extends the plot limits
-    # to encompass all values in the contour.
-    ax.set_zlim3d(zmin, zmax)
-    ax.set_xlim3d(xmin, xmax)
-    ax.set_ylim3d(ymin, ymax)
-
-    plt.show()
-
-
-if __name__ == '__main__':
-    plot_implicit(heart_3d)
\ No newline at end of file

Python/005/readme.md

-# 安装所需库
-
-```
-pip install matplotlib
-pip install numpy
-```
-
-# 运行
-
-```
-python 文件名.py
-```
-
-![](./img/eg.jpg)
-
-# 运行
-
-```
-python 文件名.py
-```
-
-![](./img/1.jpg)

Python/006/index.py

-import turtle
-import math
-
-turtle.pen()
-t = turtle
-t.up()
-t.goto(0, 150)
-t.down()
-t.color('red')
-t.begin_fill()
-t.fillcolor('red')
-t.speed(1)
-t.left(45)
-t.forward(150)
-t.right(45)
-t.forward(100)
-t.right(45)
-t.forward(100)
-t.right(45)
-t.forward(100)
-t.right(45)
-t.forward(250 + math.sqrt(2) * 100)
-t.right(90)
-t.speed(2)
-t.forward(250 + 100 * math.sqrt(2))
-t.right(45)
-t.forward(100)
-t.right(45)
-t.forward(100)
-t.right(45)
-t.forward(100)
-t.right(45)
-t.forward(150)
-t.end_fill()
-t.goto(-10, 0)
-t.pencolor('white')
-# L
-t.pensize(10)
-t.goto(-50, 0)
-t.goto(-50, 80)
-t.up()
-# I
-t.goto(-100, 0)
-t.down()
-t.goto(-160, 0)
-t.goto(-130, 0)
-t.goto(-130, 80)
-t.goto(-160, 80)
-t.goto(-100, 80)
-t.up()
-# O
-t.goto(10, 25)
-t.down()
-t.right(45)
-t.circle(25, extent=180)
-t.goto(60, 55)
-t.circle(25, extent=180)
-t.goto(10, 25)
-t.up()
-t.goto(75, 80)
-t.down()
-t.goto(100, 0)
-t.goto(125, 80)
-t.up()
-t.goto(180, 80)
-t.down()
-t.goto(140, 80)
-t.goto(140, 0)
-t.goto(180, 0)
-t.up()
-t.goto(180, 40)
-t.down()
-t.goto(140, 40)
-# U
-t.up()
-t.goto(-40, -30)
-t.down()
-t.goto(-40, -80)
-t.circle(40, extent=180)
-t.goto(40, -30)
-t.hideturtle()

Python/006/readme.md

-# 安装依赖
-
-```
-pip install turtle
-```
-
-# 运行
-
-```
-python 文件名.py
-```
-
-# 效果
-
-![](./img/1.jpg)

Python/007/main.py

-from tkinter import *
-from matplotlib import pyplot as plt
-from PIL import Image
-import random
-import math
-import numpy as np
-import os
-import colorsys
-import cv2
-from scipy.ndimage.filters import gaussian_filter
-
-canvas_width = 600
-canvas_height = 600
-world_width = 0.05
-world_heigth = 0.05
-
-# 中间心的参数
-points = None
-fixed_point_size = 20000
-fixed_scale_range = (4, 4.3)
-min_scale = np.array([1.0, 1.0, 1.0]) * 0.9
-max_scale = np.array([1.0, 1.0, 1.0]) * 0.9
-min_heart_scale = -15
-max_heart_scale = 16
-
-# 外围随机心参数
-random_point_szie = 7000
-random_scale_range = (3.5, 3.9)
-random_point_maxvar = 0.2
-
-# 心算法参数
-mid_point_ignore = 0.95
-
-# 相机参数
-camera_close_plane = 0.1
-camera_position = np.array([0.0, -2.0, 0.0])
-
-# 点的颜色
-hue = 0.92
-color_strength = 255
-
-# 常用向量缓存
-zero_scale = np.array([0.0, 0.0, 0.0])
-unit_scale = np.array([1.0, 1.0, 1.0])
-color_white = np.array([255, 255, 255])
-axis_y = np.array([0.0, 1.0, 0.0])
-
-# 渲染缓存
-render_buffer = np.empty((canvas_width, canvas_height, 3), dtype=int)
-strength_buffer = np.empty((canvas_width, canvas_height), dtype=float)
-
-# 随机点文件缓存
-points_file = "temp.txt"
-
-# 渲染结果
-total_frames = 30
-output_dir = "./output"
-
-# 格式
-image_fmt = "jpg"
-
-def color(value):
-    digit = list(map(str, range(10))) + list("ABCDEF")
-    string = '#'
-    for i in value:
-        a1 = i // 16
-        a2 = i % 16
-        string += digit[a1] + digit[a2]
-    return string
-    
-
-def heart_func(x, y, z, scale):
-    bscale = scale
-    bscale_half = bscale / 2
-    x = x * bscale - bscale_half
-    y = y * bscale - bscale_half
-    z = z * bscale - bscale_half
-    return (x**2 + 9/4*(y**2) + z**2 - 1)**3 - (x**2)*(z**3) - 9/200*(y**2)*(z**3)
-
-def lerp_vector(a, b, ratio):
-    result = a.copy()
-    for i in range(3):
-        result[i] = a[i] + (b[i] - a[i]) * ratio
-    return result
-
-def lerp_int(a, b, ratio):
-    return (int)(a + (b - a) * ratio)
-
-def lerp_float(a, b, ratio):
-    return (a + (b - a) * ratio)
-
-def distance(point):
-    return (point[0]**2 + point[1]**2 + point[2]**2) ** 0.5
-
-def dot(a, b):
-    return a[0]*b[0] + a[1]*b[1] + a[2]*b[2]
-
-def inside_rand(tense):
-    x = random.random()
-    y = -tense * math.log(x)
-    return y
-
-# 生成中间心
-def genPoints(pointCount, heartScales):
-    result = np.empty((pointCount, 3))
-    index = 0
-    while index < pointCount:
-        # 生成随机点
-        x = random.random()
-        y = random.random()
-        z = random.random()
-
-        # 扣掉心中间的点
-        mheartValue = heart_func(x, 0.5, z, heartScales[1])
-        mid_ignore = random.random()
-        if mheartValue < 0 and mid_ignore < mid_point_ignore:
-            continue
-        
-        heartValue = heart_func(x, y, z, heartScales[0])
-        z_shrink = 0.01
-        sz = z - z_shrink
-        sheartValue = heart_func(x, y, sz, heartScales[1])
-
-        # 保留在心边上的点
-        if heartValue < 0 and sheartValue > 0:
-            result[index] = [x - 0.5, y - 0.5, z - 0.5]
-
-            # 向内扩散
-            len = 0.7
-            result[index] = result[index] * (1 - len * inside_rand(0.2))
-
-            # 重新赋予深度
-            newY = random.random() - 0.5
-            rheartValue = heart_func(result[index][0] + 0.5, newY + 0.5, result[index][2] + 0.5, heartScales[0])
-            if rheartValue > 0:
-                continue
-            result[index][1] = newY
-
-            # 删掉肚脐眼
-            dist = distance(result[index])
-            if dist < 0.12:
-                continue
-            
-            index = index + 1
-            if index % 100 == 0:
-                print("{ind} generated {per}%".format(ind=index, per=((index / pointCount) * 100)))
-
-    return result
-
-# 生成随机心
-def genRandPoints(pointCount, heartScales, maxVar, ratio):
-    result = np.empty((pointCount, 3))
-    index = 0
-    while index < pointCount:
-        x = random.random()
-        y = random.random()
-        z = random.random()
-        mheartValue = heart_func(x, 0.5, z, heartScales[1])
-        mid_ignore = random.random()
-        if mheartValue < 0 and mid_ignore < mid_point_ignore:
-            continue
-
-        heartValue = heart_func(x, y, z, heartScales[0])
-        sheartValue = heart_func(x, y, z, heartScales[1])
-
-        if heartValue < 0 and sheartValue > 0:
-            result[index] = [x - 0.5, y - 0.5, z - 0.5]
-            dist = distance(result[index])
-            if dist < 0.12:
-                continue
-
-            len = 0.7
-            result[index] = result[index] * (1 - len * inside_rand(0.2))
-            index = index + 1
-
-    for i in range(pointCount):
-        var = maxVar * ratio
-        randScale = 1 + random.normalvariate(0, var)
-        result[i] = result[i] * randScale
-
-    return result
-
-# 世界坐标到相机本地坐标
-def world_2_cameraLocalSapce(world_point):
-    new_point = world_point.copy()
-    new_point[1] = new_point[1] + camera_position[1]
-    return new_point
-
-# 相机本地坐标到相机空间坐标
-def cameraLocal_2_cameraSpace(cameraLocalPoint):
-    depth = distance(cameraLocalPoint)
-    cx = cameraLocalPoint[0] * (camera_close_plane / cameraLocalPoint[1])
-    cz = -cameraLocalPoint[2] * (cx / cameraLocalPoint[0])
-    cameraLocalPoint[0] = cx
-    cameraLocalPoint[1] = cz
-    return cameraLocalPoint, depth
-
-# 相机空间坐标到屏幕坐标
-def camerSpace_2_screenSpace(cameraSpace):
-    x = cameraSpace[0]
-    y = cameraSpace[1]
-
-    # convert to view space
-    centerx = canvas_width / 2
-    centery = canvas_height / 2
-    ratiox = canvas_width / world_width
-    ratioy = canvas_height / world_heigth
-
-    viewx = centerx + x * ratiox
-    viewy = canvas_height - (centery + y * ratioy)
-
-    cameraSpace[0] = viewx
-    cameraSpace[1] = viewy
-    return cameraSpace.astype(int)
-
-
-# 绘制世界坐标下的点
-def draw_point(worldPoint):
-    cameraLocal = world_2_cameraLocalSapce(worldPoint)
-    cameraSpsace, depth = cameraLocal_2_cameraSpace(cameraLocal)
-    screeSpace = camerSpace_2_screenSpace(cameraSpsace)
-
-    draw_size = int(random.random() * 3 + 1)
-    draw_on_buffer(screeSpace, depth, draw_size)
-
-# 绘制到缓存上
-def draw_on_buffer(screenPos, depth, draw_size):
-    if draw_size == 0:
-        return
-    elif draw_size == 1:
-        draw_point_on_buffer(screenPos[0], screenPos[1], color_strength, depth)
-    elif draw_size == 2:
-        draw_point_on_buffer(screenPos[0], screenPos[1], color_strength, depth)
-        draw_point_on_buffer(screenPos[0] + 1, screenPos[1] + 1, color_strength, depth)
-    elif draw_size == 3:
-        draw_point_on_buffer(screenPos[0], screenPos[1], color_strength, depth)
-        draw_point_on_buffer(screenPos[0] + 1, screenPos[1] + 1, color_strength, depth)
-        draw_point_on_buffer(screenPos[0] + 1, screenPos[1], color_strength, depth)
-    elif draw_size == 4:
-        draw_point_on_buffer(screenPos[0], screenPos[1], color_strength, depth)
-        draw_point_on_buffer(screenPos[0] + 1, screenPos[1], color_strength, depth)
-        draw_point_on_buffer(screenPos[0], screenPos[1] + 1, color_strength, depth)
-        draw_point_on_buffer(screenPos[0] + 1, screenPos[1] + 1, color_strength, depth)
-
-
-# 根据色调和颜色强度获取颜色
-def get_color(strength):
-    result = None
-    if strength >= 1:
-        result = colorsys.hsv_to_rgb(hue, 2 - strength, 1)
-    else:
-        result = colorsys.hsv_to_rgb(hue, 1, strength)
-    r = min(result[0] * 256, 255)
-    g = min(result[1] * 256, 255)
-    b = min(result[2] * 256, 255)
-    return np.array((r, g, b), dtype=int)
-
-# 可以根据深度做一些好玩的
-def draw_point_on_buffer(x, y, color, depth):
-    if x < 0 or x >= canvas_width or y < 0 or y >= canvas_height:
-        return
-
-    # 混合
-    strength = float(color) / 255
-    strength_buffer[x, y] = strength_buffer[x, y] + strength
-
-# 绘制缓存
-def draw_buffer_on_canvas(output = None):
-    render_buffer.fill(0)
-    for i in range(render_buffer.shape[0]):
-        for j in range(render_buffer.shape[1]):
-            render_buffer[i, j] = get_color(strength_buffer[i, j])
-    im = Image.fromarray(np.uint8(render_buffer))
-    im = im.rotate(-90)
-    if output is None:
-        plt.imshow(im)
-        plt.show()
-    else:
-        im.save(output)
-
-
-def paint_heart(ratio, randratio, outputFile = None):
-    global strength_buffer
-    global render_buffer
-    global points
-
-    # 清空缓存
-    strength_buffer.fill(0)
-
-    for i in range(fixed_point_size):
-        # 缩放
-        point = points[i] * lerp_vector(min_scale, max_scale, ratio)
-
-        # 球型场
-        dist = distance(point)
-        radius = 0.4
-        sphere_scale = radius / dist
-        point = point * lerp_float(0.9, sphere_scale, ratio * 0.3)
-
-        # 绘制
-        draw_point(point)
-
-    # 生成一组随机点
-    randPoints = genRandPoints(random_point_szie, random_scale_range, random_point_maxvar, randratio)
-    for i in range(random_point_szie):
-        # 绘制
-        draw_point(randPoints[i])
-
-    # 高斯模糊
-    for i in range(1):
-        strength_buffer = gaussian_filter(strength_buffer, sigma=0.8)
-
-    # 绘制缓存
-    draw_buffer_on_canvas(outputFile)
-
-def show_images():
-    img = None
-    for i in range(total_frames):
-        save_name = "{name}.{fmt}".format(name=i, fmt=image_fmt)
-        save_path = os.path.join(output_dir, save_name)
-        img = cv2.imread(save_path, cv2.IMREAD_ANYCOLOR)
-        cv2.imshow("Img", img)
-        cv2.waitKey(25)
-
-
-def gen_images():
-    global points
-
-    if not os.path.isdir(output_dir):
-        os.mkdir(output_dir)
-    
-    # 尝试加载或生成中间心
-    if not os.path.exists(points_file):
-        print("未发现缓存点，重新生成中")
-        points = genPoints(fixed_point_size, fixed_scale_range)
-        np.savetxt(points_file, points)
-    else:
-        print("发现缓存文件，跳过生成")
-        points = np.loadtxt(points_file)
-
-    for i in range(total_frames):
-        print("正在处理图片... ", i)
-        frame_ratio = float(i) / (total_frames - 1)
-        frame_ratio = frame_ratio ** 2
-        ratio = math.sin(frame_ratio * math.pi) * 0.743144
-        randratio = math.sin(frame_ratio * math.pi * 2 + total_frames / 2)
-        save_name = "{name}.{fmt}".format(name=i, fmt=image_fmt)
-        save_path = os.path.join(output_dir, save_name)
-        paint_heart(ratio, randratio, save_path)
-        print("图片已保存至", save_path)
-
-
-if __name__ == "__main__":
-    gen_images()
-    while True:
-        show_images()

Python/008/main.py

-import random
-from math import sin, cos, pi, log
-from tkinter import *
- 
-CANVAS_WIDTH = 640  # 画布的宽
-CANVAS_HEIGHT = 640  # 画布的高
-CANVAS_CENTER_X = CANVAS_WIDTH / 2  # 画布中心的X轴坐标
-CANVAS_CENTER_Y = CANVAS_HEIGHT / 2  # 画布中心的Y轴坐标
-IMAGE_ENLARGE = 11  # 放大比例
-HEART_COLOR = "#e77c8e"  # 心的颜色#ff7171
- 
- 
-def heart_function(t, shrink_ratio: float = IMAGE_ENLARGE):
-    """
-    “爱心函数生成器”
-    :param shrink_ratio: 放大比例
-    :param t: 参数
-    :return: 坐标
-    """
-    # 基础函数
-    x = 16 * (sin(t) ** 3)
-    y = -(13 * cos(t) - 5 * cos(2 * t) - 2 * cos(3 * t) - cos(4 * t))
- 
-    # 放大
-    x *= shrink_ratio
-    y *= shrink_ratio
- 
-    # 移到画布中央
-    x += CANVAS_CENTER_X
-    y += CANVAS_CENTER_Y
- 
-    return int(x), int(y)
- 
- 
-def scatter_inside(x, y, beta=0.15):
-    """
-    随机内部扩散
-    :param x: 原x
-    :param y: 原y
-    :param beta: 强度
-    :return: 新坐标
-    """
-    ratio_x = - beta * log(random.random())
-    ratio_y = - beta * log(random.random())
- 
-    dx = ratio_x * (x - CANVAS_CENTER_X)
-    dy = ratio_y * (y - CANVAS_CENTER_Y)
- 
-    return x - dx, y - dy
- 
- 
-def shrink(x, y, ratio):
-    """
-    抖动
-    :param x: 原x
-    :param y: 原y
-    :param ratio: 比例
-    :return: 新坐标
-    """
-    force = -1 / (((x - CANVAS_CENTER_X) ** 2 + (y - CANVAS_CENTER_Y) ** 2) ** 0.6)  # 这个参数...
-    dx = ratio * force * (x - CANVAS_CENTER_X)
-    dy = ratio * force * (y - CANVAS_CENTER_Y)
-    return x - dx, y - dy
- 
- 
-def curve(p):
-    """
-    自定义曲线函数，调整跳动周期
-    :param p: 参数
-    :return: 正弦
-    """
-    # 可以尝试换其他的动态函数，达到更有力量的效果（贝塞尔？）
-    return 2 * (3 * sin(4 * p)) / (2 * pi)
- 
- 
-class Heart:
-    """
-    爱心类
-    """
- 
-    def __init__(self, generate_frame=20):
-        self._points = set()  # 原始爱心坐标集合
-        self._edge_diffusion_points = set()  # 边缘扩散效果点坐标集合
-        self._center_diffusion_points = set()  # 中心扩散效果点坐标集合
-        self.all_points = {}  # 每帧动态点坐标
-        self.build(2000)
- 
-        self.random_halo = 1000
- 
-        self.generate_frame = generate_frame
-        for frame in range(generate_frame):
-            self.calc(frame)
- 
-    def build(self, number):
-        # 爱心
-        for _ in range(number):
-            t = random.uniform(0, 2 * pi)  # 随机不到的地方造成爱心有缺口
-            x, y = heart_function(t)
-            self._points.add((x, y))
- 
-        # 爱心内扩散
-        for _x, _y in list(self._points):
-            for _ in range(3):
-                x, y = scatter_inside(_x, _y, 0.05)
-                self._edge_diffusion_points.add((x, y))
- 
-        # 爱心内再次扩散
-        point_list = list(self._points)
-        for _ in range(4000):
-            x, y = random.choice(point_list)
-            x, y = scatter_inside(x, y, 0.17)
-            self._center_diffusion_points.add((x, y))
- 
-    @staticmethod
-    def calc_position(x, y, ratio):
-        # 调整缩放比例
-        force = 1 / (((x - CANVAS_CENTER_X) ** 2 + (y - CANVAS_CENTER_Y) ** 2) ** 0.520)  # 魔法参数
- 
-        dx = ratio * force * (x - CANVAS_CENTER_X) + random.randint(-1, 1)
-        dy = ratio * force * (y - CANVAS_CENTER_Y) + random.randint(-1, 1)
- 
-        return x - dx, y - dy
- 
-    def calc(self, generate_frame):
-        ratio = 10 * curve(generate_frame / 10 * pi)  # 圆滑的周期的缩放比例
- 
-        halo_radius = int(4 + 6 * (1 + curve(generate_frame / 10 * pi)))
-        halo_number = int(3000 + 4000 * abs(curve(generate_frame / 10 * pi) ** 2))
- 
-        all_points = []
- 
-        # 光环
-        heart_halo_point = set()  # 光环的点坐标集合
-        for _ in range(halo_number):
-            t = random.uniform(0, 2 * pi)  # 随机不到的地方造成爱心有缺口
-            x, y = heart_function(t, shrink_ratio=11.6)  # 魔法参数
-            x, y = shrink(x, y, halo_radius)
-            if (x, y) not in heart_halo_point:
-                # 处理新的点
-                heart_halo_point.add((x, y))
-                x += random.randint(-14, 14)
-                y += random.randint(-14, 14)
-                size = random.choice((1, 2, 2))
-                all_points.append((x, y, size))
- 
-        # 轮廓
-        for x, y in self._points:
-            x, y = self.calc_position(x, y, ratio)
-            size = random.randint(1, 3)
-            all_points.append((x, y, size))
- 
-        # 内容
-        for x, y in self._edge_diffusion_points:
-            x, y = self.calc_position(x, y, ratio)
-            size = random.randint(1, 2)
-            all_points.append((x, y, size))
- 
-        for x, y in self._center_diffusion_points:
-            x, y = self.calc_position(x, y, ratio)
-            size = random.randint(1, 2)
-            all_points.append((x, y, size))
- 
-        self.all_points[generate_frame] = all_points
- 
-    def render(self, render_canvas, render_frame):
-        for x, y, size in self.all_points[render_frame % self.generate_frame]:
-            render_canvas.create_rectangle(x, y, x + size, y + size, width=0, fill=HEART_COLOR)
- 
- 
-def draw(main: Tk, render_canvas: Canvas, render_heart: Heart, render_frame=0):
-    render_canvas.delete('all')
-    render_heart.render(render_canvas, render_frame)
-    main.after(160, draw, main, render_canvas, render_heart, render_frame + 1)
- 
- 
-if __name__ == '__main__':
-    root = Tk()  # 一个Tk
-    root.title('宝贝爱你哟')#将标题设置为'tkinter'
-    canvas = Canvas(root, bg='black', height=CANVAS_HEIGHT, width=CANVAS_WIDTH)
-    canvas.pack()
-    heart = Heart()  # 心
-    draw(root, canvas, heart)  # 开始画
-    root.mainloop()
\ No newline at end of file