未分類 – 計測・解析ラボ

好きな山域を3Dプリントしよう

山の立体ディスプレイ
遠くから眺めてよし、頂を目指すもよし。頂から遠くの山を眺めるもよし。
思い出の山、これから登る山を3Dプリントしてみましょう。地理院地図の3Dツールを使って好きな山域のSTLファイルを作成します。STLファイルは小さな３角形のメッシュで３次元形状を表現するファイル形式で、多くのアプリケーションで共通に使われています。地理院地図の3Dツールは緯度、経度に沿った長方形の範囲の標高データをSTLファイルで出力します。そのまま3Dプリントしても良いのですが、手に持って色々な方向から山域を眺めるには外形が円形の方が便利ですから、ここでは外径115mmのディスク状にする方法を説明します。
黒部渓谷・立山・剱岳を含む山域の3Dプリントの例です。PLAクリアフィラメントでプリントしてフルカラーLEDで下から照らすと立体感が増します。

大まかな作成手順
(1) 地理院地図の3Dツールを使って立体ディスプレイを作る山域（正方形）のSTLファイルを作成します。
(2)Fusion360を使って直径115mmの円筒状にくり抜くツールとなるSTLファイルを作成します。
(3) Fusion360に山域のSTLファイルを読み込んで縮尺を調整します。
(4)くり抜きツールのSTLファイルも読み込んで、山域メッシュとツールメッシュの両方に含まれる部分を結合ツール（交差）で取り出してSTL形式で保存します。
(5) 3Dプリンターでプリントします。

(1) 地理院地図の3Dツールで山域のSTLファイルを作成
まず地理院地図で作成する山域を表示します。ここでは富士山で説明します。
「作図・ファイル」ツールでプリントしたい範囲に円を描きましょう。
3Dツールでは3D化する範囲の指定方法として「大(縦横2048px)」「小(縦横1024px)」「カスタム」がありますので、「カスタム」を選び、四隅の選択点をドラッグして先ほど描いた円に外接する正方形を指定して「OK」ボタンを押します。
ブラウザに新しいタブが作られて指定した範囲の３次元の地図が表示されます。
画面の右下端に3種類のファイルのダウンロードボタンが表示されますので、STLファイルをダウンロードします。画面の立体表示には等高線などが描かれていますが、STLファイルには描かれません。dem.stlという名前になっていますので、MtFuji.stlなど分かりやすい名前に変更すると良いでしょう。このSTLデータは一辺（長方形の場合は長辺）が150mmのデータになっています。

(2) 円筒形状のツールSTLファイルの作成（2回目以降は不要）
Fusion360のソリッドモードで、左上のブラウザの原点タブにあるXY面を右クリックしてスケッチを作成します。原点を中心として直径115mmの円を描いてスケッチを終了します。スケッチで描いた円をクリックして押し出しツールで115mmの長さの円筒を作成します。ブラウザのボディタブで円筒のボディを右クリックして「メッシュとして保存」を選択します。Cylinder.stlとして保存し、このデザインを閉じます。

(3) 山域STLファイルの読み込みと縮尺変更
新規デザインで新しいページを開きます。「挿入」メニューの「メッシュを挿入」を選択し、MtFuji.stlを開きます。この時、「上方向を反転」ボタンを押して、Y軸とZ軸の方向を反転します。「中心」ボタンを押して、メッシュデータの中心を原点に移動し、「地面に移動」ボタンを押してデータの底面がXY面と一致させます。
メッシュモードを選択し、「修正→メッシュを尺度指定」を使って一辺150mmの正方形の範囲のデータを一辺115mmに縮小します。縮小してもデータの底面がXY面になるように、「点」（縮小拡大の基準点）に底面の中心を指定します。視点を「下」にすれば選択しやすくなります。「エンティティ」は縮小拡大するメッシュを指定します。150mmを115mmに縮小したいので縮小倍率は0.77（≒115/150）と入力して「OK」ボタンを押します。データの底面がXY面にあることを確認してください。
(4) ツールSTLファイルの読み込みと結合ツール（交差）
Cylinder.stlを読み込みます。メッシュモードの結合ツール（交差）を選択し、「ターゲットボディ」にMtFujiメッシュボディ、「ツールボディ」にCylinderメッシュボディを指定します。もしも厚みが大きくてプリント時間が余分にかかりそうな場合は、結合ツールを実行する前にMtFujiメッシュボディをZ軸のマイナス側に移動します。下げすぎて底に穴を開けないように注意します。
作成したメッシュボディを右クリックし、「メッシュとして保存」を選んで、例えばMtFujiDisk.stlのようにファイル名をつけて保存します。

(5) 3Dプリンターでプリント
作成したSTLファイル（たとえば、MtFujiDisk.stl）を3Dプリンターでプリントします。

ESP32でI2Sを使ってみる

ESP32はWiFiやBluetoothが使える高機能なボードです。I2S（Inter-IC Sound）というサウンドをデジタル通信できる機能も持っていて、I2S対応のマイクモジュールやアンプモジュールを接続してArduino IDEで制御することができます。

「ESP32のプログラムガイド」や「ESP32のI2Sマイク研究」、「I2Sマイク、ESP32、Arduino、SPH0645LM4H」などの記事を参考にしながら試してみました。

結論から書くと、マイクからの入力をイヤフォンで聞くことができましたが、音質がどうこうというレベルにはなっていません。タイムラグをまったく感じないのは今後に期待できるかもしれません。

興味のある方への参考と、自分のための備忘録

ESP32は「 MH-ET LIVE MiniKit」を使用しました。マイクは「Sipeed 1マイクモジュール」２個、アンプは「GY-PCM5102」を使用しました。

ESP32のI2Sは２本あり、入力か出力どちらかを選択できます。今回はI2S_0をマイク入力、I2S_1をアンプ出力に使います。それぞれバッファメモリが使えますので、マイクの入力バッファからデータを読んで、32ビット符号付整数に変換して、処理（たんなる掛け算）して、バイト配列に戻してから、アンプ出力のバッファに書き込みます。delay()を入れずに繰り返して呼び出していますが、アンプ出力を見る限りではバッファがうまく丸く収めてくれているようです。

マイクモジュールの配線
DA — IO18(1kΩでGNDに接続)
CK — IO23
LR — L:3.3V, R:GND
WS — IO19
V — 3.3V
G — GND

アンプモジュールの配線
SCK — GND
BCK — IO16
DIN — IO21
LCK — IO17
GND — GND
VIN — 3.3V

#include <driver/i2s.h>

#define I2S_NUM_MIC                     I2S_NUM_0
#define I2S_PIN_CLK_MIC                 23
#define I2S_PIN_WS_MIC                  19
#define I2S_PIN_DOUT_MIC                -1//not use
#define I2S_PIN_DIN_MIC                 18

#define I2S_NUM_AMP                     I2S_NUM_1
#define I2S_PIN_CLK_AMP                 16
#define I2S_PIN_WS_AMP                  17
#define I2S_PIN_DOUT_AMP                21
#define I2S_PIN_DIN_AMP                 -1//not use

#define I2S_SAMPLE_RATE                 44100 //44.1kHz
#define I2S_BUFFER_COUNT                4
#define I2S_BUFFER_SIZE                 512
#define Sound_LEN                       512

uint8_t Mic_Buffer[I2S_BUFFER_SIZE];    // DMA転送バッファ
uint8_t Sound[Sound_LEN];               // サウンドデータ

int L_Data, R_Data;
int Volume = 8;

void setup() {
  Serial.begin(115200);
  i2sInit_MIC();
  i2sInit_AMP();
  Serial.println("MicModAmp02.ino");
}

void loop() {
  mic2amp();
}


void mic2amp(){
  size_t transBytes;
  i2s_read(I2S_NUM_MIC, (char*)Mic_Buffer, I2S_BUFFER_SIZE, &transBytes, portMAX_DELAY);// Mic In
  //Convert 8 bytes data to L_Data and R_Data
  for(int i=0;i<transBytes;i=i+8){
    for(int j=0;j<4;j++){
      L_Data=L_Data<<8;
      L_Data=L_Data + Mic_Buffer[i+j];
    }
    for(int j=4;j<8;j++){
      R_Data=R_Data << 8;
      R_Data=R_Data + Mic_Buffer[i+j];
    }
    //Process L_Data and R_Data
    L_Data=L_Data * Volume;
    R_Data=R_Data * Volume;
    //Convert to 8 bytes data
    for(int j=0;j<4;j++){
      Sound[i+3-j]=(byte)(L_Data & 0xFF);
      L_Data=L_Data >> 8;
    }
    for(int j=0;j<4;j++){
      Sound[i+7-j]=(byte)(R_Data & 0xFF);
      R_Data=R_Data >> 8;
    }
  }
  i2s_write(I2S_NUM_AMP, (char*)Sound, I2S_BUFFER_SIZE, &transBytes, portMAX_DELAY);// Speak Out
}


void i2sInit_MIC() {
  i2s_config_t i2s_config = {
    .mode                 = (i2s_mode_t)(I2S_MODE_MASTER | I2S_MODE_RX),
    .sample_rate          = I2S_SAMPLE_RATE,
    .bits_per_sample      = I2S_BITS_PER_SAMPLE_32BIT,
    .channel_format       = I2S_CHANNEL_FMT_RIGHT_LEFT, //stereo
    .communication_format = I2S_COMM_FORMAT_I2S,
    .intr_alloc_flags     = ESP_INTR_FLAG_LEVEL1,
    .dma_buf_count        = I2S_BUFFER_COUNT,
    .dma_buf_len          = I2S_BUFFER_SIZE,
    .use_apll             = false,
    .tx_desc_auto_clear   = false,
    .fixed_mclk           = 0
  };
  
  i2s_pin_config_t pin_config = {
    .bck_io_num           = I2S_PIN_CLK_MIC,
    .ws_io_num            = I2S_PIN_WS_MIC,
    .data_out_num         = I2S_PIN_DOUT_MIC,
    .data_in_num          = I2S_PIN_DIN_MIC,
  };
 
  i2s_driver_install(I2S_NUM_MIC, &i2s_config, 0, NULL);
  i2s_set_pin(I2S_NUM_MIC, &pin_config);
}

void i2sInit_AMP() {
  i2s_config_t i2s_config = {
    .mode                 = (i2s_mode_t)(I2S_MODE_MASTER | I2S_MODE_TX),
    .sample_rate          = I2S_SAMPLE_RATE,
    .bits_per_sample      = I2S_BITS_PER_SAMPLE_32BIT,
    .channel_format       = I2S_CHANNEL_FMT_RIGHT_LEFT, //stereo
    .communication_format = I2S_COMM_FORMAT_I2S,
    .intr_alloc_flags     = ESP_INTR_FLAG_LEVEL1,
    .dma_buf_count        = I2S_BUFFER_COUNT,
    .dma_buf_len          = I2S_BUFFER_SIZE,
    .use_apll             = false,
    .tx_desc_auto_clear   = false,
    .fixed_mclk           = 0
  };
  
  i2s_pin_config_t pin_config = {
    .bck_io_num           = I2S_PIN_CLK_AMP,
    .ws_io_num            = I2S_PIN_WS_AMP,
    .data_out_num         = I2S_PIN_DOUT_AMP,
    .data_in_num          = I2S_PIN_DIN_AMP,
  };
 
  i2s_driver_install(I2S_NUM_AMP, &i2s_config, 0, NULL);
  i2s_set_pin(I2S_NUM_AMP, &pin_config);
}

100

101

102

103

104

105

106

107

108

109

110

111

112

113

114

115

116

117

118

#include <driver/i2s.h>

#define I2S_NUM_MIC I2S_NUM_0

#define I2S_PIN_CLK_MIC 23

#define I2S_PIN_WS_MIC 19

#define I2S_PIN_DOUT_MIC -1//not use

#define I2S_PIN_DIN_MIC 18

#define I2S_NUM_AMP I2S_NUM_1

#define I2S_PIN_CLK_AMP 16

#define I2S_PIN_WS_AMP 17

#define I2S_PIN_DOUT_AMP 21

#define I2S_PIN_DIN_AMP -1//not use

#define I2S_SAMPLE_RATE 44100 //44.1kHz

#define I2S_BUFFER_COUNT 4

#define I2S_BUFFER_SIZE 512

#define Sound_LEN 512

uint8_t Mic_Buffer[I2S_BUFFER_SIZE]; // DMA転送バッファ

uint8_t Sound[Sound_LEN]; // サウンドデータ

int L_Data, R_Data;

int Volume = 8;

void setup() {

Serial.begin(115200);

i2sInit_MIC();

i2sInit_AMP();

Serial.println("MicModAmp02.ino");

}

void loop() {

mic2amp();

}

void mic2amp(){

size_t transBytes;

i2s_read(I2S_NUM_MIC, (char*)Mic_Buffer, I2S_BUFFER_SIZE, &transBytes, portMAX_DELAY);// Mic In

//Convert 8 bytes data to L_Data and R_Data

for(int i=0;i<transBytes;i=i+8){

for(int j=0;j<4;j++){

L_Data=L_Data<<8;

L_Data=L_Data + Mic_Buffer[i+j];

}

for(int j=4;j<8;j++){

R_Data=R_Data << 8;

R_Data=R_Data + Mic_Buffer[i+j];

}

//Process L_Data and R_Data

L_Data=L_Data * Volume;

R_Data=R_Data * Volume;

//Convert to 8 bytes data

for(int j=0;j<4;j++){

Sound[i+3-j]=(byte)(L_Data & 0xFF);

L_Data=L_Data >> 8;

}

for(int j=0;j<4;j++){

Sound[i+7-j]=(byte)(R_Data & 0xFF);

R_Data=R_Data >> 8;

}

i2s_write(I2S_NUM_AMP, (char*)Sound, I2S_BUFFER_SIZE, &transBytes, portMAX_DELAY);// Speak Out

}

void i2sInit_MIC() {

i2s_config_t i2s_config = {

.mode = (i2s_mode_t)(I2S_MODE_MASTER | I2S_MODE_RX),

.sample_rate = I2S_SAMPLE_RATE,

.bits_per_sample = I2S_BITS_PER_SAMPLE_32BIT,

.channel_format = I2S_CHANNEL_FMT_RIGHT_LEFT, //stereo

.communication_format = I2S_COMM_FORMAT_I2S,

.intr_alloc_flags = ESP_INTR_FLAG_LEVEL1,

.dma_buf_count = I2S_BUFFER_COUNT,

.dma_buf_len = I2S_BUFFER_SIZE,

.use_apll = false,

.tx_desc_auto_clear = false,

.fixed_mclk = 0

};

i2s_pin_config_t pin_config = {

.bck_io_num = I2S_PIN_CLK_MIC,

.ws_io_num = I2S_PIN_WS_MIC,

.data_out_num = I2S_PIN_DOUT_MIC,

.data_in_num = I2S_PIN_DIN_MIC,

};

i2s_driver_install(I2S_NUM_MIC, &i2s_config, 0, NULL);

i2s_set_pin(I2S_NUM_MIC, &pin_config);

}

void i2sInit_AMP() {

i2s_config_t i2s_config = {

.mode = (i2s_mode_t)(I2S_MODE_MASTER | I2S_MODE_TX),

.sample_rate = I2S_SAMPLE_RATE,

.bits_per_sample = I2S_BITS_PER_SAMPLE_32BIT,

.channel_format = I2S_CHANNEL_FMT_RIGHT_LEFT, //stereo

.communication_format = I2S_COMM_FORMAT_I2S,

.intr_alloc_flags = ESP_INTR_FLAG_LEVEL1,

.dma_buf_count = I2S_BUFFER_COUNT,

.dma_buf_len = I2S_BUFFER_SIZE,

.use_apll = false,

.tx_desc_auto_clear = false,

.fixed_mclk = 0

};

i2s_pin_config_t pin_config = {

.bck_io_num = I2S_PIN_CLK_AMP,

.ws_io_num = I2S_PIN_WS_AMP,

.data_out_num = I2S_PIN_DOUT_AMP,

.data_in_num = I2S_PIN_DIN_AMP,

};

i2s_driver_install(I2S_NUM_AMP, &i2s_config, 0, NULL);

i2s_set_pin(I2S_NUM_AMP, &pin_config);

}

Sendai Micro Maker Faireに出展します

仙台マイクロメイカーフェアに出展します。( #SendaiMMF )

2020年1月25日 12:00-17:00
産業見本市会館(仙台市) サンフェスタ301

いろいろな素材で作ったスクエアキューブパズルの紹介とLabVIEW 2014 Home版で書いた解探索プログラムのデモを行います。

リーフレット

LabVIEWユーザー会(2019)のお知らせ

例年NI Daysで開催されてきたLabVIEWユーザー会ですが、今年は日本でのNI Daysがなくなったので、独自開催を計画中です。
日時は2019年12月14日（土曜日) 13:00からを予定しています。場所は未定(浜松町周辺が候補)で、11月15日ぐらいまで参加者の人数を把握して決定するようです。会費は無料ですが、会場費を割り勘で負担する予定です。
予定しているプログラムなど詳細は下記URLのEventsから確認してください。
https://jlabcon.wixsite.com/jlabcon-top

私は、NIのディスカッションフォーラム「北東北ユーザー会」で紹介したIndustrial ShieldsのM-DuinoとSIEMENSのSIMATIC IOT2020を持参します。また、Arduino UNOの互換機でFPGA機能があるXLR8も、使いこなせていませんが、持参します。
また、テトロキューブ16個で4x4x4の立方体をくみ上げるパズルのLabVIEW版を作成しているのですが、自動的に解を探索するLabVIEWプログラムを（出来たところまでですが、）紹介します。画面上でパズルで遊ぶ部分はできましたのでNIのディスカッションフォーラム「北東北ユーザー会」の別トピックスで紹介します。

LabVIEWユーザー会に関する問い合わせは上記URLからお願いします。

FPGAベースのArduino互換機”XLR8″のアナログ入力について

<XLR8について>
Alorium Technologyの”XLR8″(写真下)はIntel MAX10 FPGAにATmega328互換機能を書き込んだArduino UNO互換ボードです。ユーザーがFPGAのXcelerator Blocksというエリアに機能を追加することができます。FPGAプログラムはライブラリとして提供されているものを使うこともできますし、独自のプログラムを作成することもできるそうです。

<購入先の例>
https://www.digikey.jp/product-detail/ja/alorium-technology-llc/XLR8R22M08V5U0DI/2003-1000-ND/9607423

<はじめに>
XLR8のアナログ入力は12ビットで、しかも10ビットのArduino UNOよりも高いサンプリングレートで使うことができるということなので興味を持ちました。また、Arduino UNOは高いサンプリングレートで使用した時に誤差が大きくなるということですので、サンプリングレートと同時に測定データの誤差を視覚的に把握したいと思いました。
https://www.aloriumtech.com/adc-12/

<調査方法>
アナログ入力で得られたデータをいったん配列(500要素)に保存して、データがすべて書き込まれた後でシリアル出力でPCに出力するスケッチを作成しました。受け取ったデータはArduino IDEのシリアルプロッタでグラフに表示させます。
アナログ信号としてはファンクションジェネレータの三角波を入力しました。三角波の周波数を調整してシリアルプロッタに表示されたデータが目視で１周期となる周波数を求め、概略のサンプリングレートの値を計算しました。

<分かったこと>
私なりの結論ですが、XLR8のアナログ入力はUNOに比べて分解能が4倍で、サンプリング周波数が3.7倍程度高く(最大250KHz)まで使えそうなので一考に値すると思います。

<具体的な調査方法>
Arduino UNOのアナログ電圧のサンプリングレートは、ADCクロックの分周比で設定できます。下記スケッチのADCSRAの下位３ビットに”001″から”111″を書き込みます。デフォルトの分周比は128(“111”)で9.6kHzのサンプリングレートです。分周比が小さいほどサンプリングレートが高くなりますが、誤差も大きくなるそうです。三角波を入力することで誤差の程度を視覚的に判断したいと思います。三角波はおよそ0-3Vです。

Arduinoのスケッチは要素数500個の配列に測定データを記憶させて、測定が終了してからシリアルプロッタでデータをグラフに表示することにしました。

const int COUNT=500;
int Data[COUNT];

void setup() {
  Serial.begin(9600);

  ADCSRA = ADCSRA & 0xf8;   // 分周比を決めるビット(ADPS2:0)を000へ
  //ADCSRA = ADCSRA | 0x04;   // 分周比16(100)をセット
  //ADCSRA = ADCSRA | 0x07;   // 分周比128(111)をセット

  ADCSRA = ADCSRA | 0x03;   // 分周比__(xxx)をセット

}

void loop() {
  for(int i=0;i<COUNT;i++){
    Data[i]=analogRead(A0);
  }
  for(int i=0;i<COUNT;i++){
    Serial.println(Data[i]);
  }
  delay(3000);
}

const int COUNT=500;

int Data[COUNT];

void setup() {

Serial.begin(9600);

ADCSRA = ADCSRA & 0xf8; // 分周比を決めるビット(ADPS2:0)を000へ

//ADCSRA = ADCSRA | 0x04; // 分周比16(100)をセット

//ADCSRA = ADCSRA | 0x07; // 分周比128(111)をセット

ADCSRA = ADCSRA | 0x03; // 分周比__(xxx)をセット

}

void loop() {

for(int i=0;i<COUNT;i++){

Data[i]=analogRead(A0);

}

for(int i=0;i<COUNT;i++){

Serial.println(Data[i]);

}

delay(3000);

}

XLR8のデフォルトの分解能は10ビットなのでライブラリ”XLR8ADC”を使ってsetup関数の中で12ビットに変更します。
スケッチの”XLR8ADC.begin(); XLR8ADC.set12bitMode();”が12ビットに変更している箇所です。その他はUNOで使っているスケッチと同じです。

#include <XLR8ADC.h>
const int COUNT=500;
int Data[COUNT];

void setup() {
  Serial.begin(9600);
  // Enable the XLR8 ADC. Starts in 12 bit mode

  ADCSRA = ADCSRA & 0xf8;   // 分周比を決めるビット(ADPS2:0)を000へ
  //ADCSRA = ADCSRA | 0x04;   // 分周比16(100)をセット
  //ADCSRA = ADCSRA | 0x07;   // 分周比128(111)をセット

  ADCSRA = ADCSRA | 0x01;   // 分周比__(xxx)をセット

  XLR8ADC.begin();
  XLR8ADC.set12bitMode();
}

void loop() {
  for(int i=0;i<COUNT;i++){
    Data[i]=analogRead(A0);
  }
  for(int i=0;i<COUNT;i++){
    Serial.println(Data[i]);
  }
  delay(5000);
}

#include <XLR8ADC.h>

const int COUNT=500;

int Data[COUNT];

void setup() {

Serial.begin(9600);

// Enable the XLR8 ADC. Starts in 12 bit mode

ADCSRA = ADCSRA & 0xf8; // 分周比を決めるビット(ADPS2:0)を000へ

//ADCSRA = ADCSRA | 0x04; // 分周比16(100)をセット

//ADCSRA = ADCSRA | 0x07; // 分周比128(111)をセット

ADCSRA = ADCSRA | 0x01; // 分周比__(xxx)をセット

XLR8ADC.begin();

XLR8ADC.set12bitMode();

}

void loop() {

for(int i=0;i<COUNT;i++){

Data[i]=analogRead(A0);

}

for(int i=0;i<COUNT;i++){

Serial.println(Data[i]);

}

delay(5000);

}

<UNOの結果>
Arduino UNOのデフォルトの状態(分周比128)です。シリアルプロッタのグラフには500個のデータが表示されます。データを見ながら周波数を調整して、三角波の１周期が表示されるようにしたところ18Hzでした。18Hzの１周期は55.5msecで、500データありますから、0.111msec間隔となり、サンプリング周波数は9kHzぐらいと見ることができます。グラフのピーク値は1024*3/5=614あたりになっています。

16分周

8分周から波形にギザギザが少し見え始めました。

4分周で誤差が目立ち始めました。

2分周はどこが1周期なのか判別ができません。
ということで、Arduino UNOでサンプリングレートを上げたい場合は16分周ぐらいを使うようです。

128分周—–18Hz–サンプリング周波数9KHz
64———測定せず
32———測定せず
16———133Hz–サンプリング周波数66.5KHz
8———-236Hz–サンプリング周波数118KHz
4———-385Hz–サンプリング周波数192.5KHz
2———-１周期の判断できず

<XLR8の結果>
XLR8の128分周のグラフです。18Hz、0-3Vの信号ですが、12bitなので縦軸が4096*3/5=2458付近の大きな数値になっています。

128分周—–18Hz–サンプリング周波数9KHz
64———測定せず
32———測定せず
16———133Hz–サンプリング周波数66.5KHz
8———-236Hz–サンプリング周波数118KHz
4———-385Hz–サンプリング周波数192.5KHz
2———-400Hz–サンプリング周波数200KHz

FPGAで忠実にATmega328互換の動作をしているようで、サンプリングレートは同じですが、ADCがMAX10のものなので分周比2でも綺麗な波形です。

4.8KHzの信号を入れてみても大丈夫そうです。

<XLR8倍速モードの結果>
さて、XLR8は倍速モードも持っていて、FPGAを書き換えることで16MHZのATmega328から32MHZのATmega328になります。

128分周—–36Hz–サンプリング周波数18KHz
64———測定せず
32———測定せず
16———267Hz–サンプリング周波数133.5KHz
8———-464Hz–サンプリング周波数232KHz
4———-500Hz–サンプリング周波数250KHz
2———-500Hz–サンプリング周波数250KHz

サンプリング周波数は250KHzで頭打ちですが、２分周のグラフも特に変ではありません。

<結果のまとめ>
XLR8のアナログ入力はUNOに比べて分解能が4倍で、サンプリング周波数が3.7倍程度高くまで使えそうなので一考に値すると思います。

LabVIEW Community Editionの情報 (191001)

メイカー、ホームユーザー向け(非商用)に無料で使える「LabVIEW Community Edition」が来年以降オープンされるとのことです。NI Communityでの質疑から分かった情報をまとめます。

– LabVIEW Community Editionは非商用での利用に限り無料で使うことができます。

– LabVIEW Community EditionはLabVIEW Professional Editionと同じ機能が透かしなしで使うことができます。
LabVIEWがインストールされていないPCでも実行できるプログラムを作成することができます。従来のHome/Student Editionはフロントパネルやダイアグラムにエディション名が書かれた透かしが入っていましたが、LabVIEW Community Editionには透かしが入りません。

– Raspberry Pi, BeagleBone, ArduinoでLabVIEWを使うことができるLINXの新しいバージョンが提供されます。

– ベータ版は2019年第4四半期、正式版は2020年5月にLabVIEW2020版でリリースされる予定です。
LINXの新しいバージョンのベータ版はすでに”NI Software Technology Preview”で公開されています。2019年第4四半期にリリースされるベータ版も”NI Software Technology Preview”で公開されます。

以上(191001)

「スクエア&キューブパズル」地産材でのクラフト製作の勧め

“Square & cube puzzle” :Recommendation for crafting with local wood

この記事はMFTokyo2019 (Maker Faire Tokyo 2019 )エントリーシートの予稿です。

「メイカーフェア東京2019での計測・解析ラボのプロフィール」
計測・解析ラボは盛岡市を拠点に個人事業としてLabVIEWプログラムの受託開発とArduinoなどオープンなハード・ソフトの講習会を行っています。プログラムとレーザーカッターを使った木工パズルデザインは3年目となり、地域の木材活用と木工クラフト作家との連携を模索しています。MFTokyo2019ではパズル作品を通しての交流やメイカー向けLabVIEW Home版の情報共有を行いたいと思います。

「メイカーフェア東京2019での出展内容の紹介」
正方形から立方体を組むことのできる幾何学的な美しさが特徴のスクエア&キューブパズルです。立方体4個を平面に並べた5種類のピースを16個使ったパズルで、PCやスマホで疲れた目を休めて指先を使い立体図形に親しむことができます。レーザーカッターで作ったプロトタイプを見本に岩手県内の木工家に地域材を使ったクラフト製作を勧めています。もちろん3Dプリンターでの自作も可能です。

デモビデオ@ユーチューブ

パズルのデモビデオ(ユーチューブ)

「スクエア&キューブパズル」の製作の経緯
MFTokyo 2018では「チーム岩手」から「プログラミングとレーザーカッターの素敵な関係」というタイトルで参加しました。２次元の図案をレーザーカッターで加工すると厚みのある2.5次元のパズルに成ります。次回のMFTokyo 2019ではレーザーカッターで3次元のパズルを作成したいと考えていました。

ちょうどFusion 360を使い始めて立体的な箱物を作る時のはめ込み設計に便利なことに気がついて操作の練習をしていたのですが、レーザーカッターで作ったピースを組み立てて立方体にするアイデアが浮かびました。正方形から切り出したピースが立方体になるパズルは幾何的に美しい構成です。単位立方体の個数で表現すれば次のような関係に成ります。

(正方形) (立方体) (単位立方体の個数)

8 x 8 = 4 x 4 x 4 =64

27 x 27 = 9 x 9 x 9 = 729

128 x 128= 16 x 16 x 16 = 4096

とりあえず、一番小さな組み合わせの4 x 4 x 4の立方体で考えることにしました。テトロミノというテトリスで有名な形状は正方形4個の組み合わせで5種類できるので、Fusion 360で4 x 4 x 4の立方体のモデルを作って4個づつ”結合”して全てテトロミノにすることができるか実験してみました。何度か試行錯誤して16個のテトロミノに分割できました。次はその16個で正方形ができるかやってみました。正方形になりました。

平面パズルに使っている厚さ6mmのアカマツ板では2.4cmのキューブにしかならないので、厚さ10mmを注文して作成しました。できることはできるのですが、まだ小さいので、レーザーカッターの限界と言われる15mmの板を注文して作りました。正方形のサイズは12cm、立方体のサイズは6cmになります。正方形の外枠を14cm角にすると持ち運びを考えるて手頃な大きさだろうと思います。外箱は一関の金森紙器さんに作ってもらいました。

テトロミノについて知っておこう
正方形を2個つなぎ合わせた形状をドミノと言います。3個つなぎ合わせたのがトロミノ、4個つなぎ合わせたのがテトロミノです。5個つなぎ合わせるとペントミノですが、このようにN個つなぎ合わせた図形の総称をポリオミノと言います。テトロミノは5種類ありますが、図形の特徴からI、L、T、Z、Oと呼ばれています。 I、L、Tは見た目通りですが、正方形をO、残った1個がZです。
スクエア&キューブパズルには4個の単位立方体で作ったテトロミノが16個入っています。5種類のテトロミノから16 個選ぶ方法はたくさんありますが、I型が3個、L型が4個、T型が2 個、Z型が3個、O型が4個の組み合わせにしました。

地産材でのクラフト製作の勧め
岩手に旅行で来た時に地元のクラフトマンが作ったスクエア&キューブパズルを目にする日が来るかもしれません。地域の木材を使ったパズルは、木の種類によって色や木目、重さや手触りが違うことを発見して自然を身近なものとして考え、豊かな森林への関心を高めることができるでしょう。西和賀町に広がる国内有数のブナ林、北限のツバキが咲く大船渡、遠野地方の広葉樹、久慈の南部赤松、平庭高原の白樺など広い岩手県には地域を代表する木が多くあります。天然木の手触りの良さとともに地域性が感じられるスクエア&キューブパズルが各地で購入できたら嬉しいことです。もちろん秋田の杉、青森のヒバなどを使って各地のクラフトマンが作ってくれたらもっと面白いだろうなと思います。

スクエア&キューブパズルに関連する特許
独自の発想で作ったパズルですが、幾何的にきれいなパズルなので誰か他の人も考えていただろうと思って、スクエア&キューブパズルに関係する特許について調べてみました。古い特許ですが1988年にアメリカで登録された「BLOCK PUZZLE」という特許がありました。テトロミノ16個を正方形に詰めるパズルとテトロキューブ(立方体を立体的に4個つないだ図形)16個を立方体に組むパズルを併記しています。正方形に詰めるパズルの厚みを変えるだけで立方体ができるのですが、念頭にはなかったようです。

サイコロでスクエア&キューブパズルを作ろう
とりあえずパズルが欲しい人はサイコロを瞬間接着剤で接着してパズルピースを作ってみましょう。16mmのサイコロが50個入りで1100円ぐらいで購入できます。

メイカーフェア東京2018

メイカーフェア東京2018に出展します。
木目のパズル「プログラミングとレーザーカッターの素敵な関係」というタイトルで「チームいわて」のメンバーとして参加します。
ブースはC/02-02です。
「LabVIEW Home版をメイカーな人たちに広めよう」ということで、FM音源モジュールを接続したArduinoに、LabVIEWから音源パラメータを送り込むプログラムをデモします。

タッチディスプレイでLabVIEWのUIをクリックすると、
パラメータがArduinoに送られて、音源モジュールのレジスタが書き換えられ、
音源もジュールからPCのマイク端子に音声信号が送られ、
LabVIEWで信号を受け取って波形とパワースペクトルが表示され、
PCのイヤホン端子から外部スピーカーで音が出る、
というものです。

今回のデモ用にたまたまHPの安いノートPCを購入したのですが、タッチディスプレイの便利さを実感しました。
ブースレイアウトを確認してパッキングしました。

(180803)

木目のパズル「プログラミングとレーザーカッターの素敵な関係」

Wood grain puzzle “Lovely relationship between Processing programming and laser cutter”

この記事はMFT2018 (Maker Faire Tokyo 2018 )エントリーシートの予稿です。

出展者のプロフィール
レーザー加工といえばイラストレーターで図案を描いてMDFを切り抜くものという風潮に抗って、Processingプログラムでパズル図案を描いて、地元の木材をレーザーカッターで切り抜いて木目のパズルを作っています。本業は計測の自動化と解析処理用のLabVIEWプログラムを提供していますので、LabVIEWホーム版を使ったArduino+FM音源YMF825の音作りツールも紹介する予定です。

出展内容
プログラムで描いた画像をレーザー光線がたどって切り取ってくれる。画面の中の五角形が厚みを持った五角形の木片になる。ジグソーパズルの無限の組み合わせをプログラムで描いて、地元の木材でパズルを作ってみよう。絵はいらない、木目が元の場所を教えてくれるから。久慈の南部アカマツ、遠野の桜、浄法寺のウルシ、木目は木の物語。プログラミングとレーザーカッターの素敵な関係を提案します。木目のパズルを見に来てね。

ジグソーパズルの製作
(1)Processingでピースの形状をランダムに作画します。
(2)レーザーカッターで切断します。
(3)パズルの完成です。
(4)画面の図形が厚みを持ってさわれる!!!

デモビデオ「プログラミングとレーザーカッターの素敵な関係」
https://www.facebook.com/koji.ohashi.12/videos/1696150310460940/

五角形のパズル

Processingでレーザー加工用図形をプログラムしよう！！

Processingって何？

たくさんあるプログラムツールの中でProcessingを選ぶ理由は四つ。
(1) オープンソースなので無料で使うことができる。
(2)図形をPDFで出力できる。
(3)グラフィック関連の関数が一通り揃っている。
(4)プログラムサンプルが豊富で使いやすい。

タングラムという古典的な図形パズルを例にプログラム方法を紹介します。

タングラムって何？

プログラムで作成するタングラムの図案。
図案が少し違いますが、タングラムを切断加工している様子です。
プログラムサンプル

processingでのプログラムはディスプレイに表示することを主眼に作られているのでピクセルが基本単位です。72ピクセルが1インチ(25.4mm)ですから、1mmは72/25.4ピクセルとなります。長さ100mmは100*72/25.4ピクセルとなります。つまり、欲しい寸法(mm)とプログラムで使うピクセルの変換は72/25.4を掛ければ良いだけです。

また、ドキュメントサイズを変数で指定することとPDFライブラリを使用していることに注目してください。

// タングラムの加工図の例です。
//*********************************************************
//PDFが作成されるので加工の際にイラストレーターで開きます。
//イラストレーターで図形全体を選択した状態で、
//「オブジェクトメニュー>クリッピングマスク>解除」します。
//アートボード外周に不要なパスができているので削除します。
//イラストレータ形式で保存します。
//これでレーザー加工ができます。
//*********************************************************
//テンプレートとして活用してください
//計測・解析ラボ


import processing.pdf.*;

float artboardHeight=80;  //実寸(mm)
float artboardWidth=130;  //実寸(mm)
float frameHeight=70;  //実寸(mm)
float frameWidth=120;  //実寸(mm)
float puzzleHeight=50;  //実寸(mm)
float puzzleWidth=100;  //実寸(mm)

float pixConv=72/25.4;
float artboardHeightP=artboardHeight*pixConv;
float artboardWidthP=artboardWidth*pixConv;
float frameHeightP=frameHeight*pixConv;
float frameWidthP=frameWidth*pixConv;
float puzzleHeightP=puzzleHeight*pixConv;
float puzzleWidthP=puzzleWidth*pixConv;

float frameTopP=(artboardHeightP-frameHeightP)/2;;
float frameLeftP=(artboardWidthP-frameWidthP)/2;
float puzzleTopP=frameTopP+(frameHeightP-puzzleHeightP)/2;
float puzzleLeftP=frameLeftP+(frameWidthP-puzzleWidthP)/2;

float strokeWidth = 0.001*pixConv;//( TrotecLaser LineWidth:0.001mm)


void settings(){
  int docWidth=(int)(artboardWidthP);
  int docHeight=(int)(artboardHeightP);
  size(docWidth,docHeight);
}

void setup() {
  noLoop();  // Run once and stop
  beginRecord(PDF, "Tangram.pdf");
  noFill();
  strokeWeight(strokeWidth);
  stroke(0,0,0); //black

}

void draw(){
  translate(puzzleLeftP,puzzleTopP);
  
  line(0,puzzleHeightP/2,puzzleWidthP/4,0);
  line(puzzleWidthP/4,0,puzzleWidthP/4,puzzleHeightP/2);
  line(puzzleWidthP/4,puzzleHeightP/2,puzzleWidthP/2,puzzleHeightP/2);
  
  line(0,puzzleHeightP,puzzleWidthP/4,puzzleHeightP/2);
  line(puzzleWidthP/4,puzzleHeightP/2,puzzleWidthP/2,puzzleHeightP);
  line(puzzleWidthP/2,puzzleHeightP,puzzleWidthP/2,0);
  line(puzzleWidthP/2,0,puzzleWidthP,puzzleHeightP);
  
  rect(0, 0, puzzleWidthP, puzzleHeightP);
  translate(-puzzleLeftP,-puzzleTopP);
  rect(frameLeftP, frameTopP, frameWidthP, frameHeightP);
  
  endRecord();
}

// タングラムの加工図の例です。

//*********************************************************

//PDFが作成されるので加工の際にイラストレーターで開きます。

//イラストレーターで図形全体を選択した状態で、

//「オブジェクトメニュー>クリッピングマスク>解除」します。

//アートボード外周に不要なパスができているので削除します。

//イラストレータ形式で保存します。

//これでレーザー加工ができます。

//*********************************************************

//テンプレートとして活用してください

//計測・解析ラボ

import processing.pdf.*;

float artboardHeight=80; //実寸(mm)

float artboardWidth=130; //実寸(mm)

float frameHeight=70; //実寸(mm)

float frameWidth=120; //実寸(mm)

float puzzleHeight=50; //実寸(mm)

float puzzleWidth=100; //実寸(mm)

float pixConv=72/25.4;

float artboardHeightP=artboardHeight*pixConv;

float artboardWidthP=artboardWidth*pixConv;

float frameHeightP=frameHeight*pixConv;

float frameWidthP=frameWidth*pixConv;

float puzzleHeightP=puzzleHeight*pixConv;

float puzzleWidthP=puzzleWidth*pixConv;

float frameTopP=(artboardHeightP-frameHeightP)/2;;

float frameLeftP=(artboardWidthP-frameWidthP)/2;

float puzzleTopP=frameTopP+(frameHeightP-puzzleHeightP)/2;

float puzzleLeftP=frameLeftP+(frameWidthP-puzzleWidthP)/2;

float strokeWidth = 0.001*pixConv;//( TrotecLaser LineWidth:0.001mm)

void settings(){

int docWidth=(int)(artboardWidthP);

int docHeight=(int)(artboardHeightP);

size(docWidth,docHeight);

}

void setup() {

noLoop(); // Run once and stop

beginRecord(PDF, "Tangram.pdf");

noFill();

strokeWeight(strokeWidth);

stroke(0,0,0); //black

}

void draw(){

translate(puzzleLeftP,puzzleTopP);

line(0,puzzleHeightP/2,puzzleWidthP/4,0);

line(puzzleWidthP/4,0,puzzleWidthP/4,puzzleHeightP/2);

line(puzzleWidthP/4,puzzleHeightP/2,puzzleWidthP/2,puzzleHeightP/2);

line(0,puzzleHeightP,puzzleWidthP/4,puzzleHeightP/2);

line(puzzleWidthP/4,puzzleHeightP/2,puzzleWidthP/2,puzzleHeightP);

line(puzzleWidthP/2,puzzleHeightP,puzzleWidthP/2,0);

line(puzzleWidthP/2,0,puzzleWidthP,puzzleHeightP);

rect(0, 0, puzzleWidthP, puzzleHeightP);

translate(-puzzleLeftP,-puzzleTopP);

rect(frameLeftP, frameTopP, frameWidthP, frameHeightP);

endRecord();

}

5V駆動のステッピングモーターを試す

28BYJ48 DC 5V 4-フェーズ5線式ステッピングモータとULN2003ドライバボード
これも去年の春に購入したまま眠っていたパーツです。中国直送で200円前後で購入できます。

便利なライブラリを見つけて動かしてみたという段階です。通常のstepperライブラリではなく、AccelStepperというライブラリを使うことでハーフステップで滑らかに駆動することができます。28BYJ-48モーターでの使い方は42Botsの記事に従いました。

回転数は1000pps(およそ15rpm)が上限で、それ以上はトルクが小さくなり指でシャフトをつまむと止まってしまいます。

//http://42bots.com/tutorials/28byj-48-stepper-motor-with-uln2003-driver-and-arduino-uno/

//Yes. It also helps that I've done the experiment with a fragment of mirror on the motor shaft,
// a laser pointer, and a bit of tape on the wall, confirming the 4096 and refuting the 4048.
#include <AccelStepper.h>
#define HALFSTEP 8

// Motor pin definitions
#define motorPin1  8     // IN1 on the ULN2003 driver 1
#define motorPin2  9     // IN2 on the ULN2003 driver 1
#define motorPin3  10     // IN3 on the ULN2003 driver 1
#define motorPin4  11     // IN4 on the ULN2003 driver 1

// Initialize with pin sequence IN1-IN3-IN2-IN4 for using the AccelStepper with 28BYJ-48
AccelStepper stepper1(HALFSTEP, motorPin1, motorPin3, motorPin2, motorPin4);

void setup() {
  stepper1.setMaxSpeed(1000.0);
  stepper1.setAcceleration(500.0);
  stepper1.setSpeed(2000);
  stepper1.moveTo(40760);

}//--(end setup )---

void loop() {

  //Change direction when the stepper reaches the target position
  if (stepper1.distanceToGo() == 0) {
    stepper1.moveTo(-stepper1.currentPosition());
  }
  stepper1.run();
}

//http://42bots.com/tutorials/28byj-48-stepper-motor-with-uln2003-driver-and-arduino-uno/

//Yes. It also helps that I've done the experiment with a fragment of mirror on the motor shaft,

// a laser pointer, and a bit of tape on the wall, confirming the 4096 and refuting the 4048.

#include <AccelStepper.h>

#define HALFSTEP 8

// Motor pin definitions

#define motorPin1 8 // IN1 on the ULN2003 driver 1

#define motorPin2 9 // IN2 on the ULN2003 driver 1

#define motorPin3 10 // IN3 on the ULN2003 driver 1

#define motorPin4 11 // IN4 on the ULN2003 driver 1

// Initialize with pin sequence IN1-IN3-IN2-IN4 for using the AccelStepper with 28BYJ-48

AccelStepper stepper1(HALFSTEP, motorPin1, motorPin3, motorPin2, motorPin4);

void setup() {

stepper1.setMaxSpeed(1000.0);

stepper1.setAcceleration(500.0);

stepper1.setSpeed(2000);

stepper1.moveTo(40760);

}//--(end setup )---

void loop() {

//Change direction when the stepper reaches the target position

if (stepper1.distanceToGo() == 0) {

stepper1.moveTo(-stepper1.currentPosition());

}

stepper1.run();

}