漸化式で計算を求めるプログラムを作成したいのですがわからないところが出てしまったのでご教授のほうお願いします。
levin()関数の中で漸化式の計算を行っているのですが、その関数の中身の●の部分をどう表現したらいいのかわかりません。仮引数や実引数をどのように設定するか、もう一つはこのプログラムで漸化式の『alpha』を求めることができるのか、また間違っているなら指摘していただけるとありがたいです。プログラムはめちゃくちゃだと思いますがよろしくお願いします。
#include<stdio.h>
#define SIZE 100
#define M 1
double levin(●, ●);
int main(){
double a[SIZE];
double r[SIZE];
double alpha;
int i;
int m;
a[0] = a[1] = 0.0;
m = 0;
alpha = r[0];
for(i = 0; i < M; i++){
levin(●, ●);
}
return 0;
}
double levin(●, ●){
int m = 50;
int k;
double beta = 0.0;
double km;
double a[SIZE];
double r[SIZE];
double new_a[SIZE];
double alpha;
for(k = 0; k <= m; k++){
beta += ((a[k] * (r[m + 1 - k]));
}
km = - (beta / alpha);
new_a[0] = 1;
new_a[++m] = 0;
for(k = 0; k <= m; k++){
new_a[k] = ((a[k]) +((km) * (a[m - k])));
}
alpha *= 1 - ((km) * (km));
printf("%lf\n", alpha);
a = new_a;
return 0;
}
No.1ベストアンサー
- 回答日時:
まずは、簡単なところから、インデントしたほうがいいですよ。
{}のなかをTABキーで一段さげる。そうするとずっと見やすくなります。ちなみに関数の直後に{をつけるのは一般的ではない気がします。
●、●の部分ですが、引数を利用して関数の呼び出し元に計算結果を返すかどうかで変わってきます。値を戻すのならば、ポインタを使用します。
関数は、double levin(●, ●)と宣言されていますが、
実際にreturn 0;としており、levin()関数は、必ずゼロを返しています。
XXX = levin(aaa, bbb); でXXXに計算結果を期待しているのでしたら、このプログラムは、間違えています。
ここは、私の知識不足かもしれませんが、どのような漸化式か提示いただけないとプログラムの内容(アルゴリズム)が正しいかどうか判断できないと思います。
この回答への補足
計算機で高速に線形予測を行うレビンソン-ダービンの算法を用いています。このプログラムを今回作成しようと思い合っているかどうかを投稿させていただきました。
以下に、アルゴリズムである漸化式を表記します。
(初期設定)
a0^(0)=1, α0=r0
(漸化式)
M - 1 に達するまでこれらの式を繰り返します。
β_m = Σak^(m)・rm+1-k
K_m+1 = -βm / αm
a_0^(m+1) = 1
a_k^(m+1) = ak^(m)+Km+1・am+1-k^(m) (k=1,2,...,m+1)
α_m+1 = (1-Km+1^2)αm
_のついた横のmやm+1や0やkは、a0,a1,a2,a3などの順を指します。
^(m+1)と^(m)のかっこの中身は予測係数を表します。
また、Km+1^2の2は2乗を表します。
よろしくお願いします。
ご回答ありがとうございます。
はい。ご指摘ありがとうございます。これからは直すようにします。
このlevin()関数では、return 0;は意味ないですよね。
補足のほうに、漸化式を表示しますので、よろしくお願いします。
No.2
- 回答日時:
基本的に、プログラムがあっているかどうかを人に聞くものではありません。
たとえば私が「あっていますよ」と言えば、「そうか、あってるのか、よかった」となるのでしょうか?ちがいますよね。
あっているのかを知りたければ、手計算で10個ほど計算してみてプログラムの入力に入れてみるのです。出力値が手計算結果と同じになれば合っていると証明できますよね。
これが論理的なものの考え方ではないでしょうか。
回答ありがとうございます。
確かにそのとおりですね。
自分で作成したプログラムを簡単な数字など入れてトレースすることも大切ですね。
しかし、わからないところは考えてもわからないので助言がほしかったので投稿したことはわかってくださいペコリ(o_ _)o))
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- C言語・C++・C# 10個の実数に対する降順ソート結果を出力するプログラムを作りたいのですが、以下のプログラムをどう直せ 1 2022/07/09 22:16
- C言語・C++・C# c言語でユーザ関数を利用して複素数のべき乗と絶対値の数列を計算するプログラムが作りたいです。 3 2023/01/29 22:13
- C言語・C++・C# 並列プログラミングのπ計算について 1 2022/07/16 22:30
- C言語・C++・C# C 言語の Gauss Jordan 法について 2 2022/12/28 11:16
- C言語・C++・C# Cのdoubleの浮動小数点表示について 3 2023/04/17 13:14
- C言語・C++・C# LU分解法のピボッティングについて(C言語/gcc-9) 3 2022/07/11 23:10
- C言語・C++・C# プログラミングの授業の課題です 1 2023/01/17 22:15
- C言語・C++・C# バイナリファイルをコピーするのにかかる時間を測りたいのですが実行するとFatel error:gli 2 2022/11/03 01:10
- C言語・C++・C# LU分解法のピボット選択機能実装について(C言語・gcc-9) 1 2022/07/22 15:20
- C言語・C++・C# C言語 プログラミング 4 2022/05/22 11:53
関連するカテゴリからQ&Aを探す
おすすめ情報
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
c言語で、繰り返し文の中で、0....
-
プログラムでの数字につく”f”の...
-
C言語のプログラムで#include<m...
-
C言語で直角三角形の斜辺を求め...
-
float型とdouble型の変数の違い...
-
Cで3乗根を求める方法
-
C言語初心者 構造体 課題について
-
2分法で方程式の複数の解を自...
-
至急です! マクロ定義で #defi...
-
doubleの変数にintとintの割り...
-
C言語を実行すると-infが出てき...
-
doubleは常に%lfとするべきなのか
-
浮動小数点数が表示されないん...
-
difftime()について
-
関数におけるif文とreturn文に...
-
C++で外積
-
(C,C++言語)関数の引数は自動キ...
-
C言語で表記についの質問です
-
数値を指数部と仮数部に分離したい
-
方程式を2分法を用いて解くプロ...
マンスリーランキングこのカテゴリの人気マンスリー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と出てしまうのですが...
おすすめ情報