<?xml version="1.0" encoding="utf-8" ?>
<rss version="2.0" xmlns:atom="http://www.w3.org/2005/Atom">
<channel>
<title>arata-2929のブログ</title>
<link>https://ameblo.jp/arata-2929/</link>
<atom:link href="https://rssblog.ameba.jp/arata-2929/rss20.xml" rel="self" type="application/rss+xml" />
<atom:link rel="hub" href="http://pubsubhubbub.appspot.com" />
<description>ブログの説明を入力します。</description>
<language>ja</language>
<item>
<title>Matlabを使ったWavelet変換 サンプルコード付き</title>
<description>
<![CDATA[ <p>時間周波数応答を得ることができるWavelet変換。</p><p>&nbsp;</p><p>本ブログではMatlabを使ってWavelet変換を行うサンプルプログラムをご紹介します。</p><p>&nbsp;</p><p>ということで、実施することは</p><ul><li>1つの音声ファイルに対してWavelet変換</li><li>時間周波数応答のわかりやすい図(スカログラム)を作成</li></ul><p>これを実現するために</p><ol><li>音声ファイルの読み込み</li><li>Wavelet変換</li><li>図作成</li></ol><p>という3手順で説明していきます。</p><p>&nbsp;</p><p>1. 音声ファイルの読み込み</p><p>本プログラミングで用いる音声は以下のファイルです。</p><p><a href="https://stat.ameba.jp/user_images/20191112/15/arata-2929/02/3a/p/o0773044814638061843.png"><img alt="" height="243" src="https://stat.ameba.jp/user_images/20191112/15/arata-2929/02/3a/p/o0773044814638061843.png" width="420"></a></p><p>&nbsp;</p><p>Matlabで音声ファイルを読み込むときはaudioread()という関数を用います。</p><p>&nbsp;</p><p><a href="https://stat.ameba.jp/user_images/20191112/15/arata-2929/29/fa/p/o0895007814638067396.png"><img alt="" height="37" src="https://stat.ameba.jp/user_images/20191112/15/arata-2929/29/fa/p/o0895007814638067396.png" width="420"></a></p><p>&nbsp;</p><p>これで音声ファイルを読み込むことができます。</p><p>Fsにはサンプリング周波数</p><p>yには時間に対するデータ</p><p>が入っています。</p><p>&nbsp;</p><p>ワークスペースを見てみると、</p><p><a href="https://stat.ameba.jp/user_images/20191112/15/arata-2929/47/d5/p/o0496019914638070265.png"><img alt="" contenteditable="inherit" height="80" src="https://stat.ameba.jp/user_images/20191112/15/arata-2929/47/d5/p/o0496019914638070265.png" width="200"></a></p><p>というデータが入っていますね。</p><p>つまり、この音声ファイルは</p><p>サンプリング周波数が48000で</p><p>要素数が50160</p><p>であるデータと分かります。</p><p>&nbsp;</p><p>これらを使ってWavelet変換を実施します。</p><p>&nbsp;</p><p>2.<span style="display: inline !important; float: none; background-color: rgb(255, 255, 255); color: inherit; font-family: &quot;ヒラギノ角ゴ Pro W3&quot;,&quot;Hiragino Kaku Gothic Pro&quot;,&quot;ＭＳ Ｐゴシック&quot;,&quot;MS PGothic&quot;,sans-serif,&quot;メイリオ&quot;,Meiryo; font-size: 16px; font-style: normal; font-variant: normal; font-weight: 400; letter-spacing: normal; orphans: 2; text-align: left; text-decoration: none; text-indent: 0px; text-transform: none; -webkit-text-stroke-width: 0px; white-space: normal; word-spacing: 0px;">Wavelet変換</span></p><p>さて、本題のWavelet変換です。</p><p>Wavelet変換実施にあたり、主に使用する関数はwt()という関数です。</p><p>&nbsp;</p><p><a href="https://stat.ameba.jp/user_images/20191112/16/arata-2929/23/21/p/o1108033714638081088.png"><img alt="" height="128" src="https://stat.ameba.jp/user_images/20191112/16/arata-2929/23/21/p/o1108033714638081088.png" width="420"></a></p><p>&nbsp;</p><p>これでWavelet変換が実施されました。</p><p>&nbsp;</p><p>3.<span style="display: inline !important; float: none; background-color: rgb(255, 255, 255); color: inherit; font-family: &quot;ヒラギノ角ゴ Pro W3&quot;,&quot;Hiragino Kaku Gothic Pro&quot;,&quot;ＭＳ Ｐゴシック&quot;,&quot;MS PGothic&quot;,sans-serif,&quot;メイリオ&quot;,Meiryo; font-size: 16px; font-style: normal; font-variant: normal; font-weight: 400; letter-spacing: normal; orphans: 2; text-align: left; text-decoration: none; text-indent: 0px; text-transform: none; -webkit-text-stroke-width: 0px; white-space: normal; word-spacing: 0px;">図作成</span></p><p>時間周波数応答が分かりやすいように図を作成します。</p><p>Wavelet図の作成するにあたり、描き方等を一部指定しています。</p><p>&nbsp;</p><p><a href="https://stat.ameba.jp/user_images/20191112/16/arata-2929/58/9c/p/o1029021714638082906.png"><img alt="" height="89" src="https://stat.ameba.jp/user_images/20191112/16/arata-2929/58/9c/p/o1029021714638082906.png" width="420"></a></p><p>&nbsp;</p><p>そして、</p><p><a href="https://stat.ameba.jp/user_images/20191112/16/arata-2929/8b/60/j/o0875065614638084997.jpg"><img alt="" height="315" src="https://stat.ameba.jp/user_images/20191112/16/arata-2929/8b/60/j/o0875065614638084997.jpg" width="420"></a></p><p>このように、Wavelet変換をわかりやすく図に示したもの(スカログラム)を作成することができました。</p><p>&nbsp;</p><p>&nbsp;</p><p>Matlabは有料ですが、信号処理のための関数が多数用意されています。</p><p>その後のグラフ化も非常に簡潔に記述することができます。</p><p>&nbsp;</p><div contenteditable="false" style="border:1px dotted;padding:15px;border-radius:4px;"><table border="0" cellpadding="0" cellspacing="0" style="margin:0;table-layout:fixed;" width="100%"><tbody width="100%"><tr><td aligin="center" style="vertical-align:middle;" width="95"><span style="display:block;text-align:center;"><a alt0="AmebaAffiliate" alt1="MATLAB and Simulink Student Suite R2019b" alt2="Amazon" alt3="https://images-fe.ssl-images-amazon.com/images/I/412Za25hx4L._SL160_.jpg" alt4="1" href="https://www.amazon.co.jp/MATLAB-Simulink-Student-Suite-R2019b/dp/0989614034?SubscriptionId=AKIAJLD6FH2TADXIQKDQ&amp;tag=amebablog-a2506203-22&amp;linkCode=xm2&amp;camp=2025&amp;creative=165953&amp;creativeASIN=0989614034" target="_blank"><img alt="MATLAB and Simulink Student Suite R2019b" border="0" data-img="affiliate" src="https://images-fe.ssl-images-amazon.com/images/I/412Za25hx4L._SL160_.jpg" style="max-width:95px;vertical-align:middle;margin:0;"></a></span></td><td style="line-height:1.5;padding-left:15px;vertical-align:middle;"><a alt0="AmebaAffiliate" alt1="MATLAB and Simulink Student Suite R2019b" alt2="Amazon" alt3="https://images-fe.ssl-images-amazon.com/images/I/412Za25hx4L._SL160_.jpg" alt4="1" href="https://www.amazon.co.jp/MATLAB-Simulink-Student-Suite-R2019b/dp/0989614034?SubscriptionId=AKIAJLD6FH2TADXIQKDQ&amp;tag=amebablog-a2506203-22&amp;linkCode=xm2&amp;camp=2025&amp;creative=165953&amp;creativeASIN=0989614034" target="_blank">MATLAB and Simulink Student Suite R2019b</a><div style="padding: 3px 0;">10,789円</div><div style="font-size:0.83em;">Amazon</div></td></tr></tbody></table></div><p>&nbsp;</p><p>しかし、サンプルコードが少なく(主にMatlabのホームページのみ)勉強もしにくいです。</p><p>&nbsp;</p><p>ですので、Matlabの基礎を学びたい方には、サンプルコードから勉強する前に</p><p>この本がおすすめです。この本を読んでから、</p><p>サンプルコードをいじっていくような勉強スタイルが</p><p>よろしいかと思います。</p><p>&nbsp;</p><div contenteditable="false" style="border:1px dotted;padding:15px;border-radius:4px;"><table border="0" cellpadding="0" cellspacing="0" style="margin:0;table-layout:fixed;" width="100%"><tbody width="100%"><tr><td aligin="center" style="vertical-align:middle;" width="95"><span style="display:block;text-align:center;"><a alt0="AmebaAffiliate" alt1="最新 使える!MATLAB 第2版 (KS理工学専門書)" alt2="Amazon" alt3="https://images-fe.ssl-images-amazon.com/images/I/51mk3EJnFTL._SL160_.jpg" alt4="1" href="https://www.amazon.co.jp/%E6%9C%80%E6%96%B0-%E4%BD%BF%E3%81%88%E3%82%8B-MATLAB-%E7%AC%AC2%E7%89%88-KS%E7%90%86%E5%B7%A5%E5%AD%A6%E5%B0%82%E9%96%80%E6%9B%B8/dp/4061565532?SubscriptionId=AKIAJLD6FH2TADXIQKDQ&amp;tag=amebablog-a2506203-22&amp;linkCode=xm2&amp;camp=2025&amp;creative=165953&amp;creativeASIN=4061565532" target="_blank"><img alt="最新 使える!MATLAB 第2版 (KS理工学専門書)" border="0" data-img="affiliate" src="https://images-fe.ssl-images-amazon.com/images/I/51mk3EJnFTL._SL160_.jpg" style="max-width:95px;vertical-align:middle;margin:0;"></a></span></td><td style="line-height:1.5;padding-left:15px;vertical-align:middle;"><a alt0="AmebaAffiliate" alt1="最新 使える!MATLAB 第2版 (KS理工学専門書)" alt2="Amazon" alt3="https://images-fe.ssl-images-amazon.com/images/I/51mk3EJnFTL._SL160_.jpg" alt4="1" href="https://www.amazon.co.jp/%E6%9C%80%E6%96%B0-%E4%BD%BF%E3%81%88%E3%82%8B-MATLAB-%E7%AC%AC2%E7%89%88-KS%E7%90%86%E5%B7%A5%E5%AD%A6%E5%B0%82%E9%96%80%E6%9B%B8/dp/4061565532?SubscriptionId=AKIAJLD6FH2TADXIQKDQ&amp;tag=amebablog-a2506203-22&amp;linkCode=xm2&amp;camp=2025&amp;creative=165953&amp;creativeASIN=4061565532" target="_blank">最新 使える!MATLAB 第2版 (KS理工学専門書)</a><div style="padding: 3px 0;">3,080円</div><div style="font-size:0.83em;">Amazon</div></td></tr></tbody></table></div><p>&nbsp;</p><p>少しでも皆様のためになれば幸いです。</p><p><a href="https://blog.ameba.jp/reader.do?bnm=arata-2929" title="フォローしてね"><img alt="フォローしてね" height="50" src="https://stat100.ameba.jp/blog/ucs/img/decoPeta/pc/decoPeta_22.gif" width="150"></a></p><p>&nbsp;</p><p>&nbsp;</p><p>&nbsp;</p><p>&nbsp;</p><p>&nbsp;</p>
]]>
</description>
<link>https://ameblo.jp/arata-2929/entry-12544647254.html</link>
<pubDate>Tue, 12 Nov 2019 16:23:30 +0900</pubDate>
</item>
</channel>
</rss>
