18/08/09 03:02:30.33 w0c/gBS2.net
>>47
長崎はきょうです。。。
51:１３２人目の素数さん
18/08/09 14:12:16.14 DUTF3h+v.net
雨だね。
52:１３２人目の素数さん
18/08/10 03:17:49.07 MxWQLJMW.net
アメダス （なぜか関西弁）
URLﾘﾝｸ(www.jma.go.jp)
53:１３２人目の素数さん
18/08/10 03:46:03.76 MxWQLJMW.net
〔次の問題〕
周の長さが一定の四角形（長方形以外も含む）はどうしたら面積が最大になるの？
ｽﾚﾘﾝｸ(math板)
54:イナ
18/08/10 11:43:41.24 O38KwL8o.net
;;;;;;;;人;;;;;;;;;;;
;;;;;;;(＿);;;;;;;;;;
;;;;;;(＿)_)､;;;;;;;;
;;;;;(_)＿)_))､;;;;;;;
;;;;(_(＿)_)_)_);;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;>>52四辺が同じ長さのとき。前>>50 ５月５日はこどもの日だろ。６月４日は虫歯予防デー。わかるら。
55:１３２人目の素数さん
18/08/11 05:48:14.60 OvqaCq6H.net
>イナ
こいつには証明という概念はないの？
56:１３２人目の素数さん
18/08/11 18:59:45.87 05wxJPem.net
>>52
それはわかってるんですよ。基本資料だから読んでないわけがない。そのページの２ 標準偏差（２SD）って具体的にどんな数値を入れて計算するのかを聴いているんです。
57:１３２人目の素数さん
18/08/15 08:10:57.93 ui1ogxnG.net
五角形の中心と各頂点を結ぶ放射状の線だと題意をみたさないのか？
58:１３２人目の素数さん
18/08/15 11:39:53.66 6OEqz/R6.net
URLﾘﾝｸ(i.imgur.com)
のシオマネキのハサミを覗いた部分は左右対称と仮定。
５角形の底辺の中点を原点として右側の三差路の交差点の座標を(x,y)として
コンピューターに最小値を算出させてみた。
座標 (0.3799888 0.2080354)で 最長値は3.956295と出てきた。
> optim(c(0.4,0.4),sio,method='Nelder-Mead')
$`par`
[1] 0.3799888 0.2080354
$value
[1] 3.956295
$counts
function gradient
45 NA
$convergence
[1] 0
$message
NULL
59:１３２人目の素数さん
18/08/15 11:54:36.60 6OEqz/R6.net
URLﾘﾝｸ(i.imgur.com)
カブトガニでもやってみた
> optim(c(1,0.5),kab, method = 'Nelder-Mead')
$`par`
[1] 1.0614192 0.7726595
$value
[1] 3.891157
$counts
function gradient
57 NA
$convergence
[1] 0
$message
NULL
60:１３２人目の素数さん
18/08/15 11:57:20.54 6OEqz/R6.net
>>57
最長値は3.956295と出てきた。→最小値は3.956295と出てきた。
61:１３２人目の素数さん
18/08/15 12:18:27.42 6OEqz/R6.net
>>57
Rのソースはこれ
ｽﾚﾘﾝｸ(hosp板:614番)
62:１３２人目の素数さん
18/08/15 21:24:45.40 6OEqz/R6.net
>>58
表示桁を増やしてみた。
> optim(c(1,0.5),kab, method = 'Nelder-Mead')
$`par`
[1] 1.0614192103443205 0.7726595158193953
$value
[1] 3.891156833029388
$counts
function gradient
57 NA
$convergence
[1] 0
$message
NULL
URLﾘﾝｸ(i.imgur.com)
底辺の中点を原点において
真ん中の交差点の（０，ｙ１）、右側の交差点が（０．５、ｙ２）としたとき
y1=1.0614192103443205 y2=0.7726595158193953
で 最小値が 3.891156833029388
というのがコンピュータでの結論。
等高線をグラフ化してみた。
URLﾘﾝｸ(i.imgur.com)
63:１３２人目の素数さん
18/08/15 21:30:30.90 eFYm1A33.net
左右非対称解とは結局何だったのか
64:１３２人目の素数さん
18/08/15 21:40:36.66 6OEqz/R6.net
>>62
min(カブトカニ形) ＜ min(シオマネキ形)はいえるけど
カブトカニ形より小さい解がないと言えるんだろうか？
65:１３２人目の素数さん
18/08/15 23:27:10.16 m/c8a9aB.net
言えてる。もう証明も上がってる。
66:イナ
18/08/15 23:34:05.57 0zF2NnMt.net
前>>53前々>>48カブトガニ最小=3.8911568225
分岐点2つの距離をxとおくと、
シオマネキ最小=1+{(1+√5)/2-x}×2+(1-x)×2+x
=4+√5-3x―①
斜辺1で三平方の定理。
√[1^2-{(1+√5)/2-1}^2･(1/2)^2]={(1+√5)/2-x}(√3)/2+(1-x)(√3)/2
これを解いて、
x={3√3+√15-√(10+2√5)}/4√3
②を①に代入。
シオマネキ最小=4+√5-{9+3√5-(√3)√(10+2√5)}/4
={7+√5+√(30+6√5)}/4
=3.9562341＞カブトガニ最小
67:１３２人目の素数さん
18/08/15 23:37:22.56 9MuNAPAt.net
この期に及んでまだわかってないんだな。
68:イナ
18/08/15 23:46:15.81 0zF2NnMt.net
前>>65加筆。
>>48カブトガニ最小=3.8911568225と同様、シオマネキ最小についても作図により式を立て、パソコンなどに一任することなく、中学生がわかる程度の計算によって示したいと思う。
分岐点2つの距離をxとおくと、
シオマネキ最小=1+{(1+√5)/2-x}×2+(1-x)×2+x
=4+√5-3x―①
斜辺1で三平方の定理。
√[1^2-{(1+√5)/2-1}^2･(1/2)^2]={(1+√5)/2-x}(√3)/2+(1-x)(√3)/2
これを解いて、
x={3√3+√15-√(10+2√5)}/4√3―②
②を①に代入。
シオマネキ最小=4+√5-{9+3√5-(√3)√(10+2√5)}/4
={7+√5+√(30+6√5)}/4
=3.9562341＞カブトガニ最小
69:１３２人目の素数さん
18/08/15 23:54:12.57 9MuNAPAt.net
そこじゃないんだよなぁ。
それでいいん？君の数学はそのレベルで終わりでいいん？
70:イナ
18/08/16 12:10:17.67 SQPcPMZF.net
>>68ビジュアル八面体みたいに120°じゃない極値があるってのかい？ 前>>67
カブトガニ最小=3.8911568225
シオマネキ最小=3.9562341
0.555……×7=3.888……
71:１３２人目の素数さん
18/08/16 12:18:35.11 GB9705T4.net
ホントにわかってないんだなぁ。
>>63が何を疑問に思ってるのかすら君一人だけわかってないんだなぁ。
72:イナ
18/08/16 13:12:58.98 SQPcPMZF.net
前>>69仮説。
シオマネキ最小＞カブトガニ最小＞分岐線の最小=35/9=3.888……
73:イナ
18/08/16 19:49:21.02 SQPcPMZF.net
前>>71
小さい蠏。
74:１３２人目の素数さん
18/08/16 19:49:59.26 508YkaXu.net
>>71
最小解があると言うなら、その解を図示するか、または座標を示してみないかい？
75:イナ
18/08/16 20:40:46.57 SQPcPMZF.net
>>71最小値に座標は関係ない。前後左右どうまわしても最小値になるときはなる。分岐線の長さxとおくと、最小値は7x。前>>72でももう疲れた。カブトガニが最小の蠏かもね。
分岐の角度を120°に保ちながらカブトガニより小さくならないかな。
もしなるならピタゴラス。xとyの二次式から四次式ぐらいで抑えないときついよ、ピタゴラス。寝冷えだ。
76:１３２人目の素数さん
18/08/16 20:56:43.85 508YkaXu.net
>>74
ふうん。そう。
おやすみ。おだいじに。
77:１３２人目の素数さん
18/08/16 21:06:24.90 NCnvn0cf.net
なんの意味もないスレ
78:高添沼田の親父「糞関東連合テメエらまとめてぶち殺すっ!!」
18/08/16 21:15:48.20 dZ5ratnn.net
高添沼田（葛飾区青戸6－２３－２１ハイツニュー青戸103号室）の挑発
高添沼田の親父「関東連合文句があったらいつでも孫を金属バットで殴り殺しに来やがれっ!! 関東連合の糞野郎どもは俺様がぶちのめしてやるぜっ!! 賞金をやるからいつでもかかって来いっ!!糞バエ関東連合どもっ!! 待ってるぜっ!!」 (挑戦状)
79:１３２人目の素数さん
18/08/17 09:27:54.93 VyDJkqSp.net
URLﾘﾝｸ(i.imgur.com)
のように左右対称を前提とせず、カブトガニ形の繋がりのネットワークを対象として
コンピュータで長さが最小となる座標を出して描出させてみた。
つまり変数６個で計算。
> kabu3(opt$par)
[1] 3.891156823335537
> (par=opt$par)
[1] 0.000002699127792359672 0.772642997991326407892 0.500000462950937363260
[4] 1.061317339518144553523 0.999998560151592941203 0.772643514290974220415
Rのコードはここ
ｽﾚﾘﾝｸ(hosp板:621番)
この座標で描画すると
URLﾘﾝｸ(i.imgur.com)
80:１３２人目の素数さん
18/08/17 09:58:33.93 mWFnumrr.net
>>78
この試行はやってみたことある
まずランダムに３点をとって、その３点でできる三角形の短い方から２辺と、元の正五角形の各頂点から、３点のうち最も近い点とを結んでそれら７本の長さを合計する
それから３点を少しずつ動かして極小解を探していく
結果はカブトガニ型に収束するか、４より大きな解になるか、または収束しない
シミュレーションだから最適化できているかは不明だが、カブトガニ解より小さな解はこれまで見つかっていない
81:１３２人目の素数さん
18/08/17 10:15:00.75 VyDJkqSp.net
>>79
レスありがとう。
>78は座標の取りうる値に制限をつけてないので
ネットワークが重なるようなのも対象にして算出している。
URLﾘﾝｸ(i.imgur.com)
82:１３２人目の素数さん
18/08/17 11:06:28.65 VyDJkqSp.net
>78の最小値の最適解算出の初期値は
０～１の一様分布でランダムに設定したけれど
URLﾘﾝｸ(i.imgur.com)
のような五角形外から始めても左右対称のカブトガニに収束した。
> optim(runif(6,-1,2),kabu3,method='CG')
$`par`
[1] -0.0000006244007777030901 0.7726459096103870383132 0.5000003027472189609526
[4] 1.0613226509328328450721 1.0000013550880015866795 0.7726459944754771047570
$value
[1] 3.891156823329757
83:１３２人目の素数さん
18/08/17 13:17:28.10 5QyvDwxU.net
お墓参りに来ますた。南無阿弥陀仏。
（まだ生きていたとは…）
84:１３２人目の素数さん
18/08/17 13:34:20.78 VyDJkqSp.net
>79
>まずランダムに３点をとって
５角形の頂点をp1～p5として分岐点をq1～q5の５点に増やして
p1-q1,p2-q2,p3-q3,p4-q4,p5-q5とq1-q2,q2-q3,q3-q4,q4-q5の長さの合計を
最小値とする座標を計算させてみた。
明らかに最小でない、こういうの候補と
URLﾘﾝｸ(i.imgur.com)
URLﾘﾝｸ(i.imgur.com)
最小とするパラメータ
> (opt=optim(runif(10),penta,method = 'CG'))
$`par`
[1] 1.0000004341739945346035 1.0000008344836242013542 0.4999996331498277601924
[4] 0.0000003625715489940964 0.0000001533897055299300 0.4058822346805296410466
[7] 0.7726451907010785102869 1.0613213678038868614806 0.7726454972413701050016
[10] 0.3769465135213049822305
$value
[1] 3.89115682332915
これで描画すると
URLﾘﾝｸ(i.imgur.com)
p1,q1,q2が一直線、p5,q5,q4も一直線で
分岐３点のカブトガニに収束した。
85:１３２人目の素数さん
18/08/17 13:39:00.61 VyDJkqSp.net
>>83
R の コードはここに置いた
ｽﾚﾘﾝｸ(hosp板:451番)-452
86:１３２人目の素数さん
18/08/17 14:27:36.94 VyDJkqSp.net
>>83
×こういうの候補と
〇こういうの候補を含めて
87:１３２人目の素数さん
18/08/17 15:20:06.10 Xs+I9BdE.net
>>85
→こういう候補を含めて
88:１３２人目の素数さん
18/08/18 08:51:54.24 IFsyvQ2o.net
>>83
分岐点５点設定でも３分岐点に収束するのだから
３分岐点は２分岐に収束しないのでカブトガニ最小が濃厚だな。
89:１３２人目の素数さん
18/08/21 09:55:54.41 Vvkk5uPP.net
７角形でこういうモデルのとき URLﾘﾝｸ(i.imgur.com)
の最適解をコンピューターで算出してみた。
Rのコードはこれ
rm(list=ls())
n=7
ngon <- function(n,digit=TRUE,axis=FALSE,cex=1,...){ # draw n-polygon
r=exp(2*pi/n*1i)
p=complex(n)
for(i in 1:(n+1)) p[i]= (1-r^i)/(1-r)
plot(p,bty='l',type='l',axes=axis, ann=FALSE,asp=1,...)
points(1/(1-r),pch='.')
if(digit) text(Re(p),Im(p),paste('p',1:n),cex=cex)
if(axis){axis(1) ; axis(2)}
invisible(p)} # return vertex complex
seg <- function(a,b,...){# draw segment of complex a to complex b
segments(Re(a),Im(a),Re(b),Im(b),col=2,...)}
pt <- function(x,y=NULL,...){ # draw text y at complex x
text(Re(x),Im(x), ifelse(is.null(y),'+',y), ...)}
poly_demo <- function(x=runif(n),y=runif(n)){ # draw segments
graphics.off()
p=ngon(n,axis=T,col='skyblue')
Q=complex(n)
re1=re2=0
for(i in 1:n){
Q[i]=x[i]+y[i]*1i
pt(Q[i],paste('q',i))
seg(p[i],Q[i])
re1=re1+abs(p[i]-Q[i])
}
for(i in 1:(n-1)){
seg(Q[i],Q[i+1])
re2=re2+abs(Q[i]-Q[i+1])
}
return(sum(re1)+sum(re2))}
poly_demo()
p=ngon(n,axis=T,col='skyblue')
poly <- function(par){
x=par[1:n]
y=par[(n+1):(2*n)]
Q=complex(n)
re1=re2=0
for(i in 1:n){
Q[i]=x[i]+y[i]*1i
re1=re1+abs(p[i]-Q[i])}
for(i in 1:(n-1)){
re2=re2+abs(Q[i]-Q[i+1])}
return(sum(re1)+sum(re2))}
(opt=optim(runif(2*n),poly,method = 'CG'))
par=opt$par
ngon(n,axis=T,col='skyblue')
poly_demo(par[1:n],par[(n+1):(2*n)])
90:１３２人目の素数さん
18/08/21 13:26:45.02 Vvkk5uPP.net
# 各頂点から最も近い分岐点を選んで結ぶモデル
rm(list=ls()) ; graphics.off()
# draw n-polygon
ngon <- function(n,print=TRUE,digit=TRUE,axis=FALSE,cex=1,...){
r=exp(2*pi/n*1i)
p=complex(n)
for(i in 1:(n+1)) p[i]= (1-r^i)/(1-r)
if(print){
plot(p,bty='l',type='l',axes=axis, ann=FALSE,asp=1,...)
points(1/(1-r),pch='.')
if(digit) text(Re(p),Im(p),paste('p',1:n),cex=cex)
if(axis){axis(1) ; axis(2)}}
invisible(p) # return vertex complex
}
# draw segment of complex a to complex b
seg <- function(a,b,...){
segments(Re(a),Im(a),Re(b),Im(b),col=2,...)
}
# draw text y at complex x
pt <- function(x,y=NULL,...){
text(Re(x),Im(x), ifelse(is.null(y),'+',y), ...)
}
91:１３２人目の素数さん
18/08/21 14:30:05.10 Vvkk5uPP.net
各頂点からｍ個の分岐点のうち最も近い分岐点を選んで結ぶというモデルにアルゴリズムを変更してみた。
初期値は多角形に外接する正方形の座標を一様分布で選んで、分岐線の長さが最小になる値を
コンピュータに探索させた。
ｍ＝１は警視庁 URLﾘﾝｸ(i.imgur.com)
ｍ＝２はカメムシ URLﾘﾝｸ(i.imgur.com)
ｍ＝３はカブトガニ URLﾘﾝｸ(i.imgur.com)
ｍ以上はカブトガニが最小になったが、初期値によっては時々、カメムシが出てきた。
ｍ＝４ URLﾘﾝｸ(i.imgur.com)
ｍ＝５URLﾘﾝｸ(i.imgur.com)
ｍ＝６URLﾘﾝｸ(i.imgur.com)
ｍ＝７URLﾘﾝｸ(i.imgur.com)
ｍ＝８URLﾘﾝｸ(i.imgur.com)
ちなみに６角形以上では多角形の辺を結ぶ分岐線が最小になって味気ない。
92:１３２人目の素数さん
18/08/21 14:45:15.44 hlpHfk+X.net
ｍ以上、→ ３以上、
93:１３２人目の素数さん
18/08/21 19:21:56.63 Vvkk5uPP.net
３角形から８角形までコンピュータに探索させてみた。
URLﾘﾝｸ(i.imgur.com)
94:１３２人目の素数さん
18/08/31 11:02:18.25 ww21fRDi.net
>>89
ソース全文はココ
ｽﾚﾘﾝｸ(hosp板:635-638番)
95:１３２人目の素数さん
18/09/15 14:07:09.33 XAWbYE0K.net
【新定理】辺の長さが全て整数となる直角三角形と二等辺三角形の組の中には、周の長さも面積も共に等しい組が(相似を除き)1組だけ [485983549]
ｽﾚﾘﾝｸ(poverty板)
96:１３２人目の素数さん
18/09/21 20:25:39.78 WL3LZOyH.net
シャボン液を使ったら膜が動いて極小を示せるか？
最小の証明にはならんが
97:１３２人目の素数さん
20/11/25 12:59:02.23 uNm3BuF0.net
>>94
直角⊿ (a,b,c) = (135,352,377) と
二等辺Δ (x,y,y) = (132,366,366)
は条件をみたす。（L=864, S=23760)
Denis Boris (2003/Oct/22) が a<400000 には他の例がないことを示す。
URLﾘﾝｸ(mathworld.wolfram.com)
これを IBM research の Ponder this に出題
Dan Dima (2004/Feb/12 07:07:05) が解の uni
98:quity を示す。 Michael Stoll (2004/Feb/27 16:25:23) 〃 http://rio2016.5ch.net/test/read.cgi/math/1537412180/ 慶應義塾は数論幾何学における「p進 Abel積分論」と「有理点の降下法」を応用したらしいけど そんなの使う必要あんの？
99:１３２人目の素数さん
21/02/19 03:41:37.24 45fvrIx7.net
定規とコンパスにより正五角形を作図する方法
ARを直径とする円Xを描く。
これに内接する正五角形 ABCDEA を作図しよう。
A (-1, 0)
C (cos(36), sin(36))
D (cos(36), -sin(36))
R (1, 0)
T (1/2, 0)
とする。
第二余弦定理より
CT^2 = 1 + 1/4 - cos(36) = 5/4 - φ/2 = 5/4 - 29/36 = 4/9,
CT = DT = 2/3,
直径ARの4等分点Tを中心とし、ARの1/3を半径とする円Yを描く。
円X と 円Y の交点を C および D とする。
ACの垂直2等分線と円Xの交点を B とする。
ADの垂直2等分線と円Xの交点を E とする。
弦AB = BC = CD = DE = EA, (終)
URLﾘﾝｸ(suseum.jp)
100:過去ログ ★
[過去ログ]
■ このスレッドは過去ログ倉庫に格納されています