如何通過(guò)編程對(duì)喜歡的人表白？

「來(lái)源: ｜EETOP ID：eetop-1」

來(lái)源：科研狗、EETOP

這一篇圖文推送的創(chuàng)作靈感源自于6月份的一則趣聞：

2021 年 6 月，在陜西西安，一位新郎接親時(shí)，被新娘要求當(dāng)場(chǎng)寫(xiě)代碼編程。新郎忙得滿(mǎn)頭大汗，終于在5分鐘內(nèi)完成考驗(yàn)，在電腦上做出一顆紅彤彤的愛(ài)心。新娘稱(chēng)，新郎是清華大學(xué)計(jì)算機(jī)專(zhuān)業(yè)博士，出題是想提醒他不要忘記吃飯的手藝，相信他能完成。

圖片來(lái)源：梨視頻

總所周知

，一個(gè)完美愛(ài)心

的曲面方程滿(mǎn)足下式[1]：

所以，不論是采用何種編程語(yǔ)言進(jìn)行繪圖，均是基于該方程。

為了防止以后有哪位男士娶親時(shí)又遇到該要求而著急，我覺(jué)得我得做些什么

。我這個(gè)人比較實(shí)在，我總結(jié)了 MATLAB、Python 以及 Mathematica 版本的代碼，各位有需要的可以收藏一下，以備不時(shí)之需

。

MATLAB 版本

紅色

clc;clear;close all[x,y,z]=meshgrid(linspace(-2,2));val=(x.^2+(9/4)*y.^2+z.^2-1).^3-x.^2.*z.^3-(9/80).*y.^2.*z.^3;isosurface(x,y,z,val,0)p = patch(isosurface(x,y,z,val,0));set(p,'facecolor','red','edgecolor','none');view(3);axis equal;title('\fontsize{14}\fontname{songti}你愛(ài)我','color','r')set(gcf,'unit','centimeters','position',[18513.539.03],'color','white');grid on

效果圖

：

MATLAB 版本三維心型圖（紅色）

五彩斑斕

clc;clear;close allfiguref=@(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);axis equal;title('\fontsize{14}\fontname{songti}我愛(ài)你','color','r')set(gcf,'unit','centimeters','position',[18513.539.03],'color','white');grid on

效果圖

：

MATLAB 版本三維心型圖（五彩斑斕）

Python 版本

import numpy as npfrom skimage import measureimport matplotlib.pyplot as pltplt.rcParams['font.sans-serif']=['SimHei']plt.rcParams['axes.unicode_minus'] = Falsedef fun(x, y, z):return (x**2+(9/4)*y**2+z**2-1)**3-x**2*z**3-(9/80)*y**2*z**3x, y, z = np.mgrid[-2:2:100j, -2:2:100j, -2:2:100j]vol = fun(x, y, z)iso_val=0.0verts, faces, _, _ = measure.marching_cubes(vol, iso_val, spacing=(0.1, 0.1, 0.1))fig = plt.figure()ax = fig.add_subplot(111, projection='3d')ax.plot_trisurf(verts[:, 0], verts[:,1], faces, verts[:, 2],color = 'red')ax.view_init(15, -120)plt.title('蜜雪冰城甜蜜蜜',color='red')plt.rcParams['figure.dpi'] = 1500plt.show()

效果圖

：

Python 版本三維心型圖

注意：在運(yùn)行 Python 代碼時(shí)，可能會(huì)因?yàn)閭€(gè)人軟件安裝的差異出現(xiàn)不同的錯(cuò)誤，這里給出常見(jiàn)的兩種。

ModuleNotFoundError: No module named 'skimage'遇到這種情況，需要在 cmd 中運(yùn)行：pip install scikit-image該命令是為了添加“skimage”包。ValueError: Unknown projection ‘3d’遇到這種情況，需要在原 Python 代碼的最前面添加一行代碼：from mpl_toolkits.mplot3d import Axes3D

Mathematica 版本

Mathematica 版本的代碼非常簡(jiǎn)潔：

ContourPlot3D[(x^2 + (9/4)*y^2 + z^2 - 1)^3 - x^2*z^3 - (9/80)*y^2*z^3 == 0, {x, -2, 2}, {y, -1, 1}, {z, 2, -2}, PlotPoints -> 100, BoxRatios -> Automatic, Mesh -> None, Axes -> None, Boxed -> False, ContourStyle -> {Red, Specularity[White, 5]}]

效果圖

：

Mathematica 版本三維心型圖

你學(xué)廢了嗎

？

參考資料

[1]Heart Surface: https://mathworld.wolfram.com/HeartSurface.html

拋磚引玉--芯片工程師版

以下的圖片是2019年EETOP網(wǎng)友通過(guò)仿真文件生成的：春節(jié)快樂(lè)，我覺(jué)得在此基礎(chǔ)上更改應(yīng)該也是可以做出自己獨(dú)特的心型圖片的。

方法如下：

通過(guò)Vivado直接板級(jí)調(diào)試時(shí)鐘用的PCIE IP輸出的:春節(jié)快樂(lè)時(shí)序抓取圖，純粹的硬邏輯產(chǎn)生。

實(shí)現(xiàn)原理：因?yàn)轫?xiàng)目需要，需要在板子上完成PCIE接口的通信，還要完成兩路AD，兩路DA的調(diào)試，還有四個(gè)顆粒的內(nèi)存顆粒，來(lái)完成信號(hào)的收發(fā)，F(xiàn)PGA用的是XC7V690T，還有時(shí)鐘芯片AD9512，ADI的這些芯片多是用SPI來(lái)控制寄存器，因此我在寫(xiě)寄存器的時(shí)候就想把寄存器的值換成節(jié)日祝福，應(yīng)該挺有新意的，就從新例化了一個(gè)module，每配置完一個(gè)寄存器自動(dòng)跳轉(zhuǎn)到下一個(gè)，直到所有寄存器配置完成結(jié)束，每一個(gè)寄存器位數(shù)為16位，每個(gè)漢字16個(gè)值，前后加0隔開(kāi)，一共74個(gè)，時(shí)鐘也可以根據(jù)需要分頻更改。

代碼：硬核實(shí)現(xiàn) ”春節(jié)快樂(lè)” 代碼分享