プログラミング
ライプニッツの公式で円周率を求める話を書きました。www.omoshiro-suugaku.com 中学生の頃この話を聞いて、高校生になってから手に入れた電卓で円周率を求めようと頑張りました。ポケットコンピュータというのも現れ始め、これは当時で3万円くらい。何を使…
Python、グラフも色々描けそうです。比較のためなど、いっぺんにいくつか表示することもできます。パラメータはたくさんあるようですが、ここでは本当に基本的なものを。 #======================================================== #画面をいくつかの区画…
電話ので人の声を聞くと、ザラついていると言うか、いかにも「機械を通した声」と言うか、うまい表現の方法がこれ以上見つかりませんがそんな感じに聞こえます。とにかく直に聞いたのとは違いますよね。今回は肉声を電話を通したときの声に変換してみましょ…
音楽やその他の音が鳴っているとき、大きな音なら媒質（音を伝える媒体。空気や水など）は大きく、小さな音なら小さく揺れているはずです。wavファイルに記録するデータはそれをある程度反映しているのでしょう。媒質が変位していないならデータは０、変位し…
音階の話です。音の高さを少しずつ変えて音階ができあがります。音の周波数を少しずつ増やしてだんだん高くなるようにするのです。ラの音を440Hzとしましょう。まずラの音を鳴らすwavファイルを作ってみます。コードを載せておきますが、作り方は前にブログ…
ここまでで確認した事実を組み合わせて、wavファイルのノイズを消す実験をしてみましょう。まず、以下のコードでノイズ入りのwavファイルを作ります。wavファイルの作り方は過去の記事を参照してください。 www.omoshiro-suugaku.com #=====================…
もう少しFFTで遊んでみましょう。使う言語は例によってPython です。今回はノイズの混ざったデータから周波数成分を求めます。元のデータは一見グチャグチャなのですが、FFTで正しい周波数成分が分かります。こういうのを見るとFFT（あるいは数学）は凄いな…
waveファイルを作ってみました。必要なパラメータをそろえ、サウンドデータをこしらえて保存すればOKです。作曲などはできませんので、10秒間、単純に左右で異なる高さの音が鳴るだけのものにしておきます。単純な方が原理は分かりやすいので、ちょうどよい…
16ビットで録音されたwaveファイルでは（getparams()で得られる値がsampwidth=2）、サウンドデータは－32768～32767の範囲の整数値で表され、無音は0です。で、問題です。16ビットで録音されたwaveファイルで、サウンドデータの符号を全て変えたら（1234を－…
そろそろ「複数枚の写真から立体を再現」を実験してみたいと考えています。そのために必要そうなコードを整理しておきます。立体の散布図です。与えられた3次元の点をプロットします。次に(x,y)から式を使ってｚを計算し、空間内の点を求めてプロットしてみ…
タイトル通りです。Pythonでやってみました。FFT、IFFTを使います。まず、waveファイル(コードでは'BGM.wav')をオープンしてサウンドデータを読み取り、waveファイルのパラメータを表示します。表示を見てステレオ録音であることを確認してください（コード…
Pythonのmatplotlibでグラフを描画するのは便利で易しいのですが、普通にshow()で描画するとそこでいったん止まってしまいます。時間を追って変化するグラフは描けません。調べるとpause()というのがあるようで、試してみました。plt.show(0.1)の行です。描…
前回、FFTの例をお見せしました。さらに理解を確実なものにするため、もう少し別のパラメータで確認してみましょう。タイトルは「FFTとIFFT」となっていますが、今回はFFTのみで説明します。コードは以下の通りとします。 import numpy as npimport scipy as…
一部書き換え（グラフを含む）ました（2020.1.4）。 備忘録も兼ね、とにかく分かりやすい例を挙げます。説明も十分にしますので、少し前に書いた記事を参照しながらお読みください。 www.omoshiro-suugaku.com N＝100点のサンプリングデータを相手にするとし…
続きです。グラフの縦軸の目盛りに関しては別の機会に説明するとして、今回は特に横軸の目盛りの読み方について。まずは基本的な事実を確認しておきます。みなさん分かりきっているからなのか、あまりこういう形にまとめられた資料を見ることがありません。…
少し前、このブログでフーリエ変換について書きました。 www.omoshiro-suugaku.com 量が多くなりますので何回かに分けますが、実際にFFTのコードを書くにはどうしたらいいのか、出てきた結果をどう解釈したらいいのか、など具体的に書いてみます。備忘録を兼…
PCで顔認識するのは大変です。しかし色の組み合わせなら？ 例えば1000人に5ｃｍ×10ｃｍくらいのカードを持たせます。カードは4区画に分かれており、各区画には7色から選んだ色が塗られます。つまり、カードの色がPCに見分けられれば7の4乗、2401人の見分けが…
画像処理（1）につづき、もう少し紹介しましょう。PCで写真（画像データ）を処理し、埋もれてしまっていた情報を見やすくするなどします。 ●コントラストを変える近い明るさのドットばかりが揃った画像は見づらくなります。コントラストが低い画像、と言いま…
少し前、数字認識の実験をやっていました。数字認識を含む画像認識に興味があって勉強中です。画像認識とは………数字や文字をパソコンに見せて読ませたり、写真を見せて「写っているのは犬だ」などと判断させたりする技術です。画像認識はコンピュータには難し…
大学に入ってずいぶん経って先輩から使い古しのパソコンを買いました。PC-8001というやつです。もちろん当時のパソコンは今に比べればまるっきり非力で、たいしたことはあまりできませんでした。ソフトもほとんどなし。そもそも売られていないのです。いや、…
ステガノグラフィー（steganography）の話をしましょう。ステガノグラフィーとは、データを画像など、他のデータの中に隠す技術です。近い（？）ものに暗号というのがあります。暗号は文章を読めなくする技術、ステガノグラフィーは文章の存在自体を隠す技術…
wav（ウェーブ）ファイルというのがあります。MediaPlayerなどを使ってパソコンで再生できますよね。wavファイルは構造が決まっていて、ファイルの先頭付近に音楽データの長さや、サンプリング周波数、ステレオかモノラルかの区別のための数値などが入ってい…
昔のTVの画像はアナログ信号で送信されていました。そのため送信途中で何かエラーが起こっても（ノイズが入るとか）直す方法はありませんでした。今（地デジ）はデジタル信号で、ノイズが入っても、入っていない信号を再送させたりして訂正できます（コンピ…
昔、カセットテープというのがありました。古い「ウォークマン」、ご存じでしょうか。音楽をカセットテープに録音し、それをみんな聴いていたのです。写真がカセットテープです。エヴァンゲリオンで、碇(いかり)シンジ君がよく聴いていたのがこれ。 さて、何…
2036年問題について書きましょう。Windowsでは、Windows2000以降（もちろん今のWindows10なども）、インターネット上のタイムサーバー（あちこちにある。コンピュータで動く、これもソフトウェアです）を利用してパソコンの内部時計を定期的に正しい時刻に修…
今回はプログラミングのすすめ……ということで少し書いてみます。 「ゲームをするのが好き」も、いいことでしょう。ゲームをやらない人より、人生を面白くする何かを少なくともひとつ余分に知っているわけだから。しかし音楽だって聴くだけでなく自分で演奏す…
３角形が4角形の手前にあります。3角形は一部が4角形に隠されて見えます。コンピュータで描画する場合、隠れて見えないはずの線や面を描いてはいけません。隠線処理、隠面処理と言います。方法はいくつかありますが、すぐ思いつくのは視線の通る3角形上の点…
写真の1枚目と2枚目の写真、何が違っていると思いますか？ 2枚目は1枚目より見づらい感じです。でもぼやけてるわけではない。2枚目のような画像は「コントラストが低い」と言います。コントラストとは画像の最も明るい点と暗い点の明るさの差。平たく言うと…
100均の電卓で「2＋3×4＝」と打つと「2＋3×」の段階で「5」と表示され、続けて「4＝」と打つと「20」と出ます。数式の計算としては誤っていることに注意してください。正解は14です。×、÷は＋、－に優先して計算しなければならないのでした。しかし電卓は入…
専門ではないので説明が曖昧な部分もありますが、勘弁してください。ぼくはソフトを書くのが趣味なんですが、何年か前、インタプリタ作りにはまっていました。とにかく面白かったです。 例えばCという言語では、for(i = 1, s = 0; i <= 10; i++){ s += i;}と…
