matlab中如何对傅立叶变换后的函数作图

发布网友 发布时间：2022-04-24 18:24

我来回答

共3个回答

热心网友 时间：2023-10-31 19:38

离散傅里叶变化是使用fft函数而不是fourier

你现在是处理离散函数的傅里叶变化
%by dynamic
%see also http://www.matlabsky.com
%2009.2.18
%
function matlabksy
x=-2:.1:2;
y=rect(x);
%syms x y w
f=fft(y) ;
%由于傅里叶变化后会出现复数
%我不知道你是想绘制那部分的图形
%%绘制幅值 ——一般都是这个
%plot(abs(f))
%x实部 y虚部
%plot(f)
%只是画实部
%plot(real(f))
%只是绘制虚部
%plot(img(f))

function y=rect(x);
% 矩形函数 rect
if length(size(x))>2;
error('the size of x must less than 3')
end
y=zeros(size(x));
y(abs(x)<0.5)=1;

热心网友 时间：2023-10-31 19:39

使用fourier 进行傅里叶变换后，使用如下函数进行作图处理：
function matlabksy
x=-2:.1:2;
y=rect(x);
%syms x y w
f=fft(y) ;
%由于傅里叶变化后会出现复数
%%绘制幅值
plot(abs(f))
%x实部 y虚部
plot(f)
%只是画实部
plot(real(f))
%只是绘制虚部
plot(img(f))

function y=rect(x);
% 矩形函数 rect
if length(size(x))>2;
error('the size of x must less than 3')
end
y=zeros(size(x));
y(abs(x)<0.5)=1;

热心网友 时间：2023-10-31 19:39

t = 0:0.001:0.6;
x = sin(2*pi*50*t)+sin(2*pi*120*t);
y = x + 2*randn(size(t));
plot(1000*t(1:50),y(1:50))
title('Signal Corrupted with Zero-Mean Random Noise')
xlabel('time (milliseconds)')

Y = fft(y,512);

Pyy = Y.* conj(Y) / 512;

f = 1000*(0:256)/512;
plot(f,Pyy(1:257))
title('Frequency content of y')
xlabel('frequency (Hz)')