SimulinkでCO2センサのSCD30を動かす

　Simulink Support Package for Raspberry Piをmatlabにインストールをして、ラズパイ4のI2CバスにAdafruitのSCD30センサ・ボードを接続しました。
<img width="400" alt="ファイル名" src="https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/407635/a55208f4-8974-b968-e19e-9aa23234bc10.png">


### 環境

- matlab home License R2022b update2
- Simulink Support Package for Raspberry Pi Hardware バージョン22.2.0
- Windows 10 Pro 22H2　メモリ40Gバイト
- Raspberry Pi 4 Model B 2MB

　インストールは、次の記事を参照してください。
　　[Simulink Support Package for Raspberry Pi ① インストール](https://www.denshi.club/matlab/raspberrypi/simulink-support-package-for-raspberry-pi.html)

### 動作確認
　ライブラリを使わず（存在しない？）、レジスタをたたいてセンサのデータを得ます。下記のマニュアルを参照しました。
　　[Sensirion_CO2_Sensors_SCD30_Interface_Description.pdf](https://files.seeedstudio.com/wiki/Grove-CO2-Temperature-Humidity-Sensor-SCD30/res/Sensirion_CO2_Sensors_SCD30_Interface_Description.pdf)
　動作を確認するため、CircuitPythonでアクセスしました。
　　[CircuitPython 10行プログラミング Step9 (7) CO2センサSCD30](https://www.denshi.club/pc/python/circuitpython/circuitpython-10-step9-7-co2scd30.html)

### プログラム
##### IEEE754 単精度浮動小数点数の形式：binary32 の実数への変換
　センサから読み取ったCO2、温度、湿度の値は、4バイトで構成されるbinary32の形式です。Pythonはstructのunpackで簡単に変換できます。
　simulinkのfunctionでは、
　　[matlab SimulinkのBLE functionを作ってみた　その2　備忘録](https://qiita.com/jamjam/items/fb78c0d916c3b4e121c7)
でプログラムを実装しましたが、今現在、 extrinsic declarationsがエラーになってしまいます。そこで、ブロックを組み合わせて、サブシステムを作ります。
　この実数表現は、メモリの中では32ビットの領域を使います。

- 符号ビット: 1ビット
- 指数部の幅: 8ビット
- 仮数部の幅: 23ビット

![IEEE.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/407635/67e7b0a7-2e33-e39f-6ac0-85ee9f465292.png)
　　[IEEE754 単精度浮動小数点数の形式：binary32 の変換](https://www.denshi.club/cookbook/arduino/esp32/ieee-754-binary32.html)

　この記事にあるように、共用体unionを使う方法、ポインタを使う方法がありますが、どうも、matlabでは実現できないようです。
　[メーカのプログラム](https://github.com/Sensirion/arduino-core/blob/a3ab4f0112e44072af90f891662ccd11803cdc25/src/SensirionRxFrame.cpp)ではunionを使っていました。
　そこで、オーソドックスに、符号、指数、仮数をそれぞれ計算して、最後に掛け合わせて実数に変換する方法をとりました。
　Arduinoだと、次のコードになります。iiがbinary32のデータです。

```cpp
int8_t f = pow(-1, int(bitRead(ii,31)));
double  k = 1 + (double)((((ii<<1)<<8)>>9))/pow(2,23); 
float s = (float)pow(2,(int8_t)((uint8_t)((ii<<1)>>24)-127));
Serial.print("内部形式を実数に戻した　");Serial.println(f*k*s,6);
```

　ほぼ、上記のスケッチを実装しました。
　　[ウキ「単精度浮動小数点数」](https://ja.wikipedia.org/wiki/%E5%8D%98%E7%B2%BE%E5%BA%A6%E6%B5%AE%E5%8B%95%E5%B0%8F%E6%95%B0%E7%82%B9%E6%95%B0#:~:text=%E7%8F%BE%E5%9C%A8%E3%81%AE%E6%B5%AE%E5%8B%95%E5%B0%8F%E6%95%B0%E7%82%B9%E5%BD%A2%E5%BC%8F,%EF%BC%888%E3%82%AA%E3%82%AF%E3%83%86%E3%83%83%E3%83%88%EF%BC%89%E3%81%A7%E3%81%82%E3%82%8B%E3%80%82)
も参考になります。
　動作を確認して、サブシステムにしました。
![co104a.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/407635/767f0b8d-7db7-2ad5-5cfe-b02eaeabc999.png)

##### メインのプログラム



![co104.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/407635/134e423c-0f57-1383-cc5c-2e17fedfffc7.png)

　左上から、順にパラメタを示します。

![co201.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/407635/18e5f02c-fc24-12e5-d170-da4b48e959e0.png)

![co202.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/407635/ceaed300-0679-0277-954a-a913e30979ea.png)

![co203.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/407635/f9cfbd50-227e-1171-013d-a999dc15a7da.png)

![co204.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/407635/d0b2424f-cf3c-08f9-82d4-f5a517f47f03.png)

![co205.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/407635/a223b845-4736-7b8c-c078-ec0cb6610a5b.png)

![co206.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/407635/fe5d3953-ccb6-91ee-9280-edf01e8862ef.png)

![co207.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/407635/642d494a-ef54-1b6b-61c9-87d8d73fc35f.png)



　functionです。18バイトを読み取り、2バイトごとにCRCが入っているので、そこをスキップして、4バイトごとにデータをまとめています。

```cpp
function [c,t,h]= fcn(u)

persistent c1;
persistent c2;
persistent c3;
persistent c4;

persistent t1;
persistent t2;
persistent t3;
persistent t4;

persistent h1;
persistent h2;
persistent h3;
persistent h4;

c1=u(1)*256*256*256;
c2=u(2)*256*256;
c3=u(4)*256;
c4=u(5);
c = uint32(c1+c2+c3+c4);

t1=u(7)*256*256*256;
t2=u(8)*256*256;
t3=u(10)*256;
t4=u(11);
t = uint32(t1+t2+t3+t4);

h1=u(13)*256*256*256;
h2=u(14)*256*256;
h3=u(16)*256;
h4=u(17);
h = uint32(h1+h2+h3+h4);
```

### 配置と実行順
　試行錯誤していると、配置が乱雑になって、デバッグの情報のオーバーレイ、実行順序を実行すると、変な順番に並んでいることがあります。修正する方法がわからないので、あたらしい空のモデルを作り、順番にコピペしていきます。また、プロパティの優先順位も、順番に数値を振っておきます。
でも、readコマンドを送ったのに読み取りができないとかが起こることがあります。
　I2Cバスの転送速度は100kHzですが、50kHzに変更する方法を見つけていません。SCLクロック・ストレッチは実装されているようですが、時々変なデータ（0xff）を読むのは、クロックが早いのが理由かもしれません。

### I2Cバスでの確認
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/407635/f8f5b72a-6226-8e6d-6afa-72a6efb99039.png)
　正しいデータが読めていませんが、シーケンスは正しい順番で実行されています。
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/407635/bc9d59ef-5fd2-bac8-9600-045caf771992.png)
　このデバイス、ヒータがあるのが理由かもしれませんが、1か月ほど連続使用したら、読み出した値がおかしくなっていました。
　その後、電源を入れて数分は正常な値が返えることが多くあります。壊れかけのようです。