現在、C言語でチェビシェフ窓を計算することを考えてます。
一応、エラーもなく実行できるのですが、中央の値が1になりません。
これで正しいのかどうかもわからないので、どなたか教えていただけないでしょうか。
よろしくお願いいたします。
///////////////////以下 プログラム /////////////////
#include <stdio.h>
#include <math.h>
/**
* 逆双曲余弦関数の定義
*/
double Acosh(double x);
/**
*チェビシェフ関数の定義
* n:関数の次数
*/
double Chev(int n, double x);
void main(void)
{
int i, N=128, k, M=(int)floor((N-1)/2);
///data:入力 w:出力(ウィンドウ係数) ii ,kk インデックスの倍精度型, nn サンプル数の倍精度型
double data[128], w[128], ii, kk, nn;
// s :メインローブ対サイドローブの比
double t0, s, u, sum;
////円周率の定義
const double pi=3.1415926;
nn=(double)N, s=1000.0;
//入力をすべて1に初期化
for(i=0;i<N;i++){data[i]=1.0;}
/////チェビシェフ窓の計算ルーチンと出力
t0=cosh((1.0/(nn-1))*Acosh(s));
for(i=0;i<N;i++){
ii=(double)i;
for(k=1;k<=M;k++){
kk=(double)k;
u=ii-(nn-1)/2.0;
if(k==1){
sum=Chev(N-1,t0*cos(kk*pi/nn))*cos(2.0*pi*kk*u/nn);
}else{
sum+=Chev(N-1,t0*cos(kk*pi/nn))*cos(2.0*pi*kk*u/nn);
}
}
w[i]=data[i]*(s+2.0*sum)/nn;
printf("w[%d]=%4.3f\n",i,w[i]);
}
}
double Acosh(double x)
{
double x2=x*x;
return log(x+sqrt(x2-1));
}
double Chev(int n, double x)
{
if(fabs(x)<=1){
return cos((double)n*acos(x));
}else{
return cosh((double)n*Acosh(x));
}
}
/////////////// ここまで ///////////////////
//////////////////以下出力結果/////////////
w[0]=0.628
w[1]=0.285
w[2]=0.348
w[3]=0.419
w[4]=0.499
w[5]=0.588
w[6]=0.686
w[7]=0.795
w[8]=0.914
w[9]=1.044
[中略]
w[50]=14.268
w[51]=14.559
w[52]=14.832
w[53]=15.085
w[54]=15.320
w[55]=15.533
w[56]=15.725
w[57]=15.895
w[58]=16.041
w[59]=16.164
w[60]=16.264
w[61]=16.338
w[62]=16.388
w[63]=16.413<-中央値
w[64]=16.413<-中央値
w[65]=16.388
[中略]
w[121]=0.686
w[122]=0.588
w[123]=0.499
w[124]=0.419
w[125]=0.348
w[126]=0.285
w[127]=0.628
No.1ベストアンサー
- 回答日時:
各 w[ i ] を 最大値 16.413 で割ってやらないといけない、ということじゃないんでしょうか。
http://www.dsprelated.com/showarticle/42.php
のページの w(k) ( 右から2列目 ) を出力されていると思います。
(実際に、N = 128 を N = 9 に変更してみるとわかると思います。)
参考URL:http://www.dsprelated.com/showarticle/42.php
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- C言語・C++・C# c言語でユーザ関数を利用して複素数のべき乗と絶対値の数列を計算するプログラムが作りたいです。 3 2023/01/29 22:13
- C言語・C++・C# 10個の実数に対する降順ソート結果を出力するプログラムを作りたいのですが、以下のプログラムをどう直せ 1 2022/07/09 22:16
- C言語・C++・C# C 言語の Gauss Jordan 法について 2 2022/12/28 11:16
- C言語・C++・C# プログラミングの授業の課題です 1 2023/01/17 22:15
- C言語・C++・C# Cのdoubleの浮動小数点表示について 3 2023/04/17 13:14
- C言語・C++・C# C言語のエラーについて 2 2022/07/11 13:56
- C言語・C++・C# C言語のマクローリン展開ローラン展開のコードについて 3 2022/12/15 14:45
- C言語・C++・C# バイナリファイルをコピーするのにかかる時間を測りたいのですが実行するとFatel error:gli 2 2022/11/03 01:10
- C言語・C++・C# このプログラミングの問題を教えてほしいです。 キーボードからデータ数nとn個のデータを入力し、平均値 3 2022/12/19 22:51
- C言語・C++・C# C言語初心者 構造体 課題について 1 2023/03/10 19:30
関連するカテゴリからQ&Aを探す
おすすめ情報
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
プログラムでの数字につく”f”の...
-
doubleの変数にintとintの割り...
-
C言語を実行すると-infが出てき...
-
float型とdouble型の変数の違い...
-
テイラー展開(C言語)
-
C言語 関数プロトタイプ宣言の...
-
c言語で、繰り返し文の中で、0....
-
C言語で-23乗を取り扱うには
-
c言語のコンパイルエラー canno...
-
関数におけるif文とreturn文に...
-
C 開放してるのにエラー(doubl...
-
c言語 標準体重計算のプログラ...
-
doubleは常に%lfとするべきなのか
-
long doubleの表示方法
-
C言語でdouble型の小数点の引き...
-
待ち行列シュミレーションのプ...
-
C言語(プログラミング)関連の質...
-
計算速度をはやくするには
-
処理時間計測
-
ニュートン法
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
プログラムでの数字につく”f”の...
-
doubleの変数にintとintの割り...
-
C 開放してるのにエラー(doubl...
-
Cで3乗根を求める方法
-
float型とdouble型の変数の違い...
-
至急です! マクロ定義で #defi...
-
C言語の型による処理速度の違い
-
int とdoubleの比較
-
関数におけるif文とreturn文に...
-
C言語初心者 構造体 課題について
-
c言語のコンパイルエラー canno...
-
C言語 関数プロトタイプ宣言の...
-
C言語を実行すると-infが出てき...
-
float?数字の後にLがつくもの
-
数値を指数部と仮数部に分離したい
-
difftime()について
-
浮動小数点数が表示されないん...
-
たくさんの数の平均を求める方...
-
DWORDの警告
-
-1.#IND00と出てしまうのですが...
おすすめ情報