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个数进行反傅里叶变换。
另外就是关于反傅里叶变换后的到的实部和 虚部 该怎么表达呢？
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个数进行反傅里叶变换。
另外就是关于反傅里叶变换后的到的实部和 虚部 该怎么表达呢？