果树种植最优解精美图形作法探讨 - 第13页 - 数学欣赏 - 数学研发论坛

mathe 发表于 2019-12-19 20:40:13

1. 变换矩阵最后一列不需要变为1。
2. 是矩阵的右边乘上列向量而不是行向量乘上矩阵
我计算的S就是变换阵,U1,V1,Q1,N1就是四点变换后的坐标，构成正方形.

数学星空 发表于 2019-12-19 22:04:22

经过楼上mathe的指点终于得到了想要的图形：

[-1.00000000031813, 0.384042943445230, A],
,
[-3.49395920936664, 3.49395920766090, C],
[-3.49395921045821, -3.49395920661485, D],
[-9.09783468542929, 0.999999999819806, E],
,
[-9.09783468537715, -0.999999999820944, G],
,
[-1.00000000050038, -0.384042943042200, I],
,
,
,
[-3.49395920997241, -0.384042943096800, M],
,
,
,
[-1.00000000017199, 1.00000000018635, Q],
,
,
[-9.09783468531214, -3.49395920644334, T],
,
[-1.00000000064661, -0.999999999727832, V],
[-9.09783468549431, 3.49395920717132, W],
[-3.49395920985244, 0.384042943358285, X]

数学星空 发表于 2019-12-20 19:28:45

很奇怪，用同样方法计算111#，也没得到mathe的图样，不知mathe能否指点一下？

由t^2-t-1=0,得到t=1.618033988和t=-0.6180339880

下面选择t=-0.6180339880计算得到；

[, , , , , , , , , , , , , , , , , , , , , , , ]

除掉无穷远点K, L, U, V ,画图得到：

按照楼上的计算方法，选择四个点W,U,T,X做射影变换(取变换后WUT构成三角形的顶点，边长为1，X为正三角形的中心）

即W->,U[ 0., 1., 0.]->[-1/2,-sqrt(3)/6,1] , T-> , X->

求得射影变换矩阵J1为： [[-0.6545084975, -0.5000000000, 0.5000000005], , ])

注：楼下mathe得到的变换矩阵J2为 [, [-0.6114237460, 0.1556953931, 0.1556953933], [-0.4999999974, -0.5393446628, -0.5393446636]]；

J1,J2若将每行都化为z=1,就是一样均为 [[-1.309016994, -0.9999999988, 1], [-3.927050974, 0.9999999986, 1], ]

由J1射影变换得到

, , , [-1.53262380618613, 1.97860882314602, D], , [-0.0801787283406760, 0.438476638034528, F], [-0.250000000177685, 0.144337566957350, G], , , , [-0.654508499830173, 0.844966838960361, K], [-0.809017000485346, 1.97860881867382, L], [-5.54508511075584, 5.64715142478797, M], , , , , , , , [-0.500000000000000, -0.288675134700000, U], [-0.706011333290230, 1.22284749819228, V], , [-5.94287545476067*10^(-11), -7.18818497189416*10^(-11), X]

画图显示比较杂乱，肯定是哪里出了问题？？点 [-5.54508512863216, 5.64715144008720, M]比较异类离的较远，难道是因为选择作变换的四个点（WUTX)中有一个无穷远点U造成的？？？

mathe 发表于 2019-12-20 19:49:22

你怎么计算的？显然你计算出来W,U,X,T的坐标都不符合要求

mathe 发表于 2019-12-20 19:55:38

直接对上次代码改数据(字母名称就不改了，所以和你的不匹配）
Q=;
V=;
U=;
N=;
T=matrix(12,13);
for(u=1,3, T=Q); T=0;
for(u=1,3,T=Q); T=sqrt(3)/3;
for(u=1,3,T=Q); T=1;
for(u=1,3, T=V); T=-1/2;
for(u=1,3,T=V); T=-sqrt(3)/6;
for(u=1,3,T=V); T=1;
for(u=1,3, T=U); T=1/2;
for(u=1,3,T=U); T=-sqrt(3)/6;
for(u=1,3,T=U); T=1;
for(u=1,3, T=N); T=0;
for(u=1,3,T=N); T=0;
for(u=1,3,T=N); T=1;
R=matker(T)~;
S=matrix(3,3);
for(u=1,3,S=R;S=R;S=R)
Q1=S*Q~;Q1=Q1/Q1
V1=S*V~;V1=V1/V1
U1=S*U~;U1=U1/U1
N1=S*N~;N1=N1/N1
S

结果:
(19:54) gp > Q1=S*Q~;Q1=Q1/Q1
%21 = [-1.6837458349204428006 E-39, 0.57735026918962576450914878050195745565, 1.0000000000000000000000000000000000000]~
(19:54) gp > V1=S*V~;V1=V1/V1
%22 = [-0.50000000000000000000000000000000000000, -0.28867513459481288225457439025097872783, 1.0000000000000000000000000000000000000]~
(19:54) gp > U1=S*U~;U1=U1/U1
%23 = ~
(19:54) gp > N1=S*N~;N1=N1/N1
%24 = ~
(19:54) gp > S
%25 =

[-0.61142374604378055516562730110475899584 0.15569539313815436224534312637626872883 0.15569539334936645493708909649317362345]

[-0.49999999740620259722917961798090153930 -0.53934466283938819896148377628445891984 -0.53934466357104835039158612314592208695]

mathe 发表于 2019-12-20 20:53:24

113#的一些解也可以试验一下，比如:
ADGJBEIJCDHKAFIKCEGLBFHLCJMODINODLMPAHNPGKOPBGMQFJNQAEOQEHMRBKNRCFPRILQRABCSDEFSGHITJKLTMNST
Start group:
ABCDEFGHIJKLMNOPQRST 00 01 02 03 04 05 06 07 08 09 10 11 12 13 14 15 16 17 18 19 20 21 22
GHIJKLABCDEFNMOQPRTS 00 02 01 04 03 05 07 06 12 11 13 09 08 10 15 14 17 16 20 21 18 19 22
EDFBACIHGJLKNMQROPST 01 00 05 04 03 02 12 11 15 14 17 07 06 13 09 08 16 10 19 18 20 21 22
CBADFEKJLHGINMPORQST 02 05 00 04 03 01 09 08 07 06 10 15 14 16 12 11 13 17 18 19 21 20 22
IHGJLKEDFBACMNQORPTS 01 05 00 03 04 02 11 12 06 07 13 14 15 17 08 09 10 16 20 21 19 18 22
KJLHGICBADFEMNPROQTS 02 00 05 03 04 01 08 09 14 15 16 06 07 10 11 12 17 13 21 20 18 19 22
FDEBCALJKHIGMNRQPOST 05 02 01 03 04 00 14 15 11 12 17 08 09 16 06 07 13 10 19 18 21 20 22
LJKHIGFDEBCANMRPQOTS 05 01 02 04 03 00 15 14 09 08 16 12 11 17 07 06 10 13 21 20 19 18 22
这是一个8阶置换，里面好像包含四阶置换（AIFK)，所以可以试验把这四个点依次映射成一个正方形顶点看看。

   Parameter [-1+4*t+1*t^2=0]
   A[+1 ,0 , 0]
   B
   C[-1 ,+1 , +1]
   D[+1 ,-1/2-1/2*t , 0]
   E
   F[+1 ,0 , +1]
   G[+1 ,-1/2+1/2*t , 0]
   H[+3/2+1/2*t ,-1/2-1/2*t , +1]
   I
   J
   K[+1/2+1/2*t ,0 , +1]
   L[+1/2+1/2*t ,+1/2-1/2*t , +1]
   M[-1 ,+3/2-1/2*t , +1]
   N[+1 ,-1/2-1/2*t , +1]
   O[-1 ,+1/2+1/2*t , +1]
   P[+2+1*t ,-1/2-1/2*t , +1]
   Q[+1 ,+1/2+1/2*t , +1]
   R[+2/5+1/5*t ,+3/10-1/10*t , +1]
   S[-1/2-1/2*t ,+1 , +1]
   T[+1/2+1/2*t ,-1*t , +1]

mathe 发表于 2019-12-20 20:58:14

CFGNCHIODGKQEFJPDIJREHKSJKNOGIPSFHQRFIKLGHJMDLNPEMNQCLQSCMPRDMOSELORBNRSAOPQABLM
Start group:
ABCDEFGHIJKLMNOPQRS 00 01 02 03 04 05 06 07 08 09 10 11 12 13 14 15 16 17 18 19
BAECDJFKHIGMLPSRNOQ 03 05 00 04 01 02 07 08 06 10 09 14 11 12 16 13 15 18 17 19
BACDEHIFGKJMLONSRQP 01 00 04 05 02 03 06 07 08 10 09 15 16 14 13 11 12 18 17 19
ABDECIJGKHFLMRQOPSN 04 02 03 01 05 00 08 06 07 09 10 16 14 11 15 12 13 17 18 19
ABECDKHJFGILMSPQONR 05 03 01 02 00 04 07 08 06 09 10 13 15 16 12 14 11 17 18 19
BADCEKGJIHFLMQRSNOP 02 04 00 05 01 03 08 07 06 09 10 13 12 11 15 14 16 18 17 19
ABCEDGFIHKJMLNOQPSR 00 01 03 02 05 04 06 08 07 10 09 12 11 14 13 16 15 17 18 19
BAEDCFJHKGILMPSNRQO 03 05 04 00 02 01 07 06 08 09 10 11 14 16 12 15 13 18 17 19
BACEDIHGFJKLMONRSPQ 01 00 05 04 03 02 06 08 07 09 10 16 15 13 14 12 11 18 17 19
ABDCEJIKGFHMLRQPONS 04 02 01 03 00 05 08 07 06 10 09 14 16 15 11 13 12 17 18 19
ABEDCHKFJIGMLSPOQRN 05 03 02 01 04 00 07 06 08 10 09 15 13 12 16 11 14 17 18 19
BADECGKIJFHMLQRNSPO 02 04 05 00 03 01 08 06 07 10 09 12 13 15 11 16 14 18 17 19
可以试验CDE映射为正三角形，M为中心

这个只有复数解:
   [+1*t^2+1/3]
   -1*v11-1*t
   +1*v1-3/2*v11-1/2
   +1*v0+3/2*v11-1/2
   +1*v2-3/2*v11-1/2
   +1*v3-3/2*v11-1/2
   +1*v4-1/2*v11-1/2
   +1*v5-1/2*v11-1/2
   +1*v6+3/2*v11-1/2
   +1*v7+3/2*v11-1/2
   +1*v8-1*v11-1
   +1*v9-1*v11
   +1*v10-1*v11-1
   A(v6,v2)
   B(v3,v7)
   C(v5,v4)
   D(+0,+0)
   E(+1,+1)
   F(v5,v8)
   G(v5,v9)
   H(v10,v4)
   I(v11,v4)
   J[+1,v0,0];
   K[+1,v1,0];
   L(+0,+1)
   M(+1,+0)
   N[+0,+1,0];
   O[+1,+0,0];
   P(+0,v2)
   Q(+1,v2)
   R(v3,+1)
   S(v3,+0)

数学星空 发表于 2019-12-23 18:50:19

125# mathe的计算方法，我不太确定下面理解是否正确：

设变换前的四点坐标 $Q=,V=,U=,N=$

变换后的四点坐标$Q1=,V1=,U1=,N1=$

构造12X13阶矩阵：

\[T = \begin{pmatrix}
Q_{1}&Q_{2}& Q_{3}&0&0&0&0&0&0&Q1_{1}&0&0&0\\
0&0&0&Q_{1}&Q_{2}&Q_{3}&0&0&0&Q1_{2}&0&0&0\\
0&0&0&0&0&0&Q_{1}&Q_{2}&Q_{3}&Q1_{3}&0&0&0\\
V_{1}&V_{2}&V_{3}&0&0&0&0&0&0&0&V1_{1}&0&0\\
0&0&0&V_{1}&V_{2}&V_{3}&0&0&0&0&V1_{2}&0&0\\
0&0&0&0&0&0&V_{1}&V_{2}&V_{3}&0&V1_{3}&0&0\\
U_{1}&U_{2}&U_{3}&0&0&0&0&0&0&0&0&U1_{1}&0\\
0&0&0&U_{1}&U_{2}&U_{3}&0&0&0&0&0&U1_{2}&0\\
0&0&0&0&0&0&U_{1}&U_{2}&U_{3}&0&0&U1_{3}&0\\
N_{1}&N_{2}&N_{3}&0&0&0&0&0&0&0&0&0&N1_{1}\\
0&0&0&N_{1}&N_{2}&N_{3}&0&0&0&0&0&0&N1_{2}\\
0&0&0&0&0&0&N_{1}&N_{2}&N_{3}&0&0&0&N1_{3}
\end{pmatrix}_{12 \times 13 } \]

然后求解$TR=0$的非零解$R=$

(maple的求解命令为NullSpace）

则可以得到变换矩阵S为

\[S = \begin{pmatrix}
R_1&R_2&R_3\\
R_4&R_5&R_6\\
R_7&R_8&R_9
\end{pmatrix}_{3 \times 3 }\]

注：楼上123#的结果比较

利用我得到的变换矩阵J1可以得到

利用mathe的变换矩阵J2可以得到

好像结果几乎是一样的，问题是都不是我们想要的答案！！！我也没有找到到底哪里出了问题？

mathe 发表于 2019-12-24 08:08:08

分析思路如下，
设变换矩阵$S=\begin{pmatrix}R_1&R_2&R_3\\R_4&R_5&R_6\\R_7&R_8&R_9\end{pmatrix}$将
$Q=\begin{pmatrix}Q_1\\Q_2\\Q_3\end{pmatrix},V=\begin{pmatrix}V_1\\V_2\\V_3\end{pmatrix},U=\begin{pmatrix}U_1\\U_2\\U_3\end{pmatrix},N=\begin{pmatrix}N_1\\N_2\\N_3\end{pmatrix}$变换为$Q1=\begin{pmatrix}Q1_1\\Q1_2\\Q1_3\end{pmatrix},V1=\begin{pmatrix}V1_1\\V1_2\\V1_3\end{pmatrix},U1=\begin{pmatrix}U1_1\\U1_2\\U1_3\end{pmatrix},N1=\begin{pmatrix}N1_1\\N1_2\\N1_3\end{pmatrix}$
于是存在常数$-R_{10},-R_{11},-R_{12},-R_{13}$使得
$SQ=-R_{10} Q1,SV=-R_{11} V1， SU=-R_{12} U1, SN=-R_{13} N1$
展开后就出现楼上的矩阵$T$了,即
$T\begin{pmatrix}R_1\\R_2\\R_3\\R_4\\R_5\\R_6\\R_7\\R_8\\R_9\\R_{10}\\R_{11}\\R_{12}\\R_{13}\end{pmatrix}=0$.
解这个方程得出任意的非零向量，取前9个分量构成矩阵S即可

mathe 发表于 2019-12-24 15:13:35

75#的20棵树23行可以做出一种更加对称的形式，只是直线CDJK是无穷远直线:

页: 3 4 5 6 7 8 9 10 11 12 [13] 14 15 16 17 18

数学研发论坛's Archiver