C語言如何實現繪制LoveBeat愛心曲線

這篇文章主要講解了“C語言如何實現繪制LoveBeat愛心曲線”，文中的講解內容簡單清晰，易于學習與理解，下面請大家跟著小編的思路慢慢深入，一起來研究和學習“C語言如何實現繪制LoveBeat愛心曲線”吧！

心形曲線

給出心形曲線參數方程如下：

x = sin^3(θ)

y = (13 * cos(θ) - 5 * cos(2θ) - 3 * cos(3θ) - cos(4θ)) / 16

對 x, y 同時乘以半徑 R，即可對其放大

通過上述方程可得到若干在曲線上的點，記為集合 S

曲線內點的生成

對于內部的點，我們則將S向內擴散，使其符合指數分布，得到S'

令 e ∈ [m, n] 且 e ~ E(λ)

對 P (x, y) ∈ S 作向內擴散得到點 P' ∈ S':

P' = (x, y) * e

擴散程度取決于參數 m, n, λ

曲線外點的生成

對于外部的點，我們則將S向外擴散，使其符合均勻分布，得到S''

令 u ~ U [1, 1 + b]

對 P (x, y) ∈ S 作向外擴散得到點 P'' ∈ S'':

P'' = (x, y) * u

擴散程度取決于參數 b

制作動畫

有了上述三個點集合后，就可以描述一個靜態的心形圖案

現在，我們為其添加周期動畫效果

首先引入時間參數 t 和周期函數 T(t) = sin^2(t) （可以是其他周期函數，視效果而定）

對于 P ∈ S U S'

其周期縮放程度與其到原點的距離 d 成反比，可用 R/d 來衡量 （R為心形曲線的半徑）

為其添加階數 i 來擴大距離對縮放程度的影響 (R/d)^i

我們可以得到如下函數

d' = d * (1 + a * T(t) * (R/d)^i)

外部的點在內部點收縮到最小值時，到達最大值，所以和上式相差一個相位

我們可以還做一些額外的處理：

• 對內部點和外部點添加隨機擾動

• 繪制時，隨機點的大小

• 繪制時，隨機點的顏色、亮度

示例代碼

// 環境：Visual Studio 2022 C++ 20
// EasyX版本：EasyX 2022-9-1
 
#include <random>
#include <unordered_set>
#include <graphics.h>
#include <cmath>
 
#define PINK LIGHTRED | 0x6055ff
#define LIGHTPINK LIGHTRED | 0x6655ff
#define DRAW(vecs, color) for(auto& vec : vecs)Draw(vec, color, distribution(engine))
 
using namespace std;
 
constexpr float Pi = 3.1416f;
constexpr float Rad = Pi / 180;
constexpr int ScreenWidth = 800;
constexpr int ScreenHeight = 600;
constexpr int OX = ScreenWidth / 2;
constexpr int OY = ScreenHeight / 2;
 
static default_random_engine engine;
 
struct Vec2
{
	float X = .0f;
	float Y = .0f;
 
	Vec2() {}
	Vec2(float x, float y) { X = x; Y = y; }
 
	int GetX() const { return static_cast<int>(X); }
	int GetY() const { return static_cast<int>(Y); }
	
	float Magnitude() const { return sqrt(X * X + Y * Y); }
	Vec2 operator*(float num) const{ return Vec2(X * num, Y * num); }
 
	struct VecHash
	{
		size_t operator()(const Vec2& vec) const noexcept
		{
			return std::hash<float>()(vec.X) ^ std::hash<float>()(vec.Y);
		}
	};
 
	struct VecCompare
	{
		bool operator()(const Vec2& vec1, const Vec2& vec2) const noexcept
		{
			return fabsf(vec1.X - vec2.X) < 1e-2f && fabsf(vec1.Y - vec2.Y) < 1e-2f;
		}
	};
};
using VecSet = unordered_set<Vec2, Vec2::VecHash, Vec2::VecCompare>;
 
 
float CalculateX(float t){ return powf(sin(t), 3.0f); }
float CalculateY(float t){ return -(13 * cosf(t) - 5 * cosf(2 * t) - 2 * cosf(3 * t) - cosf(4 * t)) / 16.0f; }
 
 
VecSet InitHeart(float startAngle, float endAngle, float radius, size_t count)
{
	VecSet set;
	float rad = startAngle * Rad;
	float step = (endAngle - startAngle) * Rad / count;
	float endRad = endAngle * Rad;
 
	for (; rad < endRad; rad += step)
        set.insert(Vec2(CalculateX(rad) * radius, CalculateY(rad) * radius));
 
	return set;
}
 
 
VecSet BuildInside(const VecSet& set, size_t frequency, float lambda, float range, float min)
{
	VecSet retSet;
	exponential_distribution<float> distribution(lambda);
 
	for (size_t i = 0; i < frequency; i++)
	{
		for (auto& vec : set)
		{
			float pX = distribution(engine);
			float scalarX = (pX < 1.0 ? 1.0f - pX : 1.0f) * range + min;
 
			float pY = distribution(engine);
			float scalarY = (pY < 1.0 ? 1.0f - pY : 1.0f) * range + min;
 
			retSet.insert(Vec2(vec.X * scalarX, vec.Y * scalarY));
		}
	}
 
	return retSet;
}
 
 
VecSet BuildOutside(const VecSet& set, size_t frequency, float max)
{
	VecSet retSet;
	uniform_real_distribution<float> distribution(1.0f, max);
 
	for (size_t i = 0; i < frequency; i++)
	{
		for (auto& vec : set)
            retSet.insert(Vec2(vec.X * distribution(engine), vec.Y *  distribution(engine)));
	}
 
	return retSet;
}
 
 
VecSet Undulate(const VecSet& set, float radius)
{
	VecSet retSet;
	uniform_real_distribution<float> distribution(-radius, radius);
 
	for (auto& vec : set)
		retSet.insert(Vec2(vec.X + distribution(engine), vec.Y + distribution(engine)));
 
	return retSet;
}
 
 
VecSet Zoom(const VecSet& set, float factor, float radius, float t, float idx)
{
	VecSet retSet;
 
	for (auto& vec : set)
        retSet.insert(vec * (1.0f + factor * sin(t * Pi) * powf((radius / vec.Magnitude()), idx)));
 
	return retSet;
}
 
 
void Draw(const Vec2& vec, COLORREF color, int radius)
{
	putpixel(vec.GetX() + OX, vec.GetY() + OY, color);
 
	if(radius >= 2)
		putpixel(vec.GetX() + OX + 1, vec.GetY() + OY, color);
 
	if(radius >= 3)
		putpixel(vec.GetX() + OX, vec.GetY() + OY + 1, color);
}
	
 
int main()
{
	float radius = 160.0f;
	auto border = InitHeart(0, 360, radius, 480);
	auto inside = BuildInside(border, 30, 5.0f, 0.85f, 0.15f);
 
	initgraph(ScreenWidth, ScreenHeight);
	BeginBatchDraw();
 
	float t = .0f;
	float tStep = 0.05f;
	uniform_int_distribution<int> distribution(1, 3);
	ExMessage msg{};
 
	while(!peekmessage(&msg, EX_KEY))
	{
		auto ps1 = Zoom(border, 0.1f, radius, t, 1.3f);
		auto ps2 = Undulate(Zoom(inside, 0.1f, radius, t, 1.3f), 3.0f);
		auto ps3 = Undulate(BuildOutside(border, 10 - static_cast<size_t>(sin(t) * 5), 1.35f - sin(t) * 0.15f), 3.0f);
 
		DRAW(ps1, LIGHTPINK);
		DRAW(ps2, LIGHTPINK);
		DRAW(ps3, PINK);
 
		FlushBatchDraw();
		Sleep(40);
 
		t += tStep;
		if (t > 1.0f)
			t = .0f;
 
		cleardevice();
	}
 
	EndBatchDraw();
	closegraph();
 
	return 0;
}

感謝各位的閱讀，以上就是“C語言如何實現繪制LoveBeat愛心曲線”的內容了，經過本文的學習后，相信大家對C語言如何實現繪制LoveBeat愛心曲線這一問題有了更深刻的體會，具體使用情況還需要大家實踐驗證。這里是億速云，小編將為大家推送更多相關知識點的文章，歡迎關注！