離散余弦變換，DCT for Discrete Cosine Transform。
是與傅里葉變換相關(guān)的一種變換，它類似于離散傅里葉變換(DFT for Discrete Fourier Transform),但是只使用實數(shù)。
離散余弦變換相當(dāng)于一個長度大概是它兩倍的離散傅里葉變換，這個離散傅里葉變換是對一個實偶函數(shù)進行的（因為一個實偶函數(shù)的傅里葉變換仍然是一個實偶函數(shù)），在有些變形里面需要將輸入或者輸出的位置移動半個單位(DCT有8種標(biāo)準(zhǔn)類型，其中4種是常見的)。

離散余弦變換DCT具有信號譜分量豐富、能量集中，且不需要對語音相位進行估算等優(yōu)點，能在較低的運算復(fù)雜度下取得較好的語音增強效果。

離散余弦變換，尤其是它的第二種類型，經(jīng)常被信號處理和圖像處理使用，用于對信號和圖像(包括靜止圖像和運動圖像)進行有損數(shù)據(jù)壓縮。這是由于離散余弦變換具有很強的"能量集中"特性:大多數(shù)的自然信號(包括聲音和圖像)的能量都集中在離散余弦變換后的低頻部分；

在mpeg協(xié)議或jpeg格式的編碼中包含離散余弦變換，記不太清，可以自己查；

下面來看C語言實現(xiàn)；

#include <windows.h> #include <math.h> #include "resource.h"LRESULT CALLBACK WindowProcedure (HWND, UINT, WPARAM, LPARAM);HINSTANCE hInst; TCHAR szClassName[] = TEXT("dctDemo"); #define PI 3.141593 #define N 4void dct(int x[4], HDC ); int WINAPI WinMain (HINSTANCE hThisInstance,HINSTANCE hPrevInstance,LPSTR lpszArgument,int nFunsterStil) {HWND hwnd;MSG messages;WNDCLASSEX wincl;hInst = hThisInstance;wincl.hInstance = hThisInstance;wincl.lpszClassName = szClassName;wincl.lpfnWndProc = WindowProcedure;wincl.style = CS_DBLCLKS;wincl.cbSize = sizeof (WNDCLASSEX);wincl.hIcon = LoadIcon (NULL, IDI_APPLICATION);wincl.hIconSm = LoadIcon (NULL, IDI_APPLICATION);wincl.hCursor = LoadCursor (NULL, IDC_ARROW);wincl.lpszMenuName = MAKEINTRESOURCE (IDC_DCTDEMO);wincl.cbClsExtra = 0;wincl.cbWndExtra = 0;wincl.hbrBackground = (HBRUSH)(COLOR_WINDOW+1);if (!RegisterClassEx (&wincl))return 0;hwnd = CreateWindowEx (0,szClassName,TEXT("離散余弦變換 Demo"),WS_OVERLAPPEDWINDOW,CW_USEDEFAULT,CW_USEDEFAULT,300,300,HWND_DESKTOP,NULL,hThisInstance,NULL);ShowWindow (hwnd, nFunsterStil);while (GetMessage (&messages, NULL, 0, 0)){TranslateMessage(&messages);DispatchMessage(&messages);}return messages.wParam; }LRESULT CALLBACK WindowProcedure (HWND hwnd, UINT message, WPARAM wParam, LPARAM lParam) {PAINTSTRUCT ps;HDC hdc;RECT rt; char szBuffer[100]; int x1[] = {2,3,3,2};switch (message){case WM_COMMAND:switch (LOWORD(wParam)){case IDM_dct: hdc=GetDC(hwnd);dct(x1,hdc);break;case IDM_ABOUT:MessageBox (hwnd, TEXT ("dctDemo v1.0\nCopyright (C) 2020\n by bo"),TEXT ("離散余弦變換Demo"), MB_OK | MB_ICONINFORMATION);break;case IDM_EXIT:DestroyWindow(hwnd);break;default:return DefWindowProc(hwnd, message, wParam, lParam); }break;case WM_CREATE:break;case WM_PAINT:hdc = BeginPaint(hwnd, &ps); GetClientRect(hwnd, &rt); EndPaint(hwnd, &ps);break;case WM_DESTROY:PostQuitMessage (0);break;default:return DefWindowProc (hwnd, message, wParam, lParam);}return 0; }void dct(int x[4], HDC hdc) {int row=0;char szBuffer[20];int k,n;double A,s,X[4];for(k = 0; k < N; k++){s = 0;if(k==0)A = sqrt(1.0/N); //計算k=0時的系數(shù)elseA = sqrt(2.0/N); //計算k!=0時的系數(shù)for(n = 0; n < N; n++){double tmp = x[n]*cos((PI*(2*n+1)*k)/(2*N));s = s + tmp; //累加求和}X[k] = A * s; //X[k]等于累和結(jié)果s乘以系數(shù)A}for(k = 0; k < N; k++) {sprintf(szBuffer,"%.7lf",X[k]);TextOut(hdc,10,10+row*25,szBuffer,lstrlen(szBuffer));row=row+1;}}

這是計算4個點的DCT；要計算的值序列作為數(shù)組傳給DCT函數(shù)；然后輸出；

控制臺的版本輸出如下；?

控制臺版本代碼；

?

#include<stdio.h> #include<math.h> #define PI 3.141593 #define N 4 int main() {int k, n;int x[] = {2,3,3,2};double A, s, X[4];for(k = 0; k < N; k++){s = 0;if(k==0)A = sqrt(1.0/N); //計算k=0時的系數(shù)elseA = sqrt(2.0/N); //計算k!=0時的系數(shù)for(n = 0; n < N; n++){double tmp = x[n]*cos((PI*(2*n+1)*k)/(2*N));s = s + tmp; //累加求和}X[k] = A * s; //X[k]等于累和結(jié)果s乘以系數(shù)A}for(k = 0; k < N; k++) printf("%f\n", X[k]); }

輸出的時候如果小數(shù)點格式化不對就會出現(xiàn)問題；正確的應(yīng)該是，

sprintf(szBuffer,"%.6lf",X[k]);?

工程；?

資源和頭文件；

#include "resource.h" #include <windows.h>/ // // Menu //IDC_DCTDEMO MENU BEGINPOPUP "&File"BEGINMENUITEM "離散余弦變換Demo", IDM_dctMENUITEM "E&xit", IDM_EXITENDPOPUP "&Help"BEGINMENUITEM "&About ...", IDM_ABOUTEND END #define IDM_EXIT 10001 #define IDM_ABOUT 10002#define IDC_DCTDEMO 10101 #define IDD_ABOUTBOX 10102 #define IDM_dct 40001

?

《新程序員》：云原生和全面數(shù)字化實踐50位技術(shù)專家共同創(chuàng)作，文字、視頻、音頻交互閱讀

總結(jié)

以上是生活随笔為你收集整理的离散余弦变换和C语言实现-win32版本的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。