MATLAB中关于FFT的问题
现成的FFT程序网上可以下到很多，但如何定义相似度需要考虑一下。可以考虑用相对变化率的形式（相对变化率=（数值1-数值2）/数值1 或者 相对变化率=（数值1-数值2）/数值2 ）。有一种情况跟你的需求很像：设计滤波器后说明滤波效果。这种情况，需要比较变化前后的信号进行幅频特性、相频特性曲线，以此说明滤波效果。在这一过程中就需要对两个信号分别进行FFT变换，以求得幅频曲线和相频曲线。具体过程如下：step1：将横坐标定义为t（matlab赋值语句t=[*,*,*];），将纵坐标定义为y（matlab赋值语句y=[*,*,*];），采样频率就是临近两个横坐标差值的倒数（一般横坐标都应为时间）；step2：新建一个m文件（快捷键为ctrl+N）；step3：将如下程序复制到m文件中 t=[*,*,*];%填入横轴数据 y=[*,*,*];%填入纵轴数据 N = size(t,2);%行向量时，列向量是为N = size(t,1); N = 2^(nextpow2(N)-1); Y = fft(y,N); mag_Y = abs(Y)/N*2;%各个频率点处的幅值，这个就是fft变换后的数据 f = fs/2*linspace(0,1,N/2+1);%对应的频率值 plot(f,mag_Y(1:N/2+1))%显示变换后的曲线step4：亲测实例 t = 1:0.001:10; y= 2*sin(2*pi*100*t);%频率为100Hz %=================以下同step3中代码============== N = size(t,2);%行向量时，列向量是为N = size(t,1); N = 2^(nextpow2(N)-1); Y = fft(y,N); mag_Y = abs(Y)/N*2;%各个频率点处的幅值，这个就是fft变换后的数据 f = fs/2*linspace(0,1,N/2+1);%对应的频率值 plot(f,mag_Y(1:N/2+1))%显示变换后的曲线 %=================以上同step3中代码============== 结果：由于matlab自带FFT函数有些小瑕疵，所以在100Hz处幅值虽然接近2，但还是有些偏差的，对于偏差的修正网上也有相应的方法，如果需要在留言。
matlab如何用fft
matlab自带的fft函数是快速傅里叶变换函数。主要用于降噪处理，通过使用傅里叶变换求噪声中隐藏的信号的频率分量。该函数使用方法：方法一：Y = fft(X) 用快速傅里叶变换 (FFT) 算法计算 X 的离散傅里叶变换 (DFT)。如果 X 是向量，则 fft(X) 返回该向量的傅里叶变换。如果 X 是矩阵，则 fft(X) 将 X 的各列视为向量，并返回每列的傅里叶变换。如果 X 是一个多维数组，则 fft(X) 将沿大小不等于 1 的第一个数组维度的值视为向量，并返回每个向量的傅里叶变换。方法二：Y = fft(X,n) 返回 n 点 DFT。如果未指定任何值，则 Y 的大小与 X 相同。如果 X 是向量且 X 的长度小于 n，则为 X 补上尾零以达到长度 n。如果 X 是向量且 X 的长度大于 n，则对 X 进行截断以达到长度 n。如果 X 是矩阵，则每列的处理与在向量情况下相同。如果 X 为多维数组，则大小不等于 1 的第一个数组维度的处理与在向量情况下相同。我们通过下例，来了解fft函数使用过程：第一步、指定信号的参数，采样频率为 1 kHz，信号持续时间为 1.5 秒。Fs=1000；%采样频率T=1/Fs；%采样周期L=1500；%信号长度t=（0:L-1）*T；%时间向量第二步、构造一个信号，其中包含幅值为 0.7 的 50 Hz 正弦量和幅值为 1 的 120 Hz 正弦量。S = 0.7*sin(2*pi*50*t) + sin(2*pi*120*t);第三步、用均值为零、方差为 4 的白噪声扰乱该信号。X = S + 2*randn(size(t));第四步、在时域中绘制含噪信号。通过查看信号 X(t) 很难确定频率分量。plot(1000*t(1:50),X(1:50))title('Signal Corrupted with Zero-Mean Random Noise')xlabel('t (milliseconds)')，ylabel('X(t)')第五步、计算信号的傅里叶变换。Y = fft(X);第六步、计算双侧频谱 P2， 计算单侧频谱 P1。P2 = abs(Y/L); P1 = P2(1:L/2+1);P1(2:end-1) = 2*P1(2:end-1)第七步、定义频域 f 并绘制单侧幅值频谱 P1f = Fs*(0:(L/2))/L;plot(f,P1) title('Single-Sided Amplitude Spectrum of X(t)')xlabel('f (Hz)')，ylabel('|P1(f)|')运行结果。
fft 采样频率Fs和采样点数N应该怎么确定呢 matlab
1、首先点击Matlab数据导入菜单-->选择Excel文件-->打开。2、导入-->数据类型-->数值矩阵-->导入所选数据。3、然后在Matlab工作区即可看到导入的数据是一个矩阵。4、然后新建一个simulink仿真模型-->器件库(simulink Library Browser)-->找到From warkspace拖到模型中-->找到Powergui拖到模型中-->找到示波器拖到模型中。5、在仿真模型中点击simulation-->Modle configuration parameters。6、solve->设置如下图。7、打开From Workspaces设置参数为y-->采样时间为2e-6-->打开示波器设置参数如图所示。8、最后运行仿真-->打开powergui-->Tools-->FFT Analysis-->选择信号设置参数即可-->点击Display。
