//else
//{
//printf(“%f\n”,AAD3[i]);
//fprintf(fp4,“%f\n”,AAD3[i]);
//}
}
for(i=0;i《dad3;i++)
{//if(i==0)
printf(“%s\n”,“DAD3”);
printf(“%f\n”,DAD3[i]);
fprintf(fp5,“%f\n”,DAD3[i]);
//}
//else
//{
//printf(“%f\n”,DAD3[i]);
//fprintf(fp5,“%f\n”,DAD3[i]);
//}
}
for(i=0;i《add3;i++)
{//if(i==0)
printf(“%s\n”,“ADD3”);
printf(“%f\n”,ADD3[i]);
fprintf(fp6,“%f\n”,ADD3[i]);
//}
//else
//{
//printf(“%f\n”,ADD3[i]);
//fprintf(fp6,“%f\n”,ADD3[i]);
//}
//
}
for(i=0;i《ddd3;i++)
{//if(i==0)
printf(“%s\n”,“DDD3”);
printf(“%f\n”,DDD3[i]);
fprintf(fp7,“%f\n”,DDD3[i]);
//}
//else
//{
//printf(“%f\n”,DDD3[i]);
//fprintf(fp7,“%f\n”,DDD3[i]);
//}
}
//关闭文件
fclose(fp);
fclose(fp0);
fclose(fp1);
fclose(fp2);
fclose(fp3);
fclose(fp4);
fclose(fp5);
fclose(fp6);
fclose(fp7);
}
/*runthisprogramusingtheconsolepauseroraddyourowngetch,system(“pause”)orinputloop*/
/*尝试不对
doubleA1[(LENGTH+DB_LENGTH)/2];
doubleD1[(LENGTH+DB_LENGTH)/2];
inta1=(LENGTH+DB_LENGTH)/2;
intd1=(LENGTH+DB_LENGTH)/2;
doubleAA2[(LENGTH+DB_LENGTH)/4];
doubleDA2[(LENGTH+DB_LENGTH)/4];
inta2=(LENGTH+DB_LENGTH)/4;
intd2=(LENGTH+DB_LENGTH)/4;
doubleAAA3[(LENGTH+DB_LENGTH)/8];
doubleDAA3[(LENGTH+DB_LENGTH)/8];
inta3=(LENGTH+DB_LENGTH)/8;
intd3=(LENGTH+DB_LENGTH)/8;
doubleAAAA4[(LENGTH+DB_LENGTH)/16];
doubleDAAA4[(LENGTH+DB_LENGTH)/16];
Covlution(input,h,temp,n,m);
for(i=0;i《n+m-2;i+=2)
{
A1[i/2]=temp[i];//一层分解的低频部分,进行了2抽值,即尺度空间,低频概貌部分
}
Covlution(input,g,temp,n,m);
for(i=0;i《n+m-2;i+=2)
{
D1[i/2]=temp[i];//一层分解的高频部分,已经进行了2抽取,即高频细节部分
}
Covlution(A1,h,temp,a1,m);
for(i=0;i《a1+m-2;i+=2)
{
AA2[i/2]=temp[i];//一层分解的低频部分,进行了2抽值,即尺度空间,低频概貌部分
}
Covlution(A1,g,temp,a1,m);
for(i=0;i《a1+m-2;i+=2)
{
DA2[i/2]=temp[i];//一层分解的低频部分,进行了2抽值,即尺度空间,低频概貌部分
}
Covlution(AA2,h,temp,a2,m);
for(i=0;i《a2+m-2;i+=2)
{
AAA3[i/2]=temp[i];//一层分解的低频部分,进行了2抽值,即尺度空间,低频概貌部分
}
Covlution(AA2,g,temp,a2,m);
for(i=0;i《a2+m-2;i+=2)
{
DAA3[i/2]=temp[i];//一层分解的低频部分,进行了2抽值,即尺度空间,低频概貌部分
}
Covlution(AAA3,h,temp,a3,m);
for(i=0;i《a3+m-2;i+=2)
{
AAAA4[i/2]=temp[i];//一层分解的低频部分,进行了2抽值,即尺度空间,低频概貌部分
}
Covlution(AAA3,g,temp,a3,m);
for(i=0;i《a3+m-2;i+=2)
{
DAAA4[i/2]=temp[i];//一层分解的低频部分,进行了2抽值,即尺度空间,低频概貌部分
}
for(i=0;i《4116;i++)
{
if(i《=259)
output[i]=AAAA4[i];
elseif(i《=519)
output[i]=DAAA4[i-260];
elseif(i《=1035)
output[i]=DAA3[i-520];
elseif(i《=2064)
output[i]=DA2[i-1036];
else
output[i]=DA2[i-2065];
}
*/
评论
查看更多