作者erica0707 (快快乐乐)
看板C_and_CPP
标题[问题] 模拟太多次就会出错
时间Wed Mar 18 16:49:41 2009
小的最近刚学写程式
刚刚遇到一个问题 一直无法解决 想请板上各位高手帮忙侦错一下
大约模拟超过2万次 视窗就会自己关闭 实在是很奇怪...也不知道错在哪
以下是完整程式码 不好意思有点长
先谢谢大家了
#include<stdio.h>
#include<math.h>
#include<stdlib.h>
#include<time.h>
double normal()
{
int i ;
double N;
for(i=0;i<1;i++)
{
N=double(rand())/RAND_MAX;
}
return N;
}
int main(void)
{
int i,k,j;
int n;
printf("请输入要模拟的次数:");
scanf("%d",&n);
double S[5][5];
S[0][0]=0.0006058718;
S[0][1]=0.000339878;
S[0][2]=0.0003308925;
S[0][3]=0.00004346825;
S[0][4]=0.0003617589;
S[1][0]=0.000339878;
S[1][1]=0.001853979;
S[1][2]=0.0007299478;
S[1][3]=0.001010278;
S[1][4]=0.0006364092;
S[2][0]=0.0003308925;
S[2][1]=0.0007299478;
S[2][2]=0.001150133;
S[2][3]=0.0005555975;
S[2][4]=0.0008855597;
S[3][0]=0.00004346825;
S[3][1]=0.001010278;
S[3][2]=0.0005555975;
S[3][3]=0.001122313;
S[3][4]=0.0003779205;
S[4][0]=0.0003617589;
S[4][1]=0.0006364092;
S[4][2]=0.0008855597;
S[4][3]=0.0003779205;
S[4][4]=0.0009437375;
/*
for(i=0;i<5;i++)
{
for(j=i;j<5;j++)
{
printf("请输入第%d和第%d个资产共变异度",i+1,j+1);
scanf("%lf",&S[i][j]);
S[j][i]=S[i][j];
}
}
*/
double W[n][5],sum,N[5];
srand((unsigned)time(NULL));
for(i=0;i<n;i++)
{
for(sum=0,k=0;k<5;k++)
{
N[k]=normal();
sum+=N[k];
}
for(j=0;j<5;j++)
{
N[j]=N[j]/sum ;
W[i][j]=N[j];
}
}
double Y[n][5];
for(i=0;i<n;i++)
{
for(j=0;j<5;j++)
{
for(Y[i][j]=0,k=0;k<5;k++)
{
Y[i][j]+=W[i][k]*S[k][j];
}
}
}
double var[n];
for(i=0;i<n;i++)
{
for(var[i]=0,j=0;j<5;j++)
{
var[i]+=Y[i][j]*W[i][j];
}
}
int h;
double min;
for(min=100000,i=0;i<n;i++)
{
if(var[i]<min)
{
min=var[i];
h=i;
}
}
printf("portfolio-var=%lf\n",min);
for(j=0;j<5;j++)
{
printf("第%d个资产配置为%lf\n",j+1,W[h][j]);
}
getchar();
getchar();
return 0;
}
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.114.225.29
1F：→ MOONRAKER:到command line底下执行 至少要看到结束执行时的讯息 03/18 16:53
2F：→ TConvertible:double W[n][5],double var[n],double Y[n][5] 03/18 18:01
3F：→ TConvertible:这几行,请改用 动态记忆体配置. 不然 n 值太大时, 03/18 18:02
4F：→ TConvertible:Stack 就爆了,程式就当了. 有问题可以来信/水球发问 03/18 18:03
5F：→ erica0707:T大真的很热心 我问题已经解决了 板上真的很多热心高手! 03/18 19:57