松屋で海老のチリソース定食ポテトサラダセットを食す。海老が予想以上に多く、甘み多めの甘辛のエビチリでした。

実はあまりエビが得意じゃない(特に生のを食べるとかゆくなる、軽いアレルギー）のだが、火が通っているなら大丈夫だろう、ということで珍しい松屋のエビチリを注文。ものすごくエビが多い。しかもぷりぷり。これなら私でも食べられる。

が、甘い、、、甘辛ではあるが甘みが強い。まあ万人向けにするにはこうしないと、ということか。

卓上の唐辛子をかなり多めにかけて辛さを調整していただきました。

FPUT問題(Fermi-Pasta-Ulam-Tsingou)の2次の非線形の振動子の動きを Python＋Scipyでルンゲクッタ8次のDOP853(Dormand＆Prince)を使って計算

昨日はカテナリになるのをみたけれど、ほぼ同じようなプログラムでFPUT問題（過去はFPU問題と言われていた）が計算できる。

（FPUTについてはこちら）

やっているのはこちらと同じく、

d²x_i/ dt² = (x_i+1 + x_i-1 - 2x_i) + α[(x_i+1 - x_i)²- (x_i - x_i-1)²]

の計算。通常はエネルギーのモードごとの図を描きますが、

実際の動きをGIFアニメにした。

こちら。

これだけ見てても再帰現象はわからんな、、、

鎖をそのまま置くと懸垂曲線（カテナリ―）になるというのをPython＋Scipyでルンゲクッタ8次のDOP853(Dormand＆Prince)を使って計算

Fermat Libraryでこういうのを見た。

あれ？カテナリ―（懸垂曲線）になるんじゃ？と思いますが、張力が場所によって違う効果を入れると
d²y/dx²=ρg/T √(1+(dy/dx)²)
となってカテナリ(coshの形）になります。この話はほぼまっすぐになっていて√の中身が1になる近似を使っている。 https://t.co/gk3vjTDntA

— tomo (@tonagai) September 27, 2020

ばねでつながれた質点系でカテナリになるのかどうか試してみようか、とやってみた。使うのはPython＋ScipyのDOP853.

計算する式はこちらで。

d²x_i/ dt² = -γdx_i/dt + k(x_i+1 + x_i-1 - 2x_i) -g

結果：黒いのがカテナリ。

おお、なるほどカテナリになってる。

角宮神社＠長岡京でお参り。

ようやく涼しくなったので散歩がてらお参りしてきました。

「ブルータスの心臓（新装版）」（東野圭吾さん）を読んだ。元デンソーの技術者だったことがうかがえる産業用ロボットが絡む完全犯罪もの。

とはいっても最初の方で完全犯罪は頓挫するという、、、

あらすじは「最新ロボット研究者で野心家の末永拓也は、勤務先の創業者令嬢・星子との結婚をもくろんでいた。しかし遊び相手の康子から妊娠を告げられる、、、しかし康子は拓也以外の男とも関係を持ち、それぞれを脅していた。星子の兄の直樹もその一人であり、直樹から拓也ともう一人の男に康子を殺さないか、という殺害計画を持ち掛けられる。その計画とはその3人を使った完全犯罪であった。計画は順調に見えたが、、、」というもの。

最初に産業ロボットの事故で作業員が1人死ぬ、という話がどいう絡んでくるのかと思っていたらそう来ましたか、というもの。

ブルータスは最新型の産業ロボットですが、心臓は？と思ったら最後の最後で「心臓はない」ということが効くのか。タイトルが凝っているということが最後にわかる。

高石神社（大阪）でお参り。

浜寺公園へ行ってきたときにお参りしてきたところ。ここから（続く）。

また印象操作グラフか、、、早稲田アカデミーやフジテレビの次は福島テレビか。コロナウイルスにかかるのは本人責任かというグラフ。ちゃんとしたスケールで描いてみると？あれ？これで十分伝わるじゃないか。なぜこんなことしたのか。。。

以前に、

早稲田アカデミーの合格者数グラフをまじめに考える。

早稲田アカデミーの合格者数グラフを真面目に考える（2012）

フジテレビの警察官の世代別不祥事（懲戒者数）の円グラフがひどすぎ、、、→なんとそんなグラフが簡単に描ける”Wonder Graph Generator”が！

とか印象操作グラフの話を書いてましたが、また出てきました。

これ。

番組内のグラフに「印象操作」　SNS物議で福島テレビが謝罪「体制の不備でした」

なるほど。ちゃんとしたスケールで描いてみよう。

右が今回のグラフ、左がちゃんとしたグラフ。

あれ？日本が突出しているというのは十分伝わるじゃないか。

なんでこんな印象操作グラフを作ったんだろう。。。

怪盗不思議紳士（我孫子武丸さん）を読んだ。関智一さんの劇団向けの脚本から小説が生まれたということで、場面が思い浮かぶように面白い！

我孫子さんの作品を読むのも本当に久々ですが、期待を裏切らない面白い作品。

あらすじは

「終戦間もない日本。戦災孤児の瑞樹はひょんなことから探偵、九条響太郎（空条承太郎じゃない）の助手になる。人を殺さないことで有名な伝説の怪盗不思議紳士とも渡り合う九条だが、ある事件で殺さないはずの不思議紳士が一家を惨殺する。その捜査中、九条の乗った車が大爆発し、、、」

というもの。

こういう展開かな？これはこの人が変装しているのかな？というのがことごとく外れるという楽しい展開で、最後は、、、キャッ〇〇イ？

一人だけ可哀そうな方は出てくるものの、とても面白い展開でした。これ続くんだろうか。関さんによるとこれは始まりの物語と言う。続編でるなら読んでみたい。

阪急電車のすみっコぐらしラッピングに乗った。癒される。。。

天下茶屋から京都まで阪急に乗ったのですが、たまたますみっコぐらし号でした。これは癒される、、、

浜寺公園とその中のばら庭園に行ってきた。そんなにばらは咲いてはいないがそれでも楽しい。

連休でもあまり人がいるところには行けない、ということで浜寺公園へ。とはいっても遊具があるようなところは子供が一杯だが、ばら庭園付近は人もまばらでした。

シーズンではないといっても少しずつはバラが咲いていて綺麗。

松のやでたっぷりタルタルメガチキンかつ定食を食す。とんでもないタルタルソースの量！

メガチキンが復活していたので注文。前はトマトチーズを食べたが今回はタルタルで。

ものすごいタルタルソースの量でした。どうしてもちょっと残っちゃった。

ご飯は無料の大盛にしたのにタルタルが余るとは、、、ケチってない感じがいいです。

なか卯で黄金の親子丼と酸辣湯うどんを食す。どちらもレベル高い！

まずあんかけがかかっている黄金の親子丼。これはとても美味しい。しかも生姜がついていてそれを入れるとさらにおいしさUP。

通常の親子丼も美味しいがこれはすごく好きな味。

酸辣湯うどんも、辛さは足りないが（なので卓上の唐辛子で調整）、酢と具の多さがいい。どちらもお勧め。

映画TENET/テネットをexpocityのIMAXレーザーGTで観てきた。すごいなこれは。カメラを止めるな＋ドラえもんのある映画＋ホイーラー（ファインマンの師匠）＋マクスウェルの悪魔か？2回見ないとわからないシーンが多すぎる！あと対消滅しないのが納得できん！

だいぶ前に同じくIMAXで長い予告編が流れていてめちゃくちゃ面白そうだったので期待していた。

もう6か月くらいコロナの影響で映画を観に行けてなかったけれど、最初に観に行くのはこれと決めていたのだ。

もちろんエキスポシティのIMAXの大画面で。なんと冒頭のシーンがあの長い予告編のシーンだった！

クリアファイルもらった。左側がIMAXの天井の証明が写っている。

内容ですが、もう感想として”すごいなこれは！”しか出てこない。映像がすごいのと、音楽が非常に不安にさせるものでこれもすごい。

ネタバレしてもたぶん、何言っているかわからないと思います（観ないと）。ただ観ても？？？というのも結構ある。

感想は

・前のシーンで出てきたあれ？というのが後になってわかるというのはカメラを止めるなを思い出す。

・黒幕とか、最後に助けてくれる人とかはドラえもん　のび太の大魔境（オシシ仮面かもしれないが）。

・なぜ時間が、、、というのはものすごくすっとばしているが、思わせぶりに電子が逆行すると陽電子になるという話を。

　これ、ファインマンの先生だったホイーラーが、電子が一種類しかない理由は全部同じもので、時間を逆行すると陽電子になる、

　という話を思い出す。

　しかし2人が出会って対消滅しないのが納得できーん。会うなって言ってたような気がする（保護スーツ着れば大丈夫とか言ってたか）。

　いやでもほんとだったら逆行した瞬間に何にぶつかっても空気にぶつかっても対消滅すると思う。

・時間を逆行する回転ドアはマクスウェルの悪魔のオマージュでしょうね。

とにかく、最初にぼーっとみていると後で起きたことがあれ？ということになる。2回見て、このシーンは後でこうなるというのを覚えておかないと、、、絶対もう一回観に行く。

「感電の危険性（人体電流と通電時間による症状）」の計算の自作式をカシオの高精度計算サイトkeisan.casio.jpにUP!IEC60479-1に基づくもので、低電圧（＜600V）の研修を受けてきたので。WebPlotDigitizerを利用した。

今の会社に入社以来、ほぼ高周波（ミリ波くらいまで）の開発をしていて、まあ扱う電圧と言えば5V以下程度（MEMSスイッチの内部で数十Vに昇圧というのはあるけど）しか扱ったことはなかった。

ちょっと事情があって、低電圧講習（といってもDC 750VでAC 600Vで私が扱ったことのないレベル）を受けてきた。

※ついでに、というか感電したときの救助の研修、AEDとか胸骨圧迫とかの研修もついていてそっちがメインかも。

そこで出てきた話で、人体がどれくらいの電流がどのくらいの時間流れると危ないかというのがあって、ちょっとこれは注意喚起の意味でも自作式を作ってみようと思った。

分類がIEC60479-1に記載されているそうだ（どうやって実験したんだろうというのは怖い気がするが）。
DC版とAC版（15-100kHz)があるが、ACのほうが被害大らしいのでそっちを作った。

どうやってつくったか？というと、グラフの画像はあるが関数形はまったくわからない。そこで、ダウンロード不要でグラフの数値読み取りができる

WebPlotDigitizer

https://automeris.io/WebPlotDigitizer/

を使って、図からわざわざグラフの点を読み取り、適当な関数でフィッティングした。

※y=(a0+a1*x+a2*x^2)*tanh(b0+b1*x+b2*x^2+b3*x^3)+c0+c1*x+c2*x^2 という関数です。

それがこれ。めんどくさかった、、、

関数形が分かればどの分類かは簡単にわかる。

こういう自作式をカシオの高精度計算サイトkeisan.casio.jpにアップロードしておいた。

リンクはこちら：

https://keisan.casio.jp/exec/user/1600342821

50mAが一秒流れると死の危険というのがわかる。

感電の危険性（人体電流と通電時間による症状）

低電圧でも人体に電流が流れると危険です。特に流れる電流値（人体電流）と何秒間流れるか（通電時間）で危険度が変わります。直流より交流の方が危険であり、IEC60479-1に記載されている分類を計算します(15-100kHz)

分類： AC-1：無反応 AC-2：有害な生理的影響なし AC-3：痙攣性の筋収縮や呼吸困難，一時的な心拍停止や心房細動を含んだ回復可能な心臓障害が生じる AC-4：心停止，呼吸停止，重度の熱傷が生じ死の危険が生じる 　AC-4.1：心室細動の確率およそ 5％以下． 　AC-4.2：心室細動の確率およそ 50％以下． 　AC-4.3：心室細動の確率およそ 50％以上．

高周波（RF・マイクロ波・ミリ波・5G）関連ニュース2020年9月17日 IEEE microwave magazine10月号はマイクロ波顕微鏡特集、 Microwave Journalはミリ波パッケージ。Vayyarの4Dイメージングレーダの仕組み、MIPIのA-PHY、など。

IEEE Microwave magazineの特集はマイクロ波”顕微鏡”の特集。物質の構造だけでなく、生物や、SAWのようなアコースティックデバイスの表面まで見られるとか。これは今後も進む分野だと期待してます。

https://ieeexplore.ieee.org/xpl/mostRecentIssue.jsp?punumber=6668

一方、Microwave Journalはミリ波部品とパッケージ。今後ますますパッケージは重要になりますが、AiPやFOWLPよりPolyStrata押しの記事。

https://www.microwavejournal.com/publications/1/editions/275

Vayyarの4Dイメージングレーダの仕組みとは？

Under the Hood: Vayyar’s 4D Imaging Radar

https://www.eetimes.com/under-the-hood-vayyaars-4d-imaging-radar/

MIPIのA-PHY（自動車向けインターフェース）の紹介

MIPI Alliance Releases A-PHY SerDes Interface for Automotive

https://www.mwrf.com/markets/automotive/article/21141858/mipi-alliance-releases-aphy-serdes-interface-for-automotive

Microsoft Surface Duoの分解

Microsoft Surface Duo Teardown

https://www.ifixit.com/Teardown/Microsoft+Surface+Duo+Teardown/136576

めしや食堂（大日）でトロサバの煮付けと夏野菜の白和えを食す。

もう暑くて何も食べられない、、、と思ったけどサバの煮付けなら、ということで注文。なんとか食べられた。

神津神社＠十三でお参り。

久しぶりに十三に来たので、駅近くのこちらの神社でお参り。

餃子の王将で麻婆茄子（フェアセットA＋餃子追加）を食す。これは茄子もたっぷりでレベル高い！

麻婆茄子は大好きなんですが、まあそんなに茄子は入ってないだろう、とか期待してなかったのですが、、、

これはたっぷり＋トマトまで入っていて辛さもちゃんとあってすごくレベル高い。とても美味しかった。

ジャージャー麺の方も食べてみたい。

 

すき家でお好み牛玉丼広島MIX（特盛）を食す。茶色い、、、そして一体何を食べているのか途中でわからなくなる。

元々茶色い牛丼に焼きそばが乗ってどこから見ても茶色い。しかも焼きそばだけが口に入ると？？？これ一体何食べてたんだっけ？と混乱する。ジャンク感は半端ないが、面白い味ではあります。

松屋で味噌漬けトンテキ丼生野菜セット（大盛、小鉢は とろろ）を食す。味がすごく濃くて量もたっぷり。

思っていた以上に肉がたっぷり。そして味噌漬けの味が非常に濃くてご飯にとてもあう。なので大盛りご飯にしたが、さらについでにハイボールも注文。とろろもいいアシスト。

天理スタミナラーメン＠近鉄奈良駅前をいただく。ニンニクがすごい。

彩華ラーメンは何度か食べたことがあるが、天理スタミナラーメンは初めて。これもニンニクが一口食べただけでぶわーと効いてきて、本当にスタミナが付きそうな味。スープまで完食しました。

JR奈良駅にアマビエ様とせんとくんがいた。

行きは近鉄だが帰りはJRで、ということでJR奈良駅へ。昔と全然変わって綺麗になってビエラもできてた。

駅構内にはせんとくんと、その反対側にアマビエ様がいた。でも人が少なかった。。。コロナも早く終息してほしい。

氷室神社＠奈良でお参り。

奈良国立博物館の道を挟んだ前にあります。

四捨五入殺人事件（井上ひさしさん）を読んだ。確かにすごいどんでん返しだ。東野圭吾さんのある作品のルーツになっているのかな？

これは確かに帯にあるようにだいどんでんがえしだ。あらすじは、

「大作家と新鋭作家が講演のため向かった先はテレビもない山間の温泉宿。大雨のためその温泉宿がある村に一つしかない橋が流され、孤立してしまった。その2人の前に起こる連続殺人事件。事件の背後に横たわるのは、何世代にもわたる村人の怨念なのか、、、嫌疑をかけられた大作家。そして新鋭作家が推理する。」

というもの。時代を感じさせる表現（登場人物が女性をいやらしい目で見すぎとかストリップとか）などはありますが、それはまあ仕方ない。農業問題などは今でも解決はしてないし。

そして大どんでん返しのラスト。そして実はフェアな証拠がだいぶ前の方から示されていたという、、、

東野圭吾さんのある作品（ラストが結構ハッピーエンドなやつ）を思い出した。あのルーツはここにあるのかな？と思ったり。

なかなか面白かったです。

かつやで鶏ちゃん焼チキンカツ定食（ご飯大盛）をいただく。

鶏ちゃんってなかなか食べる機会がないけれど、メニューにあれば必ず頼む。かつやで販売しているということで食べてみた。

鶏ちゃんだけでなくてチキンカツも入っている。量がかなり多くてかつ濃い味付けでご飯が進む一品でした。

奈良では鹿が赤信号も無視して道路を横断している、、、

奈良国立博物館に行ったときの写真。昔は基本、鹿は奈良公園にだけいた印象ですが、最近は人が多いところまで降りてきているんでしょうか。道路で、車がいようが赤信号だろうが関係なく歩いている。

なので人間の方が気を付けろと鹿に注意の看板が出ていた。

よみがえる正倉院宝物-再現模造にみる天平の技-@奈良国立博物館を観てきた。太古のすごい技術が現代でも作れるというのがうれしくなる。明日9/6まで。

行こう行こうと思っていたが、終了間近でようやく観に行けた。とても面白かった。

特に５弦の琵琶。材料から今では入手が困難（ウミガメの甲羅とか）なものをそろえ、科学的にも分析し現代に模造品を作るというもの。

しかももはや5弦の琵琶は現存してないので、調律すら不明だとか。

なら仏像館も同じチケットで入れるのでお得。明日までなのでぜひ。

よみがえる正倉院宝物-再現模造にみる天平の技-@奈良国立博物館を観てきた。太古のすごい技術が現代でも作れるというのがうれしくなる。明日9/6まで。

行こう行こうと思っていたが、終了間近でようやく観に行けた。とても面白かった。

特に５弦の琵琶。材料から今では入手が困難（ウミガメの甲羅とか）なものをそろえ、科学的にも分析し現代に模造品を作るというもの。

しかももはや5弦の琵琶は現存してないので、調律すら不明だとか。

なら仏像館も同じチケットで入れるのでお得。明日までなのでぜひ。

住吉公園をぶらぶら、、、だが暑すぎてすぐに退散

住吉大社の次は住吉公園へ、、、だがめちゃくちゃ暑い！自動販売機でペットボトルのスポーツドリンク飲んですぐに退散、、、、

住吉大社でお参り。暑かった、、、ただ近所のおもかる石はコロナの影響で参拝できず。

久しぶりに南海電車に乗って住吉大社駅まで。暑かった。

いつもこの橋を登るときこけそうになる、、、、

しかしここに来ると必ず行っていたおもかる石は中止。コロナの影響がこんなところにも。

阪急電車で日能研の理科の問題（スーパーボールを1mの高さから落として何回もバウンドする）、勝手に難しく考えて運動方程式と反発係数で、、、とか頭で解いてた。そんな必要全くなかった。

阪急電車に乗っているとこういう日能研の理科の問題が出ていた。

あれ？こんなの小学生解けるの？運動方程式と反発係数で、、、とか考えて1mから落とすと1/2m跳ね返るんだから

z軸を上方向として、運動方程式はスーパーボールの質量をm, 速度をvとして

m*dv/dt = -m*g

dv/dt=-g

v(t)=vo-g*t   --> v(0)=0なのでvo=0      つまりv(t)=-g*t

hoの位置から落とすとすると、これを積分して

z(t)=ho-1/2 *g*t^2

地面に落ちる時間はこれを0としてto=√(2*ho/g)

このときの速度は-g*toで、反発係数をrとすると地面で跳ね返ったときの速度は

v1=-r*(-g*to)=e*√(2*ho*g)となる。

これで打ち上げって最大の高さは速度が０になるところでその時間は

v(t1)=v1-g*t1=0         

t1=r*√(2*ho/g)となる。

このときの高さは

z(t1)=v1*t1-1/2*g*t1^2=r*√(2*ho*g)*r*√(2*ho/g)-1/2*g*r^2*2*ho/g = r^2*ho

これが最初の半分,1/2hoになるので

r^2=1/2

r=√(1/2)

になる、、、とここまで考えて、あれ？それだと次の反射で最高点は最初のhoを1/2hoに置き換えただけなんで1/4ho、、、

ああ、問題文にいつも落とした高さの1/2まで上がると書いてあった。反発係数とか運動方程式なんかまったくいらない、、、

頭が固くなっていた。

なので、2回で1/4m=25cm, 3回弾むと(1+1/2+1/2+1/4+1/4+1/8)=2+ 5/8 m = 200+12.5*5cm=262.5cm≒263cm

1cm以下になるのは

1回　50cm 2回　25㎝、3回　12.5cm　4回　6.25㎝　5回　3.125cm　6回　　1.5625㎝、、、なので7回目。

松屋で豚キムチ丼（大盛、半熟卵）を食す。辛さはちょうどよい感じ。ただちょっと塩気が強いか。

キムチが新しくなったということで食べてみた。確かにニンニクの芽とも相まって辛さ的にもちょうどいい感じになっている。豚もかなり大量ですが、どうも私には塩辛いかな、という感じ（私は唐辛子などの辛さはいいんですが塩の辛さに弱い、、、）

でも美味しいのは美味しかったですよ。半熟卵でいい感じにまろやかになる。辛いのが苦手な方は生卵でさらにまろやかにしたほうがいいかも。