求一个已知的四面体体内任意取四点所构成新四面体的体积的期望值（详见内） - C语言论坛

已结贴√ 问题点数：20 回复次数：9

求一个已知的四面体体内任意取四点所构成新四面体的体积的期望值（详见内）

计算物理作业，求一个已知的四面体,在其体内任意取四点所构成新四面体的体积大小的期望值，原四面体四点坐标为（0，0，0),(1,0,0),(0,1,0)(0,0,1)

我现在已经写出一个求四面体体积的程序，我想问下是不是设定一个次数m，任意取m次新四面体，再求所有取得的新四面体的体积的平均值就行呢？
有人能点拨我一下吗？万分感谢
Q：369020275

搜索更多相关主题的帖子: 期望值　体积　

第 2 楼

得分:0

各位有知道的吗

第 3 楼

得分:0

帮帮忙啊

第 4 楼

得分:14

没有这种说法，你想想在一个四面体中取一个新四面体有无数种可能，这不符合期望值的定义

多少恨，昨夜梦魂中。还似旧时游上苑，车如流水马如龙；花月正春风！

第 5 楼

得分:0

可是这是我们老师给的期末作业题目啊，应该没错吧......我现在做了个这个程序，要输入想取得的小四面体的个数，再取其平均值，但是我总觉得不对。。。
#include <conio.h>
#include<stdlib.h>
#include<time.h>
#include <iostream>
#include  <math.h>
using namespace std;
double  main()
{
while(1)
{double v=0;
int m;
int n;
cout<<"请输入想要取的小四面体个数:"<<endl;
cin>>n;
cout<<n<<"个小四面体"<<endl;
srand((unsigned)time(NULL));
int i=0;
double N[4][3];
double k;
for (m=0;m<n;m++)
{do
{double a=fabs((rand()%1000)*0.001),b=fabs((rand()%1000)*0.001),c=fabs((rand()%1000)*0.001);
    k=a+b+c-1;
    if(k<0)
    {N[i][0]=a;
    N[i][1]=b;
    N[i][2]=c;
    i++;}
}while(i<4);
if(i=4){double x1,y1,z1,x2,y2,z2,x3,y3,z3,x4,y4,z4,A1,A2,A3,A4;
x1=N[0][0];y1=N[0][1];z1=N[0][2];
x2=N[1][0];y2=N[1][1];z2=N[1][2];
x3=N[2][0];y3=N[2][1];z3=N[2][2];
x4=N[3][0];y4=N[3][1];z4=N[3][2];
A1=x2*(y3*z4-y4*z3)-x3*(y2*z4-y4*z2)+x4*(y2*z3-y3*z2);
A2=x1*(y3*z4-y4*z3)-x3*(y1*z4-y4*z1)+x4*(y1*z3-y3*z1);
A3=x1*(y2*z4-y4*z2)-x2*(y1*z4-y4*z1)+x4*(y1*z2-y2*z1);
A4=x1*(y2*z3-y3*z2)-x2*(y1*z3-y3*z1)+x3*(y1*z2-y2*z1);
v=fabs((A2-A1+A4-A3)/6);
v+=v;
i=0;
}
}
double E;
E=v/n;
cout<<"E="<<E<<endl;
cout<<endl;
}
return 0;
}