**%程序名:傅里叶变换的意义
%程序说明:
%本程序先感受以下为什么要学傅里叶变换。
% 首先产生两个时间上可以分开的信号(无噪声),然后在原信号上加上噪声,此时,信号时域已无法
% 观察出信号格式,没法很好的清晰的看出此信号的频率是多少,但是通过观察信号的傅里叶变化(频谱)
% ,即可以清除地观察信号分量
clear;%清除内存
close all;%关闭所有图形
x=sin(2*pi*(0:199)*0.5/10);%产生某时段的正弦信号
x(201:400)=sin(2*pi*(0:199)*2/10);%产生另一个时段的正弦信号
figure;
plot(x);%画图
grid on;%并打网格
title('时域可分信号','FontSize',14);
ylabel('x(n)','FontSize',14);
xlabel('n','FontSize',14);%标注X轴,Y轴坐标
x=x+0.8*randn(1,400);%信号加上噪声
figure;
plot(x);
grid on;
title('信号加噪声','FontSize',14);
ylabel('x(n)','FontSize',14);
xlabel('n','FontSize',14);
%标注X轴和Y轴的坐标
xf=fft(x,512);%FFT(X,N)是N点的FFT,如果X小于0,就用0填充
% 超过N个点,如果有更多的话就会被截断。
figure;
plot((0:255)*10/512,abs(xf(1:256)));
grid on;%画图
xlabel('KHZ','FontSize',14);
ylabel('信号频谱','FontSize',14);%标注XY轴坐标
title('信号频谱','FontSize',14);**
关于请教:
怎么确定这个频率就是5KHZ,不是MHZ,HZ???
第29行代码解释
根据代码的意思,是想把10当成采样率为10KHz