FINS5542-Oliver
Q1：（建议）合并+Diff (复制完检查一下，不要带页码和此行！)
//双斜线后面的是备注，对 Code 无用，这几行不要直接复制
//我列这几行，是想告诉大家，这几行【双斜线】之后的你随便改
//例如像老师教材里一样 this program generates 12-day and 65-day MV and differences…
//下列 code 你可以像老师教材一样的形式调整空格、空行等，【但不要动符号前后位置】
//ranseed 那行我给删了，这个作业里用不到，是一行没有用的随机数生成命令
#include
main()
{
decl i, pdt, mv12av, j, k, mv65av, mvdifav;
pdt=new matrix [2275][19];
pdt=loadmat("pdatau.prn");
mv12av=new matrix [2275][19];
mv65av=new matrix [2275][19];
mvdifav= new matrix [2275][19];
for(i=0;i<19;i++)
{
for (k=0;k<11;k++)
{
mv12av[k][i]=0;
}
for (j=11;j<2275;j++)
{
mv12av[j][i]=pdt[j][i];
for (k=1;k<12;k++)
{
mv12av[j][i]+=pdt[j-k][i];
}
mv12av[j][i]=(mv12av[j][i])/12;
}
for (k=0;k<64;k++)
{
mv65av[k][i]=0;
}
for (j=64;j<2275;j++)
{
mv65av[j][i]=pdt[j][i];
for (k=1;k<65;k++)
FINS5542-Oliver（注意不要复制到页眉，ps：写中文是为了让你一眼看得到）

{
mv65av[j][i]+=pdt[j-k][i];
}
mv65av[j][i]=(mv65av[j][i])/65;
}
}
for (i=0; i<19; i++)
{
mvdifav[][i]=mv12av[][i]-mv65av[][i];
}
for (i=0; i<19; i++)
{
for (k=0; k<64; k++)
{
mvdifav[k][i]=0;
}
}
for (k=0;k<2275;k++)
{
decl file=fopen("mv12ave.out","a");
fprint(file,"\n"," ", mv12av[k][0] ," ",
mv12av[k][1] ," ",
mv12av[k][2] ," ",
mv12av[k][3] ," ",
mv12av[k][4] ," ",
mv12av[k][5] ," ",
mv12av[k][6] ," ",
mv12av[k][7] ," ",
mv12av[k][8] ," ",
mv12av[k][9] ," ",
mv12av[k][10] ," ",
mv12av[k][11] ," ",
mv12av[k][12] ," ",
mv12av[k][13] ," ",
mv12av[k][14] ," ",
mv12av[k][15] ," ",
mv12av[k][16] ," ",
mv12av[k][17] ," ",
mv12av[k][18]);
fclose(file);
FINS5542-Oliver（注意不要复制到页眉，ps：写中文是为了让你一眼看得到）

}
for (k=0;k<2275;k++)
{

decl file=fopen("mv65ave.out","a");
fprint(file,"\n"," ", mv65av[k][0]," ",
mv65av[k][1]," ",
mv65av[k][2]," ",
mv65av[k][3]," ",
mv65av[k][4]," ",
mv65av[k][5]," ",
mv65av[k][6]," ",
mv65av[k][7]," ",
mv65av[k][8]," ",
mv65av[k][9]," ",
mv65av[k][10]," ",
mv65av[k][11]," ",
mv65av[k][12]," ",
mv65av[k][13]," ",
mv65av[k][14]," ",
mv65av[k][15]," ",
mv65av[k][16]," ",
mv65av[k][17]," ",
mv65av[k][18]);
fclose(file);
}

for (k=0;k<2275;k++)
{
decl file=fopen("mvdifave.out","a");
fprint(file,"\n"," ", mvdifav[k][0]," ",
mvdifav[k][1]," ",
mvdifav[k][2]," ",
mvdifav[k][3]," ",
mvdifav[k][4]," ",
mvdifav[k][5]," ",
mvdifav[k][6]," ",
mvdifav[k][7]," ",
mvdifav[k][8]," ",
mvdifav[k][9]," ",
FINS5542-Oliver（注意不要复制到页眉，ps：写中文是为了让你一眼看得到）
mvdifav[k][10]," ",
mvdifav[k][11]," ",
mvdifav[k][12]," ",
mvdifav[k][13]," ",
mvdifav[k][14]," ",
mvdifav[k][15]," ",
mvdifav[k][16]," ",
mvdifav[k][17]," ",
mvdifav[k][18]);
fclose(file);
}
}
Q2：（建议）合并+Diff
#include
main()
{
decl i, pdt, mv45av, j, k, mv250av, mvdifav;
pdt=new matrix [2275][19];
pdt=loadmat("pdatau.prn");
mv45av=new matrix [2275][19];
mv250av=new matrix [2275][19];
mvdifav= new matrix [2275][19];
for(i=0;i<19;i++)
{
for (k=0;k<44;k++)
{
mv45av[k][i]=0;
}
for (j=44;j<2275;j++)
{
mv45av[j][i]=pdt[j][i];
for (k=1;k<45;k++)
{
mv45av[j][i]+=pdt[j-k][i];
}
mv45av[j][i]=(mv45av[j][i])/45;
}
for (k=0;k<249;k++)
FINS5542-Oliver（注意不要复制到页眉，ps：写中文是为了让你一眼看得到）

{
mv250av[k][i]=0;
}
for (j=249;j<2275;j++)
{
mv250av[j][i]=pdt[j][i];
for (k=1;k<250;k++)
{
mv250av[j][i]+=pdt[j-k][i];
}
mv250av[j][i]=(mv250av[j][i])/250;
}
}
for (i=0; i<19; i++)
{
mvdifav[][i]=mv45av[][i]-mv250av[][i];
}
for (i=0; i<19; i++)
{
for (k=0; k<249; k++)
{
mvdifav[k][i]=0;
}
}
for (k=0;k<2275;k++)
{
decl file=fopen("mv45ave.out","a");
fprint(file,"\n"," ",
mv45av[k][0] ," ",
mv45av[k][1] ," ",
mv45av[k][2] ," ",
mv45av[k][3] ," ",
mv45av[k][4] ," ",
mv45av[k][5] ," ",
mv45av[k][6] ," ",
mv45av[k][7] ," ",
mv45av[k][8] ," ",
mv45av[k][9] ," ",
mv45av[k][10] ," ",
mv45av[k][11] ," ",
FINS5542-Oliver（注意不要复制到页眉，ps：写中文是为了让你一眼看得到）
mv45av[k][12] ," ",
mv45av[k][13] ," ",
mv45av[k][14] ," ",
mv45av[k][15] ," ",
mv45av[k][16] ," ",
mv45av[k][17] ," ",
mv45av[k][18]);
fclose(file);
}
for (k=0;k<2275;k++)
{
decl file=fopen("mv250ave.out","a");
fprint(file,"\n"," ",
mv250av[k][0]," ",
mv250av[k][1]," ",
mv250av[k][2]," ",
mv250av[k][3]," ",
mv250av[k][4]," ",
mv250av[k][5]," ",
mv250av[k][6]," ",
mv250av[k][7]," ",
mv250av[k][8]," ",
mv250av[k][9]," ",
mv250av[k][10]," ",
mv250av[k][11]," ",
mv250av[k][12]," ",
mv250av[k][13]," ",
mv250av[k][14]," ",
mv250av[k][15]," ",
mv250av[k][16]," ",
mv250av[k][17]," ",
mv250av[k][18]);
fclose(file);
}
for (k=0;k<2275;k++)
{
decl file=fopen("mvdifaveQ2.out","a");
fprint(file,"\n"," ",
mvdifav[k][0]," ",
FINS5542-Oliver（注意不要复制到页眉，ps：写中文是为了让你一眼看得到）
mvdifav[k][1]," ",
mvdifav[k][2]," ",
mvdifav[k][3]," ",
mvdifav[k][4]," ",
mvdifav[k][5]," ",
mvdifav[k][6]," ",
mvdifav[k][7]," ",
mvdifav[k][8]," ",
mvdifav[k][9]," ",
mvdifav[k][10]," ",
mvdifav[k][11]," ",
mvdifav[k][12]," ",
mvdifav[k][13]," ",
mvdifav[k][14]," ",
mvdifav[k][15]," ",
mvdifav[k][16]," ",
mvdifav[k][17]," ",
mvdifav[k][18]);
fclose(file);
}
}