matlabfft简单例子(matlabfft简单例子)
matlabfft简单例子(matlabfft简单例子)
N=128;
N=0:N-1;
T=n/fs;
F0=10;% sets sine signal frequency
% generates sinusoidal signals
X=sin (2*pi*f0*t);
Figure (1);
Subplot (231);
Plot (T, x)% sinusoidal waveform in time domain
Xlabel ('t');
Ylabel ('y');
Title (` sine signal, y=2*pi*10t time domain waveform ');
Grid;
% to do FFT transform and do frequency spectrum
Y=fft (x, N);% performs the FFT transform
Mag=abs (y);% for magnitude
F= (0:length (y) -1)'*fs/length (y);% corresponding frequency
conversion
Figure (1);
Subplot (232);
Plot (F, MAG):% do spectrogram
Axis ([0100,0,80]);
Xlabel (` frequency (Hz) ');
Ylabel ('amplitude');
Title (sine signal y=2*pi*10t amplitude spectrum N=128');
Grid;
MATLAB's FFT program 2009-11-08 16:30clc
X% input signal
M=nextpow2 (x);% x length corresponds to the minimum power m
of 2
N=2^M;
If length (x)
题
快递公司问题件快递公司问题件货款处理关于圆的周长面积重点题型关于解方程组的题及答案关于南海问题
(“正弦信号y = 2×π×10t功率谱”);
网格;
%求对数谱
= log(SQ);
图(1);
情节(235);
情节(f,LN);
xlabel('频率(Hz)”);
ylabel('对数谱”);
标题(“正弦信号y = 2×π×10t对数谱”);
网格;
用IFFT恢复原始信号%
xifft = IFFT(Y);
magx =房(xifft);
Ti = [ 0:长度(xifft)- 1 ] / FS;
图(1);
情节(236);
情节(Ti,magx);
xlabel(不);
ylabel(Y);
标题(“通过IFFT转换的正弦信号波形”);
网格;
****************矩形波**************** % % 2。
FS = 10;%设定采样频率
T = - 5:0.1:5;
rectpuls(t,x = 2);
X = X(1:99);
图(2);
情节(231);
情节(T(1:99),x);%作矩形波的时域波形
xlabel(不);
ylabel(Y);
标题(“矩形波时域波形”);
网格;
进行FFT变换并做频谱图%
Y = FFT(X);%进行FFT变换
MAG = ABS(Y);%求幅值
F =(0:长度(Y)- 1)×FS /长度(Y);%进行对应的频率转换
图(2);
情节(232);
图(F,MAG);%做频谱图
xlabel('频率(Hz)”);
ylabel('幅值”);
标题(“矩形波幅频谱图”);
网格;
%求均方根谱
= abs(y);
图(2);
情节(233);
情节(f,平方);
xlabel('频率(Hz)”);
ylabel('均方根谱”);
标题(“矩形波均方根谱”);
网格;
%求功率谱
功率=平方厘米2;
图(2);
情节(234);
情节(F,权力);
xlabel('频率(Hz)”);
ylabel('功率谱”);
标题(“矩形波功率谱”);
网格;
%求对数谱
= log(SQ);
图(2);
情节(235);
情节(f,LN);
xlabel('频率(Hz)”);
ylabel('对数谱”);
标题(“矩形波对数谱”);
网格;
用IFFT恢复原始信号%
xifft = IFFT(Y);
magx =房(xifft);
Ti = [ 0:长度(xifft)- 1 ] / FS;
图(2);
情节(236);
情节(Ti,magx);
xlabel(不);
ylabel(Y);
标题(“通过IFFT转换的矩形波波形”);
网格;
****************白噪声**************** % % 3。
FS = 10;%设定采样频率
T = - 5:0.1:5;
x,0(1100);
x(50),100000;
图(3);
情节(231);
情节(T(1:100),x);%作白噪声的时域波形
xlabel(不);
ylabel(Y);
标题(“白噪声时域波形”);
网格;
进行FFT变换并做频谱图%
Y = FFT(X);%进行FFT变换
MAG = ABS(Y);%求幅值
F =(0:长度(Y)- 1)×FS /长度(Y);%进行对应的频率转换
图(3);
情节(232);
图(F,MAG);%做频谱图
xlabel('频率(Hz)”);
ylabel('幅值”);
标题(“白噪声幅频谱图”);
网格;
%求均方根谱
= abs(y);
图(3);
情节(233);
情节(f,平方);
xlabel('频率(Hz)”);
ylabel('均方根谱”);
标题(“白噪声均方根谱”);
网格;
%求功率谱
功率=平方厘米2;
图(3);
情节(234);
情节(F,权力);
xlabel('频率(Hz)”);
ylabel('功率谱”);
标题(“白噪声功率谱”);
网格;
%求对数谱
= log(SQ);
图(3);
情节(235);
情节(f,LN);
xlabel('频率(Hz)”);
ylabel('对数谱”);
标题(“白噪声对数谱”);
网格;
用IFFT恢复原始信号%
ifft xifft = (y);
magx = real (xifft);
ti = [0: length (xifft) - 1] / fs;
figure (3);
subplot (236);
plot (ti, magx);
xlabel ("t");
ylabel ('y');
title ('通过ifft转换的白噪声波形');
grid;
fs = 60;
n = 0: 1: 1 359 * / fs / fs;
x = sin (2 * pi * 20 * n) + sin (2 * pi * 15 * n);
n = length (n);
fft (x = x);
f = (0: n - 1) * fs / n;
px1 = abs (x) / n,% 归一化
subplot (211);
plot (d, px1);
grid on;
title ('幅度谱');
xlabel ('-' 频率f hz);
ylabel (幅度 '');
pha = angle (x);
subplot (212);
plot (d, pha);
title ('相位谱');
xlabel ('-' 频率f hz);
ylabel ("相位 - - 弧度 ');
r = length (x);
for i: = 1 to r - 1
for j = i + 1: r
if x > x (i), (j).
temp = x (i);
x = x (i), (j);
x (j) = temp;
end
end
end
y = x;
以上为按照升序排列的, 若要降序, 则
if x < x (i), (j).
temp = x (i);
x = x (i), (j);
x (j) = temp;
冒泡%
the clc; clear;
x = rand (1,10);
for i = 1: length (x)
for j = 1: length (x) - (i)
if (x > x (j), (j + 1))
t = x (j);
x (j) = x (j + 1);
x (j + 1) = t;
end
end
end
disp (x)
图自己画
core = [1 2 3 5 7 9 13 14 15]
n = _ core core
n = length (core)
n (1) = _ core 将第一个元素移除% []
_ core (n) = n (n)% 最后一元素与mag最后一元素相同 core
_ = n _ del core core core 求二者之差%