お願いします
Re: お願いします
この式は意味がわかりません。tarou さんが書きました:f(x)=exp(-x2/2)/(2π)1/2
乗算の部分は *をつけてください。
-x2/2 ってのは何でしょう?
(2π)1/2
の部分もどう判断してよいのかわからない。
で、仮に式がわかったとして、微分や積分のやり方を知りたいのですか?
non
Re: お願いします
追伸
/ はそのまま書いてあったのでスイマセンよくわかりません、
* はつけてなかったです
f(x)=exp*(-x*x/2)/(-2π)*2分の一
よくわからないんですが、こうだとおもいます(2分の一打ち方わかりません)
/ はそのまま書いてあったのでスイマセンよくわかりません、
* はつけてなかったです
f(x)=exp*(-x*x/2)/(-2π)*2分の一
よくわからないんですが、こうだとおもいます(2分の一打ち方わかりません)
Re: お願いします
一応式はあるんですけど書き換えがよくわからなくいです><
何回もカキコんでしまってすいません
よろしくお願いします
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
extern double dif();
int main()
{
double xa, xb, x, dx, y, dy;
int i;
static int nloop = 10; /* number of steps */
xa = -1.0e+1.0; /* x_minimum */
xb = -2.0e+2.0; /* x_maximum */
dx = (xb - xa) / ((double) nloop); /* step size of integration */
y = ((double) 0);
x = xa + 0.5 * dx;
/* main part */
for(i = 0;i < nloop;i++)
{
dy = dif(x);
y += (dy * dx);
x += dx;
}
printf("The answer is %e \n", y);
exit(0);
}
double dif(double x)
{
/* an equation to be integrated */
double a;
a =exp(x * x /2)/(2* \pi) * 1/2 ;
return a;
}
何回もカキコんでしまってすいません
よろしくお願いします
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
extern double dif();
int main()
{
double xa, xb, x, dx, y, dy;
int i;
static int nloop = 10; /* number of steps */
xa = -1.0e+1.0; /* x_minimum */
xb = -2.0e+2.0; /* x_maximum */
dx = (xb - xa) / ((double) nloop); /* step size of integration */
y = ((double) 0);
x = xa + 0.5 * dx;
/* main part */
for(i = 0;i < nloop;i++)
{
dy = dif(x);
y += (dy * dx);
x += dx;
}
printf("The answer is %e \n", y);
exit(0);
}
double dif(double x)
{
/* an equation to be integrated */
double a;
a =exp(x * x /2)/(2* \pi) * 1/2 ;
return a;
}
Re: お願いします
>式がわからないといった感じです
http://detail.chiebukuro.yahoo.co.jp/qa ... 1466599879
こちらの人は式はわかっているみたい。
でも
a =exp(x * x /2)/(2* \pi) * 1/2 ;
とは違うので、別の問題かもしれませんね。
http://detail.chiebukuro.yahoo.co.jp/qa ... 1466599879
こちらの人は式はわかっているみたい。
でも
a =exp(x * x /2)/(2* \pi) * 1/2 ;
とは違うので、別の問題かもしれませんね。
Re: お願いします
同じ人だと思ったけど、別人さんでしたか。
だとすれば、
問題文 f(x)=exp(-x2/2)/(2π)1/2
を出題した先生がおかしい。
「式が意味不明なので解けません」が正解。
#同じ学校だからでしょうか、「><」を文末につけるのが流行りなのですか?><
だとすれば、
問題文 f(x)=exp(-x2/2)/(2π)1/2
を出題した先生がおかしい。
「式が意味不明なので解けません」が正解。
#同じ学校だからでしょうか、「><」を文末につけるのが流行りなのですか?><
Re: お願いします
あっちのハンドルネームはkyouto_univでしたが、まさか京都大学ってことはないよね。
式がわからないって書いてあったから、書いたのに他に何が必要なの?
>これで次はどうすればいいんですか?
>公式もうえのやり方で解けますか?
大学生の質問だとは思えない。
数学的にx=0のf(0)を求めて、それを1/2にし、そのときの解を求める。
後は、定積分でしょ。
式がわからないって書いてあったから、書いたのに他に何が必要なの?
>これで次はどうすればいいんですか?
>公式もうえのやり方で解けますか?
大学生の質問だとは思えない。
数学的にx=0のf(0)を求めて、それを1/2にし、そのときの解を求める。
後は、定積分でしょ。
non
Re: お願いします
あちらの方と私同じ人じゃないんでしりませんが・・・
学部が理系じゃないのでさっぱりなんですよね
上のプログラムでやっているんですがコンパイルしてもエラーがでてしまいます
学部が理系じゃないのでさっぱりなんですよね
上のプログラムでやっているんですがコンパイルしてもエラーがでてしまいます
Re: お願いします
あっちの人の問題をよく見たらルートの範囲が違ってた。
f(x)=exp(-x^2/2)/(2*pi)^(1/2)
だね。
グラフは、見た目はあまり変わらないけど、このようになります。
f(x)=exp(-x^2/2)/(2*pi)^(1/2)
だね。
グラフは、見た目はあまり変わらないけど、このようになります。
- 添付ファイル
-
- 無題.png (39.09 KiB) 閲覧数: 5754 回
non
Re: お願いします
>上のプログラムでやっているんですがコンパイルしてもエラーがでてしまいます
>xa = -1.0e+1.0; /* x_minimum */
>xb = -2.0e+2.0; /* x_maximum */
は
xa = -1.0e+1; /* x_minimum */
xb = -2.0e+2; /* x_maximum */
にしないとダメです。
この場合、x_maximum < x_minimum なのは別問題ですが。
もう一つ、
a =exp(x * x /2)/(2* \pi) * 1/2 ;
の中の「\pi」 はいけませんね。
加えてpi は未定義となりますので
double pi = 3.14159265359;
とか
#define pi 3.14159265359
とかしておかないとエラーになります。
ところで、上のプログラムとやらは
何をしようとしているのかはご存知でしょうか?
>xa = -1.0e+1.0; /* x_minimum */
>xb = -2.0e+2.0; /* x_maximum */
は
xa = -1.0e+1; /* x_minimum */
xb = -2.0e+2; /* x_maximum */
にしないとダメです。
この場合、x_maximum < x_minimum なのは別問題ですが。
もう一つ、
a =exp(x * x /2)/(2* \pi) * 1/2 ;
の中の「\pi」 はいけませんね。
加えてpi は未定義となりますので
double pi = 3.14159265359;
とか
#define pi 3.14159265359
とかしておかないとエラーになります。
ところで、上のプログラムとやらは
何をしようとしているのかはご存知でしょうか?
Re: お願いします
どっからどう見ても正規分布だと思うのですが。
ガウス分布、確率密度関数、この辺でググったら幸せになれると思いますよ。
何をやってるかですか?
ただの区分求積だと思いますね。
ガウス分布、確率密度関数、この辺でググったら幸せになれると思いますよ。
何をやってるかですか?
ただの区分求積だと思いますね。
- tk-xleader
- 記事: 158
- 登録日時: 14年前
- 連絡を取る:
Re: お願いします
とりあえず2問目と3問目について数学的なヒント、というよりも2問目はもはや答えか…
2問目は方程式の問題なので、xを少しずつ動かして、最も近い値を解とでもしてしまうか、もしくはニュートン法を使うのがいいと思います。
3問目は定積分の問題ですが、申し訳ないのですが、どうしてもf(x)の不定積分が求まりませんでした。そこで、f(x)のグラフを短冊みたいに無限に細かく切って、その1本1本の面積(無限に切っているので長方形とみなしても問題はないはず。)の和を求めるという方法をとりました。
そのため、定積分の最終的な答えがΣの極限の形になっています。
2問目は方程式の問題なので、xを少しずつ動かして、最も近い値を解とでもしてしまうか、もしくはニュートン法を使うのがいいと思います。
3問目は定積分の問題ですが、申し訳ないのですが、どうしてもf(x)の不定積分が求まりませんでした。そこで、f(x)のグラフを短冊みたいに無限に細かく切って、その1本1本の面積(無限に切っているので長方形とみなしても問題はないはず。)の和を求めるという方法をとりました。
そのため、定積分の最終的な答えがΣの極限の形になっています。
- 添付ファイル
-
- answer.png (30.83 KiB) 閲覧数: 5614 回
Re: お願いします
ガウス関数の積分は0→∞の範囲でしたらうまくいくんですけどね。まんまガウス積分ですけど・・・tkmakwins15 さんが書きました:どうしてもf(x)の不定積分が求まりませんでした。
D(t) = { (x,y) | x^2 + y^2 ≦ t^2, x≧0, y ≧ 0 }
E(t) = { (x,y)| 0≦x≦t, 0≦y≦t }
D(t) ⊂ E(t) ⊂ D( (√2)*t ) より
∬[ D(t) ] { e^-(x^2+y^2) } dxdy ≦ ∬[ E(t) ] { e^-(x^2+y^2) } dxdy ≦ ∬[ D( (√2)*t ) ] { e^-(x^2+y^2) } dxdy
ここで極座標形式 x = rcosΘ、 y = rsinΘ を考えるとヤコビアン ∂(x,y)/∂(r,Θ) = rより
∬[ D(t) ] { e^-(x^2+y^2) } dxdy
= ∫[0→π/2] dΘ ∫[0→t] {e^(-r^2)} rdr
= π{1- e^(-t^2) }/4
→π/4 (t→∞)
∬[ E(t) ] { e^-(x^2+y^2) } dxdy = ( ∫[0→∞] e^(-x^2)dx )^2 だから
∫[-∞→∞] e^(-x^2)dx
=2*∫[0→∞] e^(-x^2)dx
→ √π
だから質問者さんの奴は
X = x/√2 で置き換えて√2dX = dx だから
全部積分すると
「1」になるんですね。なんと美しいかな・・・
ちなみに本題ですが
「標準正規分布表」という便利なものを使えばプログラムに頼らずとも正確な答えを導出できますよ(3番)
まぁちょっとしたことするには使えますね。