嗨,哈
1。
运行脚本的一个方法是窗口fft_nor和fft_mur f_nor和f_nur碰巧的长度短于各自的信号。
而不是
图(3)
情节(f_nor、abs (fft_nor));
使用
图(3)
情节(f_nor、abs (fft_nor([1 +元素个数(fft_nor)元素个数(f_nor):结束))));
或
情节(f_nor、abs (fft_nor([1:元素个数(f_nor)))));
而不是
图(4)
情节(f_mur、abs (fft_mur));
使用
情节(f_mur、abs (fft_mur([1:元素个数(f_mur)))));
或
情节(f_mur、abs (fft_nor([1 +元素个数(fft_mur)元素个数(f_mur):结束))));
2。
另一种方法是首先定义参数的时间以及频率上你想要的数量,然后,不是之前,计算FFT,设定频率的数量你想已经定义的频率,这样的
. .
t_nor = 6;% s
N_nor =长度(海军);
fs_nor = N_nor / t_nor;%赫兹
% fs_murmur
t_mur = 8;% s
N_mur =长度(杂音);
fs_mur = N_mur / t_mur;%赫兹
% %计算频率轴
f_nor = fs_nor / 2 * linspace (1, 1, fs_nor);
f_mur = fs_mur / 2 * linspace (1, 1, fs_mur);
% %傅里叶变换
fft_nor = fft(海军研究实验室,元素个数(f_nor));
fft_mur = fft(杂音元素个数(f_mur));
% %策划频谱
%正常的声音
图(3)
情节(f_nor、abs (fft_nor));
包含(的频率(赫兹),“字形大小”12“fontweight”,“b”);
ylabel (“级”,“字形大小”12“fontweight”,“b”);
标题(“正常信号的FFT级”,“字形大小”14岁的“fontweight”,“b”);
%低语的声音
图(4)
情节(f_mur、abs (fft_mur));
包含(的频率(赫兹),“字形大小”12“fontweight”,“b”);
ylabel (“级”,“字形大小”12“fontweight”,“b”);
标题(“杂音信号的FFT级”,“字形大小”14岁的“fontweight”,“b”);
3所示。
或者你也可以oversample:增加Fs在两个信号直到你得到相同数量的频率fft()分配默认情况下。
有不同的方法,但在这种情况下,事实上,如果你样品速度由以下因素
分别6和8、信号和杂音,我看到K1 K2因素引入f_nor和f_mur的定义
然后脚本还好,添加这样的:
% %傅里叶变换
fft_nor = fft(海军);
fft_mur = fft(杂音);
% %计算fs
% fs_nrl
t_nor = 6;% s
N_nor =长度(海军);
fs_nor = N_nor / t_nor;%赫兹
% fs_murmur
t_mur = 8;% s
N_mur =长度(杂音);
fs_mur = N_mur / t_mur;%赫兹
% %计算频率轴
K1 =元素个数(fft_nor) / fs_nor
K2 =元素个数(fft_mur) / fs_mur
f_nor = fs_nor / 2 * linspace (1, 1, fs_nor * K1);
f_mur = fs_mur / 2 * linspace (1, - 1, fs_mur * K2);
% %策划频谱
%正常的声音
图(3)
情节(f_nor、abs (fft_nor));
包含(的频率(赫兹),“字形大小”12“fontweight”,“b”);
ylabel (“级”,“字形大小”12“fontweight”,“b”);
标题(“正常信号的FFT级”,“字形大小”14岁的“fontweight”,“b”);
%低语的声音
图(4)
情节(f_mur、abs (fft_mur));
包含(的频率(赫兹),“字形大小”12“fontweight”,“b”);
ylabel (“级”,“字形大小”12“fontweight”,“b”);
标题(“杂音信号的FFT级”,“字形大小”14岁的“fontweight”,“b”);
如果你发现这些线有用请是考虑将我的答案标记为接受答案?
其他读者,如果你发现这个答案的任何帮助请点击竖起大拇指投票链接,
提前感谢你的时间和精力
约翰BG
