質問をすることでしか得られない、回答やアドバイスがある。

15分調べてもわからないことは、質問しよう!

新規登録して質問してみよう
ただいま回答率
85.35%
C

C言語は、1972年にAT&Tベル研究所の、デニス・リッチーが主体となって作成したプログラミング言語です。 B言語の後継言語として開発されたことからC言語と命名。そのため、表記法などはB言語やALGOLに近いとされています。 Cの拡張版であるC++言語とともに、現在世界中でもっとも普及されているプログラミング言語です。

Q&A

解決済

1回答

435閲覧

フーリエ級数展開の展開係数について

murimurikimoi

総合スコア1

C

C言語は、1972年にAT&Tベル研究所の、デニス・リッチーが主体となって作成したプログラミング言語です。 B言語の後継言語として開発されたことからC言語と命名。そのため、表記法などはB言語やALGOLに近いとされています。 Cの拡張版であるC++言語とともに、現在世界中でもっとも普及されているプログラミング言語です。

0グッド

0クリップ

投稿2020/05/30 11:22

編集2020/05/30 11:30

前提・実現したいこと

f(t)=-1(-π<=t<=0)
= 1( 0<=t<=π)
のフーリエ級数展開をC言語で求めようとしています。しかし、anが0になるはずが0にならず困っています。下のa[n]のプログラム部分で間違っているところがあればご指摘頂けると幸いです。積分は台形公式で求めています。

発生している問題・エラーメッセージ

結果一部: a[5] が -2.16178e^(-16) 0にならない

該当のソースコード

C言語 ソースコード 一部分 ↓   #define A 10000 #define N 50 /* a[n] */ double a_1, a_2, x_1, x_2, a[A]; double f_1 = -1.0, f_2 = 1.0, f_3, f_4; double T = 2.0 * M_PI, w = 2.0 * M_PI / T; double c = -M_PI, d = M_PI; dx = (d - c) / A; for(int i = 1; i <= N; i++){ a_1 = 0.0; a_2 = 0.0; for(int j = 0; j <A; j++){ x_1 = c + dx * j; x_2 = c + dx * (j + 1); if(0 <= j &&j <A / 2){ f_3 = f_1 * cos(i * w * x_1); f_4 = f_1 * cos(i * w * x_2); a_1 += (f _3 + f _4) * dx / 2.0; } else if(A / 2 <= j &&j <A){ f_3 = f_2 * cos(i * w * x_1); f_4 = f_2 * cos(i * w * x_2); a_2 += (f _3 + f _4) * dx / 2.0; } } a[i] = (2.0 / T) * (a_1 + a_2); }

気になる質問をクリップする

クリップした質問は、後からいつでもMYページで確認できます。

またクリップした質問に回答があった際、通知やメールを受け取ることができます。

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

cateye

2020/05/30 13:59

float.hで定義されている、(マクロ)DBL_EPSILON以下の数値は判別不能です。
guest

回答1

0

ベストアンサー

数値計算においては誤差は避けられません。まず、誤差許容範囲を定めましょう。

投稿2020/05/30 11:51

HogeAnimalLover

総合スコア4830

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

murimurikimoi

2020/05/30 16:52

わかりました!ありがとうございます!
guest

あなたの回答

tips

太字

斜体

打ち消し線

見出し

引用テキストの挿入

コードの挿入

リンクの挿入

リストの挿入

番号リストの挿入

表の挿入

水平線の挿入

プレビュー

15分調べてもわからないことは
teratailで質問しよう!

ただいまの回答率
85.35%

質問をまとめることで
思考を整理して素早く解決

テンプレート機能で
簡単に質問をまとめる

質問する

関連した質問