我试图在不使用Matlab内置函数(如fft() )的情况下实现一维DFT。这是我的密码
function [Xk] = dft1(xn)
N=length(xn);
n = 0:1:N-1; % row vector for n
k = 0:1:N-1; % row vecor for k
WN = exp(-1j*2*pi/N); % Twiddle factor (w)
nk = n'*k; % creates a N by N matrix of nk values
WNnk = WN .^ nk; % DFT matrix
Xk = (WNnk*xn );
在使用以下命令后运行代码时:
I = imread('sample.jpg')
R = dft1(I)
我得到了这个特殊的错误:使用* MTIMES的错误并不完全支持整数类。至少一个输入必须是标量。若要计算元素时间,请使用时间(.*)代替.
有人能帮我解决这个问题吗?
注:我还处于学习Matlab的初级阶段,非常感谢。
发布于 2014-11-30 23:38:08
您只需将数据强制转换为,然后再次运行代码即可。基本上,错误的意思是,当应用两个变量之间的矩阵乘法时,您试图将类数据混合在一起。具体来说,您在dft1
中定义的数值向量和矩阵都是double
类型,但是您的图像可能是通过imread
读取时的uint8
类型。这就是为什么您要得到这个整数错误,因为uint8
是一个整数类,并且您试图用这个数据类型和double
数据类型执行矩阵乘法。请记住,只要一个数字是单个数字/标量,就可以混合数据类型。这也是错误所暗示的。MATLAB中不支持非浮点变量的矩阵乘法(double
,single
),因此在应用算法之前,需要确保图像数据和DFT矩阵是相同类型的。
因此,只需:
I = imread('sample.jpg');
R = dft1(double(I));
小音符
这段代码非常聪明,它(默认情况下)将一维DFT应用于图像的所有列。输出将是与I
相同大小的矩阵,其中每一列都是来自I
的每一列的一维DFT结果。
这是要考虑的问题,但是如果要将其应用于映像的所有行,则只需将I
转到dft1
中,以便行变为列,并且可以对这些新的“列”进行操作。完成之后,您只需将结果转回来,这样您就可以从dft1
成形中获得输出,从而将结果应用于每一行。因此:
I = imread('sample.jpg');
R = dft1(double(I.')).';
希望这能有所帮助!祝好运!
https://stackoverflow.com/questions/27219235
复制相似问题