Matlab自带这个函数，使用方法：
y = dct(x)
y = dct(x,n)
举例：
x = (1:100) + 50*cos((1:100)*2*pi/40);
X = dct(x);
[XX,ind] = sort(abs(X)); ind = fliplr(ind);
i = 1;
while (norm([X(ind(1:i)) zeros(1,100-i)])/norm(X)<.99)
i = i + 1;
end
i =
3
内容全部来源于网络收集，如有侵权，请联系网站删除：QQ号：24596024