通信人家园

标题: 全相位OFDM系统仿真代码 [查看完整版帖子] [打印本页]

时间: 2013-5-28 20:13

作者: niuxiaohua0123 标题: 全相位OFDM系统仿真代码

求MATLAB 大神帮忙啊基于全相位FFT的仿真代码！

时间: 2013-10-9 20:22

作者: fanyiyi

谢谢您的帮助，望通信人精神代代相传！

时间: 2015-7-10 09:47

作者: dab123

谢谢分享！！

时间: 2015-7-11 18:47

作者: zhwang554

good
谢谢分享！！

时间: 2015-7-11 19:37

作者: zhwang554

本帖最后由 zhwang554 于 2015-7-12 02:49 编辑

apFFT_Spectrum_Est

close all;clc;clear all;
N=256;
fs=256;
t=-N+1:N-1;
%t=0:2*N-1;
y= 0.0*rand(size(t))+2*cos(2*1*pi*t*100.3/fs+pi/2)+2*cos(1*2*pi*t*50.3/fs+pi/3);
y1=y(N : (2*N-1));

win=hanning(N)';
win1=win/sum(win);
y11=y1.*win1;

y11_fft=fft(y11,N);

a1=abs(y11_fft);
p1=mod(angle(y11_fft)*180/pi,360);
y2=y(1 : (2*N-1));
winn=conv(win,win);
win2=winn/sum(winn);
y22=y2.*win2;

y222=y22(N:end)+[0 y22(1:N-1)];

y2_fft=fft(y222,N);
a2=abs(y2_fft);
p2=mod(angle(y2_fft)*180/pi,360);
p21=p2(51)
p22=p2(101)
tt=0:N-1;

%figure('color','w');
figure(1)
subplot(2,1,1)
plot(tt,10*log10(a1),'k');
title('FFT amplitude');
ylim([-150,10]);
xlim([0,N/2]);
xlabel('frequency(Hz)');
ylabel('dB');
grid

subplot(2,1,2)
plot(tt,10*log10(a2),'k');
title('APFFT amplitude');
ylim([-150,10]);
xlim([0,N/2]);
xlabel('frequency(Hz)');
ylabel('dB');
grid

%figure('color','w');
figure(2)
subplot(2,1,1)
plot(tt,p1,'k.');

title('FFT phase');
ylim([0,400]);
xlim([0,N/2]);
xlabel('frequency(Hz)');
ylabel('degree');
grid

subplot(2,1,2)
plot(tt,p2,'k.');
title('APFFT phase');
ylim([0,400]);
xlim([0,N/2]);
xlabel('frequency(Hz)');
ylabel('degree');
grid

figure(3);
plot(y,'b.-');
grid

时间: 2015-7-11 19:37

作者: zhwang554

本帖最后由 zhwang554 于 2015-7-11 19:39 编辑

apFFT_Spectrum

时间: 2020-3-6 17:52

作者: Leo妮子

回帖

通信人家园 (https://www.txrjy.com/)