Myu: 2025

2025年7月2日水曜日

今年から市民農園再開

しばらくやめていた市民農園に再挑戦で、今年は小玉スイカの空中栽培も試してみた。植え付け時は、炭疽病を心配し、乗り越えたかと思えば、今度はウリハムシ、そして受粉がなかなかうまくいかないとか、いろいろありましたが、今のところ、それほど雨に悩まされることなく病気の症状もおさまっているようです。このまま梅雨明けしてくれるといいのですが。

　虫や病気の心配も少なくネットもいらず、片付けも楽、意外と手間がかからず、コンスタントに収穫できるのがつるなしインゲンです。

2025年6月26日木曜日

ロープの結び方

畑作業でもよく使うので、男結び等について　図にまとめてみました。（トラックロープの結び方は南京結びともいうようです。）　他にも、巻き結び、もやい結び、結束結び、自在結び、ひとえつぎ結び等、覚えておくと便利なようです。

GASではJSON.stringfyをつかったほうがいいようです

畑の収穫記録をGASで管理（収穫の個数を記録しその都度、合計の確認）しようかと、ChatGPTにコード頼んでみたけど、なかなか一筋縄ではいかない。AIにより、だいぶコーディングは楽になっているが、やはり、ある程度やりとりしながら、人間側が修正してあげないと動かなかった。

　最初、ChatGPTはサーバ（GAS）からのレスポンスデータを単純にオブジェクトで渡すコードを提示してくれたけれど、それではうまくいかず、やはり、JSON.stringfyを使わないとだめだった。

//***********************html.index**************************
<!DOCTYPE html>
<html>
  <head>
    <base target="_top">
    <style>
      body { font-family: sans-serif; margin: 20px; }
      label, select, input, button { margin: 5px; }
      div.result-entry { margin-top: 10px; }
    </style>
  </head>
  <body>
    <h2>畑の収穫記録</h2>
    <label>野菜：</label>
    <select id="vege"></select>

<label>メモ：</label>
    <ttextarea id="memo" name="memo" rows="5" cols="40"></ttextarea>　←ttextareaをtextareaに変更必要

<div id="result"></div>
    <h3 id="totalSum">合計金額: -- 円</h3>
<br /> <a href="https://docs.google.com/spreadsheets/d/シートID/edit?gid=0#gid=0">Sheet</a>
    <script>
      let vegeData = {};

function submitData() {
        const name = document.getElementById("vege").value;
        const count = parseInt(document.getElementById("count").value);
        const memo = document.getElementById("memo").value;

console.log("name:", name, "count:", count);

if (!name || !count || count < 1) {
          alert("野菜と個数を正しく入力してください。");
          return;
        }

google.script.run.withSuccessHandler(displayResult).addRecord({ name: name, count: count,memo: memo });
      }

function displayResult(res) {

let response = JSON.parse(res);         
 
        const r = response.record;
        
        const div = document.getElementById("result");

document.getElementById("totalSum").innerText = `合計金額: ${response.totalSum} 円`;
      }

loadVegetables();
    </script>

//*************コード.js****************
const SPREADSHEET_ID = 'シートID';
const SHEET_NAME = 'Records'; // シート名（1枚目なら "Sheet1" のままでもOK）

// 野菜マスタ（名前と単価）
const vegetablePrices = {
  "トマト": 100,
  "ナス": 80,
  "キュウリ": 50,
  "ピーマン": 60,
    "ゴーヤ": 60,
      "インゲン": 5,
        "スイカ": 300,
          "ミニトマト": 10,
            "サツマイモ": 130,
              "サトイモ": 60,
                "シソ": 5,
                  "エダマメ": 10,

};

function doGet() {
  return HtmlService.createHtmlOutputFromFile('index');
}

function getVegetables() {
  return vegetablePrices;
}

function addRecord(data) {
    const ss = SpreadsheetApp.openById(SPREADSHEET_ID);
 // const ss = SpreadsheetApp.getActiveSpreadsheet();
  const sheet = ss.getSheetByName(SHEET_NAME);

if (!sheet) {
    Logger.log('シート名「' + SHEET_NAME + '」は存在しません');
  } else {
    Logger.log('シート名「' + SHEET_NAME + '」が見つかりました');
  }

const name = data.name;
  const count = Number(data.count);
  const unitPrice = vegetablePrices[name] || 0;
  const total = unitPrice * count;
  const date = new Date();
  const memo = data.memo

sheet.appendRow([date, name, count, unitPrice, total,memo]);

const sum = getTotalSum();  // 合計金額取得

// return {
  //  record: { date: date, name: name, count: count, unitPrice: unitPrice, total: total },totalSum: sum
 // };

return JSON.stringify({record: { date: date, name: name, count: count, unitPrice: unitPrice, total: total ,memo : memo},totalSum: sum });
}

function getTotalSum() {

const ss = SpreadsheetApp.openById(SPREADSHEET_ID);
  const sheet = ss.getSheetByName(SHEET_NAME);
  const data = sheet.getRange(2, 5, sheet.getLastRow() - 1).getValues(); // 合計金額列
  return data.reduce((sum, row) => sum + (Number(row[0]) || 0), 0);
}

※　MEMO欄を追加しました。2025/6/27

2025年6月2日月曜日

ガロア理論入門（アルティン）　問１１－３（２）

P117　ｐ素数　Φp^n(x)=Φp(x^(p^(n-1)))が既約であることを示す問題の途中の変形がしばらくわからず悩む。

Φp(x^(p-1)+(pでわりきれる項）＋１）＝{x^(p^(n-1))＋（ｐでわりきれる項)}^(p-1)＋（ｐで割り切れる項）＋ｐとなる理由

∵前問で①Φp(x+1)=((x+1)^p-1)/((x+1)-1)=x^(p-1)+Σk=1~(p-2)　pCk x^(p-k-1)+p が言えていたので、これを使うとよいことにきづく。今の場合、X＝x^(p-1)+(pでわりきれる項）　としてこれを①のxをXで置き換えるとうまくいくようです。

　前問をつかっているということに、きづかなかったためしばらく悩みました。よく読むと問題にも（１）と同様にして解く　と書いてあるのを見逃していました。

2025年5月17日土曜日

アルティンのガロア理論入門　P106

途中の式変形で、すこし手こずる。x^(q^n)-α^(q^n)-(x-α)を(x-α）で割ってから、x=αを代入した式がq^n*α^(q^n-1)-1となることを示したい。ChatGPTにすぐ頼ってしまったが、いまいちすっきりしない。

　与式をf(x)とすると　　f(x)=(x-α)h(x)とみなし　h(x)=f(x)/(x-α)だから、h(x)=f'(x)とみなせるからという説明だった。たしかに、x→αでは、解析的には微分だが、、、、(自分には高度すぎる？）。

　と、少し間をおいてあらためて考えると単純に式の割り算をすると、

Σ i=0~q^(n-1) α^i*x^(q^n-i-1) 　-1 となることに気づく。これにx=αを入れると、結論が言えるようだ。自分としては、こちらのほうが納得できる。

2025年3月29日土曜日

モナドの活用

Haskell入門　を久しぶりに読み直してみた。P168にMaybeモナドについて少し復習してみた。

getItemWithMonad menu category name =do
subMenu <- lookup category menu
price <- lookup name subMenu
return (category,name,price)

途中、Nothingを記述しなくても自動的に、該当しなければ、最後はNothingになる。うまくできていると思う。
これを、手続き型言語で、実現しようとすると条件分岐に頼らざるを得ないけど、その必要もない。
このへんが、関数型言語のモナドの強みなのかもしれない。

追記＞試しに、ChatGPTで、ELMのMaybeを使うコードに変換してみた。ELMはモナドをサポートしていないため、少しコード量が増えてしまうようだ。

type alias Menu = Dict String (Dict String Float)

getItemWithMonad : Menu -> String -> String -> Result String (String, String, Float)

getItemWithMonad menu category name =

case Dict.get category menu of

Just subMenu ->

case Dict.get name subMenu of

Just price ->

Ok (category, name, price)

Nothing ->

Err "Item not found"

Nothing ->

Err "Category not found"

2025年3月19日水曜日

1960年代からの国債残高の増加率をグラフにしてみた。

　国債の増え方が話題になっているようなので、実際のところどうなのか、公式データを加工して処理してみました。この計算の仕方が、実際の増え方を表しているかどうか議論のあるところかもしれませんが、ある程度の参考にはなりそうです。

　縦軸は前年比の増加率（％）です。ちなみに、1873年　0.234億円　2024年 11053645億円ということのようなので　この期間の１年あたりの平均増加率を計算すると (11053645/0.234)^(1/(2024-1873))=1.124　　　ということで、ほぼ12%程度です。

　１９８５年あたりから１０％切るようになっていますね。それ以来ほとんど、(1998年以外は)平均に届いてないようです。

追加：ちなみに、近年のものですがアメリカも調べてみました。日本のように、５％切ることは少ないような感じです。

2025年3月14日金曜日

ルベーグ積分　理論と計算手法　P199

補題9.13(2)でまた、つまづいてしまった。仕方なく、だめもとでChatGPTに頼ってみたら、きちんと回答してくれた。AIの力も凄いものだと驚きました。ここ何か月間の間にパワーアップした？

　質問は　「ルベーグファトウの定理を使い、∫A max k≦n fk μ ≦　ν(A) から ∫ sup n∈N fn μ ≦ ν(A)をいえますか？」としました。回答に少し補足を入れてまとめてみた。

∵）　gm=max k≦n fk （これは増加列　①）として　定理を適用すると

∫A lim inf n→∞ gn μ ≦ lim inf n→∞　∫A gn μ

左辺は ∫A sup n∈N fn μ　　② ↓ここで①を使う

また　sup n∈N inf k≧n gn = lim n→∞ inf gn=lim n→∞ gn =sup n∈N fn これが②の理由

∫A gn μ　≦　ν(A) で n→∞　ならば　　lim inf n→∞　∫A gn μ≦　ν(A)

ゆえに　結論　　∫A sup n∈N fn　μ　≦　ν(A)

定理に必要な優関数がないので、質問したら、優関数はf1(x)という回答。専任の家庭教師がついているような感覚になる。

2025年2月20日木曜日

Σsinkx/k　の求め方で苦戦

ルベーグ積分　理論と計算手法　P183　例8.9(1)で、

∞
Σsinkx/k=(π-x)/2 (2π>x>0)になる理由がわからず、苦戦。

log(1-z)を微分し　-1/(1-z) = Σz^n　なので

積分して　　log(1-z)=-∫Σz^n dz =-Σ∫z^n dz = - Σ z^(n+1)/(n+z)

= - Σ z^n/n (Σが n=0からとなっていたのをn=1からとすることで変形）

これを使うと

Im(log(1-e^ix))=-ΣIm(e^ixk)/k = -Σ sin kx /k

これは arg(1-e^ix)=-(π-x)/2　(0<x<2π)　でもある

∵）log(re^ix)=logr+log(e^ix)=logr+ix

∵) 1-e^ix=(1-cosx)-isinx r=√((1-cosx)^2+sin^2x =2sin(x/2)

-sinx / 2sin(x/2)=sinΘ　なら　Θ=-(π-x)/2になる。

∵） 2sin(x/2)sin((π-x)/2) = sinx このへんは地道に三角関数の計算を確認すると

　　確かにそうなる。それにしても計算が大変だった。

フーリエ級数使うともう少し楽に計算する方法もありそうなので、そちらも試してみたい。

追記(2/22)：その後、ディリクレージョルダンの定理(P182)を使うと、計算できることが確認できた。

f(x)=-x/(2π)+1/2*sgn(x)　として bk=1/kとできるので、これを定理に適用するといいようだ。最初うまくいかなかったのは、ちょっとした計算ミスのためだった。だいぶ遠回りしてしまった。

ついでに、(2)も解いていみた。

x∈(-2π,2π)を　x'=x/2∈(-π,π)とすると　　2Σsin2kx'/2k=Σsin2kx'/k=-x'+π/2*sgnx'

（１）より 2(-x/2+π/2*sgnx)=2Σsinkx/k=2( Σsink2x/2k + Σsin(2k-1)x/(2k-1) )

=-x+π/2*sgnx + 2*(求めたい式）となるが、これを整理すると

　求めたい式＝ Σsin(2k-1)x/(2k-1)=π/4*sgnxとなる。

PC用CW発振器（ツインT型）作成

peanutでは、PCでCW交信ができるようになっているが、そのための発振器を作ってみた。最初、２SC1815と２SA1015を組み合わせた簡易的なものを試したが、周波数がうまく変えられなかったので、ツインT型を試したらうまくいった。入力はUSBへの変換アダプタを使うことにしている。

以前作成したPICによるキーヤーのコレクタ部分に上記回路のアースを接続することで、うまく動作させることができた。

追記：送信練習用に、以前作ったモールスデコーダにもつなげるようにしてみた。ただ、音量が不足していたので、こちらのサイトを参考させていただき、一石のアンプを追加した。マイクを通さず、デコーダにつなぐと、雑音の影響がなくなり、ほぼ１００％の認識率になるようです。送信練習にはとても役立ちます。

2025年1月28日火曜日

JMoocに少し挑戦してみた

　無料の大学講座JMoocに興味を持ち、少し勉強してみた。修了証もらうまで行く人が１割ということらしい（無料ということとも影響している？）。単元テストは計算力を必要とされる問題が多く時間がかかった。

　抵抗２つ、コンデンサ２つを並列や直列組み合わせた回路で、抵抗値が流れる電流に影響しなくなる条件を求めよという問題や、交流回路の複素計算など、ちょっとしたミスがあるとなかなか解けなかったりする。久しぶりに、繁分数の式を変数に置き換えたり、〇.〇×10^△の形にして、少しでも計算が楽にならないか考えたり、LTSPiceなどPCの便利な計算に慣れてしまった昨今、逆に脳にはいい刺激になるかもしれないと思った。後半の三相交流の話になってくると、電力なども含めて考えると複雑になってきた。

　電気回路に続いて、電子回路も受講、電気回路が物理に近い内容で計算中心だったが、電子回路のほうほ工学的内容だった。後半は集積回路に組み込むために容量の大きなコンデンサを使わない回路の工夫について説明があった。MosFETなどＨＦリニアアンプの製作にもいろいろと役立つ情報があった。とりあえず、２つの講座ともなんとか修了できた。

2025年1月26日日曜日

中足骨骨頭痛

　以前から、歩くと靴の中に小石があるような痛みがあり、でも確かめると何もないという症状ががあり気になっていたが、よく調べると中足骨骨頭痛という症状らしい。足のアーチが弱くなるのが原因らしい。歩き方も踵に重心が来るように歩くといいらしい（足を前に出そうとするのでなく、上にあげてそのまま真下に落とすという感覚）。この症状が出ているときは、ウオーキングもほどほどにしたほうがいいようだ。

2025年1月25日土曜日

USB電力不足

　以前から気になっていたけれど、デスクトップPCのバスパワーＵＳＢハブが電力不足のためか、マウス以外つなげない状態。試しに、ケーブルの途中に三又状態にして、余ったＵＳＢ充電アダプタとUSBプラグを使って５Ｖを追加できるようにしてみたらうまくいった（若干コネクタ部分の接触が気になる感じはあったが、ＰＣからの経路すべてに接点復活材をスプレーしたらだいぶ改善された）。

　本来であれば、セルフパワーのＵＳＢハブを購入すればいいのだろうけど、数千円はするので、余り物でなんとかできたので、少し節約はできたかもしれません。（というか、これ以上、５Ｖ電源アダプタは増やしたくなかったというのもありますが）

2025年1月24日金曜日

ChatGPTは計算が苦手？

ルベーグ積分理論と計算手法P159の計算をChatGPTに聞いてみたが、書籍の答えと一致しない。

∫(0,+∞) (y/(y^x+1) - y/(y+1)^2) λ(dy)　が確かlog2と答えた。

正しくは、たぶん

[　log(y^2+1) /2 - ( y/(y+1) - ∫1/(y+1)λ(dy) )　] (0,∞) 　かと思う。

↑部分積分

上記で計算すると確かに１になるので。

ChatGPTも、人間と同じようにまだ間違うこともあるようだ。

ChatGPTはLLMなので、どちらかというと帰納的推論？なのかもしれない。対して、wolframなどは、演繹的推論らしいけれど。。

それにしても、AIで、プログラミングもそうですが、いろんな独学もだいぶやりやすくなったと思う。

2025年1月19日日曜日

PIC12F1822(CPS)の実装

ブレッドボードでの実験はなんとかうまくいきそうなので、実装してみた。これからケースに収める予定。静電容量計測のケーブルもコネクタで取り外しできるようにしてみた。

　前回、arduinoで灯油残量センサを作ったときは、誘導雷で壊れたこともあるので、灯油タンクまでの配線が長くなるので、サージアブソーバのようなものもつけないとだめかもしれない。

追記：CPSを測定しているが、灯油タンクに、導線２本入れて試しているものの、なぜか、測定値が安定しない。原因はまだ不明。おそらく、タンクまでの配線の容量が大きく、また、温度等でRC発信周波数が安定しないためかと思われる。

2025年1月18日土曜日

PIC12F1822でCPSの実験（AQM1602Aを使用)

　静電容量の変化をCPSという機能で確認できるらしいので、実験してみた。以前、灯油タンクの残量検知を超音波センサとarduinoで試したがうまくいかなかったので、静電容量を使ってできないか、そのうち試してみたいという目的もあります。

・カウントのオーバーフローするとマイナス表示になるようなので、delay関数の値を調整する必要があった。

・RA3を入力で使えるようにしたかったので、#pragma config MCLRE = OFF #pragma config LVP = OFF とするとよいという情報がネット上にあった。（LVPの必要性はよくわからない？が）

・RA5のLED出力とRA3のスイッチ入力は、とくに機能は実装していないけれど、今後の何らかの機能が必要になったときのために、コードで反応を確認してみた。

・少しわかりにくかったのは、秋月電子のLCD（AQM1602A）のコントラストがコードで設定するようになっていたこと。以下の(1)は72の２の部分が4bitで　(2)は0101の4bitで5で　0100が4ですが、その下位2bitのみ関係しているので、00→01→10→11 で54~57の範囲になります。

　　LCD_cmd(0x72); //Contrast set 下位4bits...(1)
LCD_cmd(0x54); //Power/ICON/Contrast Control 下位2bits...(2)
　 (2)と(1)あわせて　2^6=64 stepの範囲から指定するという意味だった。（(2)のほうが上位） ※5Vだと上記ぐらいの値だけれど、3.6v(ニッケル水素充電池3本）だと、(1)0x73 (2)0x55ぐらいが必要だった。

// CONFIG1
#pragma config FOSC = INTOSC    // Oscillator Selection (INTOSC oscillator: I/O function on CLKIN pin)
#pragma config WDTE = ON        // Watchdog Timer Enable (WDT enabled)
#pragma config PWRTE = OFF      // Power-up Timer Enable (PWRT disabled)
#pragma config MCLRE = OFF      // MCLR Pin Function Select (MCLR/VPP pin function is digital input)
#pragma config CP = OFF         // Flash Program Memory Code Protection (Program memory code protection is disabled)
#pragma config CPD = OFF        // Data Memory Code Protection (Data memory code protection is disabled)
#pragma config BOREN = ON       // Brown-out Reset Enable (Brown-out Reset enabled)
#pragma config CLKOUTEN = OFF   // Clock Out Enable (CLKOUT function is disabled. I/O or oscillator function on the CLKOUT pin)
#pragma config IESO = ON        // Internal/External Switchover (Internal/External Switchover mode is enabled)
#pragma config FCMEN = ON       // Fail-Safe Clock Monitor Enable (Fail-Safe Clock Monitor is enabled)

// CONFIG2
#pragma config WRT = OFF        // Flash Memory Self-Write Protection (Write protection off)
#pragma config PLLEN = ON       // PLL Enable (4x PLL enabled)
#pragma config STVREN = ON      // Stack Overflow/Underflow Reset Enable (Stack Overflow or Underflow will cause a Reset)
#pragma config BORV = LO        // Brown-out Reset Voltage Selection (Brown-out Reset Voltage (Vbor), low trip point selected.)
#pragma config LVP = OFF         // Low-Voltage Programming Enable (Low-voltage programming enabled)  ****off for use RA3

#include <xc.h>
#include <stdio.h>
#include <stdlib.h>
#include <string.h>

// クロック周波数指定
//--------------------------------------------------------------
#define _XTAL_FREQ 32000000
#define LCD_addr 0x7C
// プロトタイプ宣言
//--------------------------------------------------------------
void I2C_init(void);
void I2C_start(void);
void I2C_stop(void);
void I2C_write(unsigned char dat);
void LCD_cmd(unsigned char cmd);
void LCD_char(unsigned char dat);
void LCD_init(void);
void LCD_clear(void);
void LCD_home(void);
void LCD_cursor(unsigned x,unsigned y);

/*　　 　　　　　    　 12F1822
 *                   ??????        
 *             VDD ?| 1      8 |? VSS
 *             RA5 ?| 2      7 |? AN0/RA0
 *    CPS**AN3/RA4 ?| 3      6 |? AN1/RA1/SCL → LCD(SCL)
 *             RA3 ?| 4      5 |? AN2/RA2/SDA → LCD(SDA)
 *                   ??????      
 */

void main(void) {
     char s[8] ;
     unsigned cap ;

OSCCON  = 0b01110010 ; // 内部クロックは8ＭＨｚとする 
   //  OSCCON  = 0b01110000;    // 内部クロック周波数：8MHz
     ANSELA  = 0b00010000 ; // アナログはAN3(RA4)を使用し、残りをすべてデジタルI/Oに割当  >>>Capacity
     TRISA   = 0b00011110 ; // AN3(RA4/CPS3)   入力その他のピンは出力に割当てる(RA3は入力)  RA5>>LED SCL[6Pin]、SDA[5Pin]を入力ピンに設定
     
     PORTA   = 0b00000000 ; // 出力ピンの初期化(全てLOWにする)
     // 容量検知モジュール(ＣＰＳＭ)の設定
     CPSCON0 = 0b00001000 ; // オシレータは中範囲(中速の発信周波数)で利用する
     CPSCON1 = 0b00000011 ; // CPS3から入力する
     // タイマー１の設定
     T1CON   = 0b11000001 ; // 容量検知オシレータでTIMER1をｶｳﾝﾄする、ﾌﾟﾘｽｹｰﾗｶｳﾝﾄ値 1:1
     TMR1H   = 0 ;          // タイマー1の初期化
     TMR1L   = 0 ;
     PEIE    = 1 ;          // 周辺装置割り込みを許可する
     GIE     = 1 ;          // 全割り込み処理を許可する 
     CPSON   = 1 ;          // 容量検知モジュール開始

//----------------------------------------------------------
    //            シリアル通信（I2C）初期設定
    //----------------------------------------------------------
    SSP1STAT = 0b10000000;  // 標準モード(100kHz)
    SSP1CON1 = 0b00101000;  // マスター（制御する側）に設定
    SSP1ADD = 0b01001111;   // 通信速度 = Fosc/(4*Clock)-1 = 32MHz/(4*100kHz)-1 = 79（10進数）     
     
 // ディスプレイの初期化
     LCD_init();       
     RA5=1; //LED On
     LCD_cmd(0x80);    // １列目に表示         
     printf("cap count:");
      // ２列目に表示
     LCD_cmd(0xC0);
     printf("Initializing..");
     
    while (1) {     
         LCD_cmd(0xC0);      
         printf("%10d",cap);                
     
         if(RA3 == 0){
             RA5=1;
         }else{
             RA5=0;
         }
         CPSON = 0 ;                 // 容量検知モジュール停止
          cap =(unsigned) (TMR1H*256) + TMR1L ; // カウント値を読み込む
         TMR1H = 0 ;                 // タイマー1の初期化
         TMR1L = 0 ;
          CPSON = 1 ; 
          // 容量検知モジュール開始
        //  RA5=1;
          __delay_ms(10) ;
          
       
      }
   return;    
}

void I2C_start() {
    SEN = 1;
    while(SEN);
}

void I2C_stop() {
    SSP1IF = 0;
    PEN = 1;
    while(PEN);
    SSP1IF = 0;
}

void I2C_write(unsigned char dat) {
    SSP1IF = 0;
    SSP1BUF = dat;
    while(!SSP1IF);
}

void LCD_init() {
    __delay_ms(40); //40ms wait
    LCD_cmd(0x38);  //8bit,2line
    LCD_cmd(0x39);  //IS=1 : extention mode set
    LCD_cmd(0x14);  //Internal OSC Frequency
    LCD_cmd(0x72);  //Contrast set                 Lower4bits...(1)
    LCD_cmd(0x54);  //Power/ICON/Contrast Control  Lower2bits...(2)    (2)(1)=64bit step
    LCD_cmd(0x6C);  //Follower control
    __delay_ms(200);//200ms wait
    LCD_cmd(0x38);  //IS=0 : extention mode cancel
    LCD_cmd(0x0C);  //Display ON
    LCD_cmd(0x01);  //Clear Display
    __delay_ms(2);  //wait more than 1.08ms
}

void LCD_cmd(unsigned char cmd) {
    I2C_start();
    I2C_write(LCD_addr);
    I2C_write(0x00);    // 0x00（コマンドの送信を示すbit）
    I2C_write(cmd);
    I2C_stop();
}

void putch(unsigned char text) {
    I2C_start();
    I2C_write(LCD_addr);
    I2C_write(0x40);    // 0x40（データの送信を示すbit）
    I2C_write(text);
    I2C_stop(); 
}

2025年1月17日金曜日

デスクトップPCをraspi経由で電源On

PCの電源Onは、USB、内部タイマー、キーボード等々いろいろあるけれど、スリープ状態にするため、夜間勝手にアップデートなどでOnなるのも気になるので、完全にraspiで操作できるようにしてみた。

　PCの電源は、導通でOnなので、例によってトランジスタ2sc1815を使った。raspiのGPIOは出力が3.3Vで、前回ブルーレイな場合と同じ回路（raspiのGPIOからトランジスタ経由でPCスイッチへ）にしてみたらうまくいった。最初、マザボピンの分岐コネクタに手持ちのピンでつないだが、接触不良起こしやすいようだったので、はんだ付けにする。PCの余分な穴にDC電源用のソケットをグルーガンで固定して、raspiとの接続ケーブルを取り外しできるようにした。（raspi側にも同様にソケットをつける。）

2025年1月14日火曜日

極限の計算のコツ

ルベーグ積分理論と計算手法P152　あたりで、ちょっとした極限の計算が分からない。最近は、どうにも解決できそうにないときChatGPTに聞くことにしている。

　1<α<3 lim x→0 (sinyx-yx/(1+x^2))/αｘ^α　が、さらっと０になることが書いてあるが、理由がわからなかった。

　こういうときは、いろいろ展開を使うといいらしい。

　　　sin yx = sin0 +yx*cos0 -(yx)^2*sin0/2! -(yx)^3*cos0/3!.....

=yx-(yx)^3/6 ①

1/(1+x^2)=1+(-x^2)+(-x^2)^2+(-x^2)^3+.....

なので yx/(1+x^2)=yx(1-x^2+x^4...)=yx-yx^3+yx^5..... ②

　　　与式=(①-②)/αx^α=(y-y^3/6)x^3... /αx^α だが

　　　　　α<3より 3-α>0　で x→0なら　　確かに　x^(3-α)　→　0がいえるようだ。（最初、ChatGPTも途中の計算間違っていたので、再度、確認したら正しい説明をしてくれた。やはり、ChatGPTにとっても少し難しくなるとミスもしやすくなるんだろうか？）

　学生時代にこういった勉強も多少はやったはずだけど、ほとんど記憶に残っていない、やはり使わないと忘れてしまうようです。

2025年1月13日月曜日

raspiで秋月のSSR操作

raspiにSSRを組み合わせてタイムスイッチつくってみた。

秋月電子のSSRキット(ゼロクロス、25A)を使った。100Vなので絶縁には細心の注意が必要。たまたま、先日、分解廃棄したDVDレコーダの中からヒートシンクを取り出していたので、早速活用。

フォトカプラを使っているようで、入力で15mmA消費とあったので、raspiのGPIOでも許容範囲であり、そのままつなげている。

ChatGPTに助けてもらい、以下のようなスクリプトを使ってうまく動作させることができた。(rapiでは、OSの中にすでに以下のようなスクリプトが使える環境がある）

*****ssr_control.sh**********

#!/bin/bash
GPIO_PIN=3 # GPIO3
if [ ! -d /sys/class/gpio/gpio$GPIO_PIN ]; then
echo "$GPIO_PIN" > /sys/class/gpio/export
echo "out" > /sys/class/gpio/gpio$GPIO_PIN/direction

fi
if [ "$1" == "on" ]; then
echo "1" > /sys/class/gpio/gpio$GPIO_PIN/value
elif [ "$1" == "off" ]; then
echo "0" > /sys/class/gpio/gpio$GPIO_PIN/value
else
echo "Usage: $0 {on|off}"
exit 1
fi

cronを設定すれば、タイムスイッチとして使える

30 5 * * * cd /home/pi && ./ssr_control.sh on

30 6 * * * cd /home/pi && ./ssr_control.sh off

2025年1月11日土曜日

PIC16F628AでUART for Dfplayer（備忘録）

古い型のPICはDfplayerとのUART通信の資料が少なく苦労した。
arduino系と違ってPICはよりハードに近いため、通信のタイミングなどより細かい調整が必要だとわかった。外部発振子が必要だという情報もあったが、内部発振でも可能なような気もしたので、少し頑張ってみた。

確認したこと：各ポートの入力、出力の設定（TRISA,TRISB)、内蔵発振４MHｚが使えること（#define _XTAL_FREQ 4000000と PCON = 0b00001000）、UART関連設定（TXSTA,RCSTA）　Config系（FOSC = INTOSCCLKで内蔵発振を指定）　SPBRG=25（４MHｚと9600Bps、非同期、ハイスピードなどの条件から計算できることが、データシートに書かれている）
ここまでは、自力でなんとかできたが、問題は、UARTの送信関数

最初は
void dfply(char cmd, char parameter) {
while (TXIF == 0);
TXREG = 0x7E;
　　...........
としたが、これだと、データが完全に送信しないうちに次を送ろうとしてだめらしい。ChatGPTから、以下の方法を教えてもらい、試したらうまくいった。
void dfply(char cmd, char parameter) {
while (!TXIF);
TXREG = 0x7E;
while (!TRMT);
　　....

　以下のようなコードでフォルダ、ファイル番号をスイッチ長押しでカウントして取得し、Mp3を指定して再生することはうまくできた。確認のため、フォルダやファイル番号をモールス信号をLEDで表示するようにした。

#define _XTAL_FREQ 4000000
#pragma config FOSC = INTOSCCLK // Oscillator Selection bits (INTOSC oscillator: CLKOUT function on RA6/OSC2/CLKOUT pin, I/O function on RA7/OSC1/CLKIN)
#pragma config WDTE = OFF
#pragma config PWRTE = OFF
#pragma config MCLRE = OFF
#pragma config BOREN = ON
#pragma config LVP = ON
#pragma config CPD = OFF
#pragma config CP = OFF

#include <xc.h>
#include <stdio.h>
#include <stdlib.h>
#include <string.h>

void dfply(char, int);
void dfply_(char, int,int) ;

void convertToMorse(int,char * ) ;
void playMorse(char* ) ;
void dashu();
void dot();
void play(int ,int);
int dfply_qry(char, int, int );

int flg=0;
int flg2=0;

static int cnum=0;  //current file number
static int cfnum=0; //current folder number 0:mp3

int main(int argc, char** argv) {
    
    TRISA = 0b00011000;  //  RA3,4 input Ok
    TRISB = 0b00000010 ; // RB1(RX)はINPUT　RB2(TX)とその他は全てOUTPUT
    PORTA = 0b00000000;
    PCON   = 0b00001000 ;//4MHz 内部クロックは４ＭＨｚとする
    CMCON  = 0b00000111 ;     // コンパレータは使用しない(RA0-RA4はデジタルピンで使用)
    VRCON  = 0b00000000 ;     // RA2はデジタルピンで使用する            
    TXSTA = 0b00100100; // 高速ボーレート、送信有効  
    RCSTA = 0b10010000; // 受信有効、シリアルポート有効
    SPBRG=25; //クロック数4Mhzなら：25、8Mhz：51、16:Mhz：103、32Mhz:207
    __delay_ms(2000);
    dfply(0x06, 0x14);  // vol 20
    __delay_ms(1000);
   while (1) {
   //RA2 LED          
        if(cnum==0 && cfnum==0 ){  //最初のみ点滅
              RA2=1;
              __delay_ms(2);
              RA2=0;
             __delay_ms(1000);
          }       
          
    while (RA3==1 && RA4==0 ){  //長押し　フォルダカウント
        cfnum++;
         dot();
         flg=1;
          __delay_ms(200);
     }

while (RA3 == 0 && RA4==1){  //長押し　ファイルカウント
          cnum++;
         dot();
         flg2=1;      
         __delay_ms(200);
      }   
 
     if(flg){  
         cnum=0;
         flg=0; 
         __delay_ms(500);
                    char morseCode[15];
            convertToMorse(cfnum,morseCode);
         playMorse(morseCode);         
      }

if(flg2)  {       
           __delay_ms(200);
           char morseCode[15];
            convertToMorse(cnum,morseCode);
           playMorse(morseCode);
              if(cfnum==0){
               dfply(0x03,cnum);
               dfply(0x0d,0);
           }else{
               play(cfnum,cnum);
           }           
           flg2=0;
        }       
      __delay_ms(300);
      cnum=0;
      //cfnum=0;
   }
   return (EXIT_SUCCESS);
}

void play(int cfnum_,int cnum_){
           dfply_(0x0f, cnum_, cfnum_); // 初期化再生
           dfply(0x0d, 0); // 再生コマンド
    }

void dfply(char cmd, int parameter){
    dfply_(cmd,parameter,0);    
}

void dfply_(char cmd, int parameter_dec,int folder_dec) {
    char parameter = (char)parameter_dec;          // 10進数を16進数に変換
    char folder = (char)folder_dec;
    
     while (!TXIF); // TXREGが空になるのを待つ
    TXREG = 0x7E;  // 送信開始デリミタ
    while (!TRMT); // データが完全に送信されるのを待つ

while (!TXIF);
    TXREG = 0xFF;
    while (!TRMT);

while (!TXIF);
    TXREG = 0x06;
    while (!TRMT);

while (!TXIF);
    TXREG = cmd;
    while (!TRMT);

while (!TXIF);
    TXREG = 0x00;
    while (!TRMT);

while (!TXIF);
    TXREG = folder;
    while (!TRMT);

while (!TXIF);
    TXREG = parameter;
    while (!TRMT);

while (!TXIF);
    TXREG = 0xEF; // 送信終了デリミタ
    while (!TRMT);
    
}

void dot() {
    RA2=1;
    __delay_ms(100);
    RA2=0;
    __delay_ms(100);
}

void dash() {
    RA2=1;
    __delay_ms(300);
    RA2=0;
    __delay_ms(100);
}

void playMorse(char* morseCode) {
int i = 0; // ループカウンタを初期化
while (morseCode[i] != '\0') { // 終端文字である '\0' をチェック
    if (morseCode[i] == '.') {
        dot(); // 点を表す関数を呼び出し
    } else if (morseCode[i] == '-') {
        dash(); // 線を表す関数を呼び出し
    } else if (morseCode[i] == ' ') {
        __delay_ms(300); // スペースの場合は遅延
    }
    i++; // 次の文字へ進む
}
}

void convertToMorse(int num, char* result) {
    // Morse code for numbers 0-9
    char* morseNumbers[] = {
        "-----",  // 0
        ".----",  // 1
        "..---",  // 2
        "...--",  // 3
        "....-",  // 4
        ".....",  // 5
        "-....",  // 6
        "--...",  // 7
        "---..",  // 8
        "----."   // 9
    };

char temp[20] = "";  // 一時的な結果を保存する配列
    result[0] = '\0';    // 結果を初期化

// 処理対象の数字が 0 の場合は特別扱い
    if (num == 0) {
        strcpy(result, morseNumbers[0]);
        return;
    }

while (num > 0) {
        int digit = num % 10;           // 現在の桁を取得
        strcpy(temp, morseNumbers[digit]); // 桁のモールスコードを取得
        strcat(temp, " ");              // 桁ごとにスペースを追加
        strcat(temp, result);           // 既存の結果に連結
        strcpy(result, temp);           // 結果を更新
        num /= 10;                      // 次の桁へ
    }
}

2025年1月8日水曜日

ルベーグ積分　理論と計算手法　P150

P150で「 t　∈　R|Z　で、Γ(t)Γ(1-t)sinπtが周期1である」の理由が省略されていたので、気になり、考えてみた。
Γ(n+t)=(n-1+t)Γ(n-1+t)=.....=(n-1+t).....tΓ(t) ①
{1-(n+t)}Γ(1-(n+t)) = Γ(1- {(n-1)+t})
[1-{(n-1)+t}]Γ(1-{(n-1)+t}) = Γ(1- {(n-2)+t})
......
{1-(1+t)}Γ(1-(1+t)) = Γ(1-t)
これより、 Γ(1-t)={1-(n+t)}.....{1-(1+t)}Γ(1-(n+t)) ②
①、②より
Γ(n+t)Γ(1-(n+t))=(n-1+t).....tΓ(t)×Γ(1-t)/{(1-(n+t)).....(1-(1+t))}
=(-1)^n Γ(t)Γ(1-t)
ここで sinπ(n+t)=(-1)^n sin πtなので
Γ(n+t)Γ(1-(n+t))sinπ(n+t)=(-1)^n Γ(t)Γ(1-t) (-1)^n sinπt
=(-1)^2n Γ(t)Γ(1-t) sin πt = Γ(t)Γ(1-t) sin πt 　n∈Zなら　これより、周期1を示しているだろうと、考えた。

2025年1月7日火曜日

室内室外温度計修理

　長く使っていた温度計が突然壊れた。最初表示されないので、電池接触かと思い、確認したら液漏れしかけた電池があったので、交換。それでもだめなので、接点復活材を使う。今度は表示されたが、操作ボタンが反応しない。操作ボタンの接触不良を疑い、裏ブタを外して、接点復活材を使う。これで、ようやく使えるようになった。電池で使うようなものは、ほとんど酸化被膜が原因の接触不良が多いので、部品交換せずに修理できるようだ。

　便利なものがあれば、それなりに故障もあるので、メンテナンスで手間はかかる。買い替えよりは、修理のほうが、環境にもやさしいかもしれない。