![\"\"](\"https:\/\/123210.net\/wordpress\/tatenaoshi\/wp-content\/uploads\/sites\/3\/2026\/07\/\u30ea\u30f3\u30af\u30a8\u30e9\u30fc\u306e\u30e1\u30c3\u30bb\u30fc\u30b8\u3068\u5b9f\u884c\u72b6\u6cc12.jpg\")
<\/a><\/p>\n\n
C++Builder2007\u30bd\u30fc\u30b9\u30b3\u30fc\u30c9<\/p>\n
<\/p>\n
\/\/*** Unit1.h
\n\/\/—————————————————————————<\/p>\n
#ifndef Unit1H
\n#define Unit1H
\n\/\/—————————————————————————
\n#include <Classes.hpp>
\n#include <Controls.hpp>
\n#include <StdCtrls.hpp>
\n#include <Forms.hpp>
\n#include “Chart.hpp”
\n#include “TeEngine.hpp”
\n#include “TeeProcs.hpp”
\n#include <ExtCtrls.hpp><\/p>\n
#include <Series.hpp> \/\/ \u2190 \u3053\u308c\u3092\u8ffd\u52a0\uff08\u6700\u91cd\u8981\uff09
\n\/\/—————————————————————————<\/p>\n
\/\/————————————————
\n\/\/ \u5c65\u6b74\u7528\u69cb\u9020\u4f53
\n\/\/————————————————
\nstruct OctoState {
\nint step;
\nint depth;
\ndouble c[8];
\ndouble norm2;
\ndouble N;
\nint lastDeriv;
\n};<\/p>\n
\/\/————————————————
\n\/\/ \u516b\u5143\u6570
\n\/\/————————————————
\nstruct Octonion {
\ndouble e[8];<\/p>\n
Octonion(double a0=0.0, double a1=0.0, double a2=0.0, double a3=0.0,
\ndouble a4=0.0, double a5=0.0, double a6=0.0, double a7=0.0);<\/p>\n
Octonion operator+(const Octonion& o) const;
\nOctonion operator-(const Octonion& o) const; \/\/ \u6e1b\u7b97
\nOctonion operator-() const; \/\/ \u5358\u9805\u30de\u30a4\u30ca\u30b9\uff08\u5ff5\u306e\u305f\u3081\uff09
\nOctonion operator*(const Octonion& o) const; \/\/ \u2190 \u5909\u66f4\uff1a\u5b8c\u5168\u4e57\u7b97
\nOctonion operator*(double scalar) const; \/\/ \u30b9\u30ab\u30e9\u30fc\u4e57\u7b97\u306f\u6b8b\u3059
\ndouble Norm() const;<\/p>\n
Octonion Conjugate() const; \/\/ \u6a19\u6e96\u7684\u306a\u5171\u5f79
\nOctonion PhaseConjugate() const; \/\/ \u4f4d\u76f8\u5171\u5f79\u7528\uff08\u5fc5\u8981\u306b\u5fdc\u3058\u3066\u62e1\u5f35\u53ef\u80fd\uff09
\nOctonion ApplyDerivation(const Octonion& a, const Octonion& b) const; \/\/ D_a,b (this)
\n};<\/p>\n
\/\/————————————————
\n\/\/ \u518d\u5e30\u516b\u5143\u6570
\n\/\/————————————————
\nclass RecursiveOcto {
\npublic:
\nint depth;
\nOctonion value;
\nRecursiveOcto* child; \/\/ \u901a\u5e38\u306e\u5b50
\nRecursiveOcto* childConj; \/\/ \u5171\u5f79\u306e\u5b50\uff08\u8ffd\u52a0\uff09
\nOctonion interference; \/\/ \u5e72\u6e09\u7d50\u679c\u3092\u4fdd\u6301\u3059\u308b\u30e1\u30f3\u30d0\u30fc\u3092\u8ffd\u52a0<\/p>\n
\/\/ \u2605 \u8ffd\u52a0
\ndouble N; \/\/ \u6570\u6f14\u7b97\u5b50\uff08\u52b1\u8d77\u5ea6\uff09
\ndouble Q; \/\/ \u96fb\u8377\u7684\u306a\u91cf = N \/ 3<\/p>\n
RecursiveOcto(int d, Octonion init);
\n~RecursiveOcto();
\nvoid Generate(int maxDepth);
\nvoid Dump(TMemo* memo);<\/p>\n
void ComputeNumberOperator(); \/\/ \u2605 \u8ffd\u52a0<\/p>\n
\/\/ \u89b3\u5bdf\uff08Observe\uff09\u6a5f\u80fd\u3068\u306e\u9023\u643a
\nint preferredDir; \/\/ \u89b3\u5bdf\u3067\u56fa\u5b9a\u3055\u308c\u305f\u65b9\u5411\uff08-1 = \u672a\u89b3\u5bdf\uff09
\nvoid Observe(int dir); \/\/ dir = 0?6\uff08e1?e7\uff09
\n};<\/p>\n
\/\/————————————————
\n\/\/ \u30d5\u30a9\u30fc\u30e0
\n\/\/————————————————
\nclass TForm1 : public TForm
\n{
\n__published: \/\/ IDE-managed Components
\nTButton *ButtonExit;
\nTMemo *Memo1;
\nTButton *Button2;
\nTChart *ChartTime;
\nTChart *ChartNorm;
\nTChart *Chart3D;
\nvoid __fastcall ButtonExitClick(TObject *Sender);
\nvoid __fastcall Button2Click(TObject *Sender);
\nvoid __fastcall FormCreate(TObject *Sender);
\nprivate: \/\/ User declarations
\nTList* History;<\/p>\n
void ClearHistory();
\nvoid CollectHistory(RecursiveOcto* node, int& counter);
\nvoid __fastcall InitTimeCharts();
\nvoid Init3DChart();
\nvoid __fastcall UpdateTimeCharts();
\nvoid __fastcall Update3DChart();
\npublic: \/\/ User declarations
\n__fastcall TForm1(TComponent* Owner);
\n};
\n\/\/—————————————————————————
\nextern PACKAGE TForm1 *Form1;
\n\/\/—————————————————————————
\n#endif<\/p>\n
\/\/*** Unit1.cpp
\n\/\/—————————————————————————<\/p>\n
#include <vcl.h>
\n#pragma hdrstop<\/p>\n
#include <math.h>
\n#include “Unit1.h”
\n\/\/—————————————————————————
\n#pragma package(smart_init)
\n#pragma link “Chart”
\n#pragma link “TeEngine”
\n#pragma link “TeeProcs”
\n#pragma resource “*.dfm”
\nTForm1 *Form1;<\/p>\n
\/\/\u30ea\u30f3\u30ab\uff5c\u51fa\u529b\u30aa\u30d7\u30b7\u30e7\u30f3\uff5c\u6700\u5927\u30b9\u30bf\u30c3\u30af\u30b5\u30a4\u30ba\u30000x64000000\uff08\u7d041.56GB\uff09\u5927\u304d\u3059\u304e\u3067\u3042\u3063\u305f\u306e\u3067\u30010x01000000\u306b\u8a2d\u5b9a
\n\/\/—————————————————————————
\n__fastcall TForm1::TForm1(TComponent* Owner)
\n: TForm(Owner)
\n{
\n}
\n\/\/—————————————————————————
\nvoid __fastcall TForm1::ButtonExitClick(TObject *Sender)
\n{
\nClose();
\n}
\n\/\/—————————————————————————
\n\/\/—————————————————————————<\/p>\n
Octonion::Octonion(double a0, double a1, double a2, double a3,
\ndouble a4, double a5, double a6, double a7) {
\ne[0]=a0; e[1]=a1; e[2]=a2; e[3]=a3;
\ne[4]=a4; e[5]=a5; e[6]=a6; e[7]=a7;
\n}<\/p>\n
Octonion Octonion::Conjugate() const {
\nreturn Octonion(
\ne[0],
\n-e[1], -e[2], -e[3], -e[4], -e[5], -e[6], -e[7]
\n);
\n}<\/p>\n
Octonion Octonion::PhaseConjugate() const {
\nreturn Conjugate();
\n}<\/p>\n
Octonion Octonion::operator+(const Octonion& o) const {
\nOctonion res;
\nfor(int i=0; i<8; i++) res.e[i] = e[i] + o.e[i];
\nreturn res;
\n}<\/p>\n
\/\/ \u4e8c\u9805\u6e1b\u7b97
\nOctonion Octonion::operator-(const Octonion& o) const {
\nOctonion res;
\nfor(int i = 0; i < 8; i++) {
\nres.e[i] = e[i] – o.e[i];
\n}
\nreturn res;
\n}<\/p>\n
\/\/ \u5358\u9805\u30de\u30a4\u30ca\u30b9
\nOctonion Octonion::operator-() const {
\nOctonion res;
\nfor(int i = 0; i < 8; i++) {
\nres.e[i] = -e[i];
\n}
\nreturn res;
\n}<\/p>\n
\/\/ \u30b9\u30ab\u30e9\u30fc\u4e57\u7b97\u306f\u6b8b\u3059
\nOctonion Octonion::operator*(double scalar) const {
\nOctonion res;
\nfor(int i=0; i<8; i++) res.e[i] = e[i] * scalar;
\nreturn res;
\n}<\/p>\n
\/\/ \u5b8c\u5168\u306a\u516b\u5143\u6570\u4e57\u7b97
\nOctonion Octonion::operator*(const Octonion& o) const {
\nOctonion res;<\/p>\n
\/\/ e0\uff08\u5b9f\u90e8\uff09
\nres.e[0] = e[0]*o.e[0] – e[1]*o.e[1] – e[2]*o.e[2] – e[3]*o.e[3]
\n– e[4]*o.e[4] – e[5]*o.e[5] – e[6]*o.e[6] – e[7]*o.e[7];<\/p>\n
\/\/ e1
\nres.e[1] = e[0]*o.e[1] + e[1]*o.e[0] + e[2]*o.e[3] – e[3]*o.e[2]
\n+ e[4]*o.e[5] – e[5]*o.e[4] – e[6]*o.e[7] + e[7]*o.e[6];<\/p>\n
\/\/ e2
\nres.e[2] = e[0]*o.e[2] – e[1]*o.e[3] + e[2]*o.e[0] + e[3]*o.e[1]
\n+ e[4]*o.e[6] + e[5]*o.e[7] – e[6]*o.e[4] – e[7]*o.e[5];<\/p>\n
\/\/ e3
\nres.e[3] = e[0]*o.e[3] + e[1]*o.e[2] – e[2]*o.e[1] + e[3]*o.e[0]
\n+ e[4]*o.e[7] – e[5]*o.e[6] + e[6]*o.e[5] – e[7]*o.e[4];<\/p>\n
\/\/ e4
\nres.e[4] = e[0]*o.e[4] – e[1]*o.e[5] – e[2]*o.e[6] – e[3]*o.e[7]
\n+ e[4]*o.e[0] + e[5]*o.e[1] + e[6]*o.e[2] + e[7]*o.e[3];<\/p>\n
\/\/ e5
\nres.e[5] = e[0]*o.e[5] + e[1]*o.e[4] – e[2]*o.e[7] + e[3]*o.e[6]
\n– e[4]*o.e[1] + e[5]*o.e[0] – e[6]*o.e[3] + e[7]*o.e[2];<\/p>\n
\/\/ e6
\nres.e[6] = e[0]*o.e[6] + e[1]*o.e[7] + e[2]*o.e[4] – e[3]*o.e[5]
\n– e[4]*o.e[2] + e[5]*o.e[3] + e[6]*o.e[0] – e[7]*o.e[1];<\/p>\n
\/\/ e7
\nres.e[7] = e[0]*o.e[7] – e[1]*o.e[6] + e[2]*o.e[5] + e[3]*o.e[4]
\n– e[4]*o.e[3] – e[5]*o.e[2] + e[6]*o.e[1] + e[7]*o.e[0];<\/p>\n
return res;
\n}<\/p>\n
double Octonion::Norm() const {
\ndouble sum = 0.0;
\nfor(int i=0; i<8; i++) sum += e[i]*e[i];
\nreturn sqrt(sum);
\n}<\/p>\n
Octonion Octonion::ApplyDerivation(const Octonion& a, const Octonion& b) const
\n{
\n\/\/ 1. \u4ea4\u63db\u5b50 [a,b]
\nOctonion ab = a * b – b * a; \/\/ [a,b]<\/p>\n
\/\/ 2. \u4ea4\u63db\u5b50\u306e\u4f5c\u7528 [[a,b], this]
\nOctonion commutator_term = ab * (*this) – (*this) * ab;<\/p>\n
\/\/ 3. \u7d50\u5408\u5b50 [a,b,this] = (a*b)*this – a*(b*this)
\nOctonion assoc = (a * b) * (*this) – a * (b * (*this));<\/p>\n
\/\/ 4. \u6a19\u6e96\u7684\u306a\u5c0e\u5206\u516c\u5f0f\uff08\u4fc2\u6570\u306f\u8abf\u6574\u53ef\u80fd\uff09
\n\/\/ D = [[a,b],x] – 3[a,b,x]
\nOctonion result = commutator_term – assoc * 3.0;<\/p>\n
\/\/ \u30b9\u30b1\u30fc\u30eb\u8abf\u6574\uff08\u6442\u52d5\u304c\u5927\u304d\u3059\u304e\u306a\u3044\u3088\u3046\u306b\uff09
\n\/\/return result * 0.5;
\n\/\/return result * 0.20; \/\/ \u6442\u52d5\u304c\u5c11\u3057\u5f37\u3081\u306b\u51fa\u3066\u3044\u308b\u3088\u3046\u306a\u306e\u3067\u3001\u30b9\u30b1\u30fc\u30eb\u3092\u5c11\u3057\u5f31\u3081\u30660.5\u21920.20\uff0860%\u6e1b\uff09\u5b89\u5b9a\u6027\u3092\u78ba\u8a8d\u3059\u308b
\n\/\/return result * 0.25; \/\/ \u6442\u52d5\u304c\u5c11\u3057\u5f37\u3081\u306b\u51fa\u3066\u3044\u308b\u3088\u3046\u306a\u306e\u3067\u3001\u30b9\u30b1\u30fc\u30eb\u3092\u5c11\u3057\u5f31\u3081\u30660.5\u21920.25\uff0850%\u6e1b\uff09\u5b89\u5b9a\u6027\u3092\u78ba\u8a8d\u3059\u308b
\n\/\/return result * 0.30; \/\/ \u6442\u52d5\u304c\u5c11\u3057\u5f37\u3081\u306b\u51fa\u3066\u3044\u308b\u3088\u3046\u306a\u306e\u3067\u3001\u30b9\u30b1\u30fc\u30eb\u3092\u5c11\u3057\u5f31\u3081\u30660.5\u21920.30\uff0840%\u6e1b\uff09\u5b89\u5b9a\u6027\u3092\u78ba\u8a8d\u3059\u308b
\nreturn result * 0.35; \/\/ \u6442\u52d5\u304c\u5c11\u3057\u5f37\u3081\u306b\u51fa\u3066\u3044\u308b\u3088\u3046\u306a\u306e\u3067\u3001\u30b9\u30b1\u30fc\u30eb\u3092\u5c11\u3057\u5f31\u3081\u30660.5\u21920.35\uff0830%\u6e1b\uff09\u5b89\u5b9a\u6027\u3092\u78ba\u8a8d\u3059\u308b
\n\/\/return result * 0.40; \/\/ \u6442\u52d5\u304c\u5c11\u3057\u5f37\u3081\u306b\u51fa\u3066\u3044\u308b\u3088\u3046\u306a\u306e\u3067\u3001\u30b9\u30b1\u30fc\u30eb\u3092\u5c11\u3057\u5f31\u3081\u30660.5\u21920.40\uff0820%\u6e1b\uff09\u5b89\u5b9a\u6027\u3092\u78ba\u8a8d\u3059\u308b
\n\/\/return result * 0.45; \/\/ \u6442\u52d5\u304c\u5c11\u3057\u5f37\u3081\u306b\u51fa\u3066\u3044\u308b\u3088\u3046\u306a\u306e\u3067\u3001\u30b9\u30b1\u30fc\u30eb\u3092\u5c11\u3057\u5f31\u3081\u30660.5\u21920.45\uff0810%\u6e1b\uff09\u5b89\u5b9a\u6027\u3092\u78ba\u8a8d\u3059\u308b
\n\/\/*** \u8a55\u4fa1\u7d50\u679c\uff1a0.30\uff1a\u63a7\u3048\u3081\u3060\u304c\u52b9\u679c\u306f\u5341\u5206\u51fa\u308b ***
\n\/\/*** \u8a55\u4fa1\u7d50\u679c\uff1a0.35\uff1a\u5c0e\u5206\u306e\u5f71\u97ff\u304c\u3088\u308a\u660e\u78ba\u306b\u73fe\u308c\u3064\u3064\u3001\u5b89\u5b9a\u6027\u3082\u9ad8\u3044 ***<\/p>\n
}<\/p>\n
\/\/—————————————————————————
\nRecursiveOcto::RecursiveOcto(int d, Octonion v)
\n: depth(d), value(v), child(NULL), childConj(NULL),
\nN(0.0), Q(0.0), preferredDir(-1)
\n{
\n}<\/p>\n
RecursiveOcto::~RecursiveOcto() {
\ndelete child;
\ndelete childConj;
\n}<\/p>\n
void RecursiveOcto::Generate(int maxDepth) {
\n\/\/if (depth >= maxDepth) return;
\n\/\/if (depth >= maxDepth || depth > 20) return; \/\/ \u5b89\u5168\u4e0a\u9650\u3092\u8ffd\u52a0<\/p>\n
\/\/ \u7d76\u5bfe\u306b\u6df1\u304f\u306a\u3089\u306a\u3044\u3088\u3046\u306b\u3059\u308b
\nif (depth >= maxDepth || depth > 12) {
\nchild = NULL;
\nchildConj = NULL;
\nreturn;
\n}<\/p>\n
\/\/ === \u65e2\u5b58\u306e\u56de\u8ee2\u5b50 ===
\nOctonion rotator(cos(depth*0.15), sin(depth*0.15), 0.05, 0,0,0,0,0);<\/p>\n
\/\/ === Fano\u5e73\u9762\u306e\u65b9\u5411\u5b9a\u7fa9 ===
\nOctonion dirs[7][2] = {
\n{ Octonion(0,1,0,0,0,0,0,0), Octonion(0,0,1,0,0,0,0,0) }, \/\/ 0
\n{ Octonion(0,1,0,0,0,0,0,0), Octonion(0,0,0,0,1,0,0,0) }, \/\/ 1
\n{ Octonion(0,1,0,0,0,0,0,0), Octonion(0,0,0,0,0,0,1,0) }, \/\/ 2
\n{ Octonion(0,0,1,0,0,0,0,0), Octonion(0,0,0,0,1,0,0,0) }, \/\/ 3
\n{ Octonion(0,0,1,0,0,0,0,0), Octonion(0,0,0,0,0,1,0,0) }, \/\/ 4
\n{ Octonion(0,0,0,1,0,0,0,0), Octonion(0,0,0,0,1,0,0,0) }, \/\/ 5
\n{ Octonion(0,0,0,1,0,0,0,0), Octonion(0,0,0,0,0,1,0,0) } \/\/ 6
\n};<\/p>\n
int dir = depth % 7;<\/p>\n
\/\/ === \u7c21\u6613\u7684\u306a\u5c40\u6240N\u63a8\u5b9a\uff08Generate\u6642\u70b9\u3067\u4f7f\u3048\u308b\u3088\u3046\u306b\uff09 ===
\ndouble imagStrength = 0.0;
\nfor(int i=1; i<=7; i++) imagStrength += fabs(value.e[i]);<\/p>\n
\/\/ \u7c21\u6613N\uff08\u91cf\u5b50\u5316\u524d\uff09
\ndouble localN = 1.0 + imagStrength * 0.8; \/\/ \u30d9\u30fc\u30b9 + \u865a\u90e8\u306e\u5f37\u3055
\nif (localN > 3.0) localN = 3.0;<\/p>\n
\/\/ N\u306b\u5fdc\u3058\u305f\u30b9\u30b1\u30fc\u30eb\u4fc2\u6570\uff08N\u304c\u5927\u304d\u3044\u307b\u3069\u5f37\u304f\u3059\u308b\uff09
\n\/\/ N=0 \u2192 \u5f31\u3044, N=3 \u2192 \u5f37\u3044
\ndouble nFactor = 0.4 + (localN \/ 3.0) * 0.9; \/\/ 0.4 ? 1.3 \u7a0b\u5ea6\u306e\u7bc4\u56f2<\/p>\n
\/\/*******************************************************************
\n\/\/ Generate \u5185\u306e\u5c0e\u5206\u90e8\u5206\u3067\u3001\u89b3\u5bdf\u6e08\u307f\u306e\u5834\u5408\u306f\u56fa\u5b9a\u65b9\u5411\u3092\u907f\u3051\u308b\uff0f\u5f31\u3081\u308b
\nif (preferredDir >= 0) {
\n\/\/ \u89b3\u5bdf\u6e08\u307f\u306a\u3089\u3001\u56fa\u5b9a\u65b9\u5411\u306b\u95a2\u4fc2\u3059\u308b\u5c0e\u5206\u3092\u5f31\u3081\u308b
\n\/\/ \uff08\u7c21\u6613\u7684\u306b nFactor \u3092\u3055\u3089\u306b\u4e0b\u3052\u308b\uff09
\nnFactor *= 0.5;
\n}
\n\/\/*******************************************************************<\/p>\n
\/\/ === \u8907\u6570\u5c0e\u5206\uff08N\u4f9d\u5b58\u30b9\u30b1\u30fc\u30eb\u3092\u639b\u3051\u308b\uff09 ===
\nOctonion g2_perturb = value.ApplyDerivation(dirs[dir][0], dirs[dir][1]) * (0.07 * nFactor);<\/p>\n
int dir2 = (dir + 1) % 7;
\ng2_perturb = g2_perturb + value.ApplyDerivation(dirs[dir2][0], dirs[dir2][1]) * (0.03 * nFactor);<\/p>\n
int dir3 = (dir + 3) % 7;
\ng2_perturb = g2_perturb + value.ApplyDerivation(dirs[dir3][0], dirs[dir3][1]) * (0.015 * nFactor);<\/p>\n
\/\/ \u6b21\u306e\u72b6\u614b
\nOctonion next = (value + g2_perturb) * rotator;
\nnext = next + Octonion(0.03, 0,0,0,0,0,0,0);<\/p>\n
\/\/ \u5171\u5f79\u7248
\nOctonion nextConj = next.PhaseConjugate();<\/p>\n
\/\/ \u5b50\u751f\u6210 (\uff12\u500b\uff1a\u5171\u5f79\u5bfe\u3092\u751f\u6210)
\nchild = new RecursiveOcto(depth+1, next);
\nchild->Generate(maxDepth);<\/p>\n
childConj = new RecursiveOcto(depth+1, nextConj);
\nchildConj->Generate(maxDepth);<\/p>\n
\/\/ \u5e72\u6e09 + \u30d5\u30a3\u30fc\u30c9\u30d0\u30c3\u30af
\ninterference = (child->value + childConj->value) * 0.5;
\ndouble feedbackStrength = 0.3;
\nvalue = value * (1.0 – feedbackStrength) + interference * feedbackStrength;
\n}<\/p>\n
\/\/ \u6570\u6f14\u7b97\u5b50\u3092\u8a08\u7b97\u3059\u308b\u30e1\u30bd\u30c3\u30c9\uff08\u65b0\u898f\u8ffd\u52a0\uff09
\nvoid RecursiveOcto::ComputeNumberOperator()
\n{
\nif (depth > 12) return; \/\/ \u8ffd\u52a0<\/p>\n
\/\/ \u57fa\u672c\u52b1\u8d77\uff1a\u81ea\u5206\u81ea\u8eab\u306e\u865a\u6570\u6210\u5206\u306e\u300c\u5f37\u3055\u300d
\n\/\/ \uff08e1?e7 \u306e\u7d76\u5bfe\u5024\u306e\u5408\u8a08\u3092\u9069\u5ea6\u306b\u30b9\u30b1\u30fc\u30eb\uff09
\ndouble imagStrength = 0.0;<\/p>\n
\/\/ \u2605 \u3053\u3053\u3092 e[1] \uff5e e[7] \u306b\u4fee\u6b63
\nimagStrength += fabs(value.e[1]);
\nimagStrength += fabs(value.e[2]);
\nimagStrength += fabs(value.e[3]);
\nimagStrength += fabs(value.e[4]);
\nimagStrength += fabs(value.e[5]);
\nimagStrength += fabs(value.e[6]);
\nimagStrength += fabs(value.e[7]);<\/p>\n
\/\/ \u5b50\u4f9b\u306e\u6709\u7121\u306b\u3088\u308b\u52b1\u8d77
\nint childCount = 0;
\nif (child) childCount++;
\nif (childConj) childCount++;<\/p>\n
\/\/ \u518d\u5e30\u7684\u306b\u5b50\u4f9b\u306e N \u3092\u52a0\u7b97\uff08\u968e\u5c64\u5168\u4f53\u306e\u52b1\u8d77\u3092\u84c4\u7a4d\uff09 \/\/ \u5b50\u4f9b\u306e N \u3092\u518d\u5e30\u7684\u306b\u96c6\u3081\u308b
\ndouble childrenN = 0.0;
\nif (child) {
\nchild->ComputeNumberOperator();
\nchildrenN += child->N;
\n}
\nif (childConj) {
\nchildConj->ComputeNumberOperator();
\nchildrenN += childConj->N;
\n}<\/p>\n
\/\/ \u6700\u7d42\u7684\u306a N \u306e\u5b9a\u7fa9\uff08\u8abf\u6574\u3057\u3084\u3059\u3044\u5f62\uff09
\n\/\/ \u5b50\u4f9b\u306e\u6570 \u00d7 1.0 + \u865a\u6570\u5f37\u5ea6 \u00d7 \u4fc2\u6570 + \u5b50\u4f9b\u305f\u3061\u306e N \u306e\u4e00\u90e8
\n\/\/N = childCount * 1.0 + imagStrength * 0.8 + childrenN * 0.3;<\/p>\n
\/\/ \u2605 \u4fc2\u6570\u3092\u8abf\u6574\u3057\u305f\u30d0\u30fc\u30b8\u30e7\u30f3
\n\/\/ \u5b50\u4f9b\u306e\u5bc4\u4e0e\u3092\u5f37\u3081\u3001\u865a\u6570\u5f37\u5ea6\u306e\u5f71\u97ff\u3092\u6291\u3048\u3081\u306b\u3059\u308b
\n\/*N = childCount * 1.5
\n+ imagStrength * 0.4
\n+ childrenN * 0.25;*\/
\n\/\/ \u2605 \u3053\u3053\u3092\u5927\u304d\u304f\u5909\u66f4
\nN = childCount * 1.0 \/\/ \u5b50\u4f9b\u306e\u6570\u306e\u5f71\u97ff\u3092\u5f31\u3081\u308b
\n+ imagStrength * 0.6
\n+ childrenN * 0.05; \/\/ \u5b50\u5b6b\u304b\u3089\u306e\u84c4\u7a4d\u3092\u304b\u306a\u308a\u5f31\u304f\u3059\u308b<\/p>\n
\/\/ \u3055\u3089\u306b\u3001\u3042\u308b\u7a0b\u5ea6\u6574\u6570\u306b\u8fd1\u3065\u3051\u308b\u305f\u3081\u306e\u4e38\u3081\uff08\u4efb\u610f\uff09
\n\/\/ N = floor(N + 0.5); \/\/ \u6574\u6570\u306b\u3057\u305f\u3044\u5834\u5408\u306f\u3053\u306e\u884c\u3092\u6709\u52b9\u306b\u3059\u308b
\n\/\/ \u2605 \u3053\u3053\u304b\u3089\u91cf\u5b50\u5316\uff08\u96e2\u6563\u5316\uff09
\n\/*if (N < 0.7) N = 0.0;
\nelse if (N < 1.7) N = 1.0;
\nelse if (N < 2.7) N = 2.0;
\nelse if (N < 3.7) N = 3.0;
\nelse N = 3.0; \/\/ \u4e0a\u9650\u30923\u306b\u5236\u9650
\n*\/
\n\/\/ \u2605 \u95be\u5024\u3092\u7de9\u3081\u305f\u91cf\u5b50\u5316
\n\/*if (N < 0.8) N = 0.0;
\nelse if (N < 1.8) N = 1.0;
\nelse if (N < 2.8) N = 2.0;
\nelse N = 3.0;*\/
\n\/\/ \u91cf\u5b50\u5316
\nif (N < 0.7) N = 0.0;
\nelse if (N < 1.5) N = 1.0;
\nelse if (N < 2.3) N = 2.0;
\nelse N = 3.0;<\/p>\n
\/\/ \u96fb\u8377\u7684\u306a\u91cf
\nQ = N \/ 3.0;
\n}<\/p>\n
void RecursiveOcto::Observe(int dir)
\n{
\nif (depth > 12) return; \/\/ \u8ffd\u52a0<\/p>\n
if (dir < 0 || dir > 6) return;<\/p>\n
preferredDir = dir;<\/p>\n
\/\/ \u6307\u5b9a\u3057\u305f\u865a\u5358\u4f4d\u306e\u6210\u5206\u3092\u5f37\u8abf\u3057\u3001\u4ed6\u3092\u76f8\u5bfe\u7684\u306b\u6291\u5236\u3059\u308b\u7c21\u6613\u6295\u5f71
\n\/\/ \uff08\u5b8c\u5168\u306a\u5c04\u5f71\u3067\u306f\u306a\u304f\u3001\u89b3\u5bdf\u306b\u3088\u308b\u300c\u56fa\u5b9a\u300d\u306e\u30a4\u30e1\u30fc\u30b8\uff09
\ndouble strength = value.e[dir + 1]; \/\/ e1?e7<\/p>\n
\/\/ \u6307\u5b9a\u65b9\u5411\u3092\u6b8b\u3057\u3064\u3064\u3001\u5168\u4f53\u3092\u5c11\u3057\u6b63\u898f\u5316
\nfor(int i = 1; i <= 7; i++) {
\nif (i == dir + 1) {
\n\/\/ \u6307\u5b9a\u65b9\u5411\u306f\u305d\u306e\u307e\u307e\uff08\u307e\u305f\u306f\u5c11\u3057\u5f37\u8abf\uff09
\nvalue.e[i] = strength * 1.1;
\n} else {
\n\/\/ \u4ed6\u306e\u65b9\u5411\u3092\u5f31\u3081\u308b
\nvalue.e[i] *= 0.6;
\n}
\n}<\/p>\n
\/\/ \u5b50\u306b\u3082\u540c\u3058\u89b3\u5bdf\u3092\u4f1d\u64ad\uff08\u518d\u5e30\uff09
\nif (child) child->Observe(dir);
\nif (childConj) childConj->Observe(dir);
\n}<\/p>\n
void RecursiveOcto::Dump(TMemo* memo)
\n{<\/p>\n
if (!memo) return;
\nAnsiString s;
\ns.sprintf(“Depth=%d Norm=%.4f e0=%.4f e1=%.4f N=%.1f Q=%.3f”,
\ndepth, value.Norm(), value.e[0], value.e[1], N, Q);
\nmemo->Lines->Add(s);<\/p>\n
\/\/ \u30aa\u30d7\u30b7\u30e7\u30f3\uff1a\u865a\u90e8\u306e\u5927\u304d\u3055\u3092\u5c11\u3057\u8a73\u3057\u304f
\ndouble imag = 0;
\nfor(int i=1;i<8;i++) imag += fabs(value.e[i]);
\ns.sprintf(” imagStrength=%.4f”, imag);
\nmemo->Lines->Add(s);<\/p>\n
if (child) { memo->Lines->Add(” [\u901a\u5e38\u306e\u5b50]”); child->Dump(memo); }
\nif (childConj) { memo->Lines->Add(” [\u5171\u5f79\u306e\u5b50]”); childConj->Dump(memo); }
\n}<\/p>\n
\/\/—————————————————————————
\n\/\/—————————————————————————<\/p>\n
void __fastcall TForm1::UpdateTimeCharts()
\n{<\/p>\n
if (!ChartTime || !ChartNorm) return;
\nif (ChartTime->SeriesCount() < 8) return;
\nif (ChartNorm->SeriesCount() < 2) return;<\/p>\n
\/\/ \u30af\u30ea\u30a2
\nfor (int i = 0; i < ChartTime->SeriesCount(); i++)
\nChartTime->Series[i]->Clear();
\nChartNorm->Series[0]->Clear();
\nChartNorm->Series[1]->Clear();<\/p>\n
for (int i = 0; i < History->Count; i++) {
\nOctoState* s = (OctoState*)History->Items[i];
\nfor (int k = 0; k < 8; k++)
\nChartTime->Series[k]->AddXY(s->step, s->c[k]);
\nChartNorm->Series[0]->AddXY(s->step, s->norm2);
\nChartNorm->Series[1]->AddXY(s->step, s->N);
\n}<\/p>\n
ChartTime->Refresh();
\nChartNorm->Refresh();
\n}<\/p>\n
void TForm1::ClearHistory()
\n{
\nfor (int i = 0; i < History->Count; i++)
\ndelete (OctoState*)History->Items[i];
\nHistory->Clear();
\n}<\/p>\n
void __fastcall TForm1::Update3DChart()
\n{<\/p>\n
if (!Chart3D) return;
\nif (Chart3D->SeriesCount() == 0) return;<\/p>\n
TPointSeries* ser = dynamic_cast<TPointSeries*>(Chart3D->Series[0]);
\nif (!ser) return;<\/p>\n
ser->Clear();<\/p>\n
for (int i = 0; i < History->Count; i++) {
\nOctoState* s = (OctoState*)History->Items[i];
\ndouble x = s->c[1];
\ndouble y = s->c[2];<\/p>\n
TColor col = clBlack;
\nif (s->N >= 2.5) col = clRed;
\nelse if (s->N >= 1.5) col = clLime;
\nelse if (s->N >= 0.5) col = clBlue;<\/p>\n
ser->AddXY(x, y, “”, col);
\n}<\/p>\n
Chart3D->Axes->Bottom->Automatic = true;
\nChart3D->Axes->Left->Automatic = true;
\nChart3D->Refresh();
\n}<\/p>\n
void TForm1::CollectHistory(RecursiveOcto* node, int& counter)
\n{
\nif (!node) return;<\/p>\n
OctoState* s = new OctoState;
\ns->step = counter++;
\ns->depth = node->depth;
\nfor (int i = 0; i < 8; i++) s->c[i] = node->value.e[i];
\ns->norm2 = node->value.Norm() * node->value.Norm();
\ns->N = node->N;
\ns->lastDeriv = node->depth % 7;<\/p>\n
History->Add(s);<\/p>\n
CollectHistory(node->child, counter);
\nCollectHistory(node->childConj, counter);<\/p>\n
<\/p>\n
}<\/p>\n
void __fastcall TForm1::Button2Click(TObject *Sender)
\n{
\nMemo1->Clear();
\nClearHistory();<\/p>\n
Octonion rootVal(1.0, 0, 0, 0, 0, 0, 0, 0);
\nRecursiveOcto* root = new RecursiveOcto(0, rootVal);
\nroot->Generate(5); \/\/ \u5f8c\u30677?8\u306b\u4e0a\u3052\u3066\u304f\u3060\u3055\u3044
\nroot->Observe(2);
\nroot->ComputeNumberOperator();<\/p>\n
int counter = 0;
\nCollectHistory(root, counter);<\/p>\n
root->Dump(Memo1);<\/p>\n
UpdateTimeCharts();
\nUpdate3DChart();<\/p>\n
delete root;
\n}
\n\/\/—————————————————————————<\/p>\n
void __fastcall TForm1::FormCreate(TObject *Sender)
\n{<\/p>\n
History = new TList();<\/p>\n
\/\/ \u30c1\u30e3\u30fc\u30c8\u304c\u3061\u3083\u3093\u3068\u30d5\u30a9\u30fc\u30e0\u306b\u7d50\u3073\u3064\u3044\u3066\u3044\u308b\u304b\u78ba\u8a8d
\nif (!ChartTime || !ChartNorm || !Chart3D) {
\nShowMessage(“\u30c1\u30e3\u30fc\u30c8\u30b3\u30f3\u30dd\u30fc\u30cd\u30f3\u30c8\u304c\u30d5\u30a9\u30fc\u30e0\u306b\u6b63\u3057\u304f\u7d50\u3073\u3064\u3044\u3066\u3044\u307e\u305b\u3093\u3002\\n”
\n“IDE\u3067\u30b3\u30f3\u30dd\u30fc\u30cd\u30f3\u30c8\u540d\u3092\u78ba\u8a8d\u3057\u3066\u304f\u3060\u3055\u3044\u3002”);
\nreturn;
\n}<\/p>\n
InitTimeCharts();
\nInit3DChart();
\n}
\n\/\/—————————————————————————<\/p>\n
void __fastcall TForm1::InitTimeCharts()
\n{
\nChartTime->View3D = false;
\nChartTime->Legend->Visible = true;
\nChartTime->Title->Text->Text = “Octonion Components”;<\/p>\n
\/\/ 8\u672c\u306e\u7dda\u30b7\u30ea\u30fc\u30ba\u3092\u4f5c\u6210
\nfor (int i = 0; i < 8; i++) {
\nTLineSeries* ser = new TLineSeries(ChartTime);
\nChartTime->AddSeries(ser);
\n\/\/ser->Title = (i == 0) ? “Re” : “e” + IntToStr(i);
\nser->Title = (i == 0) ? AnsiString(“Re”) : AnsiString(“e”) + IntToStr(i);
\nser->LinePen->Width = 1;
\n\/\/ \u8272\u306f\u304a\u597d\u307f\u3067\uff08\u4f8b\uff09
\nstatic TColor cols[8] = {clBlack, clRed, clGreen, clBlue,
\nclFuchsia, clOlive, clNavy, clMaroon};
\nser->SeriesColor = cols[i];
\n}<\/p>\n
\/\/ \u30ce\u30eb\u30e0\u7528\u30c1\u30e3\u30fc\u30c8
\nChartNorm->View3D = false;
\nTLineSeries* sNorm = new TLineSeries(ChartNorm);
\nTLineSeries* sN = new TLineSeries(ChartNorm);
\nChartNorm->AddSeries(sNorm);
\nChartNorm->AddSeries(sN);
\nsNorm->Title = “|O|2”;
\nsN->Title = “N”;
\nsNorm->SeriesColor = clBlack;
\nsN->SeriesColor = clRed;<\/p>\n
<\/p>\n
}<\/p>\n
void TForm1::Init3DChart()
\n{
\nChart3D->View3D = true;
\nChart3D->View3DOptions->Orthogonal = false;
\nChart3D->Legend->Visible = false;<\/p>\n
TPointSeries* ser = new TPointSeries(Chart3D);
\nChart3D->AddSeries(ser);
\nser->Pointer->Style = psCircle;
\nser->Pointer->HorizSize = 3;
\nser->Pointer->VertSize = 3;
\n}<\/p>\n
<\/p>\n","protected":false},"excerpt":{"rendered":"
\u30b9\u30bf\u30c3\u30af\u30aa\u30fc\u30d0\u30fc\u30d5\u30ed\u30fc\u306f\u89e3\u6d88\u3057\u307e\u3057\u305f \u5b9f\u884c\u76f4\u5f8c\u306e\u5b9f\u884c\u6642\u30a8\u30e9\u30fc\u3060\u3051\u304c\u307e\u3060\u6b8b\u3063\u3066\u3044\u307e … \u7d9a\u304d\u3092\u8aad\u3080