|
20080511～
13と7と11の倍数の論理積は13と7と11の積の倍数である。
和ァ・・・
いーつもの君の笑顔がー(ライム胡椒で食べてねっ) すごくうれしいからー(ゆヵりーあなた憑かれてるのよ) もっと笑わせちゃえー(ブラッドユイー＋＋) そんなノリでいくよー(メソ・メソポタミア号) メソ何々ナシなー？  にほんブログ村 最初見たとき、むろみさんほどではないですが、どことなく９０年代っぽい顔だなーって思ったんですよ
PR
カラー荷 クォーク 混合 レプトン ニュートリノ 振動 カイラリティ ヘリシティ バリオン数 レプトン数 右巻き 左巻き
もちろんウンコ味のカレーを選ぶ。衛生的に問題ないし、数世代か数世紀経てば慣れるだろう。 ましてやカレーやウンコといったミクロの世界では、ウンコ味というものが定義できない。 もしかするとカレー振動でフレーバーが変わって、ウンコの匂いがしないウンコ味のカレーにできるかもしれないと期待している。 しかし、ウンコ味のカレーとカレー味のウンコが重ね合わされている可能性はないのか？ その可能性はない。 というのも、ウンコにはウンコ数という量子数があり、カレーにもカレー数という量子数があって 両者は保存するからである。 また、あんまり考えたくはないが、カレー振動と類似した現象として、ウンココンゴーという現象も存在する。 ウンコにはカレーにはない「カラーチャージ」が存在する。 したがって、青いカレーを見たらウンコだと思え。 カレーには質量がないとされていたが、近年質量の存在が確認され、光速で飛ばないことが確認された。 したがって、ごくわずかながらウンコを拭くほうの手でカレーを食べるという現象が存在するかもしれない  にほんブログ村
３３３＜俺たちを、誰だと思っていやがる
目がショボショボする上に目も全身も疲れていたのでメガネを外して仕事しましたきょう。 趣味でメガネをすることはあっても 今の仕事でメガネがいるということは５割もないです。 大丈夫です。たとえ設計・開発業務にいても、検査の業務にいても いつも仕事で使う関数は対数止まりでしたから・・・ 三角関数はまったく使いませんでしたから、メガネくらいどうってことないですよ・・・ははは アレルギー、というか秋の花粉症もあるかと思いますが これは風邪、ですねえ。ちょっと衣服をパージしすぎたようです。 ちょっと明日が割りとハドーワークなんですが、まあ、いつ寝るの？今でしょ！ ってな感じで寝ればなんとかならんこともない。たぶん みなさんの応援ありがとうございます 実は結構糧にしていただいてました。最近までほとんど自覚なかった！！  にほんブログ村
あばばばわわわわ・・・どうしよう何も出てこない！＞＜
このままじゃ今週の更新率がガタ落ちに・・・ とりあえず何も出てこないネタで１つ日記稼ぎをして・・・ってこれちょっと前にやった気がしたけど何日前だっけ え、えーとじゃあ、gifの元になったExcelファイルを、満を持してアップ・・・ 準備できてねぇー！これから準備する気もまったく起きねぇー！ すんません今日はなんもない日記で.TxT(今回のコレはジト目じゃないです) ところでVBAとマクロの本を、ようやく調べることにしました ようやく調べることもできました。 Excelのバージョンが古すぎて、最近の入門書は何が書いてあるのかよくわからないのです＞＜ 僕はもうダメです、あとのことは若者に任せてこのクソみたいな世界をもうすぐ去りますよー そんな中で奇跡的に見つけた、バージョン対応してる本 大学の図書館で見つけたんですよ あーこれでようやく習える～かもしんない それと、「数学ガール」の原作１巻もちょこっとだけ読んできました。 となり町のでっかい市立図書館にもあるにはあるんですがいつも貸出中で・・・ まあネットで予約ができるので予約してから借りれば車を何度も往復させる手間も省けていいのですが どうもしょっちゅう貸出中なんですよね それに対してこの大学の図書館は全巻揃ってる！ きっとあるなんて思ってもみないんだ！ あーでも夏休みが終わったら状況一変するかも知んない・・・ もう学生じゃないですしこの大学には通ったこともないので、貸出には許可証が必要でしょうねえ 何年か前に許可証を作ってもらった過去があるだけに、再発行的な意味でも不安かつ面倒くさいので なるべく館内で読んでしまおうとか思ってます。 家から自転車で行ける距離なので、ノート持参で大事なとこだけ４回くらい書き写しておきたいですね。マクロで。  にほんブログ村
①
ただ単に溜まっていたアニメを全部見終えたとか 溜めていた仕事が全部終わったとか そういうことで鬱要因が晴れ、 無自覚なうちに自由度ががらっと広大になっている 案外、そんなちっぽけな鬱要因を過大評価していることに気がつく、またはしばらく気が付かない ② 大自然に出たら、一人の人類の悩みなんてちっぽけなもんだよねー とは、たぶんまたちょっと違う感じ。 昔から思うんですけど、まったくちっぽけに思えない不思議。 いや、だって、そこでまた人類という存在がどうして存在するのかとか考えだしたらやっぱり人類すげーってなるじゃないですか。特に数学分野とか。 その人間様が考えてる悩みってやっぱ相当宇宙相当じゃないですか だめじゃん全然肩の荷軽くなんないじゃん ③ まあ、ループして困るってんならパラドックスしないループにすりゃいいんすよ 別に過去にさかのぼっても因果関係壊れない行動だってとりえるわけですし 鬱要因が鬱要因生むんだとしたら、考えないのも手ですよね 無知の知の逆といいますか、知の無知？いや意味があまり伝わる気がしない言葉だとは思いますけど。 なんだろうタイムパラドックスに関しては負帰還が正帰還っぽく、正帰還が負帰還っぽく見えるんですよね。 ①！≒？タイトル それとも、それともちょっと、ちょっと、もちょっとちがうんのじゃではまいかろうか ∋―8軒χ10 ＝人＝ いただきます  にほんブログ村
さらばー宇宙よー
いや～すっかり1ヶ月以上放置してしまいました＞＜ 先日まで、画像も出さずにコソコソと、波動方程式の数値シミュレーションを行って やったらボソボソとロビン境界条件の条件が求まらないとかなんとかブログに書いてきましたが 我が希望的観測日記ブログ枝部が目指しているのはこんなんじゃありませんのだ！ヽ｀д´ノ 僕がなぜ、Excelでのシミュレーションにこだわっているのか！？ それは！プログラミングが怖いからだ！ 複雑さを積層していくプログラミングのあり方にウマシカがあるからなんだよ！！！ソースコードらしき書体すらアレルギーなんだよ だからこその非積層型、いつでも気軽に組み立てられる平屋のシミュレーションをやってきたはずなのであって このブログで技術を積層化してしまったら意味がないではないかー！ ＝＝＝＝＝＝＝ 以前、「反射の観点から見た直流のインピーダンスマッチング」というテーマでですね、ブログを書いたことがあったんですよ。 なぜ、電源の内部抵抗と負荷の抵抗値を同じにしなくてはいけないのか？ 理由： 電源電圧Vと、内部抵抗r、負荷Rがただ直列につながった単純な回路で Rでの消費電力を最大にするには r=Rにおいて、 そのこころは Rで消費される電力PはRの関数にすると、 Pの逆数をRで微分して0になったところが消費電力の最大値なので をRについて解くと、R=r(R=-rなんて物理的にムシムシ！) 以上。 直流の観点から見ると実に簡単です。 しかし、どうして半分の電圧相当しか供給できないのか？ ということで、直流のスイッチを入れた瞬間、つまりステップ関数をめっちゃ細かく見たらどうなるだろう？ ちゅうのがコレだったのです。 負荷を鏡に見立てて 短絡だと反射波の符号が反転するから、負荷では節(ディリクレ境界条件)が維持され 開放だと反射波の符号が反転しないから、負荷では腹(ノイマン境界条件)が維持されるわけです。 この画像を作ったときは、まだ波動方程式のマックスウェル版といいますか、 シュレディンガーじゃないほうといいますか、そういった波動方程式の数値解がまだ用意できなかったんですね で、用意したのが以下になります。 ↑ディリクレ ↓ノイマン お気づきのように、進行波と反射波に分けることがまだできていないんですよ・・・ これを満たす境界条件はインピーダンス境界条件とか、ロビン境界条件とか呼ばれてるらしいんですが びっくりするほどこの境界条件の数値解析の具体例が見当たらないんですね。 そこで、先日の日記のように自力で解こうと思ったんです。 解いてみて思ったんですが、ロビン境界条件って需要がないから書籍にされないんじゃないでしょうか・・・ どうしても反射するなら、反射する手前で止めときゃいいじゃない的な。 十分に長い領域を確保して、端まで伝搬する前に、端っこの手前までを表示すればそれで済むことなんじゃないかな って思えてきたんですよねー・・・ 微分方程式を差分として解く際の境界条件は、波の振幅をuとすると ディリクレ条件は端っこでu=0 ノイマン条件は端っこでdu/dx=0なので、x=0を端に見立てると(u1-u0)/Δx=0なので、u1=u0 たったこれだけです。 ロビン境界条件にしても ディリクレとノイマンの線形結合なので、結合係数をa、bとして au+bu'=0 差分表現だと これだけです。 でもこのaとbがたいていは関数ってwikiに書いてるところが曲者なんですよねー おそらく、周波数成分に分けてからフーリエ級数展開しろって言ってるんだと思うんですよ 周波数ごとにaとbが変わってくるんですよねー じゃあ「周波数」と「領域の長さ」と境界条件の三つ巴の戦いになるじゃないですか それを解くために伝送線路やら分布定数線路の条件を解いてたわけです。 入力インピーダンスの規格化：zi=Zin/Z0 負荷インピーダンスの規格化：z=Zl/Z0とすると(特性インピーダンスで規格化すると大文字が小文字に！) 波数β=2π/λ、波長λと周波数fの関係は波の速度をcとして、c=fλ、jは虚数単位、Lは領域の長さです この式でマッチングが取れてる条件はzi=1なので分母＝分子が条件 (純拠数部分が同じ) となって これを満たすzとβはz=1(Zl=Z0)かつβl=nπ/4 条件は2つ ・ 領域の長さが8分の1波長の整数倍 ・負荷インピーダンスZl＝特性インピーダンスZ0 ってことになります。 これで周波数に当たりをつけることができたため、 あとは領域の長さと境界条件を戦わせるだけになりました。 ディリクレとノイマンの結合比は、aとbとはなっていますが片方を固定して片方を動かせばいいので、 領域とbを固定すればaだけ動かして見ればすむはずです。 1発だけのsin波をぶち込んで試してみたんですけどねえ 差分法の精度が緩すぎるのか、なんかそれっぽいところに行けばいくほど気持ち悪く波が乱れるんですよー でも、この領域の中途半端な部分を端っこに見立ててもいいんじゃないか？ってことがわかったっぽいので まあ有意義だったんでしょう。 計算は至ってコンパクトに収まっていますよ。 ほらこのように。 詳しくはコチラからExcelファイルDLでお願いシマス！TーT（ﾄﾞﾔｧ・・・ (循環参照を使っているのでツール→オプション→計算方法→反復計算→ON 最大反復回数=1、変化の最大値=0.001にしてください) 0～20までの21個の領域で、右端と左端に境界条件を設けてます。 縦は時間で、t=0と1が初期条件、スイッチを押した瞬間はt=0と1を参照したt=2と3で計算し あとはt=2と3を参照したt=4と5で2、3と4、5を循環参照させます。 ちなみにここで、ノイマン条件かディリクレ条件かを選べるのは右端だけです。 左端からはステップ関数が入ってるので、いわばステップ関数というディリクレ条件オンリーなわけですよ だから左端は必ずu=1の節にならざるを得ないというわけです。 ↓の周期が2往復分になっているのはそのためです。 この、進行波だけを見るためには、反射する前を見ればいいので、こうなります 横幅を狭めて見てるだけです。時間も途中まで。 r=R(電源の内部抵抗＝負荷抵抗) つまり、開放したときに反射して戻ってきたのも「込み」で、”電圧1”として見ていたから マッチングが取れた「進行波だけ」だと半分に見えていただけなのです。 まあ、ステップ関数だけだとつまらないと思いますので、正弦波入力もしてみましょう 左端にゼロ抵抗のLC共振回路があるような感じだと思ってください。 反射して定在波が立っているのがわかると思います。 ディリクレ(負荷で節：電圧だと短絡、負荷抵抗＝０) ノイマン(負荷で腹、電圧だと開放、負荷抵抗＝∞) その他条件：負荷電圧が節でも腹でもない(でもマッチング取れてるまでいかない) 結局、このシミュレーションで、ロビン境界条件の係数aとbが、どのように特性インピーダンスと関連しているのかはわからずじまいでした。 そもそも、どこに特性インピーダンスを設定する項目があったのか、自分で作っていてわからないんですよ・・・だから理論的にズボシーって華麗に決めることが出来ず、aを逐次変えて模索する他なかったのです。 クーラン数(ズ) ゲート  にほんブログ村
昨日は頭痛の中、波動方程式の数値シミュレーションにおけるインピーダンスマッチングを行っていてブログ更新とかもろもろ置いてけぼりになってしまいました。
波動方程式ですので当然波になるわけです。 そうするとおのずと反射が気になるわけですよ 進行波に反射波があると定在波になってしまうので この反射波をなくすために、境界条件を特殊なものにしなければならず 端っこがいわゆる波の「腹」になるノイマン条件と、波の「節」になるディリクレ条件の間を取ってやらんといかんわけです それが、位置による0F(ゼロ階)微分(ディリクレ)と1F(1階)微分(ノイマン)を線形結合したロビン境界条件というものでして(混合境界条件とは異なります) しかしながらこの結合の割合をどうしたらマッチング取れたシミュレーションになるのかが なかなか探してもいい資料に出会えないので 分布定数線路と伝送線路から当たりをつけてみることにしました。 おそらく、伝送線路の入力インピーダンスと負荷インピーダンスが特性インピーダンスと等しくなる条件でマッチングが取れるのだと思うので、 未知数であった「領域の長さ」「振動数」「条件の割合」の三角関係のうち前者2つの関係は解決しそうな感じです。(1/8波長) 問題は、シミュレーションの定数に、基本的には波の伝搬速度しか与えてないのに この系からどうやって特性インピーダンスが算出されるのか、です。 ここがよくわからないので、 とりあえずノイマンとディリクレの結合比を変えてトライアンドエラーしつつ、 最適解が見つかれば、それにあとからもっともらしい理屈付けをしてやる と、そんな感じになればいいなと思っています。 Excelの循環参照と差分法の陽解法を使ってきましたが、マッチングが取れる界隈に近づく(たぶん)ほど波がぐちゃぐちゃになるように思えるのは必然なのか偶然なのか これは、循環参照と解除して、陰解法に手を出したほうがいいかもしれませんねえ 精度が足りないような気がします。 ところで、ロビン境界条件の使い道はどうも、波動方程式よりも熱の伝達、とりわけニュートンの冷却の法則界隈みたいですね 波動方程式は解析的に解けるから用済みといった感じなのでしょうか ニュートンの冷却の法則は以前、冷蔵庫の消費電力の測定絡みで少しだけかじったことがあるのでまだいいのですが 熱力学全般苦手です＞＜ なので熱力の言葉で書かれると何が書いてあるのかよくわかりません そういえば 今は入力信号をパルス波ではなく正弦波にして、あとからフーリエ級数展開で合成してパルス波にしようと思っているのですが 元々フーリエ級数展開が生まれたきっかけも、電気ではなく熱力学からの需要だったようですね 文字ばっかりなのもアレなんで 以前のドラクエドーナツの日記に追記した画像を上げておきますね  にほんブログ村
先日の日記にも書きましたが
無限広角式メルカトル４次元遠近法という描画方法を感覚的に作り出せないものかと思いましてね MMDで背景だけ持ってきて、カメラ回して動画キャプチャーしたのをトリミングしながらつなげてみたんですよ これをゆくゆくは、四畳半の部屋でやりたいのです。 壁からどのくらいの距離に立ち、どのくらいズームしてカメラ回しながら撮影するべきですかねえ そして部屋は自作か、他作か。 繋ぎ目をよくしようと思えば、やはり枚数もそれなりに増えてきますし ここはいっちょ、マクロとやらに手を出してみましょうかねえ とすると、トリミング範囲や個々の画像同士の間隔をもちょっと定量的にする必要がありそうですねえ 部屋をどのくらいの大きさにするべきか、とか  にほんブログ村 丸いものが描けないので長らく丸いものの描画をコンピュータに任せてきましたが、ここまできたらもう意地でもコンピュータさんに描いてもらいますよ！
左：ファミレスの上司(行方不明)、右：公務員の上司(ぬいぐるみ) サービスワーキング サーバントはどこかへ行きました fate/zero  にほんブログ村 なんてこった！鯖が「ルーシーのおっぱい」を侵略にきてやがる！これは確実に獲物を狩る野獣の目ですぜ
公務員だって、生活を守れるんだ！要出典[誰の？] 実際踊らせてみると分かるが、普通の巨大 バニー＠ジャック・バウアー  にほんブログ村
数値解析の本をいくつか読んでます。
「数値解析」と名のつく本でも結構色々バリエーションありますね 最初から偏微分方程式の数値解を扱ってるものもあれば ニュートン法や二分法、数値積分のやり方から始めている本もあって 戸惑いながらも納得のバリエーションでした。 偏微分方程式の数値解法は僕は素直に差分法しかしらなかったのですが 名前だけ聞いたことが合った有限要素法のほかに、境界要素法なんてのもあるようですね 表面だけを扱って次元を1つ落としてシミュレーションができるという強みはあるものの はたから見ると弱点ばかりのように見えなくもないです＾＾； ちょうど先日計算していたフラーレンを思い出しました。 なんというか、計算してみればあっさり終わってしまった印象があったのです。 なんというか思ったより浅かったというか。 それもそのはず、このような多面体には文字通り表面しかなかったわけです。どおりで薄っぺらいはずですよ ところで、差分法について調べている最中、どうして行列を駆使しているんだろう？ という文脈がどの本にも出てきていて驚きました。 僕の知っている差分法と違う。 と思ったら、差分法には陽解法と陰解法があって 僕はその陽解法しか知らなかったらしいんですね おそらく、おそらくですが 差分法にした時点で微分方程式が連立方程式になるわけですよ それを逐次ではなく行列方程式にしていっぺんに解いているのが陰解法なのでしょう。 まあ、当面は差分法の陽解法にしか興味出なさそうな気がします。はい  にほんブログ村
なんでかは知らないことにしよう
（ホントに知らないんです。しるか！）bps 秋月 超法規的措置 ところで7話の8話用次回予告クイズは 満を持してサバじゃねえ！！！！１というところですよ。背泳ぎしか泳げません 準備はいいですか 今言わずしていついうのですか！？e-maですよねえ！？シャンゼリオン  レッドさんとヴァンプ将軍 って、このままだと今週末のニコチャンに上がるまで話題に上がらないかもしれない・・・orz 僕は好きですよ、このアニメも長谷部も。 日本人はもっと長谷部であるべきですよ。もっとこうノーテンキラキラで。 にほんブログ村
|
カレンダー
カテゴリー
ブログランキング
ブログランキング参戦中
よかったらポチッとお願いします＾＾
最新CM
[12/30 buy steroids credit card]
[09/26 Rositawok]
[03/24 hydraTep]
[03/18 Thomaniveigo]
[03/17 Robertaverm]
最新記事
(01/01)
(01/03)
(09/23)
(09/23)
(02/11)
(05/30)
(05/28)
(05/28)
(05/27)
(08/04)
(10/24)
(06/08)
(05/22)
(01/13)
(11/04)
最新TB
プロフィール
ＨＮ：
量子きのこ
年齢：
43
ＨＰ：
性別：
男性
誕生日：
1981/04/04
職業：
WinDOS.N臣T
趣味：
妄想・計算・測定・アニメ
自己紹介：
日記タイトルの頭についてるアルファベットは日記の番号です
26進数を右から読みます 例：H→7番目、XP→15(P)×26＋23(X)=413番目。 A=0とする仕様につき一番右の桁はAにできませんのでご了承くださいｽﾞｺｰ
バーコード
ブログ内検索
アーカイブ
最古記事
(05/11)
(05/11)
(05/13)
(05/13)
(05/13)
(05/13)
(05/13)
(05/13)
(05/14)
(05/14)
(05/14)
(05/14)
(05/16)
(05/16)
(05/16)
アクセス解析
|
