做一个MATLAB程序，将时间转换为频率

i=imread('');

i=double(i);

ffti=fft(i(:)

rr=real(ffti);以傅里叶变换的实部为例。

ii=imag(ffti);以傅里叶变换的虚部为例。

a=sqrt(rr.^2+ii.^2);% 来计算光谱振幅。

figure,plot(a)

fft2=fft2(i);二维离散傅里叶变换。

sfft2=fftshift(fft2);直流分量移动到频谱的中心。

rrr=real(sfft2);以傅里叶变换的实部为例。

iii=imag(sfft2);以傅里叶变换的虚部为例。

a=sqrt(rrr.^2+iii.^2);% 来计算光谱振幅。

a=(a-min(min(a)))/(max(max(a))-min(min(a)))225;% 归一化。

figure;% 设置窗口。

imshow(a);显示原始图像的光谱。

在 MATLAB 中，有一个 FFT 快速傅里叶变换。

在 MATLAB 中，您可以使用傅里叶变换从时域变换到频域，然后使用逆傅里叶变换从频域返回到时域。 对于时域信号，可以使用FFT函数执行傅里叶变换，将其转换为频域信号。 然后，利用IFFT函数进行逆傅里叶变换，将频域信号转换回时域信号。

需要注意的是，如果信号在傅里叶变换过程中没有填零，则在逆变换后信号可能会被截断或失真，需要填零或其他技术处理。

可以使用快速傅里叶变换：fft。

在 MATLAB 命令窗口中，键入：help fft

MeshGrid 用于生成网格矩阵，这仅仅意味着生成给定 x 和 y 中的一对元素，因此每对 （x，y） 计算相应的 z，显然这会产生一个 z 表面。 但是这种说法是没有必要的，有时候我们只是想得到一条三维曲线，我们不想知道所有x、y元素成对组合的结果是什么，而且组合已经定义好了。

close all

ticn = 10;

x = 1:n;%x 坐标。

y = 1:n;%y 坐标。

%meshgrid 演示。

x, y] = meshgrid(x, y);meshgrid 函数用于生成网格矩阵。

z = x.^2 + y;

figure

mesh(x, y, z);

grid on

xlabel('x');

ylabel('y');

zlabel('z');

%不带网格网格。

z = x.^2 + y;

figure

mesh(x, y, z);当没有网格网格来生成基本矩阵时，这句话是错误的。

plot3(x, y, z);一组 （x， y） 对应于 z 值，因此 x 和 y 元素的数量必须相同。

grid on

xlabel('x');

ylabel('y');

zlabel('z');

我的意见是：

1.程序不一定能帮助你。

程序刚刚运行，但你仍然不知道它是如何工作的，所以你无法判断它是否正确。 为了理解FFT，需要对公式进行深入分析。

2.如何获得振幅为 1 的复正弦波？

使用 IFFT （u， 1024） 是正确的，它给出的时域波形为 1024 个点。 关键是 u 是多少。

如果您是 [1024,0,0,...，0]，1023 0s，则结果是一条频率为 0 的复正弦直线。

如果您是 [0,1024,0,...，0]，1023 个零点，则结果是一条频率为 1 的复正弦的直线。

如果您是 [0,0,1024,0,...，0]，1023 个零，则结果是一条频率为 2 的复正弦的直线。

等等。 3.如何获得具有指定振幅、指定相位为 45°、指定频率为 8 的复正弦波？

u=[0，…，0，1024*，0，…，0]，前面有 8 个零，后面有 1015 个零。能。

未完待续，先除法]。