第一周作业：音频分析和RGB文件的熵计算

最新推荐文章于 2024-08-04 17:34:05 发布

Shelly_77

最新推荐文章于 2024-08-04 17:34:05 发布

阅读量94

点赞数

文章标签：大数据

本文链接：https://blog.csdn.net/Shelly_77/article/details/125472558

版权

音频分析

清音S

时域

频域

浊音U

时域

频域

爆破音 K

时域

频域

总结：

清音时域波形类似高斯分布，频域类似均匀分布，较为平缓。

浊音时域分布总体来看较为均匀，具备短时周期性；
浊音的频域分布起伏较大，有多个波峰，类似正弦函数分布。

爆破音时域波形总体看具有短暂的“爆破脉冲”，放大后观察到类似正弦函数特性。频谱较为均匀，起伏较小。有一个明显波峰，对应时域的脉冲。

计算RGB文件的熵

实验过程

打开、创建所需文件
定义文件读取、存储数组
计算三通道数值的出现概率
定义并求解三通道信源熵

#include<fstream>
#include<iostream>
#include<cmath>
using namespace std;
#define size 196608//256*256*3
#define height 256
#define width 256

int main()
{
	unsigned char R[256*256]={0},G[256*256]={0},B[256*256]={0};//RGB分量
	double R_F[256]={0},G_F[256]={0},B_F[256]={0};//RGB概率分布
	double R_I=0,G_I=0,B_I=0;//RGB的熵
	FILE* picture,* red,*green,*blue;
	fopen_s(&picture,"d:\\down.rgb","rb");
	fopen_s(&red,"d:\\red.txt","w");
	fopen_s(&green,"d:\\green.txt","w");
	fopen_s(&blue,"d:\\blue.txt","w");
	if(picture==0)
		printf("Fail to open the picture!");
	else
	{
		unsigned char buffer[size]={0};
		fread(buffer,1,size,picture);//将picture的值写入数组buffer
		//将B、G、R的值分别写入B、G、R的数组中
		for(int i=0,j=0;i<size;i=i+3,j++)
		{
			//以BGRBGR……顺序排列
			B[j]=*(buffer+i);
			G[j]=*(buffer+i+1);
			R[j]=*(buffer+i+2);
		}
	    //统计RGB各数值次数(8bit,范围0~255)
		for(int i=0;i<256;i++)
		{
			for(int j=0;j<256*256;j++)
			{
				if(int(R[j]==i)) {R_F[i]++;}
				if(int(G[j]==i)) {G_F[i]++;}
				if(int(B[j]==i)) {B_F[i]++;}
			}
		}
		//计算RGB频率
		for(int i=0;i<256;i++)
		{
			R_F[i]=R_F[i]/(256*256);
			G_F[i]=G_F[i]/(256*256);
			B_F[i]=B_F[i]/(256*256);
		}
		//将概率写入文件
        for (int i = 0; i < 256; i++)
        {
            fprintf(red, "%d\t%f\n", i, R_F[i]);
            fprintf(green, "%d\t%f\n", i, G_F[i]);
            fprintf(blue, "%d\t%f\n", i, B_F[i]);
        }
		//计算RGB的熵
        for (int i = 0; i < 256; i++)
        {
            if (R_F[i] != 0) { R_I += -R_F[i] * log(R_F[i]) / log((double)2); }//以2为底的对数
            if (G_F[i] != 0) { G_I += -G_F[i] * log(G_F[i]) / log((double)2); }
            if (B_F[i] != 0) { B_I += -B_F[i] * log(B_F[i]) / log((double)2); }
        }
        printf("R的熵为%f\n", R_I);
        printf("G的熵为%f\n", G_I); 
        printf("B的熵为%f\n", B_I);
    }
    return 0;
}