matlab中如何对傅立叶变换后的函数作图
傅立叶变换是对符号变量进行的变换,我想画出变换后函数的图形,请问应当用哪一个函数呢?我用plot函数好像不可以啊比如我想画出rect函数变换后的sinc函数请给出具体例子...
傅立叶变换是对符号变量进行的变换,我想画出变换后函数的图形,请问应当用哪一个函数呢?我用plot 函数好像不可以啊
比如我想画出rect函数变换后的sinc函数
请给出具体例子,谢谢了,很着急 展开
比如我想画出rect函数变换后的sinc函数
请给出具体例子,谢谢了,很着急 展开
3个回答
展开全部
离散傅里叶变化是使用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;
你现在是处理离散函数的傅里叶变化
%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;
展开全部
使用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;
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;
已赞过
已踩过<
评论
收起
你对这个回答的评价是?
展开全部
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)')
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)')
已赞过
已踩过<
评论
收起
你对这个回答的评价是?
推荐律师服务:
若未解决您的问题,请您详细描述您的问题,通过百度律临进行免费专业咨询