タグ: 手習ひデジタル信号処理

Scilabの基礎もないところに泥縄式でいろいろやってみているので、毎度「よくわからんな」が続出。砂上の楼閣状態です。伝達関数を記述するときに「お世話になる」polyという関数もそう。z=poly(0, “z”)みたいな「おまじない」を唱えること度々。どうも多項式を定義するためのものらしいんだけども。。。

“手習ひデジタル信号処理(73) Scilab、今頃になってpolyの使い方を知る” の続きを読む

前々回、画像用のフィルタ関数があるのに「わざわざ普通のコンボルーション処理用の」conv2関数で画像をフィルタしてみました。今回は、本物の画像用の関数を使ってみたいと思います。冒頭のお猿様の画像、ちょいといつもよりボケてないっすか？ガウシアンフィルタの賜物？であります。でもさすがScilab、１つじゃないんだ関数が。

“手習ひデジタル信号処理(72) Scilab、画像用のフィルタ関数を使ってみる” の続きを読む

前回画像データを扱って驚愕したのは、Scilabの整数型から実数型への変換方法を知らなかった、という間抜けぶりです。慌てふたき無理やりな変換で誤魔化しました。今回はこの基本的な部分をちゃんと調べておきたいと思います。といって調べりゃ分かるのよ。ついでに変数のタイプをわかりやすい？文字列で出力する関数も作ってみました。
“手習ひデジタル信号処理(71) Scilab、整数型と実数型間の変換、そんなことも知らずに” の続きを読む

前回、画像入力ができるようになったので、今回はわざわざ信号処理用のコンボリューション関数を使ってそれにフィルタをかけてみます。「わざわざ」というのは画像処理用のモジュールに画像専用のフィルタ関数があるのに使わない、ということですが。「コンボリューション」するためにもいろいろある、と。 “手習ひデジタル信号処理(70)Scilab、わざわざconv2でソーベルフィルタをかけてみる” の続きを読む

前回に引き続き、Scilabの使い方編続きです。2次元（画像）の信号処理をしようとして気づきました。素のScilabには画像入出力の関数が見あたりません。しかし「充実したモジュール」群から適切なものをロードしてやれば画像入出力どころか何でもできる？っと。今回はScilabのモジュール管理 ATOMS を使ってみます。 “手習ひデジタル信号処理(69) Scilab、ATOMS、モジュール管理” の続きを読む

今回は信号処理ネタではなくScilabの使い方です。Scilabのスクリプトをファイルとして格納するとき、拡張子をsceとしたらよいのか、sciとしたらよいのか理解しておらず、「とりあえず」ということでsce一辺倒でした。しかし調べてみたらsciにする御利益があるのでした。これだから素人はイケない。 “手習ひデジタル信号処理(68) Scilab、sce拡張子とsci拡張子の使い分けの件” の続きを読む

前回Scilab様にIIRフィルタ設計をお願いできるようになったので、今回は実際にテキトーな信号をフィルタしてみるべし、と思い立ってハマりました。実際にフィルタ処理で気になっていたことがあったのだけれど、それ以前にハマる原因が潜んでました。Scilab恐るべし？って、なんだかな～。

“手習ひデジタル信号処理(67) Scilab、IIRフィルタでフィルタしようとしてハマった件” の続きを読む

前回はFIRローパスフィルタの「設計用」関数に「一皮」かぶせてみました。今回はIIRローパスフィルタと行きたいデス。今回は比較対象というか、目標を掲げました。以前、勝手に手習ひさせていただいていた三上直樹先生著の御本の付録のIIRフィルタ設計用関数です。「おんなじ」ような感じにいたしとうございます。 “手習ひデジタル信号処理(66) Scilab、IIRフィルタ設計用関数、周波数特性表示付” の続きを読む

前回はFIRローパスフィルタの係数を求め、求めた係数を使って「単位」インパルス入力信号をフィルタして時間応答を求め、さらにそれをFFTして周波数特性のプロットをいたしました。まどろっこしいデス。フイルタ係数を求めるときに「ついでに」周波数特性のプロットも見て、それからOKしたいです。知らんけど。今回は微妙な改善？ “手習ひデジタル信号処理(65) Scilab、FIRフィルタの係数求めるついでに周波数特性” の続きを読む

前回は差分方程式をループで回すという「原始的」な方法で１次のIIRフィルタの時間応答を計算してみました。今回はScilabのネイティブ機能を使って、100次のFIRフィルタの応答を「１撃」で計算してみたいと思います。100次のフィルタの係数の計算もScilabに「お任せ」。中身はよくわからないけれど楽ちん。いいのか？
“手習ひデジタル信号処理(64) Scilab、FIRフィルタの係数求めてフィルタしてみるの回” の続きを読む

前回、前々回とScilab上で伝達関数から周波数応答を求めるのを練習していました。しかし、そういえば、ということで考えてみると差分方程式が与えられているのであれば、差分方程式から直接、時間応答が求まるのを練習してないじゃん、と。今回はそこんトコロを練習してみます。

“手習ひデジタル信号処理(63) Scilab、差分方程式を直接計算、１次IIRフィルタ” の続きを読む

正規化周波数表示はカッコいいけど実際の周波数でないと実感がわかないなどとブーたれて、前回は伝達関数からゲインプロットを実周波数で表示する関数を作ってみました。今回は、位相プロットで実周波数表示を行う関数を作ってみたいと思います。素人が良く理解しないまま試行錯誤しているのだけれども大丈夫か？自分。 “手習ひデジタル信号処理(62) Scilab、伝達関数H(z)から位相プロット、Fs対応版” の続きを読む

前回は暫定版ですがFFT表示用の関数を作って多少お楽に。今回は以前から忸怩たる思い？をしてきた伝達関数からゲイン線図をプロットする方法を「改良」したいと思います。いつも正規化周波数に対するプロット。カッコイイお作法ではあるものの、具体的な周波数が入らないとイメージがわかない年寄です。実サンプリング周波数を与えたいです。

“手習ひデジタル信号処理(61) Scilab、伝達関数H(z)からゲイン線図、Fs対応版” の続きを読む

前回入力信号を生成する関数群を作成してみました。前々回に時間領域で離散時間っぽいプロットはできるようになってます。今回は、周波数領域でプロットする関数です。Scilabには部品となる関数は用意されているのですが、これ使えばいい的なものはないみたいっす。目標はAnalog Dicovery2の周波数プロット画面ですが。

“手習ひデジタル信号処理(60) Scilab、FFT表示用関数「とりあえず版」追加” の続きを読む