[MATLAB毕业设计] 关于ifft变换问题求教 - MATLAB爱好者论坛-LabFans.com

DMC3142 · 2012-04-06, 22:20

clear;
clc;
close all;
fid1=load('c:\hcn0001.txt');%读取数据
fid2=load('c:\hcnref0000.txt');%读取数据
subplot(221);
plot(fid1(:,1),fid1(:,2));%绘制hcn信号的时域图
fs=250000;%采样频率
N=16384;%FFT的数据长度
y1=fft(fid1(:,2),N);%对hcn信号进行快速傅里叶变换
f1=(0:N-1)*fs/N; %换算成实际的频率值
subplot(222);
plot(f1(1:N/2),abs(y1(1:N/2))*2/N);
title('y1的频域图. NFFT=16384');
xlabel('频率f');
ylabel('幅值A')
这里对傅里叶变换后的y1 进行反傅里叶变换但是在反傅里叶变换之前要把频率大于nyguist的置0 ，这该怎么表达。。
I=ifft1(y1(0:N/2));写可以么？只对前N/2个数进行反傅里叶变换。

另外就是关于反傅里叶变换后的到的实部和虚部该怎么表达呢？

2012-04-06, 22:20	#1
DMC3142 初级会员注册日期: 2012-04-06 年龄: 35 帖子: 1 声望力: 0	关于ifft变换问题求教 clear; clc; close all; fid1=load('c:\hcn0001.txt');%读取数据 fid2=load('c:\hcnref0000.txt');%读取数据 subplot(221); plot(fid1(:,1),fid1(:,2));%绘制hcn信号的时域图 fs=250000;%采样频率 N=16384;%FFT的数据长度 y1=fft(fid1(:,2),N);%对hcn信号进行快速傅里叶变换 f1=(0:N-1)fs/N; %换算成实际的频率值 subplot(222); plot(f1(1:N/2),abs(y1(1:N/2))2/N); title('y1的频域图. NFFT=16384'); xlabel('频率f'); ylabel('幅值A') 这里对傅里叶变换后的y1 进行反傅里叶变换但是在反傅里叶变换之前要把频率大于nyguist的置0 ，这该怎么表达。。 I=ifft1(y1(0:N/2));写可以么？只对前N/2个数进行反傅里叶变换。另外就是关于反傅里叶变换后的到的实部和虚部该怎么表达呢？