\n <\/td>\n | function<\/span> [result<\/span>] =<\/span> Analys_GaboriWT(W<\/span>, X<\/span>, a_N<\/span>, f_range<\/span>, N<\/span>, Fs<\/span>, t<\/span>, Vc<\/span>)<\/td>\n <\/tr>\n \n <\/td>\n | %<\/span> ---------------------------------------------------- %<\/span><\/td>\n <\/tr>\n \n <\/td>\n | %<\/span><\/span><\/td>\n <\/tr>\n \n <\/td>\n | %<\/span> Gablor inverse wavelet<\/span><\/td>\n <\/tr>\n \n <\/td>\n | %<\/span><\/span><\/td>\n <\/tr>\n \n <\/td>\n | %<\/span> @param<\/span><\/td>\n <\/tr>\n \n <\/td>\n | %<\/span> W: ウェーブレット解析後の結果<\/span><\/td>\n <\/tr>\n \n <\/td>\n | %<\/span> X: 解析信号<\/span><\/td>\n <\/tr>\n \n <\/td>\n | %<\/span> sigma: Gabor関数のパラメータ<\/span><\/td>\n <\/tr>\n \n <\/td>\n | %<\/span> a_N: 周波数分割数<\/span><\/td>\n <\/tr>\n \n <\/td>\n | %<\/span> f_min: 最低周波数<\/span><\/td>\n <\/tr>\n \n <\/td>\n | %<\/span> f_max: 最高周波数<\/span><\/td>\n <\/tr>\n \n <\/td>\n | %<\/span> Fs: サンプリングレート<\/span><\/td>\n <\/tr>\n \n <\/td>\n | %<\/span><\/span><\/td>\n <\/tr>\n \n <\/td>\n | %<\/span> ---------------------------------------------------- %<\/span><\/td>\n <\/tr>\n \n <\/td>\n | \n<\/td>\n <\/tr>\n | \n <\/td>\n | %<\/span> スケール配列<\/span><\/td>\n <\/tr>\n \n <\/td>\n | Sj =<\/span> sqrt(1<\/span> ./<\/span> f_range<\/span>); %<\/span>[1,500]<\/span><\/td>\n <\/tr>\n \n <\/td>\n | Sj =<\/span> 1<\/span> ./<\/span> Sj<\/span>;<\/td>\n <\/tr>\n \n <\/td>\n | \n<\/td>\n <\/tr>\n | \n <\/td>\n | %<\/span> 係数(良くわからないので適当)<\/span><\/td>\n <\/tr>\n \n <\/td>\n | %<\/span> 逆変換後のXnをスケーリングする係数<\/span><\/td>\n <\/tr>\n \n <\/td>\n | Sigma_j =<\/span> 0.001<\/span>;<\/td>\n <\/tr>\n \n <\/td>\n | Cs =<\/span> 1<\/span>;<\/td>\n <\/tr>\n \n <\/td>\n | W0 =<\/span> 1<\/span>;<\/td>\n <\/tr>\n \n <\/td>\n | Const =<\/span> Sigma_j<\/span> *<\/span> sqrt(1<\/span>/<\/span>Fs<\/span>) /<\/span> Cs<\/span> /<\/span> W0<\/span>;<\/td>\n <\/tr>\n \n <\/td>\n | \n<\/td>\n <\/tr>\n | \n <\/td>\n | %<\/span> 再構成<\/span><\/td>\n <\/tr>\n \n <\/td>\n | Xn =<\/span> Const<\/span> *<\/span> real(W<\/span>) *<\/span> Sj<\/span>'<\/span>;<\/td>\n <\/tr>\n \n <\/td>\n | \n<\/td>\n <\/tr>\n | \n <\/td>\n | %<\/span> 誤差の算出<\/span><\/td>\n <\/tr>\n \n <\/td>\n | er =<\/span> sum( (Xn<\/span>'<\/span>).^2<\/span> -<\/span> X<\/span>.^<\/span>2<\/span> );<\/td>\n <\/tr>\n \n <\/td>\n | fprintf('<\/span>逆変換誤差 = %f\\n<\/span>'<\/span><\/span>, er<\/span>);<\/td>\n <\/tr>\n \n <\/td>\n | rmsX =<\/span> rms(X<\/span>);<\/td>\n <\/tr>\n \n <\/td>\n | rmsXn =<\/span> rms(Xn<\/span>);<\/td>\n <\/tr>\n \n <\/td>\n | fprintf('<\/span>RMS X=%f<\/span>, Xn=%f\\n<\/span>'<\/span><\/span>, rmsX<\/span>, rmsXn<\/span>);<\/td>\n <\/tr>\n \n <\/td>\n | %<\/span> 返却<\/span><\/td>\n <\/tr>\n \n <\/td>\n | result =<\/span> Xn<\/span>;<\/td>\n <\/tr>\n \n <\/td>\n | \n<\/td>\n <\/tr>\n | \n <\/td>\n | end<\/span><\/td>\n <\/tr>\n \n <\/td>\n | \n<\/td>\n <\/tr>\n | \n <\/td>\n | \n<\/td>\n <\/tr>\n | \n <\/td>\n | \n<\/td>\n <\/tr>\n | \n <\/td>\n | \n<\/td>\n <\/tr>\n <\/table>\n<\/div>\n\n\n <\/div>\n\n <\/div>\n<\/div>\n\n <\/div>\n | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | |